GNRH1 and LTB4R Might Be Novel Immune-Related Prognostic Biomarkers in Clear Cell Renal Cell Carcinoma (ccRCC)

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

PDF

Research Article

GNRH1 and LTB4R Might Be Novel Immune-Related Prognostic Biomarkers in Clear Cell Renal Cell Carcinoma (ccRCC)

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

This work is licensed under a CC BY 4.0 License

Journal Publication

published 06 Jul, 2021

Read the published version in Cancer Cell International →

You are reading this latest preprint version

Background: Clear cell renal cell carcinoma (ccRCC) occupied most of renal cell carcinoma (RCC), which associated with poor prognosis. The purpose of this study is to screen novel and prognostic biomarkers for patients with ccRCC.

Methods and Results: Firstly, Gene Expression Omnibus database was used to collect microarray data for weighted gene co-expression network construction. Gene modules related to prognosis which interest us most were picked out. 90 hub genes were further chosen in the key modules, two of which including GNRH1 and LTB4R were screened and validated as immune-related prognostic biomarkers. Several methods including survival analysis, spearman correlation analysis, HPA, One-way ANOVA and ROC analysis were used for the validation of immune-related prognostic biomarkers. We further explored the relationship between immune-related prognostic biomarker expressions and immunocytes. Finally, gene set enrichment analysis (GSEA) demonstrated that the two immune-related prognostic biomarkers were significantly correlated with cell cycle.

Conclusions: Generally speaking, the present study has identified two novel prognostic biomarkers for patients with ccRCC, which showed strong correlation with prognosis of patients with ccRCC, could further be used as potential prognostic biomarkers in ccRCC.

Cancer Biology

immune-related genes

clear cell renal cell carcinoma

prognostic biomarkers

immune infiltration

WGCNA

RCC is an adenocarcinoma derived from renal tubular epithelial cells, 85% of which are ccRCC[1]. ccRCC occupied poor prognosis[1, 2]. There will be approximately 152,160 RCC cases diagnosed and 27,560 RCC-related deaths in the United States in 2021[2]. The prognosis of ccRCC is influenced by mechanisms covering complex networks of gene interactions[3]. What is worse, ccRCC is insensitivity to both radiotherapy and chemotherapy[3, 4]. Therefore, because of the poor prognosis and the lack of new therapeutic target, we aim to screen out novel and specific prognostic markers for patients with ccRCCs in the present study.

Immunotherapy is a type of cancer treatment, which helps immune system fight cancer[5]. Nowadays, more and more studies have indicated that tumors could be treated by immunotherapy effectively and safely[6, 7]. Thus, identification of immune-related prognostic biomarkers has been one of the research hotspots in cancer treatment.

Thus, in the present study, based on gene expression profiles collected from The Cancer Genome Atlas (TCGA) database, Gene Expression Omnibus (GEO) database, and AffyExpress database (including large size of gene expression matrix about cancers, which have been widely used in gene exploring), we attempted to screen out several immune-related prognostic biomarkers which might occur the potential of prognosis prediction and treatment. To the best of our knowledge, this study might be one of the first studies to identify immune-related prognostic biomarkers in ccRCC directly. Some widely used methods such as weighted gene co-expression network analysis (WGCNA) were also applied in this study to identify two immune-related modules and further screen out two novel immune-related prognostic biomarkers. What we should mention was that we creatively defined prognosis-related modules as modules associated with not only survival time but also survival status. Previous studies also defined prognosis-related modules as modules only associated with one of the two traits (survival time or survival status). But we realized that both of survival status and survival time were important indispensable yardsticks of prognosis. Therefore, two modules associated with both the two traits were preliminarily selected in the present study by WGCNA. Two immune-related genes (IRGs) in the prognosis-related modules were further identified and validated as potential prognostic biomarkers of patients with ccRCC.

To sum up, the present study indicated that two IRGs might be novel immune-related prognostic biomarkers for ccRCC. The two IRGs were validated to show strong potential for prognosis prediction and treatment of ccRCC.

Datasets and Immune-Related Genes (IRGs)

Figure 1 showed the workflow of the present study, describing the whole process of IRGs correlated to prognosis exploring. Raw expression data of ccRCC (TCGA-KIRC) were retrieved from TCGA database (https://genomecancer.ucsc.edu/). This expression data was firstly showed as count number. Standardization of TCGA-KIRC data were conducted using R package “DEseq.2”[8]. Generally, we also conducted normalization and log2 transformation via this tool. Finally, 530 ccRCCs containing necessary clinical feature information were chosen, which were further used for WGCNA. Furthermore, GSE29609[9] including 39 tumor samples was retrieved from Gene Expression Omnibus (GEO) database[10] (http://www.ncbi.nlm.nih.gov/geo/), which was used as a validation cohort. Complete survival information for this dataset was also downloaded in our study. This dataset was normalized and log2 transformed based on R package “affy”[11]. In addition, we directly collected 53 ccRCCs which stored in E-MTAB-3267 from AffyExpress database (https://www.ebi.ac.uk/arrayexpress/) to further verify our results.

Moreover, immune-related genes (IRGs) were collected from the ImmPort database (https://www.immport.org). In total, 2,499 IRGs were retrieved, 1811 of which overlapped with gene list of TCGA-KIRC expression matrix were further included in the present study.

Construction of Co-Expression Network

Firstly, we checked the expression matrix of all genes based on gsg (goodSamplesGenes) and sample network methods in R package “WGCNA”[12]. A standard of Z.Ku ≥ -2.5 was set to pick out qualified ccRCC samples for constructing a co-expression network. A method called scale free topology criterion was used to pick out approriate β (soft threshold power beta) in WGCNA. In addition, we transformed adjacency into TOM and further classified genes into gene modules. Manual (interactive) branch cutting approach, automatic single block analysis and 2 block analysis were applied for gene module identification[13]. In the three methods, we set the same cut-off criteria as follows: A relatively large minimum module size of minClusterSize = 30, and a medium sensitivity (deepSplit = 2). To merge modules with high correlation, a cut line was also set by reckoning dissimilarity of module eigengenes (MEs).

Prognosis-Related Module Identification

In the present study, both of survival status and survival time are regarded as important yardsticks of prognosis. Thus, we attempted to screen out hub modules, which associated with both the traits. Gene Significance (GS) which could quantify the association between genes and trait was firstly explored. Module significance (MS) was secondly defined based on GS. By doing these, modules which met the following criteria were regarded as prognosis-related modules: 1. Module must show significant association with both survival status (overall survival status (OS)) and survival time (OS time); 2. Module showed the highest MS of OS or OS time. We then obtained the Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis for genes in prognosis-related modules based on R package “clusterProfiler”[14]. In the present study, we thought BPs and KEGG pathway terms with P < 0.05 were functional enriched.

Novel Immune-Related Prognostic Biomarker Identification

We immediately screened out hub genes in WGCNA by using the standards of (|cor.geneModuleMembership| > 0.8 and |cor.geneTraitSignificance| > 0.2). Finally, genes overlapped between hub genes in WGCNA and IRGs were considered as immune-related prognostic biomarker, which were validated in subsequent analysis.

Internal Validation of Immune-Related Prognostic Biomarkers

Two survival types including overall survival (OS) and disease-free survival (DFS) were conducted by using GEPIA webtool[15], as a preliminary exploration of the prognostic value of immune-related prognostic biomarkers (http://gepia.cancer-pku.cn/). Furthermore, immune-related prognostic biomarker expression levels between normal samples and PAAD samples were compared as an internal validation. Statistical significance was measured by using unpaired t test. Besides, stage plots were plotted to measure the difference of hub gene expression in different stages. The statistical significance of stage plots was measured by using ANOVA.

mRNA Expression Level and Translation-Level Expression of Prognostic Biomarker Verification

After an internal validation by using GEPIA, we validated expression difference of prognostic biomarker in transcription level between ccRCC samples and normal samples based on the Oncomine database[16] (https://www.oncomine.org/), as an external validation. In this study, four analysis were included, which were collected by the Oncomine database. We further verified the protein expression level of prognostic biomarkers between normal tissues and ccRCC tissues based on The Human Protein Atlas (HPA) database[17–19] (https://www.proteinatlas.org/).

RNA extraction and qRT-PCR

The expression patterns of the GNRH1 and LTB4R genes were evaluated in matched ccRCC and para-cancerous tissues. The HiPure Total RNA Mini Kit was used to extract total RNA from the tissues and we used NanoDrop to quantify the RNA which was then reverse transcribed into cDNA by ReverTra Ace qPCR RT Kit (Toyobo, China). Finally, we performed qRT-PCR analysis of cDNA with iQTM SYBR® Green Supermix (Bio-Rad) in a final volume of 20 µl. LTB4R primer: 5'-AGCTTTGTGGTGTGGAGTATCC-3'(forward),5'-GCAACCAGCCAGTCCAAAAC-3'(reverse), GNRH1 primer: 5'-CAAAAACTCCTAGCTGGCCTT-3'(forward), 5'-CAGTTGACCAACCTCTTTGACT-3'(reverse). GAPDH primer: 5'-TGCACCACCAACTGCTTAG-3'(forward), 5'-GATGCAGGGATGATGTTC-3'(reverse).

Western blotting

The Caki1 cells and normal kidney cells were separately mixed with the lysis buffer and proteinase inhibitors to isolate total proteins. SDS - PAGE gel was used to separate proteins which was then transferred to PVDF membranes (Millipore). After two hours closed in skimmed milk, we incubated the membranes in primary and secondary antibodies respectively. The ECL western blotting detection kit (Millipore) was used to detect the resultant bands. We performed all experiments in triplicate at least. The primary antibodies included anti-GNRH1(1:1000, Abcam), anti-LTB4R (1:1000, Abcam) and anti-GAPDH (1:1000, Abcam).

Prognostic Value of Immune-Related Prognostic Biomarkers Validation

Furthermore, to better understand the prognostic value of novel biomarkers, we used 3 analysis methods in this step. Firstly, by using GSE29609 (n = 39) and E-MTAB-3267 (n = 53), ccRCCs in the two datasets were grouped into high- and low- expression groups via the most appropriate score cutoff for tissues splitting calculated using R package “maxstat”[20]. We immediately conducted three kinds of survival analyses (OS, CSS, PFS) via R package “survival”[21]. Secondly, by using TCGA-KIRC data, receiver operating characteristic (ROC) curves were drawn to see if immune-related prognostic biomarkers could distinguish ccRCC tissues from normal tissues. ROC analysis was conducted based on R package “pROC”[22]. Furthermore, to measure the resolving power, we calculated the area under curve (AUC). In this study, we considered an immune-related prognostic biomarker with AUC > 0.75 had strong prognostic value. Thirdly, spearman correlation analysis was obtained based on SPSS (version 21.0), R package “ggstatsplot”[23] was used for visualization.

Immune-Related Prognostic Biomarker Function Analysis

In this study, for understanding the lurking functions of immune-related prognostic biomarkers, the median value of each biomarker was firstly evaluated relying on TCGA-KIRC data. In total, 530 tumor tissues were split into 2 groups (high- and low- expression groups). In the present study, “c2.cp.kegg.v6.2.symbols.gmt” was chosen as the enrichment gene set. Gene set enrichment analysis (GSEA) (http://software.broadinstitute.org/gsea/index.jsp)[24] was obtained. A signaling pathway was meaningful when nominal P < 0.05, a gene size (n) ≥ 100 and a FDR < 25% in the present part.

Association between Immune-Related Prognostic Biomarker Expression and Immunocytes Exploring

Immunocytes have been proved to be independent predictors of survival in cancers. The ssGSEA was a tool which could quantify the relative infiltration of 28 tumor-infiltrating immune cell types by using an enrichment score for representation of relative abundance. Thus, in this study, we investigated the relationship between expression levels of selected biomarkers and immunocytes applying ssGSEA method by using R package “GSVA”[25].

Screening out of Two Prognosis-Related Modules

After weeding out 21 outlier samples, 509 ccRCCs were included for WGCNA (Fig. 2). The beta (β) = 5 (scale free R2 = 0.84) was further set as soft-thresholding for further adjacencies calculation (Fig. 3). Then we split IRGs into gene modules. High-related modules were further merged for avoiding over-fitting. Eventually, as shown in Fig. 4A, totally 18 modules were identified. In this method, no-significant genes were stored in grey module, which were excluded from subsequent analysis. These genes were regarded as biomarkers with no-significance. Furthermore, brown module was associated with OS time (P = 0.001, R2 = -0.14, Fig. 4B) and OS (P = 0.35, R2 = 0.35, Fig. 4B). Turquoise module was associated with not only OS time (P = 2e-05, R2 = -0.19, Fig. 4B) but also OS (P = 5e-05, R2 = 0.18, Fig. 4B). The MS (of OS time) of turquoise module was higher than MSs of any other modules meanwhile the MS (of OS) of brown module was the highest compared to others (Fig. 4D). Moreover, the associations between MM and GS for OS time (cor = 0.4, P = 2.3e-27)/ OS (cor = 0.66, P = 1.3e-120) in brown module were also significant (Fig. 4C). The same trends existed in turquoise module (OS time: cor = 0.65, P = 5.1e-142; OS: cor = 0.42, P = 1.4e-78). Both the two modules (brown, turquoise) reached the standards previous mentioned, which were considered as prognosis-related modules in the present study. We also created a network heatmap in the present study (Fig. 5A). As shown in Fig. 5B, the classical MDS plot concluded that the 18 modules were independent of each other.

Function Exploring of Genes in Prognosis-Related Modules

GO and KEGG pathway enrichment analysis were conducted to explore the potential of genes in prognosis-related modules. As shown in Fig. 6A, GO analysis indicated that genes in brown module were involved in mitotic nuclear division, nuclear division, organelle fission, chromosome segregation, mitotic sister chromatid segregation, sister chromatid segregation, microtubule cytoskeleton organization involved in mitosis, regulation of mitotic nuclear division, regulation of nuclear division and nulcear chromosome segregation (Fig. 6A). Moreover, genes in turquoise module were mainly enriched in epithelial cilium movement, cilium assembly, cilium organization and cilium movement (Fig. 6C). The KEGG enrichment analysis indicated that genes in brown module significantly enriched in cell cycle, progesterone-mediated oocyte maturation, p53 signaling pathway and oocyte meiosis (Fig. 6B) meanwhile genes in turquoise module were mainly enriched in herpes simplex virus1 infection (Fig. 6D).

Two Genes Were Regarded as Novel Immune-Related Prognostic Biomarkers

10 genes from brown module (Fig. 7A) and 80 genes from turquoise module (Fig. 7B) were selected by using the cut-off criterion of |cor.geneModuleMembership| > 0.8 and |cor.geneTraitSignificance| > 0.2. Finally, two genes including GNRH1 (gonadotropin releasing hormone 1) and LTB4R (leukotriene B4 receptor) overlapping in hub genes in WGCNA and IRGs were screened out (Fig. 7C).

Internal Validation of the Two Immune-Related Prognostic Biomarkers

After screening out four hub genes though comprehensive bioinformatics analysis, we validated these potential prognostic biomarkers. Firstly, based on the GEPIA webtool, Higher expression of GNRH1 was significantly related to worse OS (Hazard ratio [HR] = 1.9, P = 3.4e-05, Fig. 8A). Moreover, we concluded that patients in LTB4R high-expression group obviously occupied worse OS (HR = 1.9, P = 3.2e-05, Fig. 8D), accurately. Then the mRNA expression of prognostic biomarkers between tumor tissues and normal tissues was compared, the results indicated that both the genes including GNRH1 (P < 0.05, Fig. 8B), and LTB4R (P < 0.05, Fig. 8E) were significantly higher expressed in ccRCC samples compared to normal samples. Moreover, the results also suggested that high expression of GNRH1 (F = 2.63, P = 0.0497; Fig. 8C) and LTB4R (F = 3.16, P = 0.0243; Fig. 8f) were significantly related to higher tumor stage.

External Validation of the Two Immune-Related Prognostic Biomarkers

Expression levels of prognostic biomarkers between normal samples and ccRCC samples were compared based on Oncomine database. We indicated the same conclusion that GNRH1 (P = 0.049, Fig. 9A) and LTB4R (P = 0.003, Fig. 9D) had higher expression in ccRCC tissues than in normal tissues. By using HPA database, we investigated the translational level expression of immune-related prognostic biomarkers. GNRH1 obtained medium staining in ccRCC samples compared to normal samples (not detected) (Fig. 10A). Unfortunately, there was no significant difference between ccRCC and normal tissues on the level of translation (Fig. 10B). The prognostic values of immune-related prognostic biomarkers were further validated after verification of expressions. As the results suggested, GNRH1 expression was significantly associated with not only OS time (r = -0.19, P < 0.001, Fig. 11A) but also OS status (r = 0.21, P < 0.001, Fig. 11B). Moreover, LTB4R expression was significantly related to OS time (r = -0.17, P < 0.001, Fig. 11C) and OS status (r = 0.26, P < 0.001, Fig. 11D). Also, the survival analysis demonstrated that there was a trend that higher expression of GNRH1 was significantly correlated to worse OS (GSE29609), CSS (GSE29609), PFS (GSE29609), PFS (E-MTAB-3267); suggested by Fig. 12A (P = 0.310), Fig. 12B (P = 0.450), Fig. 12C (P = 0.660), Fig. 12D (P = 0.063), respectively. Meanwhile ccRCC patients with higher LTB4R expression was related to short OS time, accurately (P = 0.034, GSE29609, Fig. 12E). The trend that higher expression of LTB4R was associated with worse CSS was verified by using GSE29609 (P = 0.660, Fig. 12F). Furthermore, ccRCC patients with higher expression of LTB4R was significantly correlated to worse PFS, as Fig. 12G (P = 0.035, GSE29609) and Fig. 12H (P = 0.0028, E-MATB-3267) suggested.

In addition, we also conducted ROC analysis. The result demonstrated that GNRH1 could significantly distinguish ccRCC samples from normal samples, suggested by Fig. 13A (AUC = 0.770). Interestingly, LTB4R showed strong potential for ccRCC diagnosis (AUC = 0.829, Fig. 13B). All the results above indicated that the two immune-related prognostic biomarkers we screened were credible.

Experimental validation

We further validated the expression of GNRH1 and LTB4R in ccRCC cell line by qRT-PCR and western blotting. qRT-PCR analysis showed that GNRH1 and LTB4R mRNA levels were significantly higher in ccRCC tissues than normal kidney tissues (Fig. 14A-B). In other hand, compared with the normal renal epithelial cell line HK2, expression of GNRH1 and LTB4R is significantly increased in the ccRCC cell line Caki1 (Fig. 14C-D).

Identification of Immune-Related Prognostic Biomarker Related Pathways

The results of GSEA indicated that GNRH1 was significantly enriched in one KEGG signaling pathway including olfactory transduction (nominal P = 0.014, n = 368, FDR = 5.164%, Fig. 15C). Meanwhile, we found that LTB4R was significantly related to adipogenesis (nominal P = 0.035, n = 191, FDR = 22.685%).

Association between Immune-Related Prognostic Biomarker Expression and Immune Infiltration Level in ccRCC

Immune infiltration was reported to be associated with survival and progression of cancers. Thus, the association between prognostic biomarker and immune infiltration level was obtained by applying ssGSEA. GNRH1 expression was significantly associated with several immune cell types, such as central memory CD8 T cell, CD56dim natural killer cell, and immature dendritic cell (Fig. 15A). Moreover, LTB4R was significantly related to activated CD4 T cell, CD56dim natural killer cell, and MDSC. These results indicated that expression of the two prognostic biomarkers might affect the immune infiltration levels of some tumor-infiltrating immune cell types.

Renal cell carcinoma (RCC) is the second common malignancy after bladder cancer in the urinary system and ccRCC is the most common one representing approximately 85–90% of all renal cancer diagnosis which can occur at any age. The tumour microenvironment including various infiltrating immune cells, fibroblasts, numerous cytokines, malignant tumour cells and so on was considered a complex and dynamic ecosystem in which immune response had broad influences on tumour growth, invasion and metastasis[26]. Recent studies had shown that immune-associated genes (IAG) played an important role in the development and progression of ccRCC[27, 28]. There were some promising new treatment for ccRCC including programmed death-1 (PD-1)/programmed death-ligand 1 (PD-L1) or cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) inhibitors and so on while these immunotherapies were limited by low long-term response rates and lack of reliable prognostic factors[29]. So we determined to find some novel markers related to both immunity and prognosis for diagnosis and treatment of ccRCC.

In this study, gene co-expression network was constructed on the basis of resemblance of expression profiles during samples and we explored the interconnections between genes by mean of WGCNA analysis. Meanwhile, we collected immune-related genes from various databases. By taking the intersection of the relevant results of WGCNA analysis and the data of immune genes, we got two final hub genes (GNRH1, LTB4R) which related to poor prognosis and immunity. In past studies, these two genes were reported to be correlated with prognosis of cancer. Expression of LTB4R was essential for CD8 + T cell infiltration into tumours[30] and invasiveness of ovarian cancer cells and breast cancer cells was increased by MMP-2 pathways[31] and IL-8 pathway[32] respectively. In a word, these studies highlighted the critical nature of BLT2 in tumour survival, invasiveness, chemoresistance and metastasis, as well as its potential as a therapeutic target for some cancers. During another study, expression of GNRH1 was considered may be a prognostic factor for metastatic spread of tumor cells based on the result that the up-regulation of GNRH1 expression suggested the presence of tumour cells in the circulation of cancer patients[33].

Six clinical features (OS.time, OS, Age, Gender, Laterality, Pathologic-stage) were brought into our study, and we eventually screened the 19 modules related to the prognosis of ccRCC. During them, the turquoise module had highest positively relevance to the OS.time and the brown module was most positively correlated with the OS. The turquoise module incuded 1175 genes and the brown module included 956 genes. We ended up with a total of 90 hug genes from these two modules. On the other hand, we got the Immune-related genes from ImmPort database which included 1811 genes. We found out the real hub genes which common to the turquoise module, the brown module and the Immune related genes. Two final hub genes (GNRH1, LTB4R) related to poor prognosis of ccRCC were ensured and we also investigated the correlation between final hub genes, clinical features and immune cells. In order to further investigate the mechanism of these genes in regulating tumor genesis, GSEA analysis was performed by using TCGA-KIRC data from TCGA database. Finally, we used various methods to validate the hub genes internally and externally respectively. However, the biggest limitation for our study was that we did not validate these two genes in vivo and in vitro. We need to verify this at the cellular, tissue and animal level and so on.

In summary, this study used a variety of bioinformatics analysis methods to identify two new hub genes that might play an important role on the occurrence, development and prognosis of ccRCC. At the same time, we predicted the potential function of hub genes involved in cell cycle regulation, which are associated with staging and prognosis. The two hub genes may become novel biomarkers of ccRCC in human. However, this study only selected potential ccRCC biomarkers related to tumor development and prognosis and the effect of candidate biomarkers need to be verified by further molecular biology experiments.

ccRCC: Clear cell renal cell carcinoma；

GSEA: Gene set enrichment analysis；

ROC: Receiver operating characteristic curve;

HPA : The Human Protein Atlas;

ANOVA : Analysis of Variance;

TCGA: The Cancer Genome Atlas;

GEO: Gene Expression Omnibus;

WGCNA: Weighted gene co-expression network analysis;

IRG: Immune-related genes;

GS: Gene Significance;

MS: Module significance;

MEs: Module eigengenes;

OS: Overall survival status;

OS time: Overall survival time;

GO: Gene Ontology enrichment analysis;

KEGG: Kyoto Encyclopedia of Genes and Genomes enrichment analysis;

DFS: Disease-free survival;

PAAD: Pancreatic adenocarcinoma;

GEPIA: Gene Expression Profiling Interactive Analysis;

GSVA: Gene Set Variation Analysis

HR: Hazard ratio.

Data availability statement

The data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgements

We would like to acknowledge the TCGA and GEO database developed by the National Institutes of Health (NIH).

Funding

The present study was supported by the Medical science and technology innovation platform support project of Zhongnan Hospital of Wuhan University (PTXM2019006).

Author contributions

T.L., H. W., and X. Y. conceived and designed the study, H. W., and X. Y. performed the analysis procedures, G. D., X. Y., Z. C., K. T., X. B., H. W., and T. L. analyzed the outcomes, T. L., and H. W. contributed analysis tools, H. W., and X. Y. contributed to the writing of the manuscript. All authors reviewed the manuscript and consented for publication. Our study did not require ethical board approval because it did not involve human or animal trials.

Declaration of Conflicting Interests

The authors declare that no competing interests exist

Declaration

All authors consent for publication.

Ethical Statement

This study was approved by the Ethics Committee of Wuhan University Zhong nan Hospital.

Wettersten HI, Aboud OA, Lara PJ, Weiss RH: Metabolic reprogramming in clear cell renal cell carcinoma. NAT REV NEPHROL 2017, 13(7):410-419.
Siegel RL, Miller KD, Fuchs HE, Jemal A: Cancer Statistics, 2021. CA: A Cancer Journal for Clinicians 2021, 71(1):7-33.
Cerbone L, Cattrini C, Vallome G, Latocca MM, Boccardo F, Zanardi E: Combination therapy in metastatic renal cell carcinoma: Back to the future? SEMIN ONCOL 2020, 47(6):361-366.
Braun DA, Bakouny Z, Hirsch L, Flippot R, Van Allen EM, Wu CJ, Choueiri TK: Beyond conventional immune-checkpoint inhibition - novel immunotherapies for renal cell carcinoma. NAT REV CLIN ONCOL 2021.
Marabelle A, Tselikas L, de Baere T, Houot R: Intratumoral immunotherapy: using the tumor as the remedy. ANN ONCOL 2017, 28(suppl_12):i33-i43.
Sprooten J, Ceusters J, Coosemans A, Agostinis P, De Vleeschouwer S, Zitvogel L, Kroemer G, Galluzzi L, Garg AD: Trial watch: dendritic cell vaccination for cancer immunotherapy. ONCOIMMUNOLOGY 2019, 8(11):e1638212.
Xie F, Xu M, Lu J, Mao L, Wang S: The role of exosomal PD-L1 in tumor progression and immunotherapy. MOL CANCER 2019, 18(1):146.
Love MI, Huber W, Anders S: Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. GENOME BIOL 2014, 15(12):550.
Edeline J, Mottier S, Vigneau C, Jouan F, Perrin C, Zerrouki S, Fergelot P, Patard JJ, Rioux-Leclercq N: Description of 2 angiogenic phenotypes in clear cell renal cell carcinoma. HUM PATHOL 2012, 43(11):1982-1990.
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. NUCLEIC ACIDS RES 2002, 30(1):207-210.
Gautier L, Cope LBolstad BM, Irizarry RA: affy - analysis of Affymetrix GeneChip data at the probe level. BIOINFORMATICS 2004, 20(3):307-315.
Langfelder P, Horvath S: WGCNA: an R package for weighted correlation network analysis. BMC BIOINFORMATICS 2008, 9:559.
Ivliev AE, Pa TH, Sergeeva MG: Coexpression network analysis identifies transcriptional modules related to proastrocytic differentiation and sprouty signaling in glioma. CANCER RES 2010, 70(24):10060.
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.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z: GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. NUCLEIC ACIDS RES 2017, 45(Web Server issue).
Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM: ONCOMINE: a cancer microarray database and integrated data-mining platform. NEOPLASIA 2004, 6(1):1-6.
Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson Å, Kampf C, Sjöstedt E, Asplund A et al: Proteomics. Tissue-based map of the human proteome. SCIENCE 2015, 347(6220):1260419.
Uhlen M, Zhang C, Lee S, Sjöstedt E, Fagerberg L, Bidkhori G, Benfeitas R, Arif M, Liu Z, Edfors F et al: A pathology atlas of the human cancer transcriptome. SCIENCE 2017, 357(6352).
Thul PJ, Åkesson L, Wiking M, Mahdessian D, Geladaki A, Ait BH, Alm T, Asplund A, Björk L, Breckels LM et al: A subcellular map of the human proteome. SCIENCE 2017, 356(6340).
Hothorn T, Lausen B: On the Exact Distribution of Maximally Selected Rank Statistics. COMPUT STAT DATA AN 2002, 43:121-137.
Therneau TM: survival: Survival Analysis. TECHNOMETRICS 2015, 46(1):111-112.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Muller M: pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC BIOINFORMATICS 2011, 12:77.
Patil I: ggstatsplot: 'ggplot2' Based Plots with Statistical Details. 2018.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 2005, 102(43):15545-15550.
Hänzelmann S, Castelo R, Guinney J: GSVA: gene set variation analysis for microarray and RNA-seq data. BMC BIOINFORMATICS 2013, 14:7.
Binnewies M, Roberts EW, Kersten K, Chan V, Fearon DF, Merad M, Coussens LM, Gabrilovich DI, Ostrand-Rosenberg S, Hedrick CC et al: Understanding the tumor immune microenvironment (TIME) for effective therapy. NAT MED 2018, 24(5):541-550.
Ghatalia P, Gordetsky J, Kuo F, Dulaimi E, Cai KQ, Devarajan K, Bae S, Naik G, Chan TA, Uzzo R et al: Prognostic impact of immune gene expression signature and tumor infiltrating immune cells in localized clear cell renal cell carcinoma. J IMMUNOTHER CANCER 2019, 7(1):139.
Xu WH, Xu Y, Wang J, Wan FN, Wang HK, Cao DL, Shi GH, Qu YY, Zhang HL, Ye DW: Prognostic value and immune infiltration of novel signatures in clear cell renal cell carcinoma microenvironment. Aging (Albany NY) 2019, 11(17):6999-7020.
Atkins MB, Tannir NM: Current and emerging therapies for first-line treatment of metastatic clear cell renal cell carcinoma. CANCER TREAT REV 2018, 70:127-137.
Jala VR, Bodduluri SR, Satpathy SR, Chheda Z, Sharma RK, Haribabu B: The yin and yang of leukotriene B(4) mediated inflammation in cancer. SEMIN IMMUNOL 2017, 33:58-64.
Seo JM, Park S, Kim JH: Leukotriene B4 receptor-2 promotes invasiveness and metastasis of ovarian cancer cells through signal transducer and activator of transcription 3 (STAT3)-dependent up-regulation of matrix metalloproteinase 2. J BIOL CHEM 2012, 287(17):13840-13849.
Kim H, Choi JA, Park GS, Kim JH: BLT2 up-regulates interleukin-8 production and promotes the invasiveness of breast cancer cells. PLOS ONE 2012, 7(11):e49186.
Andrusiewicz M, Szczerba A, Wołuń-Cholewa M, Warchoł W, Nowak-Markwitz E, Gąsiorowska E, Adamska K, Jankowska A: CGB and GNRH1 expression analysis as a method of tumor cells metastatic spread detection in patients with gynecological malignances. J TRANSL MED 2011, 9:130.

Table1.xlsx
Table S1: GO biological processes of genes in brown module.

PDF

Journal Publication

published 06 Jul, 2021

Read the published version in Cancer Cell International →

Editorial decision: Major revision
06 May, 2021
Reviewer #4 agreed at journal
04 May, 2021
Review #3 received at journal
03 May, 2021
Review #2 received at journal
03 May, 2021
Reviewer #3 agreed at journal
02 May, 2021
Reviewer #2 agreed at journal
30 Apr, 2021
Reviewers invited by journal
25 Apr, 2021
Reviewer #1 agreed at journal
25 Apr, 2021
Review #1 received at journal
25 Apr, 2021
Reviews received at journal
25 Apr, 2021
Editor invited by journal
20 Apr, 2021
Editor assigned by journal
19 Apr, 2021
Submission checks completed at journal
19 Apr, 2021
First submitted to journal
13 Apr, 2021

You are reading this latest preprint version

GNRH1 and LTB4R Might Be Novel Immune-Related Prognostic Biomarkers in Clear Cell Renal Cell Carcinoma (ccRCC)

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Materials And Methods

Datasets and Immune-Related Genes (IRGs)

Construction of Co-Expression Network

Prognosis-Related Module Identification

Novel Immune-Related Prognostic Biomarker Identification

Internal Validation of Immune-Related Prognostic Biomarkers

mRNA Expression Level and Translation-Level Expression of Prognostic Biomarker Verification

RNA extraction and qRT-PCR

Western blotting

Prognostic Value of Immune-Related Prognostic Biomarkers Validation

Immune-Related Prognostic Biomarker Function Analysis

Association between Immune-Related Prognostic Biomarker Expression and Immunocytes Exploring

Results

Screening out of Two Prognosis-Related Modules

Function Exploring of Genes in Prognosis-Related Modules

Two Genes Were Regarded as Novel Immune-Related Prognostic Biomarkers

Internal Validation of the Two Immune-Related Prognostic Biomarkers

External Validation of the Two Immune-Related Prognostic Biomarkers

Experimental validation

Identification of Immune-Related Prognostic Biomarker Related Pathways

Association between Immune-Related Prognostic Biomarker Expression and Immune Infiltration Level in ccRCC

Discussion

Conclusions

Abbreviations

Declarations

Data availability statement

Acknowledgements

Funding

Author contributions

Declaration of Conflicting Interests

Declaration

Ethical Statement

References

Supplementary Files

Status:

Journal Publication

Version 1