- 1General Surgery Department, The First Affiliated Hospital of Nanchang University, Nanchang, China
- 2Obstetrics and Gynaecology Department, The First Affiliated Hospital of Nanchang University, Nanchang, China
- 3Queen Mary College, Nanchang University, Nanchang, China
- 4Obstetrics and Gynecology Department, The First Affiliated Hospital of Nanchang University, Nanchang, China
Tumor-infiltrating immune cells, associated with tumor progression, are promising prognostic biomarkers. However, the relationship between levels of gene expression and that of immune cell infiltration in cervical cancer prognosis is unknown. In this study, three cervical cancer gene expression microarrays (GSE6791, GSE63678 and GSE55940) were obtained from the GEO database. The IDO1 gene was identified by differentially expressed gene screening. The gene expression profiles of TCGA and GTEx databases along with comprehensive bioinformatics analysis identified that the IDO1 gene was upregulated in cervical cancer with significant difference in expression at different N stages. In addition, it was also upregulated in HPV16 positive sample. The pan-cancer analysis identified that IDO1 was highly expressed in most cancers. TIMER analysis revealed that the expression of IDO1 in CESC shows positive correlation with CD8+ T cells, CD4+ T cells, neutrophils, dendritic cells. IDO1 expression showed remarkable positive correlation with all immune cell markers except M1 macrophages. CD8+ T cell infiltration GSEA results showed that IDO1 was mainly associated with tumor immune-related signaling pathways.
Introduction
Cervical cancer is one of the most prevalent malignant diseases, affecting women worldwide (1). Persistent infection with high-risk human papillomaviruses (mainly HPV16 and HPV18) is the main cause of cervical cancer and its precursor lesions (2, 3). The incidence of cervical cancer has decreased due to the production of HPV vaccines, improved living conditions, and early screening (4). However, there is still a high mortality of cancer in Asia (5). Therefore, necessary measures are needed to reduce the incidence of cervical cancer. It is well known that the HPV-E6 and E7 genes encode two mucoproteins involved in the pathogenesis of cervical cancer. Cervical cancer is therefore a relatively immunogenic cancer that can use various mechanisms to evade immune attack by the host (6).
The interaction of tumor cells with the microenvironment plays a crucial role in the development of malignant tumors (7). The main component of the tumor microenvironment is the tumor-infiltrating lymphocytes, a specific population of T cells with a high specific immune response to tumor cells (8, 9). The immune microenvironment of cervical cancer not only has common features of solid tumors but also unique characteristics associated with HPV infection (10). The importance of infiltrating lymphocytes in predicting progression of different types of solid tumors has been demonstrated and T cells are the immune cells of choice for treating cancer. Tumor cell growth is attributed to develop immune resistance by evading immune surveillance mechanisms, leading to T cell dysfunction and depletion (11). Overcoming T-cell dysfunction in cancer patients is the focus of oncology treatment. Blocking antibodies against cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) and PD-1/PD-L1 immune checkpoints have shown durable clinical responses in a variety of cancers, including cervical cancer (12–14). However, many cancer patients show primary or secondary treatment tolerance to immune checkpoint therapy (15). Therefore, the search for more immune checkpoint-related genes to reverse T-cell dysfunction is essential. Indoleamine 2, 3-dioxygenase 1 (IDO1) is characterized by a rate-limiting metabolic enzyme that converts tryptophan (Trp) to downstream Kyn (Kyn) (16, 17). IDO1 is interferon-induced and has been shown to mediate powerful immunosuppression in cancer. Growing evidence indicates that IDO1 is overexpressed in the vast majority of solid tumors and associated with clinical prognosis, such as anal (18), esophageal (19), as well as cervical cancer (20). More importantly, several molecular drugs targeting IDO1 have been evaluated in multiple clinical trials with encouraging results (21, 22).
In this study, three cervical cancer gene expression microarrays were screened for differentially expressed genes and their expression was analyzed in Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC). The significance of the genes in the prognosis of CESC was determined by correlating the level of gene expression with that of immune cell infiltration and immune checkpoint-related gene expression. The Gene Set Enrichment Analysis (GSEA) identified the immune-related signaling pathways in the tumor.
Methods
Data Sources
Gene expression data microarrays for cervical cancer were obtained from the Gene Expression Omnibus (GEO) database: GSE6791 (23) (containing 84 samples), GSE55940 (24) (containing 10 samples) and GSE63678 (25) (containing 35 samples).
The human cancer gene expression data and related clinical information were obtained from the The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database. There were 10201 tumor samples in 33 tumors from TCGA database and 16871 normal samples from TCGA+GTEx database.
Differentially Expressed Gene
The”affy” and the “impute” R packages in R/Bioconductor software were used for GEO data processing, and then limma package of R software (version 3.40.2) was used for the analysis of differential gene expression. To correct false-positive results, adjusted P values were analyzed in TCGA or GTEx. Differential gene screening was performed based on adjusted P < 0.05, |logFC| > 1 and volcano plotting was performed using the ggplot2 R package. Overlapping genes were searched by VennDiagram package.
Immunocorrelation Analysis
The TIMER database (http://timer.comp-genomics.org/) was used to analyze the relation of IDO1 expression in CESC with the level of immune cell infiltration and immune cell markers. P<0.05 was considered statistically significant. Immune scores were assessed using the CIBERSORT algorithm. CIBERSORT is a versatile computational method for quantifying cell fractions from bulk tissue gene expression profiles (26). Spearman correlation test was used for the analysis of the correlation between gene and immune checkpoint associated gene. Visualization was achieved by the R (v4.0.3) package ggplot2.
Survival Analysis
Gene expression and immune cell infiltration levels were analyzed by the TIMER database in relation to overall survival (OS). Kaplan-Meier test was performed to analyze the difference between the survival of patients with high and low gene expression or high and low infiltration levels.
Protein Expression Validation
Immunohistochemical staining maps of CD8+ T cells markers for protein expression in both cervical cancer tissue and normal tissue were downloaded from The Human Protein Atlas (HPA) database.
Gene Set Enrichment Analysis
The samples were divided into two groups of high and low expression according to the median value of gene expression level, the group cutoff was set to median (Cutoff‐High and Cutoff‐Low are both 50%). GSEA was performed to investigate the functions correlated with different risk groups by GSEA 4.1.0, and the software was downloaded from the website for GSEA (http://www.gseamsigdb.org/gsea/downloads). The screening conditions were |NES|≥ 1, FDR < 0.25 and P < 0.05.
Results
Screening of Differentially Expressed Genes
Analytical screening of the GSE6791, GSE55940, and GSE63678 datasets yielded 1698, 14, and 540 differentially up-regulated genes and 236, 9, and 583 differentially down-regulated genes respectively (Figures 1A–C). Overlapping gene analysis of the differentially expressed gene identified a key gene IDO1 (Figure 1D).
Figure 1 Screening of differentially expressed genes. [(A) Volcano plot of DEGs in GSE6791; (B) Volcano plot of DEGs in GSE55940; (C) Volcano plot of DEGs in GSE63678; (D) Venn diagram showing overlapping genes of DEGs in the three datasets].
Expression of IDO1 in CESC and Multiple Cancer
IDO1 expression in CESC and paraneoplastic tissues was analyzed through TCGA database (306 tumor samples and 3 adjacent tumor samples) using Wilcox-tests in R software. No significant difference was observed in IDO1 expression (Figure 2A), probably because of the small number of adjacent tissue samples. IDO1 expression data was compared by integrating GTEx normal tissues which counts for 22 normal tissue. It was found that IDO1 expression was significantly upregulated in CESC compared to normal tissues (Figure 2B). The samples with unknown staging information were eliminated, there were 66 Nx samples, 134 N0 samples and 61 N1 samples from TCGA database. Analysis of IDO1 expression at different N stages revealed that IDO1 expression was highest in the N0 stage (Figure 2C). There were12 HPV-16 positive samples and 294 HPV-16 negative samples from TCGA database. IDO1 expression was remarkably higher in HPV16 positive than HPV16 negative samples (Figure 2D). Subsequently, tumor samples in TCGA were integrated with normal samples in GTEx (10201 tumor samples and 16871 normal samples from TCGA+GTEx database). The pan-cancer expression profile of IDO1 was analyzed using rank sum test. It was observed that IOD1 was remarkably down-regulated in Thyroid carcinoma (THCA) and up-regulated in the remaining cancers except for Acute Myeloid Leukemia (LAML), Lung squamous cell carcinoma (LUSC), Mesothelioma (MESO), Sarcoma (SARC), Thymoma (THYM), and Uveal Melanoma (UVM) (Figure 2E).
Figure 2 Expression level of IDO1. [(A) IDO1 expression level in cancerous and paraneoplastic tissues of CESC from TCGA datasets; (B) IDO1 expression level in cancerous and normal tissues of CESC from TCGA+GTEx datasets; (C) IDO1 expression level in different N stages of CESC; (D) IDO1 expression level in CESC with or without HPV type 16 infection; ns: P > 0.05, *P < 0.05, ***P < 0.001].
Immune Cell Infiltration Correlation
IDO1 expression levels in CESC were negatively correlated with tumor purity and macrophages but positively correlated with CD8+ and CD4+ T cells, neutrophils, and dendritic cells (Figure 3). The relationship between IDO1 expression and different immune cell type markers was further analyzed. By purity-adjusted correlation, it was identified that IDO1 expression levels show significant positive correlation with all immune cell markers except M1 macrophages (Table 1). The results showed that although IDO1 expression level had no correlation with B cell infiltration it had significant correlation with the type of B cell markers. This suggests that the correlation between IDO1 expression level and B cell infiltration could be verified by the expression of IDO1 and its type markers.
Immunoscore
The relative proportion of immune cell subsets in CESC samples with different IDO1 expression and different N stages in the TCGA dataset were counted using CIBERSORT. According to the median value of IDO1 expression, samples were divided into high and low expression groups. In Figure 4A, it can be seen that CD8+ T cells and M0 macrophages account for the majority of the 22 tumor-infiltrating immune cell subsets. Figure 4B shows that CD8+ T cells account for the highest proportion, revealing that the level of CD8+ T cell infiltration was significantly correlated with IDO1 expression and different N stages of CESC.
Figure 4 CIBERSORT statistics of relative proportions of immune cell subpopulations. [(A) Distribution of immune scores for different IDO1 expression level; (B) Distribution of immune scores for different N stages; “-”P > 0.05, *P < 0.05, **P < 0.01, ***P < 0.001].
Relationship Between CD8+ T-Cell Infiltration and Prognosis of CESC
Based on the median ratio of CD8+ T-cell infiltration and IDO1 expression, we divided the patient population into two groups (low and high ratios). By analyzing different CD8+ T-cell infiltration levels and CESC prognosis, it was found that patients with high CD8+ T-cell infiltration levels had a better prognosis as verified by CIBERSORT (Figures 5A, B), XCELL and QUNATISEQ (Supplementary Figure 1); low CD8+ effector memory T-cell infiltration levels were significantly associated with poor prognosis of patients. Further, the effects of both, the IDO1 expression and CD8+ T cell infiltration together were compared to patient prognosis (Figure 5C). It revealed that patients with high CD8+ T cell infiltration levels had a better prognosis in the case of high IDO1 expression. The prognosis of patients with high CD8+ T-cell infiltration levels was better with low IDO1 expression. CD8+ T-cell infiltration level can not only be regulated by IDO1 expression levels but also by other biological mechanisms in tumor microenvironment. Thus, co-influence of IDO1 expression level and CD8+ T-cell infiltration level on the prognosis of patients might be different from their individual impact.
Figure 5 Relationship between CD8+ T-cell infiltration and prognosis of CESC. [(A) CD8+ T-cell infiltration level and CESC prognosis based on CIBERSORT algorithm; (B) CD8+ T-cell infiltration level and CESC prognosis based on CIBERSORT-ABS algorithm; (C) CD8+ T-cell infiltration level and IDO1 expression level and CESC prognosis based on CIBERSORT algorithm].
Expression of Immune Checkpoint-Related Genes
IDO1 expression was classified as high or low according to the median expression value. The expression of immune checkpoint-related genes in high and low expression IDO1 were observed. The results indicate that the IDO1 expression shows significant positive correlation with the expression of immune checkpoint-related genes (Figures 6A, B). The correlation of expression level between IDO1 and other immune checkpoint-related genes in cancers are displayed in Figure 6C. Analysis of immune checkpoint-related gene expression, in multiple cancers, revealed that IDO1 itself is an immune checkpoint gene.
Figure 6 Correlation between IDO1 expression level and immune checkpoint-associated gene expression. [(A) Immune checkpoint-associated gene expression level in high and low expression level of IDO1 samples; (B) Heat map of immune checkpoint-associated gene expression; (C) Immune checkpoint-associated gene expression level in multiple cancers; *P < 0.05, **P < 0.01, ***P < 0.001].
Validation of the Protein Expression of CD8+ T Cells Markers
The protein expression of IDO1 and CD8+ T cells markers in cervical cancer tissues and normal cervical tissues was verified using the HPA online database (Figure 7). The results showed that CD8A and CD8B were not detected in both cancer tissues and normal liver tissues. IDO1 was highly detected in cervical cancer tissue and low expression in normal tissue.
Gene Set Enrichment Analysis
IDO1 expression was classified as high or low according to the median expression value. Red bar in Figures indicated low expression level while blue bar was high expression level.
KEGG pathway analysis showed that: cytokine and cytokine receptor interaction channels, antigen presentation and natural killer cell-mediated cytotoxicity channels were three significantly different pathways (Figure 8A). HALLMARK pathway analysis showed that the interferon-gamma response, allograft rejection and complementation were the three most significant pathways (Figure 8B).
Figure 8 Gene enrichment analysis of IDO1 in cancers. [(A) Highly expressed IDO1 in KEGG enriched top ten pathway; (B) Highly expressed IDO1 in HALLMARK enriched top ten pathway].
(KEGG: chemokine signaling pathway, cell adhesion molecules, hematopoietic cell lineage, Leishmania protozoa infection, autoimmune thyroid disease and TOLL-like receptor signaling pathway
HALLMARK: inflammatory response, IL6/JAK/STAT3 signaling pathway, interferon-alpha response, IL2-STAT5, KRAS signaling, TNF-alpha/NF-kb and apoptosis)
Discussion
Tumor cells are malignantly transformed normal cells. There are many mutated or aberrantly expressed proteins that may serve as tumor antigens (27). By recognizing these tumor antigens, T cells can clear tumor cells. Thus, the T cell immune response against tumor antigens is a central mechanism of the body’s anti-tumor immunity. Studies indicate that the level of tumor-infiltrating lymphocytes is positively correlated with the clinical prognosis of patients in a variety of solid tumors, thus demonstrating the existence of tumor immune surveillance (28–30).
In this study, we identified the IDO1 gene, by screening three cervical cancer gene expression microarrays for DEGs. These DEGs were upregulated in cervical cancer with significant difference in expression at different N stages and HPV-16 infection. Studies have shown that IDO1 is overexpressed in many solid tumor tissues. A study reported that IDO1 is with higher mRNA transcription and protein expression level than in normal cervix, and also in comparison to other cancers. We confirmed this conclusion by pan-cancer analysis (31). IDO1 is a rate-limiting enzyme that converts tryptophan to kynurenine. It plays a role in immunosuppression by increasing tryptophan metabolism in the tumor microenvironment (32). IDO1 expression may be induced by the secretion of IFN-γ by CD8+ T cells in the tumor microenvironment. Elevated IDO1 expression produces a series of effects such, as inhibition of: T cell function, CD4+ T cell differentiation pathway to Treg, and antigen presentation (33, 34). Tumors can utilize multiple escape mechanisms to avoid immune recognition. For antigen presentation, firstly, the antigen must be taken up by the dendritic cells and presented to CD8+ T cells. Secondly, the antigen must have direct tumor presentation in order to be recognized and killed by activated CD8+ T cells (35). These two processes are where immune escape occurs. The immune escape mechanisms include, regulation of antigen expression and alteration of antigen processing and presentation mechanisms in tumor cells (36, 37).
Tumor-infiltrating immune cells are major contributors to the tumor immune response. Their levels predict treatment outcome and survival (38). We found that in CESC, IDO1 expression levels were significantly correlated with the levels of infiltrating B cells, macrophages, CD8+ T cells, CD4+ T cells, neutrophils, and dendritic cells. Further analysis revealed that IDO1 expression level shows remarkable positive correlation with all immune cell markers except M1 macrophages. CD8+ T cells are a major component of tumor-infiltrating lymphocytes. Therefore, regulation of CD8+ T cell responses has been a focus of immunotherapy for cancer (39). It has been shown that a large infiltration of CD8+ T cells is associated with good prognosis in some tumors (40–42). In the present study, we found the highest and most significant proportion of CD8+ T cells in samples with different IDO1 expression and at different N stages by the CIBERSORT algorithm. We further observed the correlation between the CD8+ T cell infiltration and patient prognosis in CESC by different algorithms. The results showed that high CD8+ T-cell infiltration was remarkably associated with good prognosis in CESC patients. Combined with the expression level of IDO1, the prognosis of patients with low IDO1 expression was better at a high CD8+ T cell infiltration level. From this perspective, we analyzed the relationship between IDO1 expression and immune checkpoint-related genes. The results indicated that the IDO1 expression was significantly correlated with immune checkpoint-related gene expression. This suggests that IDO1 is itself an immune checkpoint gene. The expression of IDO1 creates two effects on its own microenvironments. It decreases tryptophan(Trp) and produces a series of toxic kynurenine(Kyn) metabolites (43). The toxic Kyn metabolites directly suppress the effector T cell response by favoring differentiation of Tregs (44, 45). The immune checkpoint is a key molecule in T-cell dysfunction. Hence, we can reverse T-cell dysfunction by regulating the expression of immune checkpoint molecules. This means that IDO1, which is an immune checkpoint gene is an important target for cancer immunotherapy intervention (46) (47). first reported the results of the phase I clinical trial of Epacadostat, a selective IDO1 inhibitor, for the treatment of advanced malignant solid tumors. The results show that Epacadostat was well tolerated and effectively inhibited IDO1 activity and Kyn levels. In recent years, IDO1 inhibitors have been used as a promising immunomodulatory agent for patients with advanced cancer (48–50). In practice, the infiltrated CD8+/Treg ratio may be a more precise parameter for the prognosis. Jallad et al. found that the triple immune therapy was capable of significantly enhancing the natural killer cell counts as well as the CD3+CD4+/Treg and CD3+CD8+/Treg ratios possibly enhancing the anti-tumorigenic environment (44). The present study also provides support for the application of IDO1 inhibitors in cervical cancer. Finally, to verify the possible mechanism of action of IDO1, the results of KEGG and HALLMARK enrichment analysis showed that IDO1 was mainly enriched in cytokine and cytokine receptor interaction, antigen presentation and natural killer cell-mediated cytotoxicity, and interferon γ response channels. This suggests that IDO1 is a tumor immunity and tumor escape related gene in CESC and can be used as a new target for cervical cancer therapy.
Conclusion
In summary, Indolamine-2,3-dioxygenase-1 (IDO1) is a cytosolic enzyme that catalyzes the conversion of+ essential amino acid Trp to kynurenine (Kyn). IDO1 is overexpressed in more than 50% of tumors and its overexpression increases the relative concentration of Kyn compared to Trp. Hence Kyn/Trp ratio can be used as a prognostic marker to monitor cancer invasiveness and progression. Our study mainly found that IDO1 can be a biomarker for prognosis prediction in CESC and was closely associated with infiltrating CD8+ T cells and immune checkpoint genes. This study provides ideas for the application of IDO1 inhibitors in the treatment of CESC and explores the potential value to enhance anti-tumor immunity and immunotherapy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
All authors (SZ, JW, MC, DC, JX, QC) statement: (1) SZ and JW as co-first author contributed ueqially to this work. (2) QC as corresponding author designed the whole project and provided financial support.
Funding
Funding was received from the following: General program of NSFC, 81560247, PPAR γ Molecular mechanism of regulating mat2a gene transcription in the pathogenesis of endometriosis; Scientific and technological research project of Jiangxi Provincial Department of education, GJJ180006, molecular regulation mechanism of ATAD3A/WASF3/KiSS1 signaling pathway in the pathogenesis of cervical cancer; Science and technology plan of Jiangxi Provincial Health and Family Planning Commission, 20185124, molecular mechanism of JAK2/STAT3 signaling pathway regulating WASF3 gene transcription on cervical cancer.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.720447/full#supplementary-material
Supplementary Figure 1 | Relationship between CD8+ T-cell infiltration and prognosis of CESC via different algorithm. (A) CD8+ T-cell infiltration level and CESC prognosis based on xCELL algorithm; (B) CD8+ infiltration level and CESC prognosis based on QUANTISEQ algorithm; (C) CD8+ effector memory infiltration level and CESC prognosis based on xCELL algorithm. T-cell infiltration level and CESC prognosis.
Supplementary Table 1 | Genes in Venn diagram
References
1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2019. CA: Cancer J Clin (2019) 69(1):7–34. doi: 10.3322/caac.21551
2. Cohen PA, Jhingran A, Oaknin A, Denny L. Cervical Cancer. Lancet (2019) 393(10167):169–82. doi: 10.1016/S0140-6736(18)32470-X
3. Lu MM, Li S, Yang H, Morrisey EE. Foxp4: A Novel Member of the Foxp Subfamily of Winged-Helix Genes Co-Expressed With Foxp1 and Foxp2 in Pulmonary and Gut Tissues. Mech Dev (2002) 119:S197–202. doi: 10.1016/S0925-4773(03)00116-3
4. Bhatla N, Aoki D, Sharma DN, Sankaranarayanan R. Cancer of the Cervix Uteri. Int J Gynecol. Obstetrics (2018) 143:22–36. doi: 10.1002/ijgo.12611
5. Yang L, Yang Y, Meng M, Wang W, He S, Zhao Y, et al. Identification of Prognosis-Related Genes in the Cervical Cancer Immune Microenvironment. Gene (2021) 766:145119. doi: 10.1016/j.gene.2020.145119
6. Piersma SJ. Immunosuppressive Tumor Microenvironment in Cervical Cancer Patients. Cancer Microenviron (2011) 4(3):361–75. doi: 10.1007/s12307-011-0066-7
7. Liang Y, Lü W, Zhang X, Lü B. Tumor-Infiltrating CD8+ and FOXP3+ Lymphocytes Before and After Neoadjuvant Chemotherapy in Cervical Cancer. Diagn Pathol (2018) 13(1):1–8. doi: 10.1186/s13000-018-0770-4
8. Kitamura T, Qian B-Z, Pollard JW. Immune Cell Promotion of Metastasis. Nat Rev Immunol (2015) 15(2):73–86. doi: 10.1038/nri3789
9. Şenbabaoğlu Y, Gejman RS, Winer AG, Liu M, Van Allen EM, de Velasco G, et al. Tumor Immune Microenvironment Characterization in Clear Cell Renal Cell Carcinoma Identifies Prognostic and Immunotherapeutically Relevant Messenger RNA Signatures. Genome Biol (2016) 17(1):1–25. doi: 10.1186/s13059-016-1092-z
10. Sheu BC, Chang WC, Lin HH, Chow SN, Huang SC. Immune Concept of Human Papillomaviruses and Related Antigens in Local Cancer Milieu of Human Cervical Neoplasia. J Obstetrics Gynaecol. Res (2007) 33(2):103–13. doi: 10.1111/j.1447-0756.2007.00492.x
11. Farhood B, Najafi M, Mortezaee K. CD8+ Cytotoxic T Lymphocytes in Cancer Immunotherapy: A Review. J Cell Physiol (2019) 234(6):8509–21. doi: 10.1002/jcp.27782
12. Duranti S, Pietragalla A, Daniele G, Nero C, Ciccarone F, Scambia G, et al. Role of Immune Checkpoint Inhibitors in Cervical Cancer: From Preclinical to Clinical Data. Cancers (2021) 13(9):2089. doi: 10.3390/cancers13092089
13. Ott PA, Hodi FS, Robert C. CTLA-4 and PD-1/PD-L1 Blockade: New Immunotherapeutic Modalities With Durable Clinical Benefit in Melanoma Patients. Clin Cancer Res (2013) 19(19):5300–9. doi: 10.1158/1078-0432.CCR-13-0143
14. De Velasco G, Je Y, Bossé D, Awad MM, Ott PA, Moreira RB, et al. Comprehensive Meta-Analysis of Key Immune-Related Adverse Events From CTLA-4 and PD-1/PD-L1 Inhibitors in Cancer Patients. Cancer Immunol Res (2017) 5(4):312–8. doi: 10.1158/2326-6066.CIR-16-0237
15. Pitt JM, Vétizou M, Daillère R, Roberti MP, Yamazaki T, Routy B, et al. Resistance Mechanisms to Immune-Checkpoint Blockade in Cancer: Tumor-Intrinsic and-Extrinsic Factors. Immunity (2016) 44(6):1255–69. doi: 10.1016/j.immuni.2016.06.001
16. Munn DH, Mellor AL. IDO in the Tumor Microenvironment: Inflammation, Counter-Regulation, and Tolerance. Trends Immunol (2016) 37(3):193–207. doi: 10.1016/j.it.2016.01.002
17. Uyttenhove C, Pilotte L, Théate I, Stroobant V, Colau D, Parmentier N, et al. Evidence for a Tumoral Immune Resistance Mechanism Based on Tryptophan Degradation by Indoleamine 2, 3-Dioxygenase. Nat Med (2003) 9(10):1269–74. doi: 10.1038/nm934
18. Mitra D, Horick NK, Brackett DG, Mouw KW, Hornick JL, Ferrone S, et al. High IDO1 Expression Is Associated With Poor Outcome in Patients With Anal Cancer Treated With Definitive Chemoradiotherapy. Oncol. (2019) 24(6):e275. doi: 10.1634/theoncologist.2018-0794
19. Kiyozumi Y, Baba Y, Okadome K, Yagi T, Ishimoto T, Iwatsuki M, et al. IDO1 Expression Is Associated With Immune Tolerance and Poor Prognosis in Patients With Surgically Resected Esophageal Cancer. Ann Surg (2019) 269(6):1101–8. doi: 10.1097/SLA.0000000000002754
20. Heeren AM, van Dijk I, Berry DR, Khelil M, Ferns D, Kole J, et al. Indoleamine 2, 3-Dioxygenase Expression Pattern in the Tumor Microenvironment Predicts Clinical Outcome in Early Stage Cervical Cancer. Front Immunol (2018) 9:1598. doi: 10.3389/fimmu.2018.01598
21. Zhai L, Spranger S, Binder DC, Gritsina G, Lauing KL, Giles FJ, et al. Molecular Pathways: Targeting IDO1 and Other Tryptophan Dioxygenases for Cancer Immunotherapy. Clin Cancer Res (2015) 21(24):5427–33. doi: 10.1158/1078-0432.CCR-15-0420
22. Li F, Zhang R, Li S, Liu J. IDO1: An Important Immunotherapy Target in Cancer Treatment. Int Immunopharmacol (2017) 47:70–7. doi: 10.1016/j.intimp.2017.03.024
23. Pyeon D, Newton MA, Lambert PF, Den Boon JA, Sengupta S, Marsit CJ, et al. Fundamental Differences in Cell Cycle Deregulation in Human Papillomavirus–Positive and Human Papillomavirus–Negative Head/Neck and Cervical Cancers. Cancer Res (2007) 67(10):4605–19. doi: 10.1158/0008-5472.CAN-06-3619
24. Sun N-X, Ye C, Zhao Q, Zhang Q, Xu C, Wang S-B, et al. Long Noncoding RNA-EBIC Promotes Tumor Cell Invasion by Binding to EZH2 and Repressing E-Cadherin in Cervical Cancer. PloS One (2014) 9(7):e100340. doi: 10.1371/journal.pone.0100340
25. Pappa KI, Polyzos A, Jacob-Hirsch J, Amariglio N, Vlachos GD, Loutradis D, et al. Profiling of Discrete Gynecological Cancers Reveals Novel Transcriptional Modules and Common Features Shared by Other Cancer Types and Embryonic Stem Cells. PloS One (2015) 10(11):e0142229. doi: 10.1371/journal.pone.0142229
26. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust Enumeration of Cell Subsets From Tissue Expression Profiles. Nat Methods (2015) 12(5):453–7. doi: 10.1038/nmeth.3337
27. Liu XS, Mardis ER. Applications of Immunogenomics to Cancer. Cell (2017) 168(4):600–12. doi: 10.1016/j.cell.2017.01.014
28. Thomas NE, Busam KJ, From L, Kricker A, Armstrong BK, Anton-Culver H, et al. Tumor-Infiltrating Lymphocyte Grade in Primary Melanomas Is Independently Associated With Melanoma-Specific Survival in the Population-Based Genes, Environment and Melanoma Study. J Clin Oncol (2013) 31(33):4252. doi: 10.1200/JCO.2013.51.3002
29. Chung YR, Kim HJ, Jang MH, Park SY. Prognostic Value of Tumor Infiltrating Lymphocyte Subsets in Breast Cancer Depends on Hormone Receptor Status. Breast Cancer Res Treat (2017) 161(3):409–20. doi: 10.1007/s10549-016-4072-9
30. Azimi F, Scolyer RA, Rumcheva P, Moncrieff M, Murali R, McCarthy SW, et al. Tumor-Infiltrating Lymphocyte Grade Is an Independent Predictor of Sentinel Lymph Node Status and Survival in Patients With Cutaneous Melanoma. J Clin Oncol (2012) 30(21):2678–83. doi: 10.1200/JCO.2011.37.8539
31. Qin Y, Ekmekcioglu S, Forget M-A, Szekvolgyi L, Hwu P, Grimm EA, et al. Cervical Cancer Neoantigen Landscape and Immune Activity Is Associated With Human Papillomavirus Master Regulators. Front Immunol (2017) 8:689. doi: 10.3389/fimmu.2017.00689
32. Sailer V, Sailer U, Bawden EG, Zarbl R, Wiek C, Vogt TJ, et al. DNA Methylation of Indoleamine 2, 3-Dioxygenase 1 (IDO1) in Head and Neck Squamous Cell Carcinomas Correlates With IDO1 Expression, HPV Status, Patients’ Survival, Immune Cell Infiltrates, Mutational Load, and Interferon γ Signature. EBioMedicine (2019) 48:341–52. doi: 10.1016/j.ebiom.2019.09.038
33. Maliniemi P, Laukkanen K, Väkevä L, Dettmer K, Lipsanen T, Jeskanen L, et al. Biological and Clinical Significance of Tryptophan-Catabolizing Enzymes in Cutaneous T-Cell Lymphomas. Oncoimmunology (2017) 6(3):e1273310. doi: 10.1080/2162402X.2016.1273310
34. Yentz S, Smith D. Indoleamine 2, 3-Dioxygenase (IDO) Inhibition as a Strategy to Augment Cancer Immunotherapy. BioDrugs (2018) 32(4):311–7. doi: 10.1007/s40259-018-0291-4
35. Jhunjhunwala S, Hammer C, Delamarre L. Antigen Presentation in Cancer: Insights Into Tumour Immunogenicity and Immune Evasion. Nat Rev Cancer (2021) 21(5):298–312. doi: 10.1038/s41568-021-00339-z
36. McCoach CE, Bivona TG. The Evolving Understanding of Immunoediting and the Clinical Impact of Immune Escape. J Thorac Dis (2018) 10(3):1248. doi: 10.21037/jtd.2018.03.60
37. Nijland M, Veenstra RN, Visser L, Xu C, Kushekhar K, van Imhoff GW, et al. HLA Dependent Immune Escape Mechanisms in B-Cell Lymphomas: Implications for Immune Checkpoint Inhibitor Therapy? Oncoimmunology (2017) 6(4):e1295202. doi: 10.1080/2162402X.2017.1295202
38. Zou R, Gu R, Yu X, Hu Y, Yu J, Xue X, et al. Characteristics of Infiltrating Immune Cells and a Predictive Immune Model for Cervical Cancer. J Cancer (2021) 12(12):3501. doi: 10.7150/jca.55970
39. Waldman AD, Fritz JM, Lenardo MJ. A Guide to Cancer Immunotherapy: From T Cell Basic Science to Clinical Practice. Nat Rev Immunol (2020) 20(11):651–68. doi: 10.1038/s41577-020-0306-5
40. Schatton T, Scolyer RA, Thompson JF, Mihm MC. Tumor-Infiltrating Lymphocytes and Their Significance in Melanoma Prognosis. Mol Diagnostics Melanoma (2014) 1102:287–324. doi: 10.1007/978-1-62703-727-3_16
41. Reissfelder C, Stamova S, Gossmann C, Braun M, Bonertz A, Walliczek U, et al. Tumor-Specific Cytotoxic T Lymphocyte Activity Determines Colorectal Cancer Patient Prognosis. J Clin Invest (2015) 125(2):739–51. doi: 10.1172/JCI74894
42. Seo A, Lee H, Kim E, Kim H, Jang M, Lee H, et al. Tumour-Infiltrating CD8+ Lymphocytes as an Independent Predictive Factor for Pathological Complete Response to Primary Systemic Therapy in Breast Cancer. Br J Cancer (2013) 109(10):2705–13. doi: 10.1038/bjc.2013.634
43. Thackray SJ, Mowat CG, Chapman SK. Exploring the Mechanism of Tryptophan 2, 3-Dioxygenase. Biochem Soc Trans (2008) 36(6):1120–3. doi: 10.1042/BST0361120
44. Jallad M-AN, Jurjus AR, Rahal EA, Abdelnoor AM. Triple Immunotherapy Overcomes Immune Evasion by Tumor in a Melanoma Mouse Model. Front Oncol (2020) 10:839. doi: 10.3389/fonc.2020.00839
45. Favre D, Mold J, Hunt PW, Kanwar B, Seu L, Barbour JD, et al. Tryptophan Catabolism by Indoleamine 2, 3-Dioxygenase 1 Alters the Balance of TH17 to Regulatory T Cells in HIV Disease. Sci Trans Med (2010) 2(32):32ra6–ra6. doi: 10.1126/scitranslmed.3000632
46. Pardee K, Green AA, Takahashi MK, Braff D, Lambert G, Lee JW, et al. Rapid, Low-Cost Detection of Zika Virus Using Programmable Biomolecular Components. Cell (2016) 165(5):1255–66. doi: 10.1016/j.cell.2016.04.059
47. Beatty GL, O’Dwyer PJ, Clark J, Shi JG, Bowman KJ, Scherle PA, et al. First-In-Human Phase I Study of the Oral Inhibitor of Indoleamine 2, 3-Dioxygenase-1 Epacadostat (INCB024360) in Patients With Advanced Solid Malignancies. Clin Cancer Res (2017) 23(13):3269–76. doi: 10.1158/1078-0432.CCR-16-2272
48. Jung KH, LoRusso P, Burris H, Gordon M, Bang YJ, Hellmann MD, et al. Phase I Study of the Indoleamine 2, 3-Dioxygenase 1 (IDO1) Inhibitor Navoximod (GDC-0919) Administered With PD-L1 Inhibitor (Atezolizumab) in Advanced Solid Tumors. Clin Cancer Res (2019) 25(11):3220–8. doi: 10.1158/1078-0432.CCR-18-2740
49. Mitchell TC, Hamid O, Smith DC, Bauer TM, Wasser JS, Olszanski AJ, et al. Epacadostat Plus Pembrolizumab in Patients With Advanced Solid Tumors: Phase I Results From a Multicenter, Open-Label Phase I/II Trial (ECHO-202/KEYNOTE-037). J Clin Oncol (2018) 36(32):3223. doi: 10.1200/JCO.2018.78.9602
Keywords: cervical cancer, differentially expressed genes, immune infiltration, tumor immunity, tumor-infiltrating CD8+ T cells
Citation: Zhang S, Wan J, Chen M, Cai D, Xu J and Chen Q (2021) Tumor-Infiltrating CD8+ T Cells Driven by the Immune Checkpoint-Associated Gene IDO1 Are Associated With Cervical Cancer Prognosis. Front. Oncol. 11:720447. doi: 10.3389/fonc.2021.720447
Received: 04 June 2021; Accepted: 23 September 2021;
Published: 27 October 2021.
Edited by:
Ye Wang, Qingdao University Medical College, ChinaReviewed by:
Tupa Basuroy, Massachusetts General Hospital and Harvard Medical School, United StatesTing Wang, Complete Genomics, United States
Copyright © 2021 Zhang, Wan, Chen, Cai, Xu and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qi Chen, MTc2OTI0NTUyQHFxLmNvbQ==; Junlan Xu, MjYwNzM2NTg4M0BxcS5jb20=
†These authors share first authorship