A five-gene signature predicts overall survival of patients with papillary renal cell carcinoma
A five-gene signature predicts overall survival of patients with papillary renal cell carcinoma
Ze Gao 0 1
Dong Zhang 1
Yi Duan 1
Lei Yan 0 1
Yidong Fan 0 1
Zhiqing Fang 0 1
Zhaoxu Liu 0 1
0 Department of Urology, Qilu Hospital of Shandong University , Jinan, Shandong , China , 2 Department of Breast Surgery, Qilu Hospital of Shandong University , Jinan, Shandong , China , 3 School of Nursing, Shandong University , Jinan, Shandong , China , 4 Department of Oncology and Pathology, Karolinska Institutet , Stockholm , Sweden
1 Editor: Aamir Ahmad, University of South Alabama Mitchell Cancer Institute , UNITED STATES
Data Availability Statement: The data underlying
this study are freely available from The Cancer
Genome Atlas (TCGA) via GDC data portal at
The authors did not have special access privileges.
Funding: This study was funded by Grant
2016GSF201172 from the Key Project of Research
and Development Plan of Shandong Province,
China, and by Grant 2014qlky16 from the Nanshan
Group of Qilu hospital.
Competing interests: The authors have declared
that no competing interests exist.
The present study aims to investigate the gene expression changes in papillary renal cell
carcinoma(pRCC) and screen several genes and associated pathways of papillary renal cell
The papillary renal cell carcinoma RNA sequencing (RNA-seq) data set was downloaded
from TCGA (The Cancer Genome Atlas). We identified the differentially expressed mRNAs
between cancer and normal tissues and performed annotation of differentially expressed
mRNAs to figure out the functions and pathways they were enriched in. Then, we
constructed a risk score that relied on the 5-mRNA. The optimal value for the
patients?classification risk level was identified by ROC analysis. The relationship between mRNA expression
and prognosis of papillary renal cell carcinoma was evaluated by univariate Cox regression
model. The 5-mRNA based risk score was validated in both complete set and testing set.
In general, the 5-mRNA (CCNB2, IGF2BP3, KIF18A, PTTG1, and BUB1) were identified
and validated, which can predict papillary renal cell carcinoma patient survival. This study
revealed the 5-mRNA expression profile and the potential function of a single mRNA as a
prognostic target for papillary renal cell carcinoma.
In addition, these findings may have significant implications for potential treatments options
and prognosis for patients with papillary renal cell carcinoma.
The incidence of kidney cancer approximately accounted for 2.4% of adult malignancies,
estimated 338,000 new cases in 2012[
]. Renal cell carcinoma (RCC), the seventh most common
histological type of cancer, affects men more frequently than women and comprises nearly
90% of all kidney tumors[
], with clear cell (70%), papillary (10?15%), and chromophobe
(5%) carcinoma the main histologic types. Papillary renal cell carcinoma (pRCC) is a renal
parenchyma malignant tumor with papillary or tubulopapillary architecture, including type 1
and type 2 pRCC[
]. According to the original description, type 1 is often multifocal,
characterized by single layered small cell and scanty cytoplasm whereas type 2 is heterogeneous,
characterized by large pseudostratified cells and eosinophilic cytoplasm[
]. The division of pRCCs
into type 1 and type 2 tumors has been proved to have prognostic significance. The mutation
of MET oncogene is a crucial step into the pathogenesis of hereditary pRCC forms, resulting
in constitutive activation of the tyrosine kinase domain, which leads to increasing unregulated
proliferation, invasion and metastases. However, they can be found only in a small proportion
of sporadic cases[
] and a greater number show somatic copy number gains involving
]. Although papillary renal cell carcinoma is indolent, bilateral and multifocal
in some patients, other patients have a separate lesion with an aggressive clinical course.
The pRCC type 1 has been associated with overexpression or activation mutation of MET
proto-oncogene which encodes for a hepatocyte growth factor receptor (HGFR)[
pRCC type 2 is concerned with activation of the NRF2/antioxidant response element (ARE)
]. It has been reported that the 5-year overall survival of localized pRCC rates of
78?79% and cancer-specific 5-year survival rates of 86?94%[
]. Unfortunately, there are
currently no effective forms of therapy for patients with advanced disease. Although many
patients have greatly different treatment responses and prognoses, their tumor types are
identical in histology.
With the rise of gene sequencing, a large number of messenger RNAs (mRNA) have been
discovered in the role of cancer. The mRNA is a large family of RNA molecules that convey
genetic information from DNA to the ribosome, where they specify the amino acid sequence
of the protein products of gene expression. A lot of evidences show that mRNAs are wildly
involved in fundamental cellular processes, such as cell differentiation, proliferation, growth,
mobility, and apoptosis, as well as carcinogenesis or cancer progression.
Recently, large-scale gene data sets regarding pRCC were downloaded from TCGA
database (https://cancergenome.nih.gov/). The R Programming Language (R) was utilized for
preprocessing and analysis of these data to obtain the differentially expressed genes (DEGs). In
the present study, we screened the differentially expressed mRNAs between pRCC tissues and
matched normal tissues, and found the association of different mRNAs expression with
clinical data. We developed and validated a risk score for a prognosis based mRNA signature to
demonstrate the relationship between mRNA and the prognosis of pRCC by using sample
segmentation and Cox regression analysis.
2.Materials and methods
2.1 Data source and data pre-processing
The RNA sequencing data with papillary renal cell carcinoma (pRCC) was downloaded from
TCGA dataset (June 13.2018). There were 321 pRCC samples in total, including 289 primary
pRCC and 32 normal ones. In total, we screened 285 available samples from 321 samples, 1
sample information was repeated, and 3 samples lost survival data. The inclusion criteria were
as follows:(1) a histological diagnosis of papillary renal cell carcinoma, which is a subtype of
2 / 16
renal cell carcinoma, (2) patients underwent radical nephrectomy or partial nephrectomy, (3)
patients that didn?t receive any neoadjuvant treatment, (4) patients with corresponding
RNAseq data, (5) patients with explicit clinical prognostic data. Besides, the corresponding clinical
characteristics, including sex, age, pathological stage, clinical stage, definite TNM stage and
histological type of our cohort were integrated as a clinical statistical chart in Table 1. Signal
value for all genes were transformed to the log base 2. Standardization of the level of mRNA
expression for each sample was normalized with median quantification. Since the data comes
from the TCGA database, no further approval is required from the Ethics Committee.
2.2 Screening of differentially expressed mRNA
The RNA sequencing data (level 3) of 321 pRCC samples was downloaded from the TCGA
data portal. The differentially expressed mRNAs between pRCC tissues and matched normal
aPathologic stage was not available in 29 case
bclinical stage was not available in 87 case
cTumor type was not available in 124 case
dT stage was not available in 2 case
eN stage was not available in 113 case
f M stage was not available in 68 case.
Age is reported as the median(range).
3 / 16
ones were analyzed using the limma package in R[
]. The edgeR package and DESeq2
package in R were subsequently used for the calculation of DEGs. The adjust P value <0.05 and |
log2FC|>1 were set as the cut-off criteria. The genes that presented in both edgeR and DESeq2
package analysis results were selected as the final DEGs.
In our study, we mainly used the program code written in R language to analyze and deal
with RNA-seq data.
2.3 Identification and selection of prognostic related mRNA
The patients were further randomly divided into the training set (n = 143) and the testing set
(n = 142) with caret package in R language. In the training set, the association between
mRNAs expression and the overall survival time of patients were analyzed. Univariate Cox
regression analysis was performed with Survival package in R language[
]. In the survival
analysis, We selected the time that patient?s initial pathological diagnosis to papillary renal cell
carcinoma as the start point, and the time that death or censored events occurred as the end
point. Messenger RNAs with expressing significance p values <0.05 were selected as a target
mRNA. To further ensure the reliability and feasibility of clinical prognosis based on mRNA,
we made a Robust likelihood-based survival analysis by using Rbsurv package in R[
285 samples were again randomly assigned to a training set or a testing set with the caret
package in R language. We fitted a gene to the training set of samples and obtained a corresponding
parameter estimate. Then we evaluated log likelihood with the parameter estimate and the
testing set of samples. This procedure was repeated 10 times independently, which resulted in 10
log-likelihood for each gene. The best gene, g (1), with the largest mean log likelihood was
selected. We searched the next best gene by evaluating every two-gene model and selected an
optimal one with the largest mean log likelihood. Akaike information criterions (AICs) for all
the candidate models were computed and an optimal predictive model was selected as with the
lowest AIC value. The above procedure was repeated 1000 times, thus obtained 1000 loglike s
for each gene. The mRNA with the highest frequency (freq>300) was selected as the final
2.4 Establishment and validation risk score formula
The association between DEGs and patients?overall survival time was analysed in training set.
Univariate Cox regression analysis was performed in R language by Survival package with the
threshold of p value set as 0.05. In the univariate Cox regression analysis of the training set,
these prognostic-related mRNAs were included in their estimated regression coefficients, and
a risk scoring formula was established. In the training set, the patients were divided into
highrisk score or low-risk score group based on the optimal risk score cutoff value, which
represented the point at which the Youden index (sensitivity + specificity? 1) reached a maximum
value in training set. Then, we applied the optimal risk score cutoff value to the validation set
and the whole pRCC cohort to evaluate the classification performance of the model. By using
survival ROC package in R, the receiver operating characteristic (ROC) curve was obtained.
Then, choosing the optimal cutoff point with maximum sensitivity and specificity. Based on
the optimal cut-off value, the survival difference between the low-risk and high-risk group was
evaluated. The accuracy of the risk score formula was then further verified in the testing set
and the complete set.
2.5 Functional and pathway enrichment analysis of DEGs
The ConsensusPathDB (http://cpdb.molgen.mpg.de/) database is a molecular functional
interaction database, integrating information on protein interactions, genetic interactions
4 / 16
signaling, metabolism, gene regulation, and drug-target interactions in humans. In order to
better understand the biological functions and characteristic, the present study performed
Gene ontology (GO) using clusterProfiler package in R and Kyoto Encyclopedia of Genes and
Genomes (KEGG) pathway analyses in ConsensusPathDB database. All genes from human
beings were used as reference. In this study, we analyzed the DEGs that were significantly
upand down-regulated as determined from pRCC data. The p value of the univariate survival
analysis was less than 0.05 and considered statistically significant. To correct errors following
multiple comparison analysis, the Benjamini-Hochberg step-up method was used to control
3.1 Differentially expressed mRNAs in pRCC
With a cut-off value of P < 0.05 and |log2FC| >1.0, 4776 and 4831 differentially expressed
mRNAs were identified by using the DESeq2 and edgeR package in R language, respectively
(S1 File). The genes that presented in both edgeR and DESeq2 package analysis results were
selected as the final DEGs. According to this standard, 2632 up-regulated mRNAs and 1921
down-regulaed ones were identified (Fig 1A and 1B). We displayed the distribution of all the
differentially expressed mRNAs in both the?log(FDR) and logFC dimensions through a
volcano plot in Fig 2.
3.2 Identification and selection of prognostic related mRNA
The total of 4553 mRNAs was further randomized into a training set and a testing set. A total
of 749 differential expressed mRNAs were identified by univariable cox regression with the
P <0.05 in training set (S2 File). The top of 20 genes with the lowest p value were selected and
presented in Table 2. Random data analysis was performed by using Robust likelihood-based
modelling for 1000 times. Statistic frequency analysis of the significantly altered mRNAs
indicated that all the selected mRNAs had a high frequency (Fig 3). In further analysis,
differentially expressed mRNA with a frequency above 300 was picked out. Finally, there were 5
mRNAs identified as prognostic feature (Table 3). Among these genes, four mRNAs (CCNB2,
Fig 1. (a)The Venn polt of up-regulated DEG. (b)The Venn polt of down-regulated DEG.
5 / 16
Fig 2. Volcano map with edgeR and DESeq2.
IGF2BP3, KIF18A, PTTG1) had positive coefficients which suggested that higher expression
level was associated with worse survival time and one (BUB1) had negative coefficients
suggested that higher levels of expression were related with better survival time.
3.3 The 5-mRNA signature predicts the survival of patients with pRCC
A risk score formula based on the expression level and coefficient of 5 mRNAs was created as
follows: Risk score = (-0.560 expression level of BUB1) + (0.337 expression level of CCNB2)
+ (0.424 expression level of IGF2BP3) + (0.516 expression level of KIF18A) +
(0.310 ? expression level of PTTG1). Then, the mRNAs signature based risk score for each
patient in training set were calculated, and patients in the cohort were assigned into high-risk
group (n = 51) and low risk group (n = 92). The risk score was inversely related to the
patients?survival time. The patients with lower risk score had longer survival time and death patients
had higher risk score (Fig 4). The patients in high risk group had a tendency to have higher
expression levels of IGF2BP3, KIF18A, PTTG1, BUB1 and CCNB2.
3.4 Available prognostic indicator in pRCC patients
We performed ROC analysis for the 5-mRNA risk score in the training set to further validate
the sensitivity and specificity of survival prediction. The optimal cut-off value for false negative
and false positive minimum was 0.951, and the area under the curve (AUC) was 0.82 (Fig 5A).
The Kaplan-Meier curves showed that patients in the high-risk group had significantly shorter
overall survival than those in the low-risk group (log-rank test p < 0.0001)(Fig 5B).
3.5 External validation of five-gene signature
The predictive efficiency of the five-mRNA signature in testing set with 142 patients was then
evaluated. The patients in the testing set were classified into high-risk (n = 53) and low-risk
groups (n = 89) by using the same model and criteria. Similar to the training set, overall
survival was significantly lower in the high-risk group than in the low-risk group (p < 0.0001)(Fig
6A). Risk score-based classification of the external complete set also yielded similar results as
shown in Fig 6B. Despite unbalanced samples in each group, the results of survival analysis
6 / 16
were similar to those of low-risk patients with significantly longer survival times than those
with high-risk patients (p < 0.0001). The clinical characteristics available from the TCGA
database were integrated and univariate Cox regression analysis was performed to detect the
candidate clinical prognostic parameters. To validate independent predictive power of five-gene
signature from clinicopathological factors in pRCC cohort such as age, pathological stage,
clinical stage, pathological type and so on, the stratified Cox proportional hazard analysis was
constructed and visualized in Fig 7 and Fig 8, which suggested this signature could distinguish the
high-risk subgroup from low-risk one in each clinical subtype and further convinced the
independent prognostic efficacy of the final gene signature.
3.6 Functional and pathway enrichment analysis of DEGs
After performing GO analysis of DEGs with clusterProfiler package in R language, the DEGs
were classified into three groups: molecular function group, biological process group and
cellular component group. The biological results revealed that DEGs were primarily enriched in
chromosome segregation, mitotic nuclear division, nuclear chromosome segregation,
urogenital system development, renal system development and kidney development. The cellular
component results indicated that DEGs were mainly riched in chromosomal region, proteinaceous
extracellular matrix, condensed chromosome, spindle, chromosome and centromeric region.
The molecular function results showed that DEGs were mainly enriched in ion channel
activity and substrate-specific channel activity (Fig 9A?9C). To investigate pathway enrichment,
KEGG signaling pathway analysis was used to identify the pathways, including cAMP
signaling pathway, phospholipase D signaling pathway, Hippo signaling pathway, cell cycle,
complement and coagulation cascades, TGF-beta signaling pathway, inflammatory mediator
regulation of TRP channels and melanoma (Fig 9D).
7 / 16
Fig 3. Statistic frequency analysis of the significantly altered mRNAs.
In the present study, the R language was used to analysis the gene data downloaded from
TCGA databases. The total of 4553 DEGs in pRCC compared with control samples were
identified, which included 2632 up-regulated and 1921 down-regulated genes. The DEGs were
mainly enriched in 60 GO terms, including ion channel activity, chromosome segregation and
chromosomal region. The KEGG pathway enrichment analysis result showed that the DEGs
were associated with calcium,cAMP,phosphollipase D, and hippo signaling pathway. The Ca2
+-mediated signaling pathways have been implicated either directly or indirectly related to
tumorigenesis and tumor progression[
]. In this pathway, Orai1 exaggerates cell
proliferation, migration, invasion and evasion of apoptosis[
]. The cAMP plays a complex and
context-dependent role in regulating cell migration[
]. Usually the focus has been mainly on
PKA-mediated migratory effects. Shaikh et al. suggests that hypoxia enhances
cAMP-dependent protein kinase activity by up-regulating PKA gene expression in a HIF dependent
mechanism and that PKA plays a key role in hypoxia-mediated EMT, migration, and invasion in
lung cancer cells[
]. There is increasing evidence that the Hippo pathway is dysregulated in
many human cancers, and dysregulation of the Hippo pathway exerts a significant impact on
cancer development, including liver, breast, lung, colon, ovary, and others[
phospholipase D (PLD) expression has been identified in multiple facets of complex pathological
states, including cancer and inflammatory diseases. PLD contributes to various mitogenic or
oncogenic signaling pathways[
]. Therefore, monitoring of these signaling pathway maybe
beneficial to understanding the mechanism of carcinogenesis and researching treatment of
We identified 5 mRNAs that are associated with the survival of patients with papillary renal
cell carcinoma,namely CCNB2, IGF2BP3, KIF18A, PTTG1, and BUB1 in the training set. The
prognostic related mRNAs were further selected to construct a risk score formula by Cox
Fig 4. mRNA risk score analysis of the training set. The mRNA signature risk score distribution heat-map of the mRNA
expression profiles.Rows represent mRNAs, and columns represent patients.
regression model. Besides, we used the ROC analysis to identify the optimal cutoff point, and
divided patient into high-risk and low-risk groups. In the univariate Cox regression model, the
survival time of high-risk patients was significantly shortened. Subsequently, the GO and
KEGG pathway analysis suggest that mRNA plays a crucial role in molecular pathogenesis and
progression of pRCC.
Among the 5 mRNAs which are associated with the prognosis of pRCC, some have been
reported to express in cancer or other diseases, but have not been examined in pRCC. For
example, elevated cytoplasmic CCNB2 protein levels are strongly associated with short-term
disease-specific survival of breast cancer patients[
]. Cyclin B2 (CCNB2), the member of
cyclin family proteins, regulates the activities of cyclin dependent kinases (CDKs) and different
cyclins function spatially and temporally in specific phases of the cell cycle[
9 / 16
Fig 5. (a) ROC for 5 mRNA,Method = KM. Receiver operating characteristic (ROC) analyzes the sensitivity and
specificity of the survival time by risk score based on the 5-mRNA signature. The red dot indicates the optimal cut-off
point. (b) Method Kaplan Meier. Kaplan-Meier estimated the survival time of the training set patients using risk score
based on the 5 mRNA signature. The plot was used to visualize the survival probability for the low-risk versus high-risk
group of patients based on the optimal cut-off point.
Fig 6. (a)Method Kaplan Meier. Kaplan-Meier estimated the survival time of the testing set patients using a risk score based on 5 mRNA signature. (b) Method
Kaplan Meier. Kaplan-Meier estimated the survival time of the complete set patients from using a risk score based on 5 mRNA signature. The plot was used to
visualize the survival probability for the low-risk versus high-risk group of patients based on the optimal cut-off point.
10 / 16
Fig 7. Stratified Cox hazard analysis of five gene signature on clinicopathological characteristics in pRCC cohort.
commonly triggers the G2/M transfer by activating CDK1 kinase[
]. CCNB2 is generally
located in the Golgi apparatus during both interphase and mitosis. Many evidences
Fig 8. The comparative gene expression level of CCNB2, IGF2BP3, KIF18A, PTTG1, and BUB1 in normal tissue and pRCC tissue from TCGA database,
11 / 16
Fig 9. (a) biological process group of GO analysis. (b) cellular component group of GO analysis. (c) molecular function group of GO analysis.(d) The KEGG
analysis. Biological function and KEGG pathway analysis of target genes. The overlapping target genes were predicted using the R language. (a-c) The enriched
GO biological processes of target genes. (d) The enriched KEGG pathways of target genes.
suggested that CCNB2 expression was increased in a variety of human cancers, such as
nonsmall cell lung cancer[
], breast carcinoma[
], gastric cancer[
], pituitary adenoma[
] and adrenocortical carcinoma[
]. Therefore, CCNB2 is
important in pRCC and may be used as a prognostic indicator. It has been reported that
insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) expression correlates with
]. Schaeffer et al. suggested that IGF2BP3 was found to be selectively overexpressed in
pancreatic ductal adenocarcinoma tissues but not in benign pancreatic tissues[
]. In ovarian
cancer, high expression of IGF2BP3 was associated with poor survival, and women diagnosed
at advanced stages with elevated IGF2BP3 was at higher risk of developing chemoresistance
]. Lochhead et al. revealed that normal colorectal epithelium was negative for IGF2BP3 in
patients of normal mucosa adjacent to carcinoma, and IGF2BP3 was associated with poor
differentiation, stage III?IV disease, BRAF mutation, and LINE-1 hypomethylation[
]. Lin et al.
demonstrated that IGF2BP3 enhances cell invasion ability and tumorigenicity in human
OSCC in vitro and in vivo[
]. In the present study, it was demonstrated that the pRCC
patients with IGF2BP3 alterations exhibited a poorer survival rate compared with those
without the genetic alterations. This result suggested that the mutation in IGF2BP3 reduces the
survival rate of patients with pRCC. Kif18A, a member of the kinesin super family of molecular
motor proteins, regulates chromosome congregation and suppresses kinetochore movements
to control mitotic chromosome alignment in the pre-anaphase state of the mammalian cell
]. Notably, Nagahara et al. demonstrated that colorectal cancer cells transfected with
Kif18A cDNA demonstrated significant enhanced migration and invasion compared to
]. KIF18A might be a biomarker for HCC diagnosis and an independent
predictor of DFS and OS after surgical resection[
]. The results of the present study were
12 / 16
consistent with the results of previous studies. KIF18A may have an important role in
progression of pRCC, however this requires further study, in order to verify the specific molecular
marker role of KIF18A in the diagnosis of patients with pRCC. Pituitary tumour transforming
gene 1 (PTTG1) is over-expressed in a vast array of malignancies including pituitary [
] and lung[
] cancer. The high level of PTTG1 is commonly
associated with an enhanced proliferative capacity, increased tumour grade and high invasive
potential. Current evidences regarding the biological role of BUB1 in cancer is contradictory.
Mutations in BUB1, some of which are functional, occur in some cancers, including those that
originate in the lung, colon, and are reported to be associated with chromosomal instability
and lymph node metastasis, suggesting that silencing of this kinase may mediate aggressive
]. Moreover, Aurora b hyperactivation caused by overexpression of BUB1
leads to misaggregation of chromosomes, thereby inducing aneuploidy[
Taken all together, by performing a comprehensive analysis for differentially expressed mRNA
profiles and corresponding clinical information, our study demonstrated that five-mRNA
signature was a potential diagnostic marker in pRCC, and was an independent prognostic factor
in pRCC patients. This signature has a lot of potential prognostic and therapeutic implications
for the pRCC patient management. However, further research is needed to validate our
findings and establish molecular mechanisms for mRNAs interactions and papillary renal cell
S1 File. Details of the result of univariable cox regression analysis.
S2 File. The mRNA with p-value less than 0.05 in univariable cox regression analysis.
This study was funded by Grant 2016GSF201172 from the Key Project of Research and
Development Plan of Shandong Province, China, and by Grant 2014qlky16 from the Nanshan
Group of Qilu hospital.
This study is based upon data generated by the TCGA Research Network: http://
Conceptualization: Zhaoxu Liu.
Data curation: Zhiqing Fang.
Investigation: Yi Duan.
Methodology: Yidong Fan.
Resources: Lei Yan.
Software: Dong Zhang.
Writing ? original draft: Ze Gao.
13 / 16
14 / 16
15 / 16
ovarian cancer. British journal of cancer. 2015; 113(3):414?24. https://doi.org/10.1038/bjc.2015.254
PMID: 26158423; PubMed Central PMCID: PMC4522643.
1. Ferlay J , Soerjomataram I , Dikshit R , Eser S , Mathers C , Rebelo M , et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012 . International journal of cancer . 2015 ; 136 ( 5 ): E359 - 86 . https://doi.org/10.1002/ijc.29210 PMID: 25220842 .
2. Badowska-Kozakiewicz AM , Budzik MP , Koczkodaj P , Przybylski J . Selected tumor markers in the routine diagnosis of chromophobe renal cell carcinoma . Archives of medical science: AMS . 2016 ; 12 ( 4 ): 856 - 63 . https://doi.org/10.5114/aoms. 2015 .51188 PMID: 27478468; PubMed Central PMCID : PMC4947610 .
3. Jemal A , Siegel R , Ward E , Hao Y , Xu J , Murray T , et al. Cancer statistics , 2008 . CA: a cancer journal for clinicians . 2008 ; 58 ( 2 ): 71 - 96 . https://doi.org/10.3322/CA. 2007 .0010 PMID: 18287387 .
4. Znaor A , Lortet-Tieulent J , Laversanne M , Jemal A , Bray F . International variations and trends in renal cell carcinoma incidence and mortality . European urology . 2015 ; 67 ( 3 ): 519 - 30 . https://doi.org/10.1016/ j.eururo. 2014 . 10 .002 PMID: 25449206 .
5. Montironi R , Cheng L , Scarpelli M , Lopez-Beltran A . Pathology and Genetics: Tumours of the Urinary System and Male Genital System: Clinical Implications of the 4th Edition of the WHO Classification and Beyond . European urology. 2016 ; 70 ( 1 ): 120 - 3 . https://doi.org/10.1016/j.eururo. 2016 . 03 .011 PMID: 26996660 .
6. Delahunt B , Eble JN . Papillary renal cell carcinoma: a clinicopathologic and immunohistochemical study of 105 tumors. Modern pathology: an official journal of the United States and Canadian Academy of Pathology, Inc . 1997 ; 10 ( 6 ): 537 - 44 . PMID: 9195569 .
7. Courthod G , Tucci M , Di Maio M , Scagliotti GV . Papillary renal cell carcinoma: A review of the current therapeutic landscape . Critical reviews in oncology/hematology . 2015 ; 96 ( 1 ): 100 - 12 . https://doi.org/10. 1016/j.critrevonc. 2015 . 05 .008 PMID: 26052049 .
8. Ridge CA , Pua BB , Madoff DC . Epidemiology and staging of renal cell carcinoma . Seminars in interventional radiology . 2014 ; 31 ( 1 ): 3 - 8 . https://doi.org/10.1055/s-0033-1363837 PMID: 24596434; PubMed Central PMCID : PMC3930658 .
9. Schmidt L , Duh FM , Chen F , Kishida T , Glenn G , Choyke P , et al. Germline and somatic mutations in the tyrosine kinase domain of the MET proto-oncogene in papillary renal carcinomas . Nature genetics . 1997 ; 16 ( 1 ): 68 - 73 . https://doi.org/10.1038/ng0597-68 PMID: 9140397.
10. Lee BH . Commentary on: "Comprehensive molecular characterization of papillary renal-cell carcinoma." Cancer Genome Atlas Research Network .: N Engl J Med . 2016 Jan 14 ; 374 ( 2 ): 135 - 45 . Urologic oncology. 2017 ; 35 ( 9 ): 578 - 9 . https://doi.org/10.1056/NEJMoa1505917 PMID: 26536169 .
11. Albiges L , Guegan J , Le Formal A , Verkarre V , Rioux-Leclercq N , Sibony M , et al. MET is a potential target across all papillary renal cell carcinomas: result from a large molecular study of pRCC with CGH array and matching gene expression array . Clinical cancer research: an official journal of the American Association for Cancer Research . 2014 ; 20 ( 13 ): 3411 - 21 . https://doi.org/10.1158/ 1078 - 0432 .CCR- 13 - 2173 PMID: 24658158 .
12. Ooi A , Wong JC , Petillo D , Roossien D , Perrier-Trudova V , Whitten D , et al. An antioxidant response phenotype shared between hereditary and sporadic type 2 papillary renal cell carcinoma . Cancer cell . 2011 ; 20 ( 4 ): 511 - 23 . https://doi.org/10.1016/j.ccr. 2011 . 08 .024 PMID: 22014576 .
13. Steffens S , Janssen M , Roos FC , Becker F , Schumacher S , Seidel C , et al. Incidence and long-term prognosis of papillary compared to clear cell renal cell carcinoma-a multicentre study . European journal of cancer . 2012 ; 48 ( 15 ): 2347 - 52 . https://doi.org/10.1016/j.ejca. 2012 . 05 .002 PMID: 22698386 .
14. Patard JJ , Leray E , Rioux-Leclercq N , Cindolo L , Ficarra V , Zisman A , et al. Prognostic value of histologic subtypes in renal cell carcinoma: a multicenter experience . Journal of clinical oncology: official journal of the American Society of Clinical Oncology . 2005 ; 23 ( 12 ): 2763 - 71 . https://doi.org/10.1200/ JCO. 2005 . 07 .055 PMID: 15837991 .
15. Ritchie ME , Phipson B , Wu D , Hu Y , Law CW , Shi W , et al. limma powers differential expression analyses for RNA-sequencing and microarray studies . Nucleic acids research . 2015 ; 43 ( 7 ):e47. https://doi. org/10.1093/nar/gkv007 PMID: 25605792; PubMed Central PMCID : PMC4402510 .
16. O 'Quigley J , Moreau T. Cox's regression model: computing a goodness of fit statistic . Computer methods and programs in biomedicine. 1986 ; 22 ( 3 ): 253 - 6 . PMID: 3524984 .
Wang JY , Tai JJ . Robust quantitative trait association tests in the parent-offspring triad design: conditional likelihood-based approaches . Annals of human genetics . 2009 ; 73 ( 2 ): 231 - 44 . https://doi.org/10.
1111/j.1469- 1809 . 2008 . 00502 . x PMID : 19183344 .
18. Renaud G , Stenzel U , Maricic T , Wiebe V , Kelso J. deML: robust demultiplexing of Illumina sequences using a likelihood-based approach . Bioinformatics . 2015 ; 31 ( 5 ): 770 - 2 . https://doi.org/10.1093/ bioinformatics/btu719 PMID: 25359895; PubMed Central PMCID : PMC4341068 .
19. Monteith GR , McAndrew D , Faddy HM , Roberts-Thomson SJ . Calcium and cancer: targeting Ca2+ transport . Nature reviews Cancer . 2007 ; 7 ( 7 ): 519 - 30 . https://doi.org/10.1038/nrc2171 PMID: 17585332 .
20. Prevarskaya N , Skryma R , Shuba Y. Calcium in tumour metastasis: new roles for known actors . Nature reviews Cancer . 2011 ; 11 ( 8 ): 609 - 18 . https://doi.org/10.1038/nrc3105 PMID: 21779011 .
21. Motiani RK , Hyzinski-Garcia MC , Zhang X , Henkel MM , Abdullaev IF , Kuo YH , et al. STIM1 and Orai1 mediate CRAC channel activity and are essential for human glioblastoma invasion . Pflugers Archiv: European journal of physiology . 2013 ; 465 ( 9 ): 1249 - 60 . https://doi.org/10.1007/s00424-013-1254-8 PMID: 23515871; PubMed Central PMCID : PMC3748246 .
22. Howe AK . Cross-talk between calcium and protein kinase A in the regulation of cell migration . Current opinion in cell biology . 2011 ; 23 ( 5 ): 554 - 61 . https://doi.org/10.1016/j.ceb. 2011 . 05 .006 PMID: 21665456; PubMed Central PMCID : PMC3175278 .
23. Shaikh D , Zhou Q , Chen T , Ibe JC , Raj JU , Zhou G. cAMP-dependent protein kinase is essential for hypoxia-mediated epithelial-mesenchymal transition, migration, and invasion in lung cancer cells . Cellular signalling . 2012 ; 24 ( 12 ): 2396 - 406 . https://doi.org/10.1016/j.cellsig. 2012 . 08 .007 PMID: 22954688 .
24. Yu FX , Zhao B , Guan KL . Hippo Pathway in Organ Size Control, Tissue Homeostasis, and Cancer . Cell . 2015 ; 163 ( 4 ): 811 - 28 . https://doi.org/10.1016/j.cell. 2015 . 10 .044 PMID: 26544935; PubMed Central PMCID : PMC4638384 .
25. Kang DW , Choi KY , Min do S. Functional regulation of phospholipase D expression in cancer and inflammation . The Journal of biological chemistry . 2014 ; 289 ( 33 ): 22575 - 82 . https://doi.org/10.1074/jbc. R114.569822 PMID: 24990948; PubMed Central PMCID : PMC4132765 .
26. Shubbar E , Kovacs A , Hajizadeh S , Parris TZ , Nemes S , Gunnarsdottir K , et al. Elevated cyclin B2 expression in invasive breast carcinoma is associated with unfavorable clinical outcome . BMC cancer . 2013 ; 13 :1. https://doi.org/10.1186/ 1471 -2407-13-1 PMID: 23282137; PubMed Central PMCID : PMC3545739 .
27. Vermeulen K , Van Bockstaele DR , Berneman ZN . The cell cycle: a review of regulation, deregulation and therapeutic targets in cancer . Cell proliferation . 2003 ; 36 ( 3 ): 131 - 49 . PMID: 12814430 .
28. Petri ET , Errico A , Escobedo L , Hunt T , Basavappa R. The crystal structure of human cyclin B . Cell cycle . 2007 ; 6 ( 11 ): 1342 - 9 . https://doi.org/10.4161/cc.6.11.4297 PMID: 17495533 .
29. Wu T , Zhang X , Huang X , Yang Y , Hua X . Regulation of cyclin B2 expression and cell cycle G2/m transition by menin . The Journal of biological chemistry . 2010 ; 285 ( 24 ): 18291 - 300 . https://doi.org/10.1074/ jbc.M110.106575 PMID: 20404349; PubMed Central PMCID : PMC2881754 .
30. Jackman M , Firth M , Pines J . Human cyclins B1 and B2 are localized to strikingly different structures: B1 to microtubules, B2 primarily to the Golgi apparatus . The EMBO journal . 1995 ; 14 ( 8 ): 1646 - 54 . PMID: 7737117; PubMed Central PMCID : PMC398257 .
31. Qian X , Song X , He Y , Yang Z , Sun T , Wang J , et al. CCNB2 overexpression is a poor prognostic biomarker in Chinese NSCLC patients . Biomedicine & pharmacotherapy = Biomedecine & pharmacotherapie . 2015 ; 74 : 222 - 7 . https://doi.org/10.1016/j.biopha. 2015 . 08 .004 PMID: 26349989 .
32. Wang DG , Chen G , Wen XY , Wang D , Cheng ZH , Sun SQ . Identification of biomarkers for diagnosis of gastric cancer by bioinformatics . Asian Pacific journal of cancer prevention: APJCP . 2015 ; 16 ( 4 ): 1361 - 5 . PMID: 25743799 .
33. Park SH , Yu GR , Kim WH , Moon WS , Kim JH , Kim DG . NF-Y-dependent cyclin B2 expression in colorectal adenocarcinoma . Clinical cancer research: an official journal of the American Association for Cancer Research . 2007 ; 13 ( 3 ): 858 - 67 . https://doi.org/10.1158/ 1078 - 0432 .CCR- 06 -1461 PMID: 17289878 .
34. De Martino I , Visone R , Wierinckx A , Palmieri D , Ferraro A , Cappabianca P , et al. HMGA proteins upregulate CCNB2 gene in mouse and human pituitary adenomas . Cancer research . 2009 ; 69 ( 5 ): 1844 - 50 . https://doi.org/10.1158/ 0008 - 5472 .CAN- 08 -4133 PMID: 19223528 .
35. Fernandez-Ranvier GG , Weng J , Yeh RF , Khanafshar E , Suh I , Barker C , et al. Identification of biomarkers of adrenocortical carcinoma using genomewide gene expression profiling . Archives of surgery . 2008 ; 143 ( 9 ): 841 - 6 ; discussion 6. https://doi.org/10.1001/archsurg.143.9.841 PMID: 18794420 .
36. Ennajdaoui H , Howard JM , Sterne-Weiler T , Jahanbani F , Coyne DJ , Uren PJ , et al. IGF2BP3 Modulates the Interaction of Invasion-Associated Transcripts with RISC . Cell reports . 2016 ; 15 ( 9 ): 1876 - 83 . https://doi.org/10.1016/j.celrep. 2016 . 04 .083 PMID: 27210763; PubMed Central PMCID : PMC4889463 .
37. Schaeffer DF , Owen DR , Lim HJ , Buczkowski AK , Chung SW , Scudamore CH , et al. Insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) overexpression in pancreatic ductal adenocarcinoma correlates with poor survival . BMC cancer . 2010 ; 10 : 59 . https://doi.org/10.1186/ 1471 -2407-10-59 PMID: 20178612; PubMed Central PMCID : PMC2837867 .
38. Hsu KF , Shen MR , Huang YF , Cheng YM , Lin SH , Chow NH , et al. Overexpression of the RNA-binding proteins Lin28B and IGF2BP3 (IMP3) is associated with chemoresistance and poor disease outcome in
39. Lochhead P , Imamura Y , Morikawa T , Kuchiba A , Yamauchi M , Liao X , et al. Insulin-like growth factor 2 messenger RNA binding protein 3 (IGF2BP3) is a marker of unfavourable prognosis in colorectal cancer . European journal of cancer . 2012 ; 48 ( 18 ): 3405 - 13 . https://doi.org/10.1016/j.ejca. 2012 . 06 .021 PMID: 22840368; PubMed Central PMCID : PMC3613860 .
40. Lin CY , Chen ST , Jeng YM , Yeh CC , Chou HY , Deng YT , et al. Insulin-like growth factor II mRNA-binding protein 3 expression promotes tumor formation and invasion and predicts poor prognosis in oral squamous cell carcinoma . Journal of oral pathology & medicine: official publication of the International Association of Oral Pathologists and the American Academy of Oral Pathology . 2011 ; 40 ( 9 ): 699 - 705 . https://doi.org/10.1111/j.1600- 0714 . 2011 . 01019 . x PMID : 23647076 .
41. Malumbres M , Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm . Nature reviews Cancer . 2009 ; 9 ( 3 ): 153 - 66 . https://doi.org/10.1038/nrc2602 PMID: 19238148 .
42. Nagahara M , Nishida N , Iwatsuki M , Ishimaru S , Mimori K , Tanaka F , et al. Kinesin 18A expression: clinical relevance to colorectal cancer progression . International journal of cancer . 2011 ; 129 ( 11 ): 2543 - 52 . https://doi.org/10.1002/ijc.25916 PMID: 21213216 .
43. Liao W , Huang G , Liao Y , Yang J , Chen Q , Xiao S , et al. High KIF18A expression correlates with unfavorable prognosis in primary hepatocellular carcinoma . Oncotarget . 2014 ; 5 ( 21 ): 10271 - 9 . https://doi. org/10.18632/oncotarget.2082 PMID: 25431949; PubMed Central PMCID : PMC4279371 .
44. Saez C , Japon MA , Ramos-Morales F , Romero F , Segura DI , Tortolero M , et al. hpttg is over-expressed in pituitary adenomas and other primary epithelial neoplasias . Oncogene . 1999 ; 18 ( 39 ): 5473 - 6 . https:// doi.org/10.1038/sj.onc.1202914 PMID: 10498901 .
45. Zhang X , Horwitz GA , Heaney AP , Nakashima M , Prezant TR , Bronstein MD , et al. Pituitary tumor transforming gene (PTTG) expression in pituitary adenomas . The Journal of clinical endocrinology and metabolism . 1999 ; 84 ( 2 ): 761 - 7 . https://doi.org/10.1210/jcem.84.2.5432 PMID: 10022450 .
46. Boelaert K , McCabe CJ , Tannahill LA , Gittoes NJ , Holder RL , Watkinson JC , et al. Pituitary tumor transforming gene and fibroblast growth factor-2 expression: potential prognostic indicators in differentiated thyroid cancer . The Journal of clinical endocrinology and metabolism . 2003 ; 88 ( 5 ): 2341 - 7 . https://doi. org/10.1210/jc.2002-021113 PMID: 12727994 .
47. Heaney AP , Singson R , McCabe CJ , Nelson V , Nakashima M , Melmed S . Expression of pituitarytumour transforming gene in colorectal tumours . Lancet . 2000 ; 355 ( 9205 ): 716 - 9 . https://doi.org/10. 1016/S0140- 6736 ( 99 ) 10238 - 1 PMID: 10703804 .
48. Rehfeld N , Geddert H , Atamna A , Rohrbeck A , Garcia G , Kliszewski S , et al. The influence of the pituitary tumor transforming gene-1 (PTTG-1) on survival of patients with small cell lung cancer and nonsmall cell lung cancer . Journal of carcinogenesis . 2006 ; 5:4 . https://doi.org/10.1186/ 1477 -3163-5-4 PMID: 16426442; PubMed Central PMCID : PMC1360069 .
49. Davidson B , Nymoen DA , Elgaaen BV , Staff AC , Trope CG , Kaern J , et al. BUB1 mRNA is significantly co-expressed with AURKA and AURKB mRNA in advanced-stage ovarian serous carcinoma . Virchows Archiv: an international journal of pathology . 2014 ; 464 ( 6 ): 701 - 7 . https://doi.org/10.1007/s00428-014- 1577-7 PMID: 24756216 .
50. Ricke RM , van Deursen JM. Aurora B hyperactivation by Bub1 overexpression promotes chromosome missegregation . Cell cycle . 2011 ; 10 ( 21 ): 3645 - 51 . https://doi.org/10.4161/cc.10.21.18156 PMID: 22033440; PubMed Central PMCID : PMC3266005 .