- Research article
- Open Access
Immune-related genes in tumor-specific CD4+ and CD8+ T cells in colon cancer
BMC Cancer volume 20, Article number: 585 (2020)
Immune escape is an immunological mechanism underlying tumorigenesis, and T cells play an important role in this process. In this study, immune-related genes were evaluated in tumor-infiltrating CD4+ and CD8+ T cells in colon cancer.
ESTIMATE was used to calculate stromal and immune scores for tumor datasets downloaded from The Cancer Genome Atlas–Colon Cancer (COAD). Differentially expressed genes (DEGs) between samples with high and low stromal and immune scores were screened, followed by a functional enrichment analysis of the overlapping DEGs. The DEGs related to CD4+ and the CD8+ T cells were then screened. Predicted miRNA–mRNA and lncRNA–miRNA pairs were used to construct a competing endogenous RNA (ceRNA) network. Furthermore, chemical–gene interactions were predicted for genes in the ceRNA network. Kaplan–Meier survival curves were also plotted.
In total, 83 stromal-related DEGs (5 up-regulated and 78 down-regulated) and 1270 immune-related DEGs (807 up-regulated and 293 down-regulated genes) were detected. The 79 overlapping DEGs were enriched for 39 biological process terms. Furthermore, 79 CD4+ T cell-related genes and 8 CD8+ T cell-related genes, such as ELK3, were screened. Additionally, ADAD1 and DLG3, related to CD4+ T cells, were significantly associated with the prognosis of patients with colon cancer. The chr22-38_28785274–29,006,793.1–miR-106a-5p-DDHD1 and chr22-38_28785274–29,006,793.1–miR-4319-GRHL1 axes obtained from CD4+ and CD8+ T cell-related ceRNAs were identified as candidates for further studies.
ELK3 is a candidate immune-related gene in colon cancer. The chr22-38_28785274–29,006,793.1–miR-106a-5p-DDHD1 and chr22-38_28785274–29,006,793.1–miR-4319-GRHL1 axes may be related to CD4+ and CD8+ T cell infiltration in colon cancer.
ELK3 is a candidate immune-related gene in colon cancer
The chr22-38_28785274–29,006,793.1–miR-106a-5p–DDHD1 axis was identified in an analysis of CD4+ T cell-related ceRNA networks
The chr22-38_28785274–29,006,793.1–miR-4319–GRHL1 axis was identified in an analysis of CD8+ T cell-related ceRNA networks
The CD4+ T cell-related genes ADAD1 and DLG3 were associated with prognosis
A total of 175 chemical–target pairs in CD4+ T cells and 9 in CD8+ T cells were obtained
Colon cancer is among the worst cancers with respect to mortality and incidence worldwide, especially in Asia . Despite considerable advances in surgical and adjuvant therapy for colon cancer, rates of recurrence for patients with stages I–III and stage IV cancer are 30 and 65%, respectively . Various issues, including resistance, relapse, and metastasis occurring after traditional radiotherapy, chemotherapy, and new targeted drug treatments, have revealed that tumors are a systemic disease and not simply a result of mutations in oncogenes and the inactivation of tumor suppressor genes [3, 4]. Tumor escape from immune surveillance mechanisms has prompted the advent of tumor immunotherapy [5, 6]. Immune evasion is one of the most important characteristics of tumors . It mainly occurs by the modification of tumor cells and changes in the tumor microenvironment. As such, an understanding of the mechanism underlying tumor immune escape can provide new strategies for immunotherapy . The T cell-mediated immune response against tumors is the basis for cancer immunotherapy and is correlated with favorable disease outcomes [9, 10].
Tumor immune escape is related to a decline in T cell responses, mainly manifested by immune tolerance to CD4+ T cells and the inhibition of CD8+ T cell activation . CD4+ T helper cells can assist in the activation of naive CD8+ T cells and can help to eliminate major histocompatibility antigen class II (MHC-II)-negative tumor cells. However, tumor cells can induce specific immune tolerance in CD4+ T cells . During the antitumor immune response, CD8+ T cells play a major role in directly killing tumor cells by recognizing tumor antigens. However, the local tumor microenvironment contains a large number of cytokines that individually or synergistically affect the activation of cytotoxic T lymphocytes (CTL) and the sensitivity of tumor cells to CTL activity . Via interactions among various immune elements, cancer cells may enter a dormant state or exhibit immune evasion, which may directly promote tumor development and progression. Therefore, we analyzed the molecular mechanisms associated with CD4+ and CD8+ T cell activity in colon cancer to explore the potential benefits of immunotherapy.
To provide a basis for the investigation of colon cancer-induced T cell-mediated immune escape, miRNA target prediction and competing endogenous RNA (ceRNA) network construction were performed. Furthermore, prognostic targets were evaluated, and drug screening was performed to screen prognostic indicators and strategies for colon cancer immunotherapy, which may provide guidance for clinical decision-making. In our study, datasets for colon cancer samples were downloaded from The Cancer Genome Atlas (TCGA) public database. Moreover, we performed an enrichment analysis of the differentially expressed genes (DEGs) correlated with both stromal and immune scores. CD4+ and CD8+ T cell-related DEGs were obtained and used for analyses of protein–protein interactions (PPI). Additionally, the ceRNA networks for CD4+ and CD8+ T cells were analyzed and small chemical molecules related to the DEGs in the ceRNA network were predicted. Finally, a survival analysis of the CD4+ and CD8+ T cell-related DEGs was performed.
The GDC TCGA Colon Cancer (COAD) dataset (version 07-19-2019) was downloaded from TCGA (https://xenabrowser.net/), including RNAseq FPKM data and corresponding clinical phenotypes. The dataset included data for 448 colon cancer samples. The data were analyzed according to the workflow illustrated in Fig. 1.
Analysis of stromal and immune scores
For the mRNA expression profile analysis, gene expression values were calculated by mapping probes (obtained from the microarray dataset and the annotation files for the chip platform) to gene symbols. The average mRNA expression level was obtained when multiple probes matched to one symbol.
To analyze stromal and immune cell infiltration in tumor tissues, Estimation of Stromal and Immune cells in Malignant Tumours using Expression data (ESTIMATE, version 1.0.13)  in R was used to predict tumor purity. Thus, the stromal and immune scores for each tumor sample were obtained.
wAnalysis of differentially expressed genes
Based on the median stromal score, the tumor samples were divided into groups with high and stromal scores. Likewise, the tumor samples were divided into groups with high and low immune scores according to the median immune score.
The typical Bayesian test in the limma package (Version 3.10.3)  was used to analyze differentially expressed mRNAs (dif-mRNAs) between the two stromal score groups and two immune score groups. mRNAs with a fold change value of more than 0.263 and a p-value of less than 0.05 were selected as dif-mRNAs.
To identify potential regulatory genes associated with both stromal and immune cell contents, the overlapping DEGs were selected and represented by a Venn diagram.
Utilizing the clusterProfiler tool  in R (Version 3.2.11), the co-expressed DEGs were evaluated by a Kyoto Encyclopedia of Genes and Genomes (KEGG)  pathway enrichment analysis and Gene Ontology  Biological Process (BP) enrichment analysis. The significantly enriched terms with a p-value of less than 0.05 and involving no less than two DEGs were selected.
Screening CD4+ and CD8+ T cell-related DEGs
Based on the DEG data for the RNA-seq expression profiles, the abundance of infiltrating immune cells in tumor samples was estimated using the Cibersort algorithm . The infiltration of six types of immune cells (CD4+ T cells, CD8+ T cells, B cells, neutrophils, dendritic cells, and lymphocytes) in tumor tissues was evaluated.
The Pearson correlation coefficients for the relationships between the expression levels of DEGs and the abundance of infiltrated CD4+ and CD8+ T cells were calculated. The immune-related DEGs with an absolute value of r > 0.15 were selected.
PPI network construction
Interactions between proteins encoded by immune-related DEGs in CD4+ T cells were retrieved from the STRING database  (version 11.0) with a PPI score setting of 0.15 (low confidence) and the species set to human. Based on the retrieved PPIs, the CD4+ T cell-related PPI network was visualized using Cytoscape  (version 3.2.0). The CD8+ T cell-related PPI network was constructed following the same method.
ceRNA network construction
miRWalk 3.0  was used to predict target–miRNA regulatory relationships. The miRNAs correlated with DEGs related to CD4+ T cells were predicted, and the species was set to human. The regulatory relationships with a score exceeding 0.95 appearing in both the TargetScan and miRDB databases were selected. The miRNAs correlated with DEGs related to CD8+ T cells were predicted using the same method.
The lncRNA–miRNA relationships for CD4+ T cells were predicted using the DIANA-LncBase database v.2 . The lncRNA–miRNA regulatory relationships with a score of 1 were selected. The lncRNA–miRNA and target–miRNA data were integrated to construct the ceRNA network of the lncRNA–miRNA–target relationships related to CD4+ T cells. A CD8+ T cell-related ceRNA network was constructed following the same method.
Small chemical molecule–target network construction
The Comparative Toxicogenomics Database (CTD)  was used to search for colon cancer-related genes and small chemical molecules. Overlapping genes associated with colon cancer and belonging to the CD4+T cell-related ceRNA network were used to screen chemical–target pairs. A CD4+ T cell-related small chemical molecule–target network was constructed using Cytoscape. A CD8+ T cell-related small chemical molecule–target network was also constructed using the same method.
Clinical phenotype data related to prognosis in TCGA were collected, including overall survival (OS). The CD4+ T cell-related genes and were divided into high/low expression groups based on the median gene expression value and the same grouping strategy was applied to CD8+ T cell-related genes. Differences were evaluated by the log-rank test and genes with a p-value of less than 0.05 were considered significantly correlated with prognosis. Furthermore, Kaplan–Meier (K-M) survival curves were plotted.
Differences in gene expression between colon cancer with high and low stromal and immune scores
There were 83 DEGs (5 up-regulated and 78 down-regulated) between the groups with high and low stromal scores. A total of 1270 DEGs (including 807 up-regulated and 463 down-regulated genes) were detected between the groups with high and low immune scores. A volcano map of DEGs is shown in Fig. 2. Furthermore, 5 up-regulated and 74 down-regulated DEGs were identified in both the stromal and immune groups, as shown in a Venn diagram in Fig. 2c and in Supplementary Table 1.
Enrichment analysis of the overlapping DEGs
In an enrichment analysis of co-expressed DEGs, up-regulated DEGs were enriched for 22 BP terms, including response to copper ions and response to electrical stimuli (Fig. 3a). The down-regulated DEGs were enriched for 17 BP terms, including immune response to tumor cells and stabilization of the membrane potential (Fig. 3b). The co-expressed DEGs were not enriched for any of the investigated KEGG pathways.
CD4+ and CD8+ T cell-related DEGs
Based on the Cibersort algorithm, the abundances of infiltrating immune cells (CD4+ T cells, CD8+ T cells, B cells, neutrophils, dendritic cells, and lymphocytes) in the tumor samples were estimated (Fig. 4). Based on Pearson correlation coefficients, 79 CD4+ T cell-related DEGs and 8 CD8+ T cell-related DEGs were identified.
PPI networks of CD4+ and CD8+ T cell-related DEGs
We generated a PPI network with 59 nodes and 77 interacting pairs for CD4+ T cell-related genes, including the ETS transcription factor ELK3, as shown in Fig. 5. There were no PPIs involving DEGs related to CD8+ T cells under the threshold used in this study.
ceRNA network of CD4+ and CD8+ T cell-related DEGs
We obtained 69 miRNA–mRNA relationships (55 miRNAs and 19 target genes) in an analysis of CD4+ T cells (Fig. 6a) and only two miRNA–mRNA relationships in an analysis CD8+ T cells, including two miRNAs and one target gene (Fig. 6b).
In total, 11 lncRNAs were predicted to interact with miRNAs in CD4+ T cells. After integrating these data and data for miRNA–target relationships, 110 ceRNA regulatory relationships were obtained for CD4+ T cells, such as chr22-38_28785274–29,006,793.1–miR-106a-5p–DDHD1, which included 45 miRNAs and 11 mRNAs (Fig. 7a).
Four predicted lncRNAs were obtained for CD8+ T cells. After integration with the miRNA–target relationships, five ceRNA regulatory relationships were obtained for CD8+ T cells, such as chr22-38_28785274–29,006,793.1–miR-4319–GRHL1, which included one miRNA and one mRNA (Fig. 7b).
Colon cancer-related small chemical molecules target mRNAs in the ceRNA network
For CD4+ T cells, 175 small chemical molecule–target interactions were obtained, including 2 lncRNAs, 10 genes, and 64 types of small chemical molecules (Fig. 8a). For CD8+ T cells, nine small chemical molecule–target interactions were obtained, including one gene and nine types of small chemical molecules (Fig. 8b).
CD4+ and CD8+ T cell-related genes associated with colon cancer prognosis.
After screening for DEGs associated with survival, adenosine deaminase domain containing 1 (ADAD1) and discs large MAGUK scaffold protein 3 (DLG3) were selected among the 77 CD4+ T cell-related DEGs. However, there were no survival-related DEGs among the eight CD8+ T cell-related genes. The K-M survival curve for ADAD1 is represented in Fig. 9a and the K-M survival curve for DLG3 is represented in Fig. 9b.
ELK3 is a candidate immune-related gene in colon cancer
In this study, by comparing samples from TCGA with high or low stromal and immune scores, we obtained 83 DEGs (5 up-regulated and 78 down-regulated genes) between the two stromal score groups and 1270 DEGs (807 up-regulated and 293 down-regulated genes) between the two immune score groups. In total, 79 co-expressed DEGs were obtained (i.e., genes associated with both stromal and immune scores). The 5 up-regulated co-expressed DEGs were enriched for 22 BP terms (response to copper ions and response to electrical stimuli) and the 74 down-regulated co-expressed DEGs were enriched for 17 BP terms (immune response to tumor cells and stabilization of membrane potential). Moreover, 79 CD4+ T cell-related DEGs were detected, with 77 edges and 59 nodes in the PPI network, compared with only eight CD8+ T cell-related DEGs. Furthermore, there were 110 ceRNA relationships in the network of CD4+ T cell-related DEGs and five ceRNA relationships in the network of CD8+ T cells. Finally, to clarify the clinical applications of our findings as well as to provide a potential strategy for the prediction of clinical prognosis and targeted drug selection, prognostic genes were analyzed and small-molecule drugs were screened. In an analysis of small chemical molecule–gene interactions, 175 pairs were obtained for CD4+ T cells and 9 pairs for CD8+ T cells. In a survival analysis, ADAD1 and DLG3 in CD4+ T cells were most strongly correlated with prognosis. As such, analyses of the molecular mechanisms underlying the roles of CD4+ and CD8+ T cells in colon cancer may provide novel targets for immunotherapy.
ELK3 is a transcription factor belonging to the E26 transformation-specific (ETS) family . The PI3K/Akt/mammalian target of rapamycin (mTOR) and ERK signaling pathways activate ELK3 . Notably, studies have shown that ELK3 regulates cell migration and invasion in hepatoma cells and breast cancer . In an analysis of the regulation of colorectal cancer stemness, Wang et al. demonstrated that ELK3 is involved. Moreover, ELK3 was identified as a potential target of miR-507, and its expression is restored by the abrogation of LINC00525-knockdown . Many important biological processes are regulated by the ETS protein family, such as immune cell functions [29, 30]. ELK3 regulates the expression of heme oxygenase-1 (HO-1) as a transcriptional repressor. Moreover, inflammatory mediators tend to affect the expression of ELK3, which is down-regulated by bacterial endotoxins. Tsoyi et al. have shown that the ETS protein family plays a role in the immune response; during the inflammatory response to infection, ELK3 and HO-1 are important for macrophage function . In our study of CD4+ T cell-related DEGs, ELK3 was downregulated.
The chr22-38_28785274–29,006,793.1–miR-106a-5p-DDHD1 and chr22-38_28785274–29,006,793.1–miR-4319-GRHL1 axes may be related to CD4+ and CD8+ T cell infiltration in colon cancer
We identified the chr22-38_28785274–29,006,793.1–miR-106a-5p–DDHD1 axis from the CD4+ T-cell related ceRNA network as a candidate for further analyses. The phospholipase A1 (PLA1) family members are classified as extracellular and intracellular and are implicated in different intracellular mechanisms. As a phosphatidic acid (PA)-preferring PLA1 (PA-PLA1), intracellular DDHD1 has been studied extensively owing to its implications for cancer development. DDHD1 is involved in the synthesis of lysophosphatidylinositol (LPI) . LPI activity is correlated with tumor growth and aggressiveness via its interaction with G protein-coupled receptor 55 (GPR-55) [33,34,35,36]. Moreover, DDHD1 supports the proliferation and survival of colon cancer cells. Studies have also demonstrated that the inhibition of the MAPK/ERK and PI3K/Akt signaling pathways reduce the viability of colon cancer cells in vitro, and apoptotic cell death is increased by the silencing of DDHD1 via small interfering RNA . Our results were consistent with those of previous studies, which supported the effects of the lysophospholipid mediator DDHD1 on tumor processes. In neoplastic cells, by interacting with GPR-55, LPI induces ERK and the Akt signaling . MiR-106a-5p belongs to the miR-17 family. According to the consensus seed region, there are three clusters in the miR-17 family. MiR-106a-5p is located on Xq26.2, which belongs to the miR-106a-363 cluster. MiR-106a-5p is highly expressed in gastric [38,39,40,41,42], breast [43, 44], colorectal , and non-small cell lung cancers . In squamous cell carcinomas , colon cancers , and gliomas , miR-106a-5p is expressed at relatively low levels. Studies have also shown that in colorectal cancer, the inhibition of the anti-metastatic gene transforming growth factor-b receptor 2 (TGFBR2) increases cell migration and invasion via miR-106a-5p . In our study, DDHD1 was identified as a potential target of miR-106a-5p in colon cancer cells, which influences disease progression.
In our study, the chr22-38_28785274–29,006,793.1–miR-4319–GRHL1 axis from the CD8+ T-cell related ceRNA network was identified as a candidate for further analyses. Grainyhead-like 1 (GRHL1) belongs to the GRHL transcription factor family, which comprises GRHL1, GRHL2, and GRHL3 . Studies have suggested that the Grainyhead family genes exhibit homologous DNA-binding immunoglobulin folding to the DNA-binding domain of the key tumor suppressor p53 and that these genes participate in wound healing, embryonic neural tube closure, and epidermal integrity [51,52,53]. Recent studies have shown that these transcription factors are involved in various cancers, such as skin squamous cell carcinoma, gastric cancer, breast cancer, colorectal cancer, and cervical cancer . Moreover, GRHL2 knockdown in colorectal cancer cells inhibits cell proliferation by targeting ZEB1 . Huang et al. revealed that the expression of miR-4319 is inversely related to patient survival in colorectal cancer. Moreover, the overexpression of miR-4319 markedly reduces colorectal cancer cell proliferation by infecting ankyrin repeat and BTB domain-containing 1 (ABTB1) and alters the cell cycle distribution . Thus, we hypothesized that the chr22-38_28785274–29,006,793.1–miR-4319–GRHL1 axis is correlated with CD8+ T cell functions and with the pathogenesis of colon cancer.
A series of bioinformatics analyses were conducted to identify and characterize DEGs related to CD4+ and CD8+ T cells in colon cancer tissues. ELK3, which was down-regulated in cancer tissues, may be correlated with colon cancer and CD4+ T cells. Moreover, the chr22-38_28785274–29,006,793.1–miR-106a-5p–DDHD1 axis from the CD4+ T cell-related ceRNA network as well as the chr22-38_28785274–29,006,793.1–miR-4319–GRHL1 axis from the CD8+ T cell-related ceRNA network were identified as candidates for further analyses. In a survival analysis of CD4+ T cell-related DEGs, ADAD1 and DLG3 were strongly correlated with prognosis. Furthermore, 175 small chemical molecule–gene interaction pairs in CD4+ T cells and 9 in CD8+ T cells were screened. Our identification of T cell-related RNAs, ceRNA network construction, and miRNA target prediction may provide a basis for further studies of colon cancer-induced T cell-mediated immune escape. The prediction of small chemical molecule drugs and survival differences based on differential RNA expression may provide a novel direction for clinical decision-making with respect to treatment and the evaluation of prognosis in colon cancer from the perspective of immunity.
Major histocompatibility antigen class II
Genomic Data Commons
The Cancer Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data
Comparative Toxicogenomics Database
Adenosine deaminase domain containing 1
Discs Large MAGUK scaffold protein 3
mammalian target of rapamycin
G protein-coupled receptor 5
Transforming growth factor-b receptor 2
Arnold M, et al. Global patterns and trends in colorectal cancer incidence and mortality. Gut. 2017;66(4):683–91.
van der Stok EP, et al. Surveillance after curative treatment for colorectal cancer. Nat Rev Clin Oncol. 2017;14(5):297.
Vasan N, Baselga J, Hyman DM. A view on drug resistance in cancer. Nature. 2019;575(7782):299–309.
Zaidi N, Jaffee EM. Immunotherapy transforms cancer treatment. J Clin Invest. 2019;129(1):46–7.
Canning C, et al. Liver immunity and tumour surveillance. Immunol Lett. 2006;107(2):83–8.
Vilgelm AE, Richmond A. Chemokines modulate immune surveillance in tumorigenesis, Metastasis and Response to Immunotherapy. Front Immunol. 2019;10:333.
Jiang X, et al. Role of the tumor microenvironment in PD-L1/PD-1-mediated tumor immune escape. Mol Cancer. 2019;18(1):10.
De Matteis S, et al. Advances in molecular mechanisms and immunotherapy involving the immune cell-promoted epithelial-to-Mesenchymal transition in lung Cancer. J Oncol. 2019;2019.
Angell H, Galon J. From the immune contexture to the Immunoscore: the role of prognostic and predictive immune markers in cancer. Curr Opin Immunol. 2013;25(2):261–7.
Ribas A, Wolchok JD. Cancer immunotherapy using checkpoint blockade. Science. 2018;359(6382):1350–5.
Dmitry G. Mechanisms and functional significance of tumour-induced dendritic-cell defects. Nat Rev Immunol. 2004;4(12):941–52.
Liu YM, Luo J. Recognition of CD4~+ T cells to MHC-II molecule and its biological significance. Immunological Journal; 2006. p. S1.
Zhang S, Zhang H, Zhao J. The role of CD4 T cell help for CD8 CTL activation. Biochem Biophys Res Commun. 2009;384(4):405–8.
Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
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.
Yu G, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Ashburner M, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25.
Newman AM, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453.
Szklarczyk D, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2014;43(D1):D447–52.
Shannon P, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Dweep H, Gretz N. miRWalk2. 0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.
Paraskevopoulou MD, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2015;44(D1):D231–8.
Davis AP, et al. The comparative toxicogenomics database: update 2019. Nucleic Acids Res. 2018;47(D1):D948–54.
Gross C, Dubois-Pot H, Wasylyk B. The ternary complex factor net/Elk-3 participates in the transcriptional response to hypoxia and regulates HIF-1α. Oncogene. 2008;27(9):1333.
Wasylyk C, et al. Inhibition of the Ras-net (Elk-3) pathway by a novel pyrazole that affects microtubules. Cancer Res. 2008;68(5):1275–83.
Lee JH, et al. ELK3 promotes the migration and invasion of liver cancer stem cells by targeting HIF-1α. Oncol Rep. 2017;37(2):813–22.
Wang S, Li J, Yang X. Long non-coding RNA LINC00525 promotes the Stemness and Chemoresistance of colorectal Cancer by targeting miR-507/ELK3 Axis. Int J Stem Cells. 2019;12(2):347.
Dittmer J, Nordheim A. Ets transcription factors and human disease. Biochim Biophys Acta. 1998;1377(2):F1–11.
Garrett-Sinha LA. Review of Ets1 structure, function, and roles in immunity. Cell Mol Life Sci. 2013;70(18):3375–90.
Tsoyi K, et al. Elk-3 is a KLF4-regulated gene that modulates the phagocytosis of bacteria by macrophages. J Leukoc Biol. 2015;97(1):171–80.
Higgs HN, et al. Cloning of a phosphatidic acid-preferring phospholipase A1 from bovine testis. J Biol Chem. 1998;273(10):5468–77.
Sutphen R, et al. Lysophospholipids are potential biomarkers of ovarian cancer. Cancer Epidemiol Prev Biomarkers. 2004;13(7):1185–91.
Pineiro R, Maffucci T, Falasca M. The putative cannabinoid receptor GPR55 defines a novel autocrine loop in cancer cell proliferation. Oncogene. 2011;30(2):142.
Andradas C, et al. The orphan G protein-coupled receptor GPR55 promotes cancer cell proliferation via ERK. Oncogene. 2011;30(2):245.
Oka S, et al. Identification of GPR55 as a lysophosphatidylinositol receptor. Biochem Biophys Res Commun. 2007;362(4):928–34.
Raimondo S, et al. The phospholipase DDHD1 as a new target in colorectal cancer therapy. J Exp Clin Cancer Res. 2018;37(1):82.
Kim BH, et al. Prognostic implications for high expression of oncogenic microRNAs in advanced gastric carcinoma. J Surg Oncol. 2013;107(5):505–10.
Wang Z, et al. miR-106a is frequently upregulated in gastric cancer and inhibits the extrinsic apoptotic pathway by targeting FAS. Mol Carcinog. 2013;52(8):634–46.
Yao Y, et al. MicroRNA profiling of human gastric cancer. Mol Med Rep. 2009;2(6):963–70.
Xiao B, et al. Detection of miR-106a in gastric carcinoma and its clinical significance. Clin Chim Acta. 2009;400(1–2):97–102.
Guo J, et al. Differential expression of microRNA species in human gastric cancer versus non-tumorous tissues. J Gastroenterol Hepatol. 2009;24(4):652–7.
Kim K, et al. Identification of oncogenic microRNA-17-92/ZBTB4/specificity protein axis in breast cancer. Oncogene. 2012;31(8):1034.
Wang F, et al. Correlation and quantitation of microRNA aberrant expression in tissues and sera from patients with breast tumor. Gynecol Oncol. 2010;119(3):586–93.
Feng B, et al. Colorectal cancer migration and invasion initiated by microRNA-106a. PLoS One. 2012;7(8):e43452.
Donnem T, et al. MicroRNA signatures in tumor tissue related to angiogenesis in non-small cell lung cancer. PLoS One. 2012;7(1):e29671.
Hummel R, et al. MiRNAs and their association with locoregional staging and survival following surgery for esophageal carcinoma. Ann Surg Oncol. 2011;18(1):253–60.
Díaz R, et al. Deregulated expression of miR-106a predicts survival in human colon cancer patients. Genes Chromosom Cancer. 2008;47(9):794–802.
Yang G, et al. MiR-106a inhibits glioma cell growth by targeting E2F1 independent of p53 status. J Mol Med. 2011;89(10):1037.
Quan Y, et al. Grainyhead-like 2 promotes tumor growth and is associated with poor prognosis in colorectal Cancer. J Cancer. 2015;6(4):342–50.
Kokoszyńska K, et al. The fold recognition of CP2 transcription factors gives new insights into the function and evolution of tumor suppressor protein p53. Cell Cycle. 2008;7(18):2907–15.
Narasimha M, et al. Grainy head promotes expression of septate junction proteins and influences epithelial morphogenesis. J Cell Sci. 2008;121(6):747–52.
Pyrgaki C, Liu A, Niswander L. Grainyhead-like 2 regulates neural tube closure and adhesion molecule expression during neural fold fusion. Dev Biol. 2011;353(1):38–49.
Mlacki M, et al. Recent discoveries concerning the involvement of transcription factors from the Grainyhead-like family in cancer. Exp Biol Med (Maywood). 2015;240(11):1396–401.
Quan Y, et al. Downregulation of GRHL2 inhibits the proliferation of colorectal cancer cells by targeting ZEB1. Cancer Biol Ther. 2014;15(7):878–87.
Huang L, et al. MiR-4319 suppresses colorectal cancer progression by targeting ABTB1. United European Gastroenterol J. 2019;7(4):517–28.
The authors gratefully acknowledge the multiple databases, which made the data available.
This work was supported by the Medical and Health Projects of Zhejiang Province (No. 2020KY301) and Public Welfare Technology Application Research Program of Huzhou (No. 2020ZDT2012). Both funders provided funding for the articles and proofread the literature.
Ethics approval and consent to participate
Consent for publication
The authors declare that no conflicts of interest exist.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yang, X., Wu, W., Pan, Y. et al. Immune-related genes in tumor-specific CD4+ and CD8+ T cells in colon cancer. BMC Cancer 20, 585 (2020). https://doi.org/10.1186/s12885-020-07075-x
- Colon cancer
- Competing endogenous RNAs
- CD4+ T cells
- CD8+ T cells