A 7-gene signature predicts the prognosis of patients with bladder cancer

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

Download PDF

Research Article

A 7-gene signature predicts the prognosis of patients with bladder cancer

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

This work is licensed under a CC BY 4.0 License

Journal Publication

published 28 Jan, 2022

Read the published version in BMC Urology →

You are reading this latest preprint version

Background: The biomarkers have an important guiding role in prognosis and treatment of patients with bladder cancer(BC). The aim of the present study was to identify and evaluate a prognostic gene signature in BC patients.

Methods: The gene expression profiles of BC samples and the corresponding clinicopathological data were downloaded from GEO and TCGA. The differentially expressed genes (DEGs) were identified by R software. Univariate Cox regression and the least absolute shrinkage and selection operator(LASSO) Cox regression were applied to construct the prognostic score model. A nomogram was established with the identified prognostic factors to predict the overall survival rates of BC patients. The discriminatory and predictive capacity of the nomogram was evaluated based on the concordance index (C-index), calibration curves and decision curve analysis (DCA).

Results: A 7-gene signature(KLRB1, PLAC9, SETBP1, NR2F1, GRHL2, ANXA1 and APOL1) was identified from 285 DEGs by univariate and LASSO Cox regression analyses. Univariate and multivariate Cox regression analyses showed that age, lymphovascular invasion, lymphatic metastasis, metastasis and the 7-gene signature risk score was an independent predictor of BC patient prognosis. A nomogram that integrated these independent prognostic factors was constructed. The C-index (0.73, CI 95%, 0.693–0.767) and calibration curve demonstrated the good performance of the nomogram. DCA of the nomogram further showed that this model exhibited good net benefit.

Conclusions: The combined 7-gene signature could serve as a biomarker for predicting BC prognosis. The nomogram built by risk score and other clinical factors could be an effective tool for predicting the prognosis of patients with BC.

Urology & Nephrology

bladder cancer

TCGA

GEO

prognosis

nomograms

Bladder cancer (BC) is a common tumour of the urinary system. In recent years, its incidence has shown an increasing trend, and BC patients have high recurrence and morbidity rates ^1,2. It is estimated that there were more than 500,000 new cases of bladder cancer and approximately 200,000 BC-related deaths worldwide in 2018 ³. According to the invasion of the muscle layer, BC can be divided into nonmuscle-invasive BC and muscle-invasive BC. Nonmuscle-invasive BC accounts for approximately 75% of bladder cancers. The 5-year survival rate of nonmuscle-invasive BC patients is high, but these patients are prone to relapse ^4,5. Muscle-invasive BC accounts for approximately 25% of BCs. Among these muscle-invasive BC patients, the degree of progression is higher and the rate of tumour metastasis is greater than that for nonmuscle-invasive BC; in addition, the 5-year survival rate is only 5%, and a more thorough treatment strategy is needed ^{6, 7}. Therefore, it is of great clinical significance to investigate factors involved in the development and prognosis of BC. At present, the TNM staging system can predict the prognosis of patients and suggest treatment strategies. However, TNM staging is determined via clinical pathology, and it does not take into account the biological heterogeneity of tumours. Therefore, a more accurate and effective prediction model is needed to improve the risk prediction of patients, which could be beneficial for guiding individualized treatment.

In recent years, with the rapid development of tumour molecular biology, gene biomarkers have gradually gained attention. With the development of new technology and directions for screening and diagnosis of tumours, the search for tumour signatures with high sensitivity and specificity has been increasingly anticipated. However, the current gene biomarkers are mostly single gene biomarkers with limited accuracy, and the tumour prediction models developed from multiple signatures may provide more reliable and accurate evaluation than single clinical signatures ⁸. In this study, the differentially expressed genes screened from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database were applied to establish a prognostic evaluation model for BC patients. This model integrates a gene signature and clinicopathological risk factors to provide an essential reference for prognostic evaluation and treatment in BC patients.

2.1 Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases

Gene expression profiles GSE7476 and GSE37815 were obtained from the GEO database (http://www.ncbi.nlm.nih.gov/geo/). GSE7476 contained 9 BC and 3 normal tissue samples. GSE37815 had 18 BC and 5 normal tissue samples. In addition, RNA-seq data and clinicopathologic characteristic data of BC patients and normal controls were obtained from the TCGA database (http://cancergenome.nih.gov/), including 414 BC and 19 normal tissue samples. The clinicopathologic characteristics and follow-up data included the patient's age, sex, diagnosis subtype, lymphovascular invasion, histologic grade, lymphatic metastasis, metastasis, and survival state and time. The exclusion criteria were as follows: (a) other cancer history; (b) lack of gene sequencing data or lack of clinical data; (c) missing important clinicopathologic characteristic data.\

2.2 Data preprocessing and screening of differentially expressed genes

The Limma package in R software (version 3.5.3) was used to screen the differentially expressed genes (DEGs) in the gene expression profiles (GSE7476 and GSE37815). The edgeR package in R software was used to analyse the DEGs in the BC RNA-seq data from the TCGA. The threshold for DEGs was set to | log2 (fold change) | > 1, and an adjusted P value < 0.05 was considered statistically significant. These DEG results were combined and plotted by the Venn Diagram package in R software.

2.3 Construction and assessment of the gene signature risk score model

The correlation between gene expression level and overall survival (OS) of BC patients was investigated by univariate Cox regression analysis using the survival package in R software, and P < 0.05 was considered statistically significant. Further, the above results were analysed by the least absolute shrinkage and selection operator (LASSO) and multivariate Cox methods using the glmmet package in R software. A gene signature risk score model for BC patients was constructed based on LASSO and multivariate Cox regression analysis results and was used to calculate the risk score for each BC patient. The risk score, which was calculated from the prognostic score model, divided all the patients into high- and low-risk groups by the median value. The Kaplan-Meier survival curve was performed for two groups by the survival package in R software. A time-dependent receiver operating characteristic (ROC) curve was performed for the risk score by the survivalROC package in R software, and the area under the curve (AUC) was calculated to judge the prediction results. The higher the AUC value is, the better the predictive performance of the model. Generally, AUC > 0.9 indicates a high accuracy, 0.7 < AUC ≤ 0.9 indicates a moderate accuracy, 0.5 < AUC ≤ 0.7 indicates a low accuracy, and AUC ≤ 0.5 indicates a lack of discriminating power.

2.4 Construction and assessment of the predictive nomogram

Univariate Cox and multivariate Cox regression analyses were performed using the survival package in R software. A nomogram was established based on the above results to predict the 1-, 3-, 5-, and 8-year OS of BC patients by using the R software rms package. The calibration plots were used to evaluate the accuracy of the nomogram. The 45° line represented the best prediction value. The closer the curve was, the more ideal the result was. Decision curve analysis (DCA) was used to test the clinical value of this model ⁹.

3.1 Screening of differentially expressed genes in bladder cancer

The R difference analysis results showed that 314 genes were upregulated and 859 genes were downregulated in GSE7476, 240 genes were upregulated and 552 genes were downregulated in GSE37815, and 2769 genes were upregulated and 1963 genes were downregulated in the TCGA dataset. A total of 68 genes were upregulated and 217 genes were downregulated after merging these three results (Fig. 1).

3.2 Identification of the 7-gene signature for survival prediction

A total of 126 prognostic-related genes were screened by univariate Cox regression analysis, and 7 key genes, including grainyhead-like 2 (GRHL2), annexin A1 (ANXA1), apolipoprotein L1 (APOL1), SET binding protein 1 (SETBP1), nuclear receptor subfamily 2 group F member 1 (NR2F1), killer cell lectin-like receptor subfamily B (KLRB1) and placenta-specific 9 (PLAC9), were identified using LASSO and multivariate Cox analysis (Fig. 2). The gene expression levels of PLAC9, SETBP1, NR2F1 and ANXA1 were positively correlated with the OS of BC patients, while KLRB1, GRHL2 and APOL1 were negatively correlated with the OS of BC patients. Based on the results of LASSO and multivariate Cox regression analysis, the risk score model was constructed. The risk score formula for this model was: (-0.183 × expression value of KLRB1) + (0.106 × expression value of PLAC9) + (0.117 × expression value of SETBP1 + (0.115 × expression value of NR2F1) + (-0.104 × expression value of GRHL2) + (0.220 × expression value of ANXA1) + (-0.139 × expression value of APOL1). The 7-gene signature prognostic model was used to calculate the prognostic risk score of each patient. BC patients were divided into high-risk and low-risk groups according to the median score. The distribution of risk score, survival status and expression pattern in high-risk and low-risk groups are shown in Fig. 3. Mortality was significantly higher in the high-risk group than in the low-risk group. The predictive performance of the 7-gene signature prognostic model was assessed by using Kaplan-Meier curves and ROC curve AUC values. The Kaplan-Meier curve results showed that the prognosis of the two groups was significantly different (p < 0.001) and that the prognosis of the high-risk group was significantly lower than that of the low-risk group (Fig. 4A). The AUC values for this model (0.711, 0.714, 0.711, and 0.761 for the 1-, 3‐, 5-, and 8‐year OS) demonstrated that the 7-gene signature prognostic model had good performance in predicting prognosis in patients with BC (Fig. 4B).

3.3 Establishment and evaluation of the nomogram

To assess whether the risk scoring model consisting of a 7-gene signature was independent of other clinical factors, univariate Cox and multivariate Cox regression analyses were performed according to clinicopathologic data and risk scores of BC patients. Univariate Cox and multivariate Cox results showed that age, diagnosis subtype, lymphovascular invasion, lymphatic metastasis, metastasis, and risk scores in BC patients were powerful and independent prognostic factors (Table 1). The nomogram was a statistical tool that could predict the overall probability of OS for BC patients with specific outcomes. Based on the results of multivariate Cox regression analysis, a nomogram was established to predict 1-, 3-, 5-, and 8-year OS (Fig. 5). The C-index of the nomogram for OS was 0.73 (CI 95%, 0.693–0.767). The calibration curve and DCA results showed that the nomogram predictions for 1-, 3-, 5-, and 8-year OS of BC patients had great accuracy and net benefit (Fig. 6, Fig. 7).

Table 1

Univariate and Multivariate cox regression analysis of Overall survival in the bladder cancer patients
Variables	N(%)	Univariate cox regression analysis		Multivariate cox regression analysis
Variables	N(%)	HR(95%CI)	P	HR (95%CI)	P
Sex
Female	106(26.6%)	Reference		Reference
Male	292(73.4%)	0.853(0.617–1.181)	0.340	Reference
Age	68.08 ± 10.653	1.034(1.019–1.050)	< 0.001	1.032(1.015–1.048)	< 0.001
Diagnosis subtype
Non-Papillary	270(67.8%)	Reference		Reference
Papillary	128(32.2%)	0.682(0.480–0.968)	0.032	0.864(0.603–1.237)	0.424
History of neoadjuvant treatment
N0	388(97.5%)	Reference		Reference
Yes	10(2.5%)	1.140(0.422–3.077)	0.796	Reference
Lymphovascular invasion
No	130(32.7%)	Reference		Reference
Yes	143(35.9%)	2.195(1.512–3.187)		1.714(1.120–2.624)	0.013
Unkown	125(31.4%)	1.629(1.094–2.426)		1.630(1.085–2.449)	0.019
Histologic grade
High Grade	378(95.0%)	Reference		Reference
Low Grade	20(5.0%)	2.890(0.715–11.683)	0.137	Reference
Lymphatic metastasis
No	130(32.7%)	Reference		Reference
Yes	143(35.9%)	2.293(1.670–3.149)	< 0.001	1.544(1.062–2.244)	0.023
Unkown	125(31.4%)	1.692(1.041–2.751)	0.034	1.937(1.179–3.183)	0.009
Metastasis
No	190(47.7%)	Reference		Reference
Yes	11(2.8%)	3.299(1.582–6.881)	0.001	1.653(0.771–3.545)	0.196
Unkown	197(49.5%)	1.551(1.140–2.111)	0.005	1.441(1.048–1.980)	0.024
Risk
Low	199(50.0%)	Reference		Reference
High	199(50.0%)	2.617(1.909–3.588)	< 0.001	2.180(1.567–3.033)	< 0.001

BC is one of the most common tumours in the urinary system. The global incidence rate of BC is increasing every year. It is one of the 10 most common cancers in men ¹⁰ and is a non-negligible hazard to human health. BC prognosis-related biomarkers have an important guiding role in diagnosis, prognosis and corresponding treatment methods. Therefore, it is one of the hot spots to find effective BC prognostic biomarkers. This study aimed to screen differentially expressed genes and construct a prognostic evaluation model of BC patients. The DEGs were obtained by R analysis of the GEO and TCGA databases in our study. Univariate Cox analysis was used for preliminary screening, and the results showed that 126 genes had significant differences that were related to the prognosis of BC. Seven genes (GRHL2, ANXA1, APOL1, SETBP1, NR2F1, KLRB1 and PLAC9) were identified to be independent prognostic factors for the prognosis of BC patients via LASSO and multivariate Cox regression analyses. These genes were used to establish a prognostic risk score model for BC patients to evaluate the OS of BC patients. The Kaplan-Meier curve and ROC curve demonstrated that the 7-gene signature prognostic model had a good prediction performance. The 7-gene signature risk score was confirmed to be independent of other clinicopathological features of BC patients by univariate Cox and multivariate Cox regression analyses. Further, the risk score and clinicopathological risk factors (age, lymphovascular invasion, lymphatic metastasis, metastasis) were used to establish a prognostic nomogram to predict the prognosis of patients at 1, 3, 5 and 8 years. The results of the calibration curve and DCA confirmed the reliability of the nomogram, which showed that the nomogram could help clinicians predict patient prognosis and provide guidance for patient rehabilitation assessment and treatment decision-making.

At present, although many genes have been identified as effective indicators for predicting the prognosis of BC patients ^11–14, their predictive performance is relatively poor due to the limited number of screened genes as well as the lack of correlation with clinicopathological data. Compared with a single clinical gene biomarker, the integration of multiple related genes into a prognostic evaluation model can effectively improve the accuracy of prediction ^15–17. Similarly, the prognostic evaluation model of the 7-gene signature has displayed satisfactory predictive performance in this study. Among the 7 genes in the signature, PLAC9, SETBP1, ANXA1 and NR2F1 were negatively correlated with OS in BC patients. PLAC9 is a placental protein that is upregulated in embryonic expression (weeks 8 to 9) ¹⁸. Ouyang C et al. reported that overexpression of PLAC9 induces G2 and M phase arrest in cell division and inhibits cell growth ¹⁹, but it has not been reported in bladder cancer. It might be necessary to further explore the specific role of PLAC9 in bladder cancer in the future. SETBP1 has been established as an important diagnostic marker for bone marrow malignancies, and its expression is closely related to the prognosis of patients. Studies have shown that it plays a key role in tumour invasion and rapid evolution ²⁰. Additionally, SETBP1 inhibits the expression of tumour protein phosphatase 2A (PP2A) ²¹ and regulates cell proliferation ²², suggesting that SETBP1 may also play a regulatory role in tumorigenesis and progression. ANXA1 functions to regulate a variety of cell biological behaviours, including membrane aggregation, phagocytosis, cell proliferation and tumorigenesis ²³. Reports have shown that ANXA1 is related to the recurrence and drug resistance of BC, and its expression level was positively correlated with T phase ²⁴. ANXA1 was shown to play an important role in the development of high-grade BC ²⁵. Our findings also demonstrated that high expression of ANXA1 is related to the occurrence and development of bladder cancer and affects the prognosis of BC patients. NR2F1 regulates the progression of cell differentiation, cancer progression, and central and peripheral neurogenesis ²⁶. Kikuchi-Koike R et al. indicated that NR2F1 is highly connected with the high-risk recurrence of breast cancer ²⁷. High NR2F1 expression has a strong effect on increasing lung metastatic potential. Additionally, the overexpression of NR2F1 was reported to enhance invasion and metastasis in salivary adenoid cystic carcinoma ²⁸. At present, the role of NR2F1 in BC has not been reported. In this study, we showed for the first time that NR2F1 is highly expressed in BC and is closely related to the prognosis of BC patients. However, the specific regulation of NR2F1 in BC still needs further exploration. The other three genes, GRHL2, APOL1 and KLRB1, were positively correlated with the OS of BC patients. GRHL2 enhances the effect of carcinogenesis and promotes tumour cell proliferation. It can also promote the growth of oral squamous cell carcinoma, colorectal cancer and hepatocellular carcinoma ²⁹. However, more reports have shown that GRHL2 can inhibit epithelial-mesenchymal transition to inhibit cell migration ³⁰ and cancer metastasis ³¹. Additionally, the overexpression of GRHL2 can effectively inhibit tumour invasion and migration in gastric and breast cancer ^{30, 32}. In this study, we confirmed that GRHL2 is significantly expressed at low levels in BC and has a positive regulatory effect on the prognosis of BC patients. GRHL2 may be an effective prognostic marker for BC. APOL1 is mainly active in the kidney ³³, and its high expression can promote the occurrence of kidney disease ³⁴. In the study of Gutiérrez OM et al, the mortality rate of patients with a high-risk APOL1 genotype was significantly lower than that of patients with low-risk genotypes in the prognostic review cohort of chronic kidney disease ³⁵. However, there are few reports about the effect of APOL1 on BC, and further research is still needed. KLRB1 is involved in the coding of NKRP1A/LLT1, which triggers the activation of T cells and B cells and can be used in the treatment of cancer ³⁶, and its interaction with LLT1 can inhibit natural killing (NK) cell cytotoxicity ³⁷. Many studies have shown that a high level of LLT1 expression can significantly restrain NK cells and enhance tumour immune escape, as knockdown of LLT1 enhances NK cell-mediated lysis of tumour tissues ^38–40. These results combined with our results demonstrate that KLRB1 is expressed at low levels in BC and is positively correlated with OS in BC patients. Because KLRB1 is capable of eliciting spontaneous antitumour immune responses, it may be a promising potential target for BC immunotherapy. Further experiments are needed to verify our hypothesis.

To better apply this model for predicting the OS of BC patients, the risk score was combined with other independent prognostic factors (age, lymphovascular invasion, lymphatic metastasis, metastasis) to establish a nomogram to predict the OS of BC patients at 1 year, 3 years, 5 years, and 8 years. The C-index, calibration curve and DCA results showed the good performance and clinically usefulness of the nomogram. These results showed that the nomogram established by combining these independent prognostic factors was an effective and accurate tool for predicting the OS of BC patients. Additionally, this nomogram played an important role in both short-term and long-term observation studies, so it might be helpful for patient consultations, clinical decision-making and determining follow-up arrangements for BC patients.

However, our research still has some limitations. First, most of the patients in the TCGA database were white or black. The racial diversity was limited, so the results may be biased towards these groups. Second, most of the selected genes were not previously reported in BC or lacked research in BC; therefore, further experimental validation was the direction of our research. Finally, the nomogram for predicting the OS of BC patients, which was based on the TCGA database, was only established and validated in a single data set, lacking further validation with external BC data. Therefore, corrections according to future verification studies in multiple different data sets are required.

In summary, using univariate Cox, LASSO and multivariate Cox analyses, we identified seven key genes (GRHL2, ANXA1, APOL1, SETBP1, NR2F1, KLRB1, and PLAC9) and constructed a 7-gene signature prognostic model for BC patients. We also combined the 7-gene signature risk score with age, lymphovascular invasion, lymphatic metastasis and metastasis to construct a nomogram for prognostic evaluation. This nomogram could provide an important reference for prognostic evaluation and treatment decision-making for BC patients. However, the 7-gene signature in the BC patient prognosis assessment model still needs more clinical validation before it can truly be applied in clinical practice.

BC: bladder cancer

GEO: Gene Expression Omnibus

TCGA: The Cancer Genome Atlas

DEGs: differentially expressed genes

LASSO: least absolute shrinkage and selection operator

ROC: receiver operating characteristic

C‐index: concordance index

DCA: decision curve analysis

OS: overall survival

AUC: area under the curve

NK cell: natural killing cell

Ethics approval and consent to participate

In the present study, the data were downloaded from the GEO and TCGA databases; therefore, no studies with human participants or animals were performed by any of the authors. Thus, no ethical approval or patient consent was required.

Consent for publication

Not applicable.

Availability of data and materials

The datasets of gene expression profiles for bladder cancer (GSE7476 and GSE37815) are available in the GEO on the NCBI website (http://www.ncbi.nlm.nih.gov/geo). The mRNA-Seq data and clinical follow-up data associated with the BC and para-bladder cancer tissues samples were downloaded from the TCGA (https://cancergenome.nih.gov/).

Conflicts of interest

All authors declare that they have no conflict of interest to state.

Funding

The present study was supported by National Key Researchand Development Program of China (ontract No. 2018YFA0902801), National Natural Science Foundation of China (contract No. 81803576), Guangdong Basic and Applied Basic Research Foundation (ontract No. 2020A1515010152), Public health research project in Futian District, Shenzhen (grant no. FTWS2020027), and the research start-up fee for the Eighth Affiliated Hospital, Sun Yat-sen University (grant no. zdbykyqdf005).

Authors’ Contributions

HZ and TF conceived this study; TF, LZ and LY conducted the work and analyzed the data; LH, HC and LZ collected and collated the data; HZ, TF, LZ and LY wrote manuscript. All the authors listed have read and approved the manuscript.

Acknowledgments

Not applicable.

Lerner SP, Bajorin DF, Dinney CP, et al. Summary and Recommendations from the National Carcinoma Institute's Clinical Trials Planning Meeting on Novel Therapeutics for Non-Muscle Invasive Bladder Carcinoma. Bladder Carcinoma. 2016;2(2):165-202.
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global carcinoma statistics, 2012. CA Carcinoma J Clin. 2015;65(2):87-108.
Bray F, Ferlay J, Soerjomataram I, et al: Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394- 424.
Cheng L, Davison DD, Adams J, et al. Biomarkers in bladder carcinoma: translational and clinical implications. Crit Rev Oncol Hematol. 2014;89(1):73-111.
Elsen S, Lerut E, Van Cleynenbreugel B, van der Aa F, van Poppel H, de Witte PA. Biodistribution of Evans blue in an orthotopic AY-27 rat bladder urothelial cell carcinoma model: implication for the improved diagnosis of non-muscle-invasive bladder carcinoma (NMIBC) using dye-guided white-light cystoscopy. BJU Int. 2015;116(3):468-77.
Christensen E, Birkenkamp-Demtröder K, Nordentoft I, et al. Liquid Biopsy Analysis of FGFR3 and PIK3CA Hotspot Mutations for Disease Surveillance in Bladder Carcinoma. Eur Urol. 2017;71(6):961-969.
Xylinas E, Kluth LA, Rieken M, Karakiewicz PI, Lotan Y, Shariat SF. Urine markers for detection and surveillance of bladder carcinoma. Urol Oncol. 2014;32(3):222-9.
Qian Z, Li Y, Fan X, et al. Prognostic value of a microRNA signature as a novel biomarker in patients with lower‐grade gliomas. J Neurooncol. 2017;137:127‐137.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Makin 2006;26(6):565-74.
Chen W, Zheng R, Baade PD, et al. Carcinoma statistics in China, 2015. CA Carcinoma J Clin. 2016;66(2):115-32.
Li S, Liu X, Liu T, et al. Identification of Biomarkers Correlated with the TNM Staging and Overall Survival of Patients with Bladder Carcinoma. Front Physiol. 2017;8,947.
Yan X, Guo Z, Liu X, et al. Four novel biomarkers for bladder carcinoma identified by weighted gene coexpression network analysis. J Cell Physiol. 2019;234(10), 19073-19087.
Yan L, Li S, Xu C, et al. Target protein for Xklp2 (TPX2), a microtubule-related protein, contributes to malignant phenotype in bladder carcinoma. Tumour Biol. 2013;34(6):4089-100.
Bi Y, Chen S, Jiang J, et al. CDCA8 expression and its clinical relevance in patients with bladder carcinoma. Medicine (Baltimore). 2018;97(34):e11899.
Cheng P. A prognostic 3-long noncoding RNA signature for patients with gastric carcinoma. J Cell Biochem. 2018;119(11):9261-9269.
Lai J, Wang H, Pan Z, Su F. A novel six-microRNA-based model to improve prognosis prediction of breast carcinoma. Aging (Albany NY). 2019;11(2):649-662.
Li X, Shi Y, Yin Z, Xue X, Zhou B. An eight-miRNA signature as a potential biomarker for predicting survival in lung adenocarcinoma. J Transl Med. 2014;12,159.
Galaviz-Hernandez C, Stagg C, de Ridder G, et al. Plac8 and Plac9, novel placental-enriched genes identified through microarray analysis. Gene. 2003;309(2):81-9.
Ouyang C, Pu YZ, Qin XH, et al. Placenta-specific 9, a putative secretory protein, induces G2/M arrest and inhibits the proliferation of human embryonic hepatic cells. Biosci Rep. 2018;38(6).
Winkelmann N, Schäfer V, Rinke J, et al. Only SETBP1 hotspot mutations are associated with refractory disease in myeloid malignancies. J Cancer Res Clin Oncol. 2017;143(12):2511-2519.
Ammatuna E, Eefting M, van Lom K, Kavelaars FG, Valk PJ, Touw IP. Atypical chronic myeloid leukemia with concomitant CSF3R T618I and SETBP1 mutations unresponsive to the JAK inhibitor ruxolitinib. Ann Hematol. 2015;94(5):879-80.
Cristóbal I, Blanco FJ, Garcia-Orti L, et al. SETBP1 overexpression is a novel leukemogenic mechanism that predicts adverse outcome in elderly patients with acute myeloid leukemia. Blood. 2010;115(3):615-25.
Guo C, Liu S, Sun MZ. Potential role of Anxa1 in carcinoma. Future Oncol. 2013;9(11):1773-93.
Yu S, Meng Q, Hu H, Zhang M. Correlation of ANXA1 expression with drug resistance and relapse in bladder carcinoma. Int J Clin Exp Pathol. 2014;7(9):5538-48.
Li CF, Shen KH, Huang LC, Huang HY, Wang YH, Wu TF. Annexin-I overexpression is associated with tumour progression and independently predicts inferior disease-specific and metastasis-free survival in urinary bladder urothelial carcinoma. Pathology. 2010;42(1):43-9.
Sosa MS, Parikh F, Maia AG, et al. NR2F1 controls tumour cell dormancy via SOX9- and RARβ-driven quiescence programmes. Nat Commun. 2015;6,6170.
Kikuchi-Koike R, Nagasaka K, Tsuda H, et al. Array comparative genomic hybridization analysis discloses chromosome copy number alterations as indicators of patient outcome in lymph node-negative breast cancer. BMC Cancer. 2019;19(1):521.
Gao XL, Zheng M, Wang HF, et al. NR2F1 contributes to cancer cell dormancy, invasion and metastasis of salivary adenoid cystic carcinoma by activating CXCL12/CXCR4 pathway. BMC Cancer. 2019;19(1):743.
Ma L, Yan H, Zhao H, Sun J. Grainyhead-like 2 in development and cancer. Tumour Biol. 2017;39(5):1010428317698375.
Cieply B, Riley P, Pifer PM, et al. Suppression of the epithelial-mesenchymal transition by Grainyhead-like-2. Cancer Res. 2012;72(9):2440-53.
Cieply B, Farris J, Denvir J, Ford HL, Frisch SM. Epithelial-mesenchymal transition and tumor suppression are controlled by a reciprocal feedback loop between ZEB1 and Grainyhead-like-2. Cancer Res. 2013;73(20):6299-309.
Xiang J, Fu X, Ran W, Wang Z. Grhl2 reduces invasion and migration through inhibition of TGFβ-induced EMT in gastric cancer. Oncogenesis. 2017;6(1):e284.
Dummer PD, Limou S, Rosenberg AZ, Heymann, J, Nelson G, Winkler CA, et al. APOL1 Kidney Disease Risk Variants: An Evolving Landscape. Semin Nephrol. 2015;35(3):222-36.
Fu Y, Zhu JY, Richman A, et al. APOL1-G1 in Nephrocytes Induces Hypertrophy and Accelerates Cell Death. J Am Soc Nephrol. 2017;28(4):1106-1116.
Gutiérrez OM, Irvin MR, Zakai NA, et al. APOL1 Nephropathy Risk Alleles and Mortality in African American Adults: A Cohort Study. Am J Kidney Dis. 2020;75(1), 54-60.
Bialoszewska A, Malejczyk J. Biological and Clinical Significance of Human NKRP1A/LLT1 Receptor/Ligand Interactions. Crit Rev Immunol. 2018;38(6):479-489.
Lin YL, Lin SC. Analysis of the CD161-expressing cell quantities and CD161 expression levels in peripheral blood natural killer and T cells of systemic lupus erythematosus patients. Clin Exp Med. 2017;17(1):101-109.
Marrufo AM, Mathew SO, Chaudhary P, Malaer JD, Vishwanatha JK, Mathew PA. Blocking LLT1 (CLEC2D, OCIL)-NKRP1A (CD161) interaction enhances natural killer cell-mediated lysis of triple-negative breast cancer cells. Am J Cancer Res. 2018;8(6):1050-1063.
Braud VM, Biton J, Becht E, et al. Expression of LLT1 and its receptor CD161 in lung cancer is associated with better clinical outcome. Oncoimmunology. 2018;7(5):e1423184.
Llibre A, López-Macías C, Marafioti T, et al. LLT1 and CD161 Expression in Human Germinal Centers Promotes B Cell Activation and CXCR4 Downregulation. J Immunol. 2016;196(5):2085-94.

No competing interests reported.

Download PDF

Journal Publication

published 28 Jan, 2022

Read the published version in BMC Urology →

Editorial decision: Major revision
08 Jul, 2021
Reviews received at journal
26 May, 2021
Reviewers agreed at journal
09 May, 2021
Reviews received at journal
26 Apr, 2021
Reviewers agreed at journal
12 Apr, 2021
Reviewers invited by journal
11 Apr, 2021
Editor assigned by journal
11 Apr, 2021
Editor invited by journal
07 Apr, 2021
Submission checks completed at journal
07 Apr, 2021
First submitted to journal
03 Apr, 2021

You are reading this latest preprint version

A 7-gene signature predicts the prognosis of patients with bladder cancer

Status:

Journal Publication

Version 1

Abstract

Figures

1. Introduction

2. Materials And Methods

3. Results

4. Discussion

5. Conclusion

Abbreviations

Declarations

References

Additional Declarations

Status:

Journal Publication

Version 1