Transcriptome of Tumor-Infiltrating T Cells in Colorectal Cancer Patients Uncovered a Unique Gene Signature in CD4+ T Cells Associated with Poor Disease-Specific Survival.

Colorectal cancer (CRC) is influenced by infiltration of immune cell populations in the tumor microenvironment. While elevated levels of cytotoxic T cells are associated with improved prognosis, limited studies have reported associations between CD4+ T cells and disease outcomes. We recently performed transcriptomic profiling and comparative analyses of sorted CD4+ and CD8+ tumor-infiltrating lymphocytes (TILs) from bulk tumors of CRC patients with varying disease stages. In this study, we compared the transcriptomes of CD4+ with CD8+ TILs. Functional annotation pathway analyses revealed the downregulation of inflammatory response-related genes, while T cell activation and angiogenesis-related genes were upregulated in CD4+ TILs. The top 200 deregulated genes in CD4+ TILs were aligned with the cancer genome atlas (TCGA) CRC dataset to identify a unique gene signature associated with poor prognosis. Moreover, 69 upregulated and 20 downregulated genes showed similar trends of up/downregulation in the TCGA dataset and were used to calculate "poor prognosis score" (ppScore), which was significantly associated with disease-specific survival. High ppScore patients showed lower expression of Treg-, Th1-, and Th17-related genes, and higher expression of Th2-related genes. Our data highlight the significance of T cells within the TME and identify a unique candidate prognostic gene signature for CD4+ TILs in CRC patients.

Epigenetic regulation of immune checkpoints and T cell exhaustion markers in tumor-infiltrating T cells of colorectal cancer patients.

Aim: To elucidate the epigenetic alterations behind the upregulation of immune checkpoints and T cell exhaustion markers in colorectal cancer (CRC) patients. Materials & methods: mRNA expressions of different immune checkpoint/exhaustion markers were analyzed by quantitative real-time reverse transcriptase PCR and epigenetic investigations were performed using bisulfite sequencing and chromatin immunoprecipitation quantitative PCR. Results: mRNA expressions of PD-1, TIM-3, CTLA-4, PD-L1 and TOX2 were significantly upregulated in CD4+ and CD8+ tumor-infiltrating lymphocytes and bulk CRC tumor tissues. Histone 3 lysine 9 trimethylation was downregulated and histone 3 lysine 4 trimethylation was upregulated in PD-L1 and TOX2 promoters in tumor tissues, suggesting that PD-L1 and TOX2 upregulation in CRC tumors could be mediated by activating histone 3 lysine 4 trimethylation. Conclusion: Epigenetic modifications in promoters of immune checkpoint and T cell exhaustion genes could induce their upregulation, and potentially implicate the use of epigenetic modifiers to enhance antitumor immunity in CRC patients.

Differential gene expression of tumor-infiltrating CD4+ T cells in advanced versus early stage colorectal cancer and identification of a gene signature of poor prognosis.

Tumor-infiltrating lymphocytes (TILs) play indispensable roles in the progression and response to treatment of solid tumors. However, the prognostic significance of CD4+ TILs is not fully disclosed in cancers generally and in CRC in particular, mainly due to the existence of different functional subsets of CD4+ T cells. We performed transcriptomic profiling of CD4+ TILs isolated from CRC patients in order to identify differentially expressed genes and their functional pathways in early versus advanced disease stages. We found that in advanced stages, genes related to immune and inflammatory responses, in particular Th1-mediated immune response and cytotoxicity-mediated genes, were downregulated; while epigenetic-mediated silencing genes were upregulated. Interestingly, we identified genes, which were steadily upregulated or downregulated in CD4+ TILs with CRC progression from stage I to IV. Additionally, of the top 200 deregulated genes, 43 upregulated and 64 downregulated genes showed similar deregulation trends in the cancer genome atlas CRC dataset. From these 97 deregulated genes, we identified a "poor prognosis CD4 gene signature (ppCD4sig)". Patients with high ppCD4sig score showed shorter disease-specific survival (DSS) and progression-free interval (PFI). The ppCD4sig was an independent prognostic indicator for DSS (HR = 1.73, 95% CI 1.32-2.27, P = 0.0001) and PFI (HR = 1.75, 95% CI 1.3-2.35, P = 0.0016). Additionally, patients at advanced stages and at a younger age (<55 years) were more likely to have a high ppCD4sig score. Altogether, our data provide novel insights and a unique prognostic gene signature of CD4+ TILs in the CRC microenvironment.

Differential expression of TIM-3 in circulation and tumor microenvironment of colorectal cancer patients.

T-cell immunoglobulin and mucin-domain containing-3 (TIM-3) is an inhibitory immune checkpoint, which suppresses anti-tumor immune responses. TIM-3 expression on different immune cells in periphery and tumor microenvironment (TME) of colorectal cancer (CRC) patients has not been fully investigated. We found that TIM-3 was mainly expressed on monocytic myeloid cells (MMCs) and antigen-presenting cells (APCs) in circulation but was mainly expressed on T cells and APCs in the TME. Additionally, TIM-3- T cells co-expressed higher levels of PD-1 than TIM-3+ T cells in normal tissue. In contrast, TIM-3+ T cells in the TME showed significantly higher PD-1 expression. Interestingly, there was a trend towards increased levels of TIM-3+ APCs with disease stages; however, levels of TIM-3+ T cells were decreased with disease stages in the TME. This study shows the differential expression of TIM-3 on different immune cells in circulation and TME of CRC patients, and their associations with disease stages.

Transcriptomic Profiling of Tumor-Infiltrating CD4+TIM-3+ T Cells Reveals Their Suppressive, Exhausted, and Metastatic Characteristics in Colorectal Cancer Patients.

T cell immunoglobulin mucin-3 (TIM-3) is an immune checkpoint identified as one of the key players in regulating T-cell responses. Studies have shown that TIM-3 is upregulated in the tumor microenvironment (TME). However, the precise role of TIM-3 in colorectal cancer (CRC) TME is yet to be elucidated. We performed phenotypic and molecular characterization of TIM-3+ T cells in the TME and circulation of CRC patients by analyzing tumor tissues (TT, TILs), normal tissues (NT, NILs), and peripheral blood mononuclear cells (PBMC). TIM-3 was upregulated on both CD4+ and CD3+CD4- (CD8+) TILs. CD4+TIM-3+ TILs expressed higher levels of T regulatory cell (Tregs)-signature genes, including FoxP3 and Helios, compared with their TIM-3- counterparts. Transcriptomic and ingenuity pathway analyses showed that TIM-3 potentially activates inflammatory and tumor metastatic pathways. Moreover, NF-κB-mediated transcription factors were upregulated in CD4+TIM-3+ TILs, which could favor proliferation/invasion and induce inflammatory and T-cell exhaustion pathways. In addition, we found that CD4+TIM-3+ TILs potentially support tumor invasion and metastasis, compared with conventional CD4+CD25+ Tregs in the CRC TME. However, functional studies are warranted to support these findings. In conclusion, this study discloses some of the functional pathways of TIM-3+ TILs, which could improve their targeting in more specific therapeutic approaches in CRC patients.

Transcriptomic profiling disclosed the role of DNA methylation and histone modifications in tumor-infiltrating myeloid-derived suppressor cell subsets in colorectal cancer.

BACKGROUND: Increased numbers of myeloid-derived suppressor cells (MDSCs) are positively correlated with poor prognosis and reduced survivals of cancer patients. They play central roles in tumor immune evasion and tumor metastasis. However, limited data are available on phenotypic/transcriptomic characteristics of the different MDSCs subsets in cancer. These cells include immature (I-MDSCs), monocytic (M-MDSCs), and polymorphonuclear/granulocytic (PMN-MDSCs). METHODS: Phenotypic characterization of myeloid subsets from 27 colorectal cancer (CRC) patients was assessed by flow cytometric analyses. RNA-sequencing of sorted I-MDSCs, PMN-MDSCs, and antigen-presenting cells (APCs) was also performed. RESULTS: We found that the levels of I-MDSCs and PMN-MDSCs were increased in tumor tissues (TT), compared with normal tissues (NT) in colorectal cancer. Our functional annotation analyses showed that genes associated with histone deacetylase (HDAC) activation- and DNA methylation-mediated transcriptional silencing were upregulated, and histone acetyl transferase (HAT)-related genes were downregulated in tumor-infiltrating I-MDSCs. Moreover, pathways implicated in cell trafficking and immune suppression, including Wnt, interleukin-6 (IL-6), and mitogen-activated protein kinase (MAPK) signaling, were upregulated in I-MDSCs. Notably, PMN-MDSCs showed downregulation in genes related to DNA methylation and HDAC binding. Using an ex vivo model, we found that inhibition of HDAC activation or neutralization of IL-6 in CRC tumor tissues downregulates the expression of genes associated with immunosuppression and myeloid cell chemotaxis, confirming the importance of HDAC activation and IL-6 signaling pathway in MDSC function and chemotaxis. CONCLUSIONS: This study provides novel insights into the epigenetic regulations and other molecular pathways in different myeloid cell subsets within the CRC tumor microenvironment (TME), giving opportunities to potential targets for therapeutic benefits.

Review of the Integrated Groundwater and Surface-Water Model (IGSM).

Development of the finite-element-based Integrated Groundwater and Surface-Water Model (IGSM) began in the 1970s. Its popularity grew in the early 1990s with its application to California's Central Valley Groundwater Surface-Water Model in support of the Central Valley Project Improvement Act. Since that time, IGSM has been applied by federal, state, and local agencies to model a number of major basins in California. Our review of the recently released version 5.0 of IGSM reveals a solution methodology that deviates from established solution techniques, potentially compromising its reliability under many circumstances. One difficulty occurs because of the semi-explicit time discretization used. Combined with the fixed monthly time step of IGSM, this approach can prevent applications from accurately converging when using parameter values typically found in nature. Additionally, IGSM fails to properly couple and simultaneously solve ground water and surface water models with appropriate mass balance and head convergence under the reasonable conditions considered herein. As a result, IGSM-predicted streamflow is error prone, and errors could exceed 100%. IGSM does not inform the user that there may be a convergence problem with the solution, but instead generally reports good mass balance. Although our review touches on only a few aspects of the code, which exceeds 17,000 lines, our experience is that similar problems arise in other parts of IGSM. Review and examples demonstrate the potential consequences of using the solution methods in IGSM for the prediction, planning, and management of water resources, and provide perspective on the roles of standards and code validation in ground water modeling.