Analyzing and validating the prognostic value and immune microenvironment of clear cell renal cell carcinoma

doi:10.21203/rs.3.rs-274762/v1

Download PDF

Research

Analyzing and validating the prognostic value and immune microenvironment of clear cell renal cell carcinoma

https://doi.org/10.21203/rs.3.rs-274762/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 29 Mar, 2022

Read the published version in Animal Cells and Systems →

Version 1

posted

You are reading this latest preprint version

Background: There is still controversy regarding immunotherapy biomarkers. Therefore, we aimed to identify prognostic biomarkers related to immunotherapy for clear cell renal cell carcinoma (ccRCC).

Methods: Fragments Per Kilobase Million (FPKM) data and clinical characteristics for ccRCC patients from The Cancer Genome Atlas (TCGA) database were downloaded. Unsupervised consensus clustering analysis was performed to divide patients into different immune subgroups according to their single-sample gene set enrichment analysis (ssGSEA) scores. Then, we validated the differences in immune cell infiltration, prognosis, clinical characteristics and expression levels of HLA and immune checkpoint genes between different immune subgroups. Weighted gene coexpression network analysis (WGCNA) was used to identify the significant modules and hub genes that were related to the immune subgroups. A nomogram was established to predict the overall survival (OS) outcomes after independent prognostic factors were identified by least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analyses.

Results: Five clusters (immune subgroups) were identified. There was no significant difference in age, sex or N stage. And there were significant differences in race, T stage, M stage, grade, prognosis and tumor microenvironment. HLA gene families and CTLA4 showed significant differences between the five clusters, while PD1 and PDL1 did not. The red module was significant, and 14 hub genes were obtained. In addition, the nomogram containing LAG3 and GZMK accurately predicted OS outcomes of ccRCC patients.

Conclusion: LAG3 and GZMK are strongly related to immunity and may provide suggestions for ccRCC immunotherapy.

Allergy & Immune Disorders

Clinical Pharmacology

Environmental Policy

Clear cell renal cell carcinoma

microenvironment

prognosis

nomogram

bioinformatics

Renal parenchymal carcinoma is adenocarcinoma derived from renal tubular epithelial cells; 85% of these tumors are clear cell carcinomas, while others are granular cell carcinomas and mixed cell carcinomas. Hemorrhage, necrosis, cystic change and calcification are common in cancer ¹. Clear cell renal cell carcinoma (ccRCC) has the highest mortality rate among urinary system cancers. Its incidence rate and mortality rate account for 2% of all tumors. Approximately 100 thousand patients worldwide die of this disease ². China's incidence rate of renal cancer is approximately 3%, and the incidence rate of renal cancer in cities is significantly higher than that in rural areas ³. Approximately 30% of ccRCC is advanced when diagnosed, and approximately 10%-20% of early ccRCC will recur and metastasize after treatment. The 5-year survival rate is approximately 11.7% ⁴.

Compared with other tumors, ccRCC is less sensitive to chemotherapy and radiotherapy, and surgery is still the main treatment for nonadvanced ccRCC patients. Many years ago, the drug treatment of advanced ccRCC was limited to interleukin-2 and interferon alpha. However, due to its serious adverse reactions and extremely low efficiency, its clinical application is not widely used ⁵. The development and application of immune checkpoint inhibitors (ICIs) have improved the prognosis of patients with RCC. The CTLA-4 pathway plays an important role in the early stage of immune system activation, and the PD-1/PD-L1 pathway plays an important role in the tumor microenvironment of the immune system in the effect stage ⁶. The first-line treatment of advanced ccRCC with ravulizumab combined with ipilimumab reported in 2018 confirmed the first-line treatment status of this regimen in advanced renal clear cell carcinoma. A total of 1096 patients with advanced renal clear cell carcinoma were enrolled in this study. They were randomly divided into two groups: the ravulizumab combined with ipilimumab group and the sunitinib group. The objective remission rate of the combined group was significantly higher than that of the sunitinib group (42% vs. 27%, P < 0.001), and the complete remission rate of the combined group was also higher than that of the sunitinib group (9% vs. 1%). After a median follow-up of 25.2 months, the survival rates of the combined group and sunitinib group were 75% and 60%, respectively ⁷.

Although the success of immunotherapy is encouraging, countless patients have achieved remarkable results with its intervention. However, some patients did not respond to immunotherapy. The complexity and diversity of the tumor immune microenvironment (TIME) may determine its important impact on immunotherapy ^8,9. Immune cells and stromal cells are two important components of the extracellular matrix in the microenvironment. They have a wide range of interactions with tumor cells, including early tumor recruitment and the activation of stromal cells to form primitive precancerous stroma. These cells can change the cancer cell phenotype to a malignant phenotype. They also establish complex cell-cell interaction networks that help to improve and maintain the immunosuppressive microenvironment, promote immune escape, and ultimately promote the development of cancer. Therefore, the TIME is a multidirectional, dynamic and complex interaction network among immune cells, stromal cells and tumor cells ¹⁰. Therefore, in this study, we aimed to identify possible biomarkers for ccRCC through the tumor immune microenvironment.

Data Collection and Processing

Level 3 fragments per kilobase of transcript per million mapped reads (HTSeq-FPKM) RNA-seq data of ccRCC patients were downloaded from the UCSC Xena database (https://xenabrowser.net/datapages/) ¹¹, and clinical information for each patient was obtained. Samples with any incomplete information were excluded, and patients who survived for less than 1 day were also excluded. We definded overall survival (OS) as the period from the diagnosis to death for any reason and progression-free survival (PFS) as the period from the diagnosis to death or progression of the tumor.

Evaluation of the Immune Cell Infiltration Level

Single‑sample gene set enrichment analysis (ssGSEA) ¹² was performed to obtain the 29 immune signatures, including immune cells or immune function, of each sample, and the activities of the immune pathways were assessed using the GSVA package ¹³. In the tumor microenvironment, immune cells and stromal cells are two main types of nontumor components and have been proposed for the diagnostic and prognostic evaluations of tumors. The ESTIMATE algorithm ¹⁴ was used to evaluate the immune cell infiltration level (immune score), stromal content (stromal score) and comprehensive environmental score (ESTIMATE score).

Identification of Molecular Subtypes

Consistent clustering is a method used to provide quantitative evidence for determining the number and members of possible clusters in a dataset (e.g., a microarray gene expression dataset). The consistency clustering method involves subsampling from a group of items and determining clusters with a specific cluster count (k). Then, for the consensus value, two items with the same clustering in the same subsample are calculated and stored in the symmetric consistent matrix of each k value. The ConsensusClusterPlus package ¹⁵ was used to perform consistent clustering analysis of each ssGSEA score of each ccRCC sample.

Weighted Gene Coexpression Network Analysis (WGCNA)

The WGCNA package ¹⁶ was used to perform this analysis. We first selected the genes that varied greatly in different samples (i.e., genes whose variance was greater than all variance quartiles). The goodsamplesgenes function was used to evaluate whether the matrix data were qualified, that is, whether there was a missing value. If so, the gene was excluded. Then, the samples were clustered, and outliers were removed. After removing the outliers, we reconstructed the sample clustering tree and selected an appropriate soft threshold (β) to ensure scale-free distribution; that is, a small number of genes showed an absolute advantage (high expression), and most showed a disadvantage (low expression). By means of the power index, a certain power index was used to determine the correlation coefficient between genes, and a suitable soft threshold (β) was determined. According to the β value, the proximity matrix and topological matrix were obtained, and the genes were clustered by dissimilarity between genes. Then, the tree was cut into different modules by the dynamic cutting method (the minimum number of genes in the module was 50). Combining the correlation between each module and sample character, we analyzed the relationship between the module and phenotype. The module with the highest correlation coefficient with clinical characteristics was considered an important module. In addition, genetic significance (GS) and module membership (MM) were calculated. GS was defined as the level of correlation between gene expression and clinical characteristics. MM was defined as the correlation between the module’s own genes and gene expression profiles. For the important modules, the inclusion criteria for the hub genes were as follows: GS >0.4 and MM > 0.7.

Functional Enrichment Analysis

For important modules and hub genes, we analyzed their biological functions and possible pathways through Gene Ontology (GO) database and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The appeal operation was operated through the clusterProfiler ¹⁷, org.Hs.eg.db (http://www.bioconductor.org/packages/release/data/annotation/html/org.Hs.eg.db.html) and enrichplot (http://www.bioconductor.org/packages/release/bioc/html/enrichplot.html) packages. Fisher exact test and Benjamin Hochberg (B-H) multiple test correction method were used to correct the occurrence of false positive. The adjustied P < 0.05 was taken as the cut-off standard.

Construction and Validation of the Nomogram

To further verify the impact of the hub genes on the prognosis of patients, we combined the clinical information and prognostic information to construct the relevant nomogram. Least absolute shrinkage and selection operator (LASSO) regression analysis was first used to filter variables and prevent overfitting using the glmnet package¹⁸. After the influencing factors were screened by LASSO regression, multivariate Cox regression analysis was used to identify independent prognostic factors using the survival package ¹⁹. The nomogram was used to show the predictive results of these independent prognostic factors for the 1-, 3-, and 5-year OS rates by multivariate Cox regression analysis using the rms package ²⁰. Calibration curves were generated to test the predictive power of the nomogram. In addition, a time-dependent receiver operating characteristic (ROC) curve was used to judge the predictive power of the nomogram.

Statistical Analysis

All the above calculations were carried out with R software (version 4.0.2, 64 bit). Corrected P < 0.05 was considered statistically significant. When comparing whether there was a significant difference between the two groups of data, we used the rank-sum test. To determine whether there was a significant difference between multiple groups of data, we used the one-way analysis of variance (ANOVA) test. The Kaplan-Meier method was used to draw the patient survival curves, and the log-rank test was used to compare whether there was any difference in prognosis between the groups.

Immune Microenvironment Landscape of ccRCC

After using the ssGSEA method to measure the immune function of The Cancer Genome Atlas (TCGA) ccRCC transcriptome, we obtained the corresponding ssGSEA score of each sample according to the 29 immune-related pathways, and immune cell infusion was used to evaluate the immune function of ccRCC, as shown in Supplementary Table 1.

Consensus Clustering According to Molecular Subtypes

The results of consensus clustering of the ssGSEA scores were visualized using an empirical cumulative distribution function (CDF) plot and a delta area plot (Figure 1A and C). K=5 was regarded as the best value to cut these clusters, and the consensus matrix heat map is shown in Figure 1B. Therefore, the total ccRCC dataset was divided into 5 clusters: C1, C2, C3, C4, and C5. By using the ESTIMATE algorithm, we obtained three scores: the immune score, stromal score and ESTIMATE score (shown in Supplementary Table 2). A heatmap of the 5 clusters with their ssGSEA and ESTIMATE scores is shown in Figure 1D. Cluster C3 showed obvious differences compared with the other clusters.

For different clusters, there was no statistical difference in age (Figure 2A, P=0.1181), sex (Figure 2B, P=0.3233) or N stage (Figure 2E, P=0.3277). However, there were significant differences in race (Figure 2C, P=0.0069), T stage (Figure 2D, P=0.0001), M stage (Figure 2F, P=0.0157) and grade (Figure 2H, P<0.0001). Regarding the prognosis of patients, there were significant differences between the five clusters. Cluster 4 had the best prognosis for both OS (Figure 2G) and PFS (Figure 2I). Moreover, the tumor microenvironment using the ESTIMATE algorithm was statistically significant between the five clusters (Figure 3A-C). However, concerning the stromal score, there was no significant difference between the C4 and C5 clusters (Figure 3A). For the three ICIs and PD1 and PDL1, only some clusters showed significant differences (Figure 3D and E). Regarding CTLA4, except for the expression of CTLA4 between cluster 1 and cluster 4, there were significant differences in the other clusters (Figure 3F). Notably, almost all HLA genes showed significant differences between the five clusters (Figure 3G).

WGCNA

A total of 4555 genes with great variation in different samples were subjected to WGCNA. After it was clear that there was no missing value, an outlier was deleted (Figure 4A). β = 3 was used to establish a proximity matrix to make our gene distribution conform to the scale-free network (Figure 4B). After clustering genes and using the dynamic cutting method to cut the tree into different modules and merge similar modules, we obtained a total of 8 modules (Figure 4C). Among them, the red module was positively correlated with our immune cluster classification (Figure 5D). The red module showed the highest GS and MM based on an intramodular analysis (Figure 5E) and was regarded as a significant module. Based on the cutoff value for hub genes (MM > 0.7 and cor GS > 0.4), we ultimately obtained 14 hub genes: CD79A, FCRLA, GPR174, GZMK, JCHAIN, LAG3, MZB1, PDCD1, PLA2G2D and POU2AF1 (supplementary table 3).

Functional Enrichment Analysis

In the GO analysis, the significant modules were mostly enriched in cell killing, the external side of the plasma membrane and receptor ligand activity (Figure 5A). In the KEGG analysis, the modules were mostly enriched in cytokine−cytokine receptor interactions (Figure 5B). The main enriched functions and pathways of the 14 hub genes were related to immunity. Examples included the regulation of regulatory T cell differentiation, the regulation of the T cell apoptotic process and the negative regulation of T cell activation (Figure 5C).

Identification of the prognostic factors of nomogram

Eleven prognostic factors, including three hub genes, MZB1, GZMK and LAG3 were identified by LASSO regression analysis(supplementary figure 1).Multivariate Cox regression analysis showed that the final independent prognostic factors were LAG3, GZMK, radiotherapy, age, T stage, M stage and grade. The results of the multivariate Cox regression analysis are shown in supplementary figure 2, and the nomograms are shown in Figure 6A. The areas under the curve (AUCs) of the time-dependent ROC curves were 0.839, 0.802 and 0.769 for 1-, 3-, and 5-year OS, respectively (Figure 6B). The calibration curves of the 1-, 3-, and 5-year survival rates also showed that the nomogram had good prediction ability (Figure 6C-E).

In our study, we first obtained RNA-seq and clinical data of ccRCC patients from the TCGA database. Moreover, ssGSEA and the ESTIMATE algorithm were used to obtain the corresponding immune microenvironment scores, and through consistency clustering analysis, we obtained the relevant immune clusters. The immune microenvironments of these clusters were significantly different, and there were also significant differences between the immunosuppressive sites of the CTLA4 and HLA gene families. The related hub genes were found after performing WGCNA among different clusters, and the prognostic signature constructed by these hub genes and clinical information has great predictive ability. These biomarkers may be useful for RCC immunotherapy.

In recent years, the drug treatment of ccRCC has changed from the era of cytokines to the era of targeted drugs and even the new era of immunotherapy. Immunotherapy has become a choice for ccRCC treatment. ccRCC is characterized by multidrug resistance and biological heterogeneity. It is not sensitive to traditional radiotherapy and chemotherapy. Because there are a large number of inflammatory cells, such as T cells, natural killer (NK) cells, dendritic cells and macrophages, around RCCs, immunotherapy may be an effective treatment ²¹. With the expeditious development of high-throughput sequencing technology, there are an increasing number of studies on tumor biomarkers using bioinformatics technology by using the RNA-aeq data from lots of databases ^22-24. Therefore, in this study, we explored immune-related biomarkers of ccRCC.

At present, many clinical and basic experimental studies have shown that a variety of factors constitute a unique tumor microenvironment and promote the immune tolerance of tumor cells. One of the most prominent factors is the expression of a variety of negative costimulatory molecules, such as CTLA-4 ²⁵, TIM3 ²⁶, and LAG3 ²⁷. Due to the poor clinical effect of CTLA-4 and PD-1 in the treatment of KIRC, many researchers are actively seeking new negative costimulatory molecules combined with blocking therapy strategies. Remarkably, LAG3 is considered another promising target for cancer therapy. In addition to its inhibitory effect on T cell activation, it also has a synergistic inhibitory effect with PD-1 in various diseases, such as viral infection, parasites, hematological tumors and solid tumors. For example, the depletion of bovine T cells is caused by cryptosporidium infection ²⁸, chronic tuberculosis infectious diseases ²⁹, viral infectious diseases ³⁰, ovarian cancer ³¹ and chronic lymphoma ³².

Although LAG3 is another potential target of immunotherapy, a study on pan-cancer from the TCGA found that although the CD8 + T cell marker CD8a and LAG-3 are strongly coexpressed in most cancers, there are still three obvious exceptions: HPV + head and neck squamous cell carcinoma, RCC and glioblastoma ³³. Zelba et al. ³⁴ isolated tumor-infiltrating lymphocytes and autologous peripheral blood mononuclear cells (PBMCs) from patients with primary RCC. After staining and examining the cells with polychromatic flow cytometry, it was found that blocking PD-1 could lead to the upregulation of LAG-3; moreover, the double blocking of PD-1 and LAG-3 and the dual blocking of PD-1 and TIM-3 resulted in an increase in IFNγ release in vitro. It has also been proven that the double blocking of PD-1 and LAG-3 is a promising checkpoint blocking combination for ccRCC.

Another hub gene that we found to have a significant impact on prognosis was GZMK. There are five granzyme genes in humans (GZMA, GZMB, GZMH, GZMK and GZMM). GZMK is an apoptotic serine protease secreted by granules. Studies have shown that GZMK has a significant proinflammatory effect. It can induce the secretion of a variety of inflammatory factors, such as interleukin (IL)-1β, IL-6, and IL-8 and monocyte chemoattractant protein-1 (MCP-1). In addition, GZMK activated protease activated receptor 1 and induced cell proliferation ³⁵. In a phase I multicenter dose escalation study for advanced solid malignancies, in which 36 RCC patients were included and the safety and therapeutic efficacy of the anti-PD 1 antibody MEDI0680 were evaluated, the degree of tumor invasion of CD8+ T cells and the expression of the GZMK gene were also increased after using MEDI0680 ³⁶. Moreover, GZMK may be involved in the pathological mechanism of many kinds of tumors, such as colorectal cancer ³⁷, lung cancer ³⁸ and breast cancer ³⁹.

Our research also has some shortcomings. First, Our research is mainly based on bioinformatics analysis, all the data are from the open database and not from our own experience; thus, in-depth functional verification is needed. Second, there may be some deviations in the data obtained from public databases. Finally, all the samples examined in this study were from patients living in the United States, most of whom were Caucasian, so our results may not be representative of the global population. Therefore, more well-designed experiments with large sample sizes need to be used to further verify our conclusions.

In our study, we found two immune treatment-related biomarkers: LAG3 and GZMK. We also found that the prognosis of ccRCC patients could be accurately determined by combining their clinical information and the expression levels of LAG3 and GZMK.

Conflict of interests statement

The Authors declares that there is no conflict of interest.

Funding

None

Acknowledgments

We thank Yi Liao for their assistance with data calculation.

Funding

Not applicable

Ethics approval

Not applicable

Data availability statement

The data generated within the study is shown in this manuscript. Any raw data or

analysis would be available from the corresponding author upon request.

Authors’ Contributions

Xin Liu and Jie Chen designed this study. Jingwei Ke performed statistical analysis. Xin Liu and Jie Chen wrote the paper. All authors read and approved the final version of the paper.

Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394-424.
Torre LA, Siegel RL, Ward EM, Jemal A. Global cancer incidence and mortality rates and trends--an update. Cancer Epidemiol Biomarkers Prev. 2016;25(1):16-27.
Pang C, Guan Y, Li H, Chen W, Zhu G. Urologic cancer in China. Jpn J Clin Oncol. 2016;46(6):497-501.
La Rochelle J, Klatte T, Dastane A, et al. Chromosome 9p deletions identify an aggressive phenotype of clear cell renal cell carcinoma. Cancer. 2010;116(20):4696-4702.
Escudier B, Pluzanska A, Koralewski P, et al. Bevacizumab plus interferon alfa-2a for treatment of metastatic renal cell carcinoma: a randomised, double-blind phase III trial. Lancet. 2007;370(9605):2103-2111.
Zhou X, Hou W, Gao L, Shui L, Yi C, Zhu H. Synergies of antiangiogenic therapy and immune checkpoint blockade in renal cell carcinoma: from theoretical background to clinical reality. Front Oncol. 2020;10:1321.
Motzer RJ, Tannir NM, McDermott DF, et al. Nivolumab plus Ipilimumab versus sunitinib in advanced renal-cell carcinoma. N Engl J Med. 2018;378(14):1277-1290.
Yang J, Zhang C. Regulation of cancer-immunity cycle and tumor microenvironment by nanobiomaterials to enhance tumor immunotherapy. WIREs Nanomed Nanobiotechnology. 2020;12(4):e1612.
Shang J, Song Q, Yang Z, et al. Analysis of PD-1 related immune transcriptional profile in different cancer types. Cancer Cell Int. 2018;18:218.
Hanahan D, Coussens LM. Accessories to the crime: functions of cells recruited to the tumor microenvironment. Cancer Cell. 2012;21(3):309-322.
Goldman M, Craft B, Swatloski T, et al. The UCSC cancer genomics browser: update 2015. Nucleic Acids Res. 2015;43(Database issue):D812-D817.
Barbie DA, Tamayo P, Boehm JS, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462(7269):108-112.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Yoshihara K, Shahmoradgoli M, Martínez E, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572-1573.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Yu G, Wang LG, Han Y, He QY. ClusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284-287.
Engebretsen S, Bohlin J. Statistical predictions with glmnet. Clin Epigenetics. 2019;11(1):123.
Rizvi AA, Karaesmen E, Morgan M, et al. Gwasurvivr: an R package for genome-wide survival analysis. Bioinformatics. 2019;35(11):1968-1970.
Huang C, Liu Z, Xiao L, et al. Clinical significance of serum CA125, CA19-9, CA72-4, and fibrinogen-to-lymphocyte ratio in gastric cancer with peritoneal dissemination. Front Oncol. 2019;9:1159.
Noessner E, Brech D, Mendler AN, Masouris I, Schlenker R, Prinz PU. Intratumoral alterations of dendritic-cell differentiation and CD8⁺ T-cell anergy are immune escape mechanisms of clear cell renal cell carcinoma. Oncoimmunology. 2012;1(8):1451-1453.
Liao Y, Xiao H, Cheng M, Fan X. Bioinformatics analysis reveals biomarkers with cancer stem cell characteristics in lung squamous cell carcinoma. Front Genet. 2020;11:427.
Liao Y, Wang Y, Cheng M, Huang C, Fan X. Weighted gene coexpression network analysis of features that control cancer stem cells reveals prognostic biomarkers in lung adenocarcinoma. Front Genet. 2020;11:311.
Liao Y, Yin G, Wang X, Zhong P, Fan X, Huang C. Identification of candidate genes associated with the pathogenesis of small cell lung cancer via integrated bioinformatics analysis. Oncol Lett. 2019;18(4):3723-3733.
Simonaggio A, Epaillard N, Elaidi R, et al. Impact des signatures moléculaires sur le choix du traitement systémique du cancer du rein métastatique [Impact of molecular signatures on the choice of systemic treatment for metastatic kidney cancer]. Bull Cancer. 2020;107(5s):S24-S34.
Li B, Severson E, Pignon JC, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174.
Alme AK, Karir BS, Faltas BM, Drake CG. Blocking immune checkpoints in prostate, kidney, and urothelial cancer: an overview. Urol Oncol. 2016;34(4):171-181.
Ivanova DL, Krempels R, Denton SL, et al. NK cells negatively regulate CD8 T cells to promote immune exhaustion and chronic toxoplasma gondii infection. Front Cell Infect Microbiol. 2020;10:313.
Boer MC, van Meijgaarden KE, Bastid J, Ottenhoff TH, Joosten SA. CD39 is involved in mediating suppression by Mycobacterium bovis BCG-activated human CD8⁺ CD39⁺ regulatory T cells. Eur J Immunol. 2013;43(7):1925-1932.
Chen H, Moussa M, Catalfamo M. The role of immunomodulatory receptors in the pathogenesis of HIV infection: a therapeutic opportunity for HIV cure? Front Immunol. 2020;11:1223.
Imai Y, Hasegawa K, Matsushita H, et al. Expression of multiple immune checkpoint molecules on T cells in malignant ascites from epithelial ovarian carcinoma. Oncol Lett. 2018;15(5):6457-6468.
Saleh RR, Peinado P, Fuentes-Antrás J, et al. Prognostic value of lymphocyte-activation gene 3 (LAG3) in cancer: a meta-analysis. Front Oncol. 2019;9:1040.
Panda A, Rosenfeld JA, Singer EA, Bhanot G, Ganesan S. Genomic and immunologic correlates of LAG-3 expression in cancer. Oncoimmunology. 2020;9(1):1756116.
Zelba H, Bedke J, Hennenlotter J, et al. PD-1 and LAG-3 dominate checkpoint receptor-mediated T-cell inhibition in renal cell carcinoma. Cancer Immunol Res. 2019;7(11):1891-1899.
Joeckel LT, Allison CC, Pellegrini M, Bird CH, Bird PI. Granzyme K-deficient mice show no evidence of impaired antiviral immunity. Immunol Cell Biol. 2017;95(8):676-683.
Naing A, Infante J, Goel S, et al. Anti-PD-1 monoclonal antibody MEDI0680 in a phase I study of patients with advanced solid malignancies. J Immunother Cancer. 2019;7(1):225.
Zhang L, Yu X, Zheng L, et al. Lineage tracking reveals dynamic relationships of T cells in colorectal cancer. Nature. 2018;564(7735):268-272.
Weng TY, Wang CY, Hung YH, Chen WC, Chen YL, Lai MD. Differential expression pattern of THBS1 and THBS2 in lung cancer: clinical outcome and a systematic-analysis of microarray databases. PLoS One. 2016;11(8):e0161007.
Lee HJ, Lee JJ, Song IH, et al. Prognostic and predictive value of nanostring-based immune-related gene signatures in a neoadjuvant setting of triple-negative breast cancer: relationship to tumor-infiltrating lymphocytes. Breast Cancer Res Treat. 2015;151(3):619-627.

S1.jpg
The results of LASSO regression analysis
S2.jpg
The results of multivariate Cox regression analysis.
supplementarytable1.xls
The results of ssGSEA score.
supplementarytable2.xls
The results of the immune score, stromal score and ESTIMATE score.
supplementarytable3.xls
The results of GS and MM value from WGCNA.

Download PDF

Journal Publication

published 29 Mar, 2022

Read the published version in Animal Cells and Systems →

Version 1

posted

You are reading this latest preprint version

Analyzing and validating the prognostic value and immune microenvironment of clear cell renal cell carcinoma

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Materials And Methods

Data Collection and Processing

Evaluation of the Immune Cell Infiltration Level

Identification of Molecular Subtypes

Weighted Gene Coexpression Network Analysis (WGCNA)

Functional Enrichment Analysis

Construction and Validation of the Nomogram

Statistical Analysis

Results

Immune Microenvironment Landscape of ccRCC

Consensus Clustering According to Molecular Subtypes

WGCNA

Functional Enrichment Analysis

Identification of the prognostic factors of nomogram

Discussion

Conclusion

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1