Abstract
Background: Traditional clinicopathological features (TNM, pathology grade) are often insufficient in predictive prognosis accuracy of clear cell renal cell carcinoma (ccRCC). The IL6-JAK-STAT3 pathway is aberrantly hyperactivated in many cancer types, and such hyperactivation is generally associated with a poor clinical prognosis implying that it can be used as a promising prognosis indicator. The relation between the IL6-JAK-STAT3 pathway and ccRCC remains unknown.
Methods: We evaluated the levels of various cancer hallmarks and filtered out the promising risk hallmarks in ccRCC. Subsequently, a prognosis model based on these hallmark-related genes was established via weighted correlation network analysis and Cox regression analysis. Besides, we constructed a nomogram based on the previous model with traditional clinicopathological features to improve the predictive power and accuracy.
Results: The IL6-JAK-STAT3 pathway was identified as the promising risk hallmarks in ccRCC, and the pathway-related prognosis model based on five genes was built. Also, the nomogram we developed demonstrated the strongest and most stable survival predictive ability.
Conclusion: Our study would provide new insights for guiding individualized treatment of ccRCC patients.
Introduction
Renal cell carcinoma (RCC), 1 of the 10 most common cancers in men and women with an increasing incidence since the 1990s, was expected to cause approximately 76,080 new cases and 13,780 deaths in 20211. Among RCC, clear cell renal cell carcinoma (ccRCC) accounts for almost 75% of detected cases (Hsieh et al., 2017), imposing a significant burden on patients, families, and society. Primary tumor resection is the standard treatment (Subramanian and Haas, 2018) for clinically localized ccRCC. For advanced ccRCC, molecularly targeted drugs and mammalian rapamycin inhibition are the first-line treatments (Powles et al., 2018), but the therapeutic impact is not promising due to the tumor microenvironment being highly dynamic, adaptive, and heterogeneous (Heidegger et al., 2019). Patients greatly suffered as a result of this embarrassing situation. Therefore, in addition to traditional indicators (TNM, pathology grade, etc.), a more stable and robust improved prognosis method is desperately required.
The IL6-JAK-STAT3 pathway is always aberrantly hyperactivated in many cancer types (Wang and Sun, 2014; Zhang et al., 2016; Chiang et al., 2020; Pan et al., 2020), and this phenomenon is often associated with poor outcomes (Kusaba et al., 2005; Macha et al., 2011; Chen et al., 2013). This may be linked to a pathway that promotes tumor growth, affects tumor microenvironment, and induces antitumor immunity (Johnson et al., 2018). Given that immune checkpoint inhibition has significantly enhanced the treatment of patients with advanced cancer, numerous studies (Austin et al., 2014; Atsaves et al., 2017; Bu et al., 2017) have explored the relationship between them. The results demonstrated that signaling via the IL6-JAK-STAT3 pathway induces PD-1 and/or PD-L1 expression. However, no IL6-JAK-STAT3 pathway-related method has been published to date to predict ccRCC prognosis.
This study first identified the IL6-JAK-STAT3 pathway among various cancer hallmarks as a promising risk hallmark in ccRCC. Secondly, a scare-free network was constructed, and the genes significantly correlated with the IL6-JAK-STAT3 pathway in the network were picked out. Subsequently, IL6-JAK-STAT3 pathway-related gene signature was established. Besides, we explored the correlation of alternative polyadenylation of these genes and ccRCC prognosis as well as the IL6-JAK-STAT3 pathway and tumor microenvironment. Finally, a nomogram based on risk gene signature with traditional indicators was developed to improve the predictive power and accuracy. The main flowchart of the study was shown in Figure 1.
FIGURE 1
Materials and Methods
Dataset Preparation and Data Processing
This research included 640 ccRCC tumors at various stages and 72 normal kidney tissues with associated clinical data. RNA-seq data of 539 ccRCC and 72 normal kidney tissue used as training and inside validation cohort were obtained from The Cancer Genome Atlas (TCGA). The microarray dataset GSE40435 was downloaded from NCBI GEO2 and was used for the testing set. Microarray and RNA-seq data included in our study were normalized, and log2 was transformed.
Promising Risk Hallmark Selection and Validation
To evaluate the levels of various hallmarks in ccRCC, TCGA gene expression matrix was transformed into a hallmark enrichment score matrix via “gsva,” an R package that contained a single-sample gene set enrichment analysis (ssGSEA) algorithm (Barbie et al., 2009). The hallmark gene sets were downloaded from the Molecular Signatures Database3. Univariate and multivariate Cox regression analysis was performed to select the promising risk hallmarks, and the hallmark that satisfied previous analysis with the lowest p-value was chosen. To validate that the risk hallmark could be a risk factor, we compared ssGSEA score of promising risk hallmarks between normal kidney tissues and tumor tissues. The previous comparison was also done between alive and dead patients in ccRCC. Moreover, we conducted Gene Set Enrichment Analysis (GSEA) between normal kidney tissues and tumor tissues. Furthermore, GSEA and survival analysis were performed in different groups, grouping by median hallmarks enrichment score.
Genes of Promising Risk Hallmark Further Selection
Based on TCGA gene expression matrix, a scale-free co-expression network was identified, and genes with similar expression patterns were clustered into the same module. Subsequently, the module with the highest correlation to risk hallmarks was identified. To pick out the hub genes from the previous module, gene significance for the module was set as 0.5, and module membership was set as 0.8. Genes that met the previous critical limitation were chosen for further study as promising candidates.
Prognosis Model Establishment and Validation
The promising gene expression matrix was divided into training sets (80%) and validation sets (20%). The least absolute shrinkage and selection operator (LASSO) Cox regression model was performed to identify the most robust prognostic genes from promising genes. Here, the IL6-JAK-STAT3 pathway-related risk score (IRS) was calculated as follows: IRS = Sum (Expression of the most robust prognostic genes∗ their LASSO Cox coefficients). Subsequently, the patients were divided into high-IRS and low-IRS groups based on median IRS. The log-rank test was used to plot the Kaplan–Meier survival curve between these two groups. To determine gene signature accuracy, the receiver operating characteristic (ROC) curve was generated via R package “timeROC.” Besides, dataset GSE 40435 from another platform was used to re-evaluate module accuracy. The IRS of dataset was calculated, and it was split in the same manner as previously done. For inaccessible follow-up information, we performed GSEA in different IRS groups and explored the relationship between IRS and pathology grade.
Further Bioinformatic and Statistical Analyses
Recently, increased focus has been paid to tumor microenvironment and alternative polyadenylation, despite the fact that their functions in ccRCC remain unknown. As a result, we herein estimated the abundance of infiltrating immune cells and predicted tumor purity in different levels of ssGSEA score. To examine the relationship between PD-1 and ssGSEA score, we compared PD-1 expression in high-ssGSEA score and low-ssGSEA score groups. Alternative polyadenylation (APA) is emerging as a new paradigm of post-transcriptional regulation for > 70% of human genes (Elkon et al., 2013). Recently, research has revealed that APA plays a significant role in various cancers (Kang et al., 2020; Zhang et al., 2021). The Cancer 3′ UTR Atlas (Feng et al., 2018) is an online database that, in conjunction with other genomics profiles in cBioPortal, may be used to investigate the functional consequences of APA events in human cancers. We explored APA events of promising genes and their correlation with prognosis in ccRCC.
Established Nomogram of Combined Signature With Clinical Characteristics and APA
To further enhance prognostic accuracy, signature with clinical characteristics (age, gender, grade, stage, and T) and external risk factor APA of FMNL1 were combined for established nomogram. This process was accomplished by R package “rms.” Here, time-dependent receiver operating characteristic analysis was used to measure the predictive power of different risk factors via R package “timeROC,” and the results were compared via R package “tidyverse.”
Results
The IL6-JAK-STAT3 Pathway Was Identified as a Promising Risk Hallmark
After performing univariate analysis and multivariate Cox regression analysis in the hallmark enrichment score matrix, the results are shown in Figure 2A. Four hallmarks were the risk factors, while six hallmarks were on the contrary. Here, the IL6-JAK-STAT3 pathway was identified as the promising risk hallmark with the lowest p-value. Figure 2B reveals that patients with a high Z-score (hallmarks enrichment score) had a worse prognosis than those with a low Z-score. More interestingly, the IL6-JAK-STAT3 pathway score in dead patients was higher when compared with alive ones, as displayed in Figure 2C. This phenomenon was also observed between normal kidney tissues and tumor tissues (Figure 2C and Supplementary Figure 1). The GSEA results also demonstrated that the IL6-JAK-STAT3 pathway was functional in a higher Z-score group (Figure 2D). All these results demonstrate that the IL6-JAK-STAT3 pathway was the suitable promising risk hallmark.
FIGURE 2
A Part of Genes in the Cyan Module Was Selected for Further Analysis
After performing WGCNA, a scale-free co-expression network was identified. The genes were clustered into seven modules with a different color (Figure 3A). Among these modules, the cyan module exhibits the highest correlation score with the IL6-JAK-STAT3 pathway at 0.8. The score of the p-value was 4e–114 (Figure 3B). In addition, module membership vs. gene significance results are shown in Figure 3C, indicating that the cyan module was closely related to the IL6-JAK-STAT3 pathway. A total of 142 genes in the cyan module met the previous limitation and were chosen for further analysis. Subsequently, the least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed and results were shown in Figures 3D,E.
FIGURE 3
Five IL6-JAK-STAT3 Pathway-Related Genes Prognosis Model Establishment
Subsequently, ccRCC samples from TCGA were randomly divided into train cohort and validation cohort via R package “caret” in ratios of 8:2. Using R package “survival and survminer,” the signature was established and the IL6-JAK-STAT3 pathway-related risk score (IRS) was calculated as follows: IRS = (0.11∗ LAIR1) + (–0.20∗ C16orf54) + (0.16∗ FMNL1) + (0.42∗ NOD2) + (–0.97∗ WDFY4). As Figure 4A displays, patients with a high IRS had a worse outcome than those with low IRS in the training cohort. The predictive performance results of IRS are illustrated in Figure 4B, and the area under the curve reached 0.778 at 1 year, 0.711 at 3 years, and 0.729 at 5 years. Consistently, the results of the Kaplan–Meier curve and time-dependent ROC curves in the inside validation cohort indicated a good performance of the signature model (Figures 4E,F). Previously, it was tested on a different platform, cohort GSE 40435. Due to a lack of follow-up data, we compared IRS in different grades, and the results presented in Figure 4C demonstrated that patients had a higher IRS as the grade elevated. The GSEA analysis indicated that the IL6-JAK-STAT3 pathway functioned on a higher IRS group as Figure 4G shown.
FIGURE 4
Univariate and multivariate Cox regression analyses showed that IRS was an independent prognostic predictor with a hazard ratio of 1.158. More intriguingly, pathology grade and tumor stage behave as independent risk factors among various factors while T did not (Figures 4D,H).
APA of Pathway-Related Genes Was Related to Prognosis in ccRCC
The boxplot in Figures 5A,C revealed the difference between immune cells and immune-related functions or pathways in different IRS groups. The score of diverse immune cell subpopulations was higher in high IRS groups except for Mast cells, while immune-related functions or pathways were all more active. Similarly, CD274 expression and immune score were upregulated in higher IRS groups (Figures 5B,D). Figures 6A,D show that APA of FMNL1 could be an independent risk factor, implying that a higher APA of FMNL1 resulted in a worse outcome.
FIGURE 5
FIGURE 6
The Nomogram Improves the Predictive Accuracy of Prognostic
As Figure 6E shows, a nomogram was built with HRS and other clinical characteristics. To evaluate nomogram accuracy, calibration analysis was performed, and the results in Figure 6C suggest high correctness. Compared with other characteristics, the nomogram (Figure 6B) has the strongest and most stable survival predictive ability, and the average AUC is higher than 0.8, much better than traditional characteristics (pathological grade and tumor stage). This means that the nomogram can significantly enhance the ability to predict prognosis.
Discussion
The IL6-JAK-STAT3 pathway had a vast influence on tumors with various biological progress, such as migration, invasion, and angiogenesis (Ni et al., 2020; Pan et al., 2020; Siersbaek et al., 2020; Tung et al., 2020). Therefore, targeting the IL6-JAK-STAT3 pathway was a promising treatment correlated with multifaceted mechanisms. Numerous new IL-6/JAK/STAT3 pathway inhibitors have been identified to date and are currently undergoing preclinical and/or clinical evaluation (Oh et al., 2015; Liu et al., 2016), whereas few studies were reported on IL6-JAK-STAT3 pathway signature in ccRCC. Considering the IL6-JAK-STAT3 pathway as a complex gene network, a scale-free co-expression network was constructed, and genes in the network that significantly correlated with the targeted pathway were used to establish the model, rather than relying on any individually identified “IL6-JAK-STAT3 pathway-related genes” from literature. In addition, the traditional prognostic system was combined with the model to optimize the clinical routine.
This study translated the gene expression matrix into a hallmark enrichment score matrix using ssGSEA algorithm. Subsequently, the IL6-JAK-STAT3 pathway was identified as promising risk hallmarks via univariate and multivariate Cox regression analysis while the performed WGCNA selected the genes closely related to hallmarks. Besides, a LASSO Cox regression model was constructed, and IRS was calculated concurrently. An inside validation cohort and test cohort from other platforms identified the prognosis value of the model. Considering APA and tumor microenvironment are tightly closed with tumor, we preliminary explored its relation with the IL6-JAK-STAT3 pathway. Besides, improving the predictive power and accuracy of prognosis in ccRCC was accomplished by establishing a nomogram including IRS with other available traditional clinicopathological features. Calibration analysis revealed that prediction was extremely close to actual survival. More interesting, the AUC(t) of the nomogram fluctuated around 0.8, much better than other predictive models at different time points during follow-up.
A part of genes we identified was studied in many reports. Mey et al. (2019) revealed that NOD2 expression is upregulated in ccRCC with a mouse xenograft model. Xu et al. (2017) demonstrated that NO2 might be a biomarker for the survival of kidney cancer patients. This phenomenon was also observed in LAIR1 (Jingushi et al., 2019). Besides, Peng et al. (2020) claimed that collagen promotes anti-PD-1/PD-L1 resistance in cancer through LAIR1-dependent CD8+ T cell exhaustion, similar to Xu et al.’s findings (Xu et al., 2020). The reports about FMNL1 extensively focused on its role in aggressiveness of many tumors, such as nasopharyngeal carcinoma, leukemia, and ccRCC (Favaro et al., 2013; Chen et al., 2018; Zhang et al., 2020), while WDFY4 extensively focused on autoimmune diseases (Kochi et al., 2018; Yuan et al., 2018; Zhu et al., 2020). Unlike previous genes, the reports about c16 or f54 were limited and still need further investigation.
Despite the good performance of the nomogram, some limitations in our study remain, requiring further investigation. First, the study was based on public datasets, requiring further verification in our own clinical cohort. Second, as a retrospective study, it needs further validation in a larger prospective cohort. Finally, further experimental studies are needed to explore the roles that five pathway-related genes play in ccRCC.
Conclusion
We established an IL6-JAK-STAT3 pathway-related gene prognosis model to predict the outcomes in ccRCC patients. Based on this model, we combined traditional clinicopathological features in advance and constructed a nomogram with more stable and accurate survival predictive ability. Moreover, this study would provide new insights for guiding individualized treatment of ccRCC patients.
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.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: TCGA and NCBI database.
Author contributions
CZ and ZW conceived and designed the study, collected the data, performed the analysis, and wrote the manuscript. XY, ZL, YR, GW, HZ, MQ, and CX contributed to data and analysis tools. JC and JL helped to revise the manuscript. CS was guidance and advice. All authors contributed to the article and approved the submitted version.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.686907/full#supplementary-material
Footnotes
References
1
AtsavesV.TsesmetzisN.ChioureasD.KisL.LeventakiV.DrakosE.et al (2017). PD-L1 is commonly expressed and transcriptionally regulated by STAT3 and MYC in ALK-negative anaplastic large-cell lymphoma.Leukemia311633–1637. 10.1038/leu.2017.103
2
AustinJ. W.LuP.MajumderP.AhmedR.BossJ. M. (2014). STAT3, STAT4, NFATc1, and CTCF regulate PD-1 through multiple novel regulatory regions in murine T cells.J. Immunol.1924876–4886. 10.4049/jimmunol.1302750
3
BarbieD. A.TamayoP.BoehmJ. S.KimS. Y.MoodyS. E.DunnI. F.et al (2009). Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1.Nature462108–112. 10.1038/nature08460
4
BuL. L.YuG. T.WuL.MaoL.DengW. W.LiuJ. F.et al (2017). STAT3 Induces immunosuppression by upregulating PD-1/PD-L1 in HNSCC.J. Dent. Res.961027–1034. 10.1177/0022034517712435
5
ChenW. H.CaiM. Y.ZhangJ. X.WangF. W.TangL. Q.LiaoY. J.et al (2018). FMNL1 mediates nasopharyngeal carcinoma cell aggressiveness by epigenetically upregulating MTA1.Oncogene376243–6258. 10.1038/s41388-018-0351-8
6
ChenY.WangJ.WangX.LiuX.LiH.LvQ.et al (2013). STAT3, a poor survival predicator, is associated with lymph node metastasis from breast cancer.J. Breast Cancer1640–49. 10.4048/jbc.2013.16.1.40
7
ChiangK. C.ChangK. S.HsuS. Y.SungH. C.FengT. H.ChaoM.et al (2020). Human Heme oxygenase-1 induced by interleukin-6 via JAK/STAT3 pathways is a tumor suppressor gene in hepatoma cells.Antioxidants (Basel)9:251. 10.3390/antiox9030251
8
ElkonR.UgaldeA. P.AgamiR. (2013). Alternative cleavage and polyadenylation: extent, regulation and function.Nat. Rev. Genet.14496–506. 10.1038/nrg3482
9
FavaroP.TrainaF.Machado-NetoJ. A.LazariniM.LopesM. R.PereiraJ. K.et al (2013). FMNL1 promotes proliferation and migration of leukemia cells.J. Leukoc. Biol.94503–512. 10.1189/jlb.0113057
10
FengX.LiL.WagnerE. J.LiW. (2018). TC3A: the cancer 3’.UTR Atlas. Nucleic Acids Res.46D1027–D1030. 10.1093/nar/gkx892
11
HeideggerI.PircherA.PichlerR. (2019). Targeting the tumor microenvironment in renal cell cancer biology and therapy.Front. Oncol.9:490. 10.3389/fonc.2019.00490
12
HsiehJ. J.PurdueM. P.SignorettiS.SwantonC.AlbigesL.SchmidingerM.et al (2017). Renal cell carcinoma.Nat. Rev. Dis. Primers3:17009. 10.1038/nrdp.2017.9
13
JingushiK.UemuraM.NakanoK.HayashiY.WangC.IshizuyaY.et al (2019). Leukocyteassociated immunoglobulinlike receptor 1 promotes tumorigenesis in RCC.Oncol. Rep.411293–1303. 10.3892/or.2018.6875
14
JohnsonD. E.O’KeefeR. A.GrandisJ. R. (2018). Targeting the IL-6/JAK/STAT3 signalling axis in cancer.Nat. Rev. Clin. Oncol.15234–248. 10.1038/nrclinonc.2018.8
15
KangD.LeeY.LeeJ. S. (2020). RNA-binding proteins in cancer: functional and therapeutic perspectives.Cancers (Basel)12:2699. 10.3390/cancers12092699
16
KochiY.KamataniY.KondoY.SuzukiA.KawakamiE.HiwaR.et al (2018). Splicing variant of WDFY4 augments MDA5 signalling and the risk of clinically amyopathic dermatomyositis.Ann. Rheum. Dis.77602–611. 10.1136/annrheumdis-2017-212149
17
KusabaT.NakayamaT.YamazumiK.YakataY.YoshizakiA.NagayasuT.et al (2005). Expression of p-STAT3 in human colorectal adenocarcinoma and adenoma; correlation with clinicopathological factors.J. Clin. Pathol.58833–838. 10.1136/jcp.2004.023416
18
LiuX.JonesG. W.ChoyE. H.JonesS. A. (2016). The biology behind interleukin-6 targeted interventions.Curr. Opin. Rheumatol.28152–160. 10.1097/BOR.0000000000000255
19
MachaM. A.MattaA.KaurJ.ChauhanS. S.ThakarA.ShuklaN. K.et al (2011). Prognostic significance of nuclear pSTAT3 in oral cancer.Head Neck33482–489. 10.1002/hed.21468
20
MeyL.JungM.RoosF.BlahetaR.HegeleA.KinscherfR.et al (2019). NOD1 and NOD2 of the innate immune system is differently expressed in human clear cell renal cell carcinoma, corresponding healthy renal tissue, its vasculature and primary isolated renal tubular epithelial cells.J. Cancer Res. Clin. Oncol.1451405–1416. 10.1007/s00432-019-02901-7
21
NiJ. S.ZhengH.OuY. L.TaoY. P.WangZ. G.SongL. H.et al (2020). miR-515-5p suppresses HCC migration and invasion via targeting IL6/JAK/STAT3 pathway.Surg. Oncol.34113–120. 10.1016/j.suronc.2020.03.003
22
OhD. Y.LeeS. H.HanS. W.KimM. J.KimT. M.KimT. Y.et al (2015). Phase i study of OPB-31121, an oral STAT3 inhibitor, in patients with advanced solid tumors.Cancer Res. Treat.47607–615. 10.4143/crt.2014.249
23
PanM. S.WangH.AnsariK. H.LiX. P.SunW.FanY. Z. (2020). Gallbladder cancer-associated fibroblasts promote vasculogenic mimicry formation and tumor growth in gallbladder cancer via upregulating the expression of NOX4, a poor prognosis factor, through IL-6-JAK-STAT3 signal pathway.J. Exp. Clin. Cancer Res.39:234. 10.1186/s13046-020-01742-4
24
PengD. H.RodriguezB. L.DiaoL.ChenL.WangJ.ByersL. A.et al (2020). Collagen promotes anti-PD-1/PD-L1 resistance in cancer through LAIR1-dependent CD8(+) T cell exhaustion.Nat. Commun.11:4520. 10.1038/s41467-020-18298-8
25
PowlesT.NecchiA.RosenG.HariharanS.ApoloA. B. (2018). Anti-programmed cell death 1/Ligand 1 (PD-1/PD-L1) antibodies for the treatment of urothelial carcinoma: state of the art and future development.Clin. Genitourin. Cancer.16117–129. 10.1016/j.clgc.2017.11.002
26
SiersbaekR.ScabiaV.NagarajanS.ChernukhinI.PapachristouE. K.BroomeR.et al (2020). IL6/STAT3 signaling hijacks estrogen receptor alpha enhancers to drive breast cancer metastasis.Cancer Cell38:e419. 10.1016/j.ccell.2020.06.007
27
SubramanianP.HaasN. B. (2018). Recent advances in localized RCC: a focus on VEGF and immuno-oncology therapies.Urol. Oncol.3623–30. 10.1016/j.urolonc.2017.09.008
28
TungK. L.WuY. T.LiuC.LinS. C.WuC. H.WuS. Y.et al (2020). EBV Rta-induced IL-6 promotes migration of bystander tumor cells through IL-6R/JAK/STAT3 pathway in vitro.Anticancer Res.403255–3264.
29
WangS. W.SunY. M. (2014). The IL-6/JAK/STAT3 pathway: potential therapeutic strategies in treating colorectal cancer (Review).Int. J. Oncol.441032–1040. 10.3892/ijo.2014.2259
30
XuD.ZhangS.ZhangS.LiuH.YuanT. J. O. (2017). NOD2 maybe a biomarker for the survival of kidney cancer patients.Oncotarget8101489–101499.
31
XuL.WangS.LiJ.LiJ.LiB. (2020). Cancer immunotherapy based on blocking immune suppression mediated by an immune modulator LAIR-1.Oncoimmunology9:1740477. 10.1080/2162402X.2020.1740477
32
YuanQ.LiY.LiJ.BianX.LongF.DuanR.et al (2018). WDFY4 is involved in symptoms of systemic lupus erythematosus by modulating B cell fate via noncanonical autophagy.J. Immunol.2012570–2578. 10.4049/jimmunol.1800399
33
ZhangM. F.LiQ. L.YangY. F.CaoY.ZhangC. Z. (2020). FMNL1 exhibits pro-metastatic activity via Cxcr2 in clear cell renal cell carcinoma.Front. Oncol.10:564614. 10.3389/fonc.2020.564614
34
ZhangN.ZengY.DuW.ZhuJ.ShenD.LiuZ.et al (2016). The EGFR pathway is involved in the regulation of PD-L1 expression via the IL-6/JAK/STAT3 signaling pathway in EGFR-mutated non-small cell lung cancer.Int. J. Oncol.491360–1368. 10.3892/ijo.2016.3632
35
ZhangY.LiuL.QiuQ.ZhouQ.DingJ.LuY.et al (2021). Alternative polyadenylation: methods, mechanism, function, and role in cancer.J Exp. Clin. Cancer Res.40:51. 10.1186/s13046-021-01852-7
36
ZhuJ.HaoS.ZhangX.QiuJ.XuanQ.YeL. (2020). Integrated bioinformatics analysis exhibits pivotal exercise-induced genes and corresponding pathways in malignant melanoma.Front. Genet.11:637320. 10.3389/fgene.2020.637320
Summary
Keywords
prognosis, IL6-JAK-STAT3, nomogram, WGCNA, renal cell carcinoma, immune
Citation
Zhan C, Xu C, Chen J, Shen C, Li J, Wang Z, Ying X, Luo Z, Ren Y, Wu G, Zhang H and Qian M (2021) Development and Validation of an IL6/JAK/STAT3-Related Gene Signature to Predict Overall Survival in Clear Cell Renal Cell Carcinoma. Front. Cell Dev. Biol. 9:686907. doi: 10.3389/fcell.2021.686907
Received
28 March 2021
Accepted
10 May 2021
Published
29 September 2021
Volume
9 - 2021
Edited by
Cathy Tournier, The University of Manchester, United Kingdom
Reviewed by
Iram Fatima, University of Nebraska Medical Center, United States; Xiaobai Li, Zhejiang Academy of Agricultural Sciences, China
Updates
Copyright
© 2021 Zhan, Xu, Chen, Shen, Li, Wang, Ying, Luo, Ren, Wu, Zhang and Qian.
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: Chong Shen, shclmy@163.com
†These authors have contributed equally to this work
This article was submitted to Signaling, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
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.