Oncotarget

Research Papers:

Meta-analysis of gene expression and integrin-associated signaling pathways in papillary renal cell carcinoma subtypes

PDF |  HTML  |  Supplementary Files  |  How to cite

Oncotarget. 2016; 7:84178-84189. https://doi.org/10.18632/oncotarget.12390

Metrics: PDF 2373 views  |   HTML 3466 views  |   ?  

Kai Zhang, Hang-Mao Lee, Gong-Hong Wei and Aki Manninen _

Abstract

Kai Zhang1, Hang-Mao Lee1, Gong-Hong Wei1, Aki Manninen1

1Biocenter Oulu, Oulu Center for Cell-Matrix Research, Faculty of Biochemistry and Molecular Medicine, University of Oulu, Oulu 90220, Finland

Correspondence to:

Aki Manninen, email: [email protected]

Keywords: meta-analysis, papillary renal cell carcinoma, gene expression, integrin, gender

Received: March 24, 2016     Accepted: September 24, 2016     Published: October 01, 2016

ABSTRACT

Papillary renal cell carcinoma (PRCC) is the second most common renal cell carcinoma (RCC) that can be further subdivided into type 1 (PRCC1) and type 2 (PRCC2) RCCs based on histological and genetic features. PRCC2 is often more aggressive than PRCC1. While integrin-associated protein complexes mediate tumorigenesis and metastases in many types of cancers it is not known whether integrin-mediated signaling impacts PRCC and differs between PRCC1 and PRCC2. In this study, we combined the analysis of five PRCC gene expression datasets derived from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) by using integrative bioinformatics pipelines. We found 1475 differentially expressed genes among which 37 genes were associated with integrin pathways. In comparison with PRCC1, PRCC2 cases showed upregulated expression of α5-integrin (ITGA5) whereas the expression of α6- (ITGA6) and β8-integrins (ITGB8) was downregulated. Because PRCC2 occurs more frequently in men, the meta-analysis was extended to explore the gender effects. This analysis revealed 8 genes but none of them was related to integrin pathways suggesting that other mechanisms than integrin-mediated signaling underlie the observed gender differences in the pathogenicity of PRCC2.


INTRODUCTION

The second most common histological subtype of renal cell carcinomas (RCC) is the papillary RCC that accounts for 10–20 percent of all renal cancer cases [1]. Papillary renal cell carcinoma (PRCC) has two subtypes defined by different histological features. PRCC1 shows both papillae and tubular structures covered by small cells with scanty cytoplasm and small oval nuclei. PRCC2 indicates only papillary structures covered by large cells with abundant eosinophilic cytoplasm and large, spherical nuclei with prominent nucleoli [2]. PRCC2 is often a more aggressive disease that is associated with less differentiated histology phenotype, high number of nodal and distant metastases and worse survival rates comparing with PRCC1 [34]. PRCC1 and PRCC2 have distinct genetic backgrounds [5]. While overexpression or activating mutations of MET proto-oncogene encoding for a hepatocyte growth factor receptor (HGFR) are common in PRCC1, PRCC2 has been associated with activation of the NRF2-ARE pathway and CDKN2A silencing [68]. Interestingly, differential expression of components regulating cell-extracellular matrix (ECM) interactions has also been implicated in PRCC1 [9]. However, how cell-ECM interactions are modulated in PRCC2 remains unknown.

Integrins are heterodimeric cell surface receptors that mediate interactions between cells and the ECM [10]. In solid tumors integrins regulate cancer initiation, stemness, drug resistance and metastasis by regulating the assembly of large multiprotein complexes [11]. These integrin-mediated complexes not only facilitate cells to adhere to the ECM but also convey various signals between cells and their immediate microenvironment to regulate cell behavior [1011]. Jones et al. reported that abnormal chemokine receptor signaling modulates the activity of α3-, α5-, β1-, β3- and β4-integrins thereby regulating the adhesive properties of clear cell renal cell carcinoma (ccRCC) cells in vitro [12]. However, possible contributions of modified integrin pathways underlying the histological differences between PRCC1 and PRCC2 have not been thoroughly addressed. Intriguingly, PRCC2 occurs more frequently and is more aggressive in male patients when compared with female patients, but the mechanism remains largely unknown [1315]. In order to address these issues, we combined available RNA-seq and microarray datasets of PRCC samples from GEO and TCGA and compared the mRNA expression profiles of PRCC1 and PRCC2. To get a comprehensive insight into molecular mechanisms and differences between PRCC1 and PRCC2, we performed a meta-analysis of gene expression focusing on the integrin-associated pathways using these large GEO and TCGA datasets. Our analysis highlighted differential regulation of integrin pathways between PRCC1 and PRCC2 but found no significant gender-associated changes in integrin pathway genes between male and female patients with PRCC2.

RESULTS

Meta-analysis of differentially expressed genes between PRCC1 and PRCC2

Four microarray datasets of PRCC patient samples for which matched clinical information was available were obtained from GEO by using GEOquery [16]. In addition, one RNA-seq dataset of PRCC samples was obtained from TCGA-kidney renal papillary cell carcinoma (KIRP) database by using TCGA-assembler [17] (Figure 1A and Table 1). After removing non-PRCC1/PRCC2 (or unidentified) samples and PRCC cases with mixed subtyping, a total of 138 PRCC1 and 135 PRCC2 samples were selected for further analysis (Figure 1A1C). PCA biplots of the six quality control criteria and the five PRCC datasets supported the inclusion of all of the datasets for meta-analysis (Figure 2A and Table 2). Three main meta-analysis methods by combining p-value in MetaDE package were employed: 1) Maximum p-value (maxP), 2) Minimum p-value (minP) and 3) r-th ordered p-value (roP) (Figure 2B) [18]. 1758, 1558 and 1976 differentially expressed (DE) genes were detected by maxP, minP and roP evaluation criteria, respectively, using detection competency curves and false discovery rate (FDR) cut-off less than 0.05 (Figure 2B, 2C and Table 3). All of the three analyses highlighted integrin pathways with significant overlap such that 37% of the DE genes within integrin pathways were shared by all three analyses and 55% were shared between two out of three analyses (Supplementary Figure S1A). Meta-analysis performed with maxP criteria, however, identified more integrin pathway-related DE genes than the other two and was thus selected for further analysis. All the 1758 DE genes obtained with maxP criteria were analyzed by using meta-analysis for pathway enrichment (MAPE) within the MetaPath package to reveal cellular pathways differentially regulated between PRCC1 and PRCC2 [19]. MAPE analysis identified 115 enriched pathways when the MAPE_I method (threshold set to 0.2) which integrates gene and pathway information from the DE genes was employed (Figure 3) [19]. Cell adhesion was one of the pathways appearing from this analysis but not among the most prominently enriched pathways. To further refine the list of DE genes between the PRCC subtypes we employed a random-effect model (REM) meta-analysis approach that estimates the differences in gene expression across all the different datasets by combining the individual effects sizes into a meta-effect size (ES) [1820]. To this end MetaDE.ES was applied for the raw data in the five datasets (Table 1) resulting in 2610 DE genes when the FDR threshold was set to 0.05 (Supplementary Figure S1B and Table 3). Finally, the DE genes obtained from each of these two meta-analyses were plotted into a Venn diagram that revealed 1475 DE genes common to both maxP meta-analysis method that combines p-values and REM meta-analysis method that combines effect sizes (Figure 4A).

Overview of the meta-analysis and the PRCC datasets.

Figure 1: Overview of the meta-analysis and the PRCC datasets. (A) A scheme of the meta-analysis workflow. (B) The number of histologically classified PRCC1 and PRCC2 samples in the five datasets. (C) The number of female (blue) and male (red) PRCC2 patients in the three datasets where this information was available.

Table 1: The source data for each of the datasets used in this study

Source

Platform

GEO Accession

Ooi A et al 2011

GPL11433

GSE26574

Furge KA et al 2011

GPL4866

GSE7023

Kort EJ et al 2008

GPL6671

GSE11024

Yang XJ et al 2013

GPL570

GSE2748

TCGA-KIRP 2016

IlluminaHiseq RNASeqV2

NA

GEO = Gene Expression Omnibus. KIRP= Kidney renal papillary cell carcinoma.

Meta-analysis of differentially expressed genes between PRCC1 and PRCC2 by combining p-values.

Figure 2: Meta-analysis of differentially expressed genes between PRCC1 and PRCC2 by combining p-values. (A) PCA biplot of quality control measures in five PRCC studies. (B) The number of differentially expressed genes plotted as a function of false discovery rate FDR in the analysis of five different datasets and the three different meta-analysis algorithms (maxP, minP and roP). (C) A heat map representation of DE genes in PRCC1 and PRCC2 samples subjected to maxP meta-analysis (FDR = 0.05).

Table 2: The quality control analysis for PRCC1 and PRCC2 in the different datasets

Number

Study

IQC

EQC

CQCg

CQCp

AQCg

AQCp

Rank

1

GSE26574

7.39

3.7

307.65

307.65

307.65

188.44

2

2

GSE2748

4.91

3.82

307.65

307.65

270.54

161.59

2.42

3

GSE7023

2.85

0.06*

307.65

307.65

307.65

182.73

2.83

4

GSE11024

3.79

3.82

307.65

224.15

175.48

68.07

3.25

5

TCGA_KIRP

0.61*

3.82

33.09

70.3

19.52

31.94

4.5

*p-value not significant after Bonferroni correction.

A heatmap of enriched pathways from maxP meta-analysis by combining p-values.

Figure 3: A heatmap of enriched pathways from maxP meta-analysis by combining p-values. A heatmap showing 115 enriched pathways detected by MAPE_I maxP under q-value = 0.2 threshold with the 1758 DE genes between PRCC1 and PRCC2 cases. Cell adhesion pathway is highlighted in red color.

Table 3: The number of differentially expressed genes in the datasets

Cut-off

GSE26574

GSE7023

GSE11024

TCGA_KIRP

GSE2748

minP

roP

maxP REM

p = 0.01

1371

642

891

2021

648

1389

1758

1587 2225

p = 0.05

2388

1402

1794

3126

1469

2322

2438

2162 2941

FDR = 0.01

546

78

106

1250

44

744

1347

1168 1837

FDR = 0.05

1314

377

609

2343

294

1558

1976

1758 2610

minimum P-value (minP), maximum P-value (maxP), rth ordered P-value (roP).

Comparative meta-analysis of differentially expressed genes between PRCC1 and PRCC2 samples.

Figure 4: Comparative meta-analysis of differentially expressed genes between PRCC1 and PRCC2 samples. (A) A Venn Diagram showing the overlap of DE genes detected by two meta-analyses, one combining the p-values (maxP) and another combining the effect sizes (REM). (B) PANTHER GO-Slim biological process analysis revealing 11 major functional categories associated with the 1475 DE genes. Biological adhesion is one of the major categories. (C) PANTHER pathway ontology analysis highlighting integrin pathway as the second-most enriched pathway. (D) A functional annotation clustering 2D report highlighting the cell adhesion related pathways associated with the integrin-related DE genes from the meta-analysis. (E) Protein-protein functional interaction network including 37 integrin pathway-associated genes extracted from the meta-analysis.

Differential regulation of the integrin pathway signaling between PRCC1 and PRCC2

The 1475 DE genes obtained from the overlay of p-value and ES combined meta-analysis were further examined using PANTHER gene ontology-slim biological process analysis. This analysis revealed that metabolic process and cellular process were among the top enriched categories (Figure 4B). PANTHER pathway ontology of the 1475 DE genes demonstrated 37 genes related to integrin pathway which was the second most enriched pathway (Figure 4C). These results were corroborated by similarly overlapping the DE genes from each of the two remaining meta-analyses based on p-values (minP and roP) with the DE genes obtained from REM meta-analysis. As expected, PANTHER pathway ontology analysis found integrin signaling pathway to be enriched within shared DE genes between REM and minP (1265 genes) and REM and roP (1300 genes) meta-analyses (Supplementary Figure S1C and S1E).

A functional annotation clustering 2D view report of the integrin pathway associated genes from the maxP/REM overlapping dataset was visualized by using DAVID functional annotation tool (Figure 4D). These genes included α5 (ITGA5)- and β8 (ITGB8)-integrins, both of which are RGD-binding receptors associated with epithelial-to-mesenchymal transition (EMT) [2122][23]. Moreover, α6-integrin (ITGA6), a laminin receptor that regulates epithelial cell polarity and growth was found to be differentially regulated between PRCC1 and PRCC2 (Table 4) [2425]. Finally, STRING interaction network analysis was performed with these 37 integrin pathway-related genes. This analysis revealed highly connected functional protein-protein interaction network that included all but one of the 37 DE genes (Figure 4E). To relate the differential expression of the three integrin genes in PRCC patients to healthy controls we extracted the expression data from two studies, the TCGA-KIRP study based on RNA-seq data which was also included in our meta-analysis and an independent study by Jones et al. (GSE15641) [7, 26]. This analysis revealed that ITGA5 was significantly downregulated in PRCC1 patients when compared with healthy controls (Figure 5A). In contrast, a tendency for modestly elevated levels of ITGA5 expression was noted for PRCC2 although this was not statistically significant (Figure 5A). ITGA5 expression levels in unclassified PRCC patients in the GSE15641 dataset tended to be higher than in healthy controls but this difference was not statistically significant (Figure 5B). Both studies were consistent with significant downregulation of ITGA6 expression that was particularly evident in PRCC2 patients (Figure 5C). The GSE15641 dataset also displayed a robust downregulation of ITGA6 in the PRCC samples (Figure 5D). In contrast to reduced ITGA6 levels, upregulation of ITGB8 expression was observed in PRCC patients with the highest expression levels seen in PRCC1 patients (Figure 5E). The GSE15641 dataset similarly showed significant upregulation of ITGB8 in PRCC patient samples (Figure 5F). Taken together, the combined meta-analysis extracted three candidate integrin genes whose differential regulation may underlie the different pathogenic properties of the two PRCC subtypes. Comparing with TCGA-KIRP study, the expression pattern of these three integrins predicts that the majority of the cancer samples in the GSE15641 dataset, where matched clinical information was not available, belong to the PRCC subtype 2.

Table 4: Differentially expressed genes of the integrin pathway in PRCC1 and PRCC2

PRCC2 > PRCC1

PTK2 SHC1 ACTN4 CAV1 COL1A1 COL4A1 FLNA ITGA5 MAP2K2 PARVA RHOD RRAS

PRCC2 < PRCC1

PTPN12 RHOA RHOB RAC1 ARF6 AKAP9 ACTR1B CSK BCAR1 FLNB KRAS RAF1 FN1 ITGA6 ITGB8 MAP2K3 MAP2K4 MAP3K1 MAP3K5 NTN4 LAMA5 LAMB1 LAMC2 MAPK10 MAPK7

PRCC1 = Papillary renal cell carcinoma type1, PRCC2 = Papillary renal cell carcinoma type2.

Analysis of the expression levels of ITGA5, ITGA6 and ITGB8 in healthy controls and in PRCC1 and PRCC2 patient samples.

Figure 5: Analysis of the expression levels of ITGA5, ITGA6 and ITGB8 in healthy controls and in PRCC1 and PRCC2 patient samples. (A) Relative expression levels of ITGA5 in the kidney from healthy control, PRCC1 and PRCC2 in the TCGA-KIRP dataset and (B) in healthy controls and PRCC patients from GSE15641 dataset where matched clinical information was not available [26]. (C) Relative expression levels of ITGA6 in the TCGA-KIRP and (D) the GSE15641 datasets. (E) Relative expression levels of ITGB8 in the TCGA-KIRP and (F) the GSE15641 datasets. Mann-Whitney U tests were performed to assess statistical significance of the observed differences between different sample groups in A–F. The horizontal lines represent median values.

Analysis of the gender specific differentially expressed genes in PRCC2 patients

Clinical data shows that PRCC2 is more frequent and tends to be more aggressive in male patients [27]. To address the possible underlying genetic factors we subjected male and female PRCC2 patient data for meta-analysis to study the DE genes. Given that information on gender was not available for two GEO datasets, the remaining three datasets (GSE11024, GSE2748 and TCGA-KIRP) were included for further analysis (Figure 1C). After removing PRCC1 cases the selected datasets consisted of 34 female and 73 male PRCC2 samples. To provide objective quality control and inclusion/exclusion criteria of the filtered datasets for meta-analysis, we performed MetaQC and runQC packages in R. All of the three datasets passed the quality control test and were included for further analysis (Figure 6A and Table 5). maxP meta-analysis detected 8 genes that were differentially regulated between female and male PRCC2 samples using a FDR cut-off under 0.05 (Figure 6B). The expression patterns of the 8 genes in the patient samples were visualized in a heat map (Figure 6C). The PANTHER GO-Slim biological process analysis revealed that metabolic process was the major functional category associated with the 8 genes (Figure 6D). STRING Interaction Network analysis of the extracted genes showed that RPS4Y1 interacted with translation elongation factor EEF1A2 (Figure 6E). However, RPS4Y1 is a Y-chromosome-specific gene and when a similar meta-analysis addressing the gender-dependent differences in gene expression was performed for PRCC1 patients where no gender bias has been reported, RPS4Y1 was also detected among the DE genes (Supplementary Figure 2). Therefore, it is unlikely that elevated RPS4Y1 expression is associated with the more aggressive pathogenesis of PRCC2 in male patients. In addition to RSP4Y1, the analysis of gender-specific DE genes in PRCC2 patients also highlighted interaction of two members of the membrane-spanning 4-domain family, subfamily A (MS4A) transmembrane (TMEM)-176A and 176B that have been found to be deregulated in multiple cancer types [28]. However, none of the genes were associated with cell adhesion or integrin pathways.

Comparative meta-analysis of DE genes in female and male PRCC2 patients.

Figure 6: Comparative meta-analysis of DE genes in female and male PRCC2 patients. (A) PCA biplot of quality control measures in the three PRCC studies after removing non-PRCC2 samples. (B) The detection competency curves of DE genes plotted against false-discovery rate in analyses of the three individual datasets and the maxP meta-analyses approaches. (C) A heat map representation of DE genes (FDR = 0.05) between female and male PRCC2 patients. (D) Molecular function analysis of the 8 DE genes revealed by the meta-analysis. (E) Protein-protein interaction network analysis of the 8 DE genes in female and male PRCC2 patients.

Table 5: Quality Control Result for PRCC2 based on gender difference

Number

Study

IQC

EQC

CQCg

CQCp

AQCg

AQCp

Rank

1

GSE2748g

20

3.52

7.83

24.43

5.16

4.59

1.5

2

GSE11024g

4

3.3

9.84

16.25

2.56

8.9

1.83

3

TCGA_KIRPg

2.67

3.7

0.11*

0.88*

0.16*

0.26*

2.67

*p-value not significant after Bonferroni correction.

DISCUSSION

High-throughput genomic datasets of PRCC patient material have accumulated in the public databases such as GEO and TCGA. Meta-analysis combining these datasets helps to increase the statistical power in mining and explaining the underlying mechanisms driving pathogenesis of PRCC subtypes. Here a meta-analysis was performed to extract the differentially expressed genes between PRCC1 and PRCC2 with a particular focus on the integrin-related pathways. The combined meta-analysis extracted 1475 DE genes, 37 of which were associated with integrin pathway. The more aggressive PRCC2 was associated with modest upregulation of ITGA5, a fibronectin receptor known to promote tumor progression in ccRCC [29]. Cell biological data supports a positive role for α5β1-integrin in promoting mesenchymal cell migration and cancer cell invasion [3031]. Ligation of α5β1-integrins have also been reported to support cell survival in suboptimal growth conditions [3233]. It is likely that upregulated ITGA5 expression in part contributes to the metastatic properties of PRCC2. ITGA6 expression was significantly downregulated in PRCC samples, especially in PRCC2, when compared with normal kidneys (Figure 5C and 5D). α6-integrin is a laminin receptor which delineates the basal membrane of polarized epithelial cells and synergizes with growth signals to support cell proliferation and migration [34]. Twist-mediated transformation of renal cancer cells led to reduced expression of ITGA6 [35]. Reduced expression of α6β4-integrin may enhance cancer cell dissemination from primary tumors and has been reported to positively correlate with prostate cancer progression [36]. However, increased expression of ITGA6 was reported for metastatic cells of many other solid cancer types [37][38]. Interestingly, α6-integrin is highly expressed in stem cells including cancer stem cells [39]. Our finding on decreased ITGA6 levels raise the possibility that enhanced dissemination of kidney cancer cells could in part explain the aggressive pathogenesis of PRCC2. However, further functional studies are required to confirm these findings. ITGB8 was the most robustly upregulated integrin in both PRCC datasets. Especially high levels of ITGB8 were seen in PRCC1 samples. ITGB8 forms an αVβ8-integrin heterodimer that plays a critical role in activation of latent TGF-β [40]. TGFβ-activation in turn drives epithelial-to-mesenchymal transition which may directly contribute to cancer cell migration and growth. TGFβ can function as both tumor suppressor and oncogene depending on other coinciding signals [41]. How elevated TGFβ-signaling might regulate PRCC pathogenesis is not known but it has been reported that TGFβ activation modulates metastatic properties of some RCCs [4243]. TGFβ activation appears to suppress growth of ccRCC tumors but the effect on PRCC remains poorly understood [44].

A meta-analysis was performed to identify gender-specific genes which may explain the different incidence and aggressiveness of PRCC2 between male and female patients. RPS4Y1, EEF1A2, TMEM176A and TMEM176B were preferentially up-regulated in male patients. RPS4Y1, a Y chromosome-linked gene [45], was found to functionally interact with a translation factor EEF1A2. While these proteins were linked to metabolic processes EEF1A2 may indirectly associate with cell adhesion related pathways as it has been shown to regulate MMP-9 expression and thereby influence migration and metastatic properties of pancreatic cancer cells [46]. Our meta-analysis approach revealed three key integrin subunits that were differentially associated with PRCC subtypes. An important limitation of the current study is that integrin activation and stability is frequently regulated at post-translational level. Therefore inclusion of proteomic data as well as functional cell biological studies are required for future studies to confirm and validate the specific roles of α5-, α6- and β8-integrins highlighted here by genetic meta-analysis.

MATERIALS AND METHODS

Data description

173 datasets were found in GEO series (http://www.ncbi.nlm.nih.gov/gds) by search term “papillary renal cell carcinoma” (Searched on 10th of February 2016). After manually removing datasets in which PRCC1 and PRCC2 had not been separately classified, four datasets were obtained (Table 1). For the GSE26574 dataset, 22 PRCC1 and 12 PRCC2 samples were analyzed by expression profiling microarray using a GPL11433 platform without gender information [47]. In GSE7023 dataset, 14 PRCC1 and 16 PRCC2 samples were analyzed by expression profiling microarray using GPL4866 platform without gender information [48]. In GSE11024 and GSE2748 datasets, gender information was available for a total of 25 PRCC1 and 21 PRCC2 samples which were analyzed by expression profiling microarray using GPL6671 and GPL570 platforms, respectively [4950]. In addition to microarrays datasets, one large PRCC dataset was available via TCGA-KIRP and consisted 77 PRCC1 and 86 PRCC2 with gender and clinical information [51]. This dataset had been analyzed by RNA-seq using Illumina Hiseq RNASeqV2-platform (http://cancergenome.nih.gov/).

Data processing

Capacity computing environment (Finland CSC Taito-shell application server: https://research.csc.fi) was used for running integrative bioinformatics pipelines. For the four GEO datasets, the raw gene expression data and clinical information was downloaded from GEO by using GEOquery package in R (BioConductor; https://www.bioconductor.org/). GEOquery is an open-source and open-development software project [52, 16]. TCGA-assembler was applied for downloading and processing the TCGA-KIRP gene expression raw data and clinical information [53, 17]. After removing non-PRCC samples from the datasets, the remaining data was saved to a comma-delimited text-file by using Excel as described previously [54]. Objective quality control analysis was performed by importing the processed raw datasets in log-format into R by using MetaDE [54]. Largest interquartile range (IQR) method was used to calculate gene-wise expression. Non-expressed (30%) and non-informative (30%) genes were filtered out. MetaQC was applied to include or exclude the processed datasets as described previously [55, 18, 56].

Data analysis

The MetaQC processed datasets were subjected to meta-analysis by combining p-values and effect sizes [54]. For the meta-analysis, the five datasets were combined and subjected to quality control and inclusion/exclusion criteria. These criteria consisted of internal quality control (IQC) index evaluating the homogeneity of coexpression and external quality control (EQC) index supervised by external pathway information, accuracy quality control indexes for genes (AQCg) and pathways (AQCp) and consistency of differential expression quality control (CQCg and CQCp) indexes that collectively depict the reproducibility and consistency of the data between the different individual the studies and results from the combined meta-analysis [57]. Three different meta-analysis methods for combining p-value; maxP, minP and roP were applied resulting in 1558, 1976 and 1758 differentially expressed (DE) genes, respectively (FDR = 0.05, Table 3) [20, 18]. For visualization, a heatmap plot of the DE genes (FDR = 0.05) was created by using MetaDE. A heatmap for pathway enrichment was provided by MetaPath packages, and MAPE_I was under q-value = 0.2 threshold [19]. Additionally, 2610 genes were obtained from a “REM” meta-analysis of the same raw data by combining effect sizes using MetaDE (FDR = 0.05) [58]. A Venn diagram was created to display an overlay of 1475 DE genes using Venny 2.1 (http://bioinfogp.cnb.csic.es/tools/venny/) [59]. The 1475 DE genes extracted from the overlay meta-analysis were further analyzed with DAVID 6.7 and PANTHER gene ontology to classify gene pathway, functions and interactions (https://david.ncifcrf.gov/) [6061]. Protein-protein interactions of the identified DE genes were also analyzed using STRING v10 online tool that visualizes known and predicted protein-protein interactions (http://string.embl.de/) [6263]. The gender comparison of PRCC2 samples were processed with MetaQC and MetaDE by combining the p-values of the raw data after PRCC1 cases were manually removed. PANTHER GO.Slim Molecular function analyses were performed as described earlier (http://pantherdb.org/) [64][65]. Expression analysis of ITGA5, ITGA6 and ITGB8 in Renal Cell Carcinomas and normal Kidney was obtained from Oncomine [26] and TCGA-KIRP.

ACKNOWLEDGMENTS

We thank Dr. Xingbin Wang and Dr. Dongwan D. Kang, Department of Human Genetics, University of Pittsburgh, for suggestions and instructions with the MetaOmics R packages.

CONFLICTS OF INTEREST

We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.

GRANT SUPPORT

This work was funded by Academy of Finland (140974, 263770, and 135560) and Biocenter Oulu.

REFERENCES

1. Eble JN, Sauter G, Epstein JI. SIA. (Eds.): WorldWorld Health Organization Classification of Tumours. Pathology and Genetics of Tumours of the Urinary System and Male Genital Organs. IARC Press. 2004.

2. Delahunt B, Eble JN. Papillary renal cell carcinoma: a clinicopathologic and immunohistochemical study of 105 tumors. Mod Pathol. 1997; 10:537–44.

3. Klatte T, Pantuck AJ, Said JW, Seligson DB, Rao NP, LaRochelle JC, Shuch B, Zisman A, Kabbinavar FF, Belldegrun AS. Cytogenetic and molecular tumor profiling for type 1 and type 2 papillary renal cell carcinoma. Clin Cancer Res. 2009; 15:1162–9.

4. Leroy X, Zini L, Leteurtre E, Zerimech F, Porchet N, Aubert J-P, Gosselin B, Copin M-C. Morphologic subtyping of papillary renal cell carcinoma: correlation with prognosis and differential expression of MUC1 between the two subtypes. Mod Pathol. 2002; 15:1126–30.

5. Sanders ME, Mick R, Tomaszewski JE, Barr FG. Unique patterns of allelic imbalance distinguish type 1 from type 2 sporadic papillary renal cell carcinoma. Am J Pathol. 2002; 161:997–1005.

6. Albiges L, Guegan J, Le Formal A, Verkarre V, Rioux-Leclercq N, Sibony M, Bernhard J-C, Camparo P, Merabet Z, Molinie V, Allory Y, Orear C, Couve S, 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. Clin Cancer Res. 2014; 20:3411–21.

7. Linehan WM, Spellman PT, Ricketts CJ, Creighton CJ, Fei SS, Davis C, Wheeler DA, Murray BA, Schmidt L, Vocke CD, Peto M, Al Mamun AAM, Shinbrot E, et al. Comprehensive Molecular Characterization of Papillary Renal-Cell Carcinoma. N Engl J Med. 2015; 374:135–45.

8. Lubensky IA, Schmidt L, Zhuang Z, Weirich G, Pack S, Zambrano N, Walther MM, Choyke P, Linehan WM, Zbar B. Hereditary and sporadic papillary renal carcinomas with c-met mutations share a distinct morphological phenotype. Am J Pathol. 1999; 155:517–26.

9. Wala SJ, Karamchandani JR, Saleeb R, Evans A, Ding Q, Ibrahim R, Jewett M, Pasic M, Finelli A, Pace K, Lianidou E, Yousef GM. An integrated genomic analysis of papillary renal cell carcinoma type 1 uncovers the role of focal adhesion and extracellular matrix pathways. Mol Oncol. 2015; 9:1667–77.

10. Manninen A. Epithelial polarity—generating and integrating signals from the ECM with integrins. Exp Cell Res. 2015; 334:337–49.

11. Seguin L, Desgrosellier JS, Weis SM, Cheresh DA. Integrins and cancer: regulators of cancer stemness, metastasis, and drug resistance. Trends Cell Biol. 2015; 25:234–40.

12. Jones J, Marian D, Weich E, Engl T, Wedel S, Relja B, Jonas D, Blaheta RA. CXCR4 chemokine receptor engagement modifies integrin dependent adhesion of renal carcinoma cells. Exp Cell Res. 2007; 313:4051–65.

13. Schrader AJ, Sevinc S, Olbert PJ, Hegele A, Varga Z, Hofmann R. [Gender-specific characteristics and survival of renal cell carcinoma]. Urologe A. 2008; 47:1182, 1184–6.

14. Sakano S, Hinoda Y, Okayama N, Kawai Y, Ito H, Nagao K, Hara T, Matsuyama H. Gender-specific association of methylenetetrahydrofolate reductase genotype and haplotype with the aggressiveness and prognosis of clear cell renal cell carcinoma in Japanese patients. BJU Int. 2010; 106:424–30.

15. Brannon AR, Haake SM, Hacker KE, Pruthi RS, Wallen EM, Nielsen ME, Rathmell WK. Meta-analysis of clear cell renal cell carcinoma gene expression defines a variant subgroup and identifies gender influences on tumor biology. Eur Urol. 2012; 61:258–68.

16. Davis S, Meltzer PS. GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics. 2007; 23:1846–7.

17. Zhu Y, Qiu P, Ji Y. TCGA-assembler: open-source software for retrieving and processing TCGA data. Nat Methods. 2014; 11:599–600.

18. Chang L-C, Lin H-M, Sibille E, Tseng GC. Meta-analysis methods for combining multiple expression profiles: comparisons, statistical characterization and an application guideline. BMC Bioinformatics. 2013; 14:368.

19. Shen K, Tseng GC. Meta-analysis for pathway enrichment analysis when combining multiple genomic studies. Bioinformatics. 2010; 26:1316–23.

20. Tseng GC, Ghosh D, Feingold E. Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 2012; 40:3785–99.

21. Mu D, Cambier S, Fjellbirkeland L, Baron JL, Munger JS, Kawakatsu H, Sheppard D, Broaddus VC, Nishimura SL. The integrin alpha(v)beta8 mediates epithelial homeostasis through MT1-MMP-dependent activation of TGF-beta1. J Cell Biol. 2002; 157:493–507.

22. Huang RY-J, Kuay KT, Tan TZ, Asad M, Tang HM, Ng AHC, Ye J, Chung VY, Thiery JP. Functional relevance of a six mesenchymal gene signature in epithelial-mesenchymal transition (EMT) reversal by the triple angiokinase inhibitor, nintedanib (BIBF1120). Oncotarget. 2015; 6:22098–113. doi: 10.18632/oncotarget.4300.

23. Schaffner F, Ray A, Dontenwill M. Integrin α5β1, the Fibronectin Receptor, as a Pertinent Therapeutic Target in Solid Tumors. Cancers (Basel). 2013; 5:27–47.

24. Kadoya Y, Kadoya K, Durbeej M, Holmvall K, Sorokin L, Ekblom P. Antibodies against domain E3 of laminin-1 and integrin alpha 6 subunit perturb branching epithelial morphogenesis of submandibular gland, but by different modes. J Cell Biol. 1995; 129:521–34.

25. Tennenbaum T, Belanger AJ, Glick AB, Tamura R, Quaranta V, Yuspa SH. A splice variant of alpha 6 integrin is associated with malignant conversion in mouse skin tumorigenesis. Proc Natl Acad Sci USA. 1995; 92:7041–5.

26. Jones J, Otu H, Spentzos D, Kolia S, Inan M, Beecken WD, Fellbaum C, Gu X, Joseph M, Pantuck AJ, Jonas D, Libermann TA. Gene signatures of progression and metastasis in renal cell cancer. Clin Cancer Res. 2005; 11:5730–9.

27. Courthod G, Tucci M, Di Maio M, Scagliotti G V. Papillary renal cell carcinoma: A review of the current therapeutic landscape. Crit Rev Oncol Hematol. 2015; 96:100–12.

28. Cuajungco MP, Podevin W, Valluri VK, Bui Q, Nguyen VH, Taylor K. Abnormal accumulation of human transmembrane (TMEM)-176A and 176B proteins is associated with cancer pathology. Acta Histochem. 2012; 114:705–12.

29. Hase H, Jingushi K, Ueda Y, Kitae K, Egawa H, Ohshio I, Kawakami R, Kashiwagi Y, Tsukada Y, Kobayashi T, Nakata W, Fujita K, Uemura M, et al. LOXL2 status correlates with tumor stage and regulates integrin levels to promote tumor progression in ccRCC. Mol Cancer Res. 2014; 12:1807–17.

30. Mierke CT, Frey B, Fellner M, Herrmann M, Fabry B. Integrin α5β1 facilitates cancer cell invasion through enhanced contractile forces. J Cell Sci. 2011; 124:369–83.

31. Ganguly KK, Pal S, Moulik S, Chatterjee A. Integrins and metastasis. Cell Adh Migr. 7:251–61.

32. Lee JW, Juliano RL. α5β1 integrin protects intestinal epithelial cells from apoptosis through a phosphatidylinositol 3-kinase and protein kinase B-dependent pathway. Mol Biol Cell. 2000;11:1973–87.

33. Matter ML, Ruoslahti E. A signaling pathway from the α5β1 and α(v)β3 integrins that elevates bcl-2 transcription. J Biol Chem. 2001; 276:27757–63.

34. Yoon S-O, Shin S, Lipscomb EA. A novel mechanism for integrin-mediated ras activation in breast carcinoma cells: the α6β4 integrin regulates ErbB2 translation and transactivates epidermal growth factor receptor/ErbB2 signaling. Cancer Res. 2006; 66:2732–9.

35. Zhu T, Zhang H, Tao J, Sheng L, Hu X, Rong R, Xu M. Twist2 promotes kidney cancer cell proliferation and invasion by regulating ITGA6 and CD44 expression in the ECM-receptor interaction pathway. Onco Targets Ther. 2016; 9:1801.

36. Schmelz M, Cress AE, Scott KM, Bürger F, Cui H, Sallam K, McDaniel KM, Dalkin BL, Nagle RB. Different phenotypes in human prostate cancer: α6 or α3 integrin in cell-extracellular adhesion sites. Neoplasia. 4:243–54.

37. Rabinovitz I, Nagle RB, Cress AE. Integrin alpha 6 expression in human prostate carcinoma cells is associated with a migratory and invasive phenotype in vitro and in vivo. Clin Exp Metastasis. 1995; 13:481–91.

38. Jones JL, Royall JE, Critchley DR, Walker RA. Modulation of myoepithelial-associated alpha6beta4 integrin in a breast cancer cell line alters invasive potential. Exp Cell Res. 1997; 235:325–33.

39. Goel HL, Gritsko T, Pursell B, Chang C, Shultz LD, Greiner DL, Norum JH, Toftgard R, Shaw LM, Mercurio AM. Regulated splicing of the α6 integrin cytoplasmic domain determines the fate of breast cancer stem cells. Cell Rep. 2014; 7:747–61.

40. Pozzi A, Zent R. TGF-β sequestration by mesangial cell integrin αvβ8: A novel mechanism of glomerular endothelial cell regulation. Am J Pathol. 2011; 178:485–9.

41. Principe DR, Doll JA, Bauer J, Jung B, Munshi HG, Bartholin L, Pasche B, Lee C, Grippo PJ. TGF-β: duality of function between tumor prevention and carcinogenesis. J Natl Cancer Inst. 2014; 106:djt369.

42. Kominsky SL, Doucet M, Brady K, Weber KL. TGF-beta promotes the establishment of renal cell carcinoma bone metastasis. J Bone Miner Res. 2007; 22:37–44.

43. Copland JA, Luxon BA, Ajani L, Maity T, Campagnaro E, Guo H, LeGrand SN, Tamboli P, Wood CG. Genomic profiling identifies alterations in TGFbeta signaling through loss of TGFbeta receptor expression in human renal cell carcinogenesis and progression. Oncogene. 2003; 22:8053–62.

44. Boström A-K, Lindgren D, Johansson ME, Axelson H. Effects of TGF-β signaling in clear cell renal cell carcinoma cells. Biochem Biophys Res Commun. 2013; 435:126–33.

45. Bergen AW, Pratt M, Mehlman PT, Goldman D. Evolution of RPS4Y. Mol Biol Evol. 1998; 15:1412–9.

46. Xu C, Hu D, Zhu Q. eEF1A2 promotes cell migration, invasion and metastasis in pancreatic cancer by upregulating MMP-9 expression through Akt activation. Clin Exp Metastasis. 2013; 30:933–44.

47. Ooi A, Wong J-C, Petillo D, Roossien D, Perrier-Trudova V, Whitten D, Min BWH, Tan M-H, Zhang Z, Yang XJ, Zhou M, Gardie B, Molinié V, et al. An antioxidant response phenotype shared between hereditary and sporadic type 2 papillary renal cell carcinoma. Cancer Cell. 2011; 20:511–23.

48. Furge KA, Chen J, Koeman J, Swiatek P, Dykema K, Lucin K, Kahnoski R, Yang XJ, Teh BT. Detection of DNA copy number changes and oncogenic signaling abnormalities from gene expression data reveals MYC activation in high-grade papillary renal cell carcinoma. Cancer Res. 2007; 67:3171–6.

49. Kort EJ, Farber L, Tretiakova M, Petillo D, Furge KA, Yang XJ, Cornelius A, Teh BT. The E2F3-Oncomir-1 axis is activated in Wilms’ tumor. Cancer Res. 2008; 68:4034–8.

50. Yang XJ, Tan M-H, Kim HL, Ditlev JA, Betten MW, Png CE, Kort EJ, Futami K, Furge KA, Takahashi M, Kanayama H-O, Tan PH, Teh BS, et al. A molecular classification of papillary renal cell carcinoma. Cancer Res. 2005; 65:5628–37.

51. TCGA. Comprehensive Molecular Characterization of Papillary Renal-Cell Carcinoma. N Engl J Med. 2016; 374:135–45.

52. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004; 5:R80.

53. Zhu Y, Qiu P, Ji Y. TCGA-assembler: open-source software for retrieving and processing TCGA data. Nat Methods. 2014; 11:599–600.

54. Wang X, Kang DD, Shen K, Song C, Lu S, Chang L-C, Liao SG, Huo Z, Tang S, Ding Y, Kaminski N, Sibille E, Lin Y, et al. An R package suite for microarray meta-analysis in quality control, differentially expressed gene analysis and pathway enrichment detection. Bioinformatics. 2012; 28:2534–6.

55. Li J, Tseng GC. An adaptively weighted statistic for detecting differential gene expression when combining multiple transcriptomic studies. Ann Appl Stat. 2011; 5:994–1019.

56. Song C, Tseng GC. Hypothesis setting and order statistic for robust genomic meta-analysis. Ann Appl Stat. 2014; 8:777–800.

57. Kang DD, Sibille E, Kaminski N, Tseng GC. MetaQC: objective quality control and inclusion/exclusion criteria for genomic meta-analysis. Nucleic Acids Res. 2012; 40:e15.

58. Choi JK, Yu U, Kim S, Yoo OJ. Combining multiple microarray studies and modeling interstudy variation. Bioinformatics. 2003; 19.

59. Oliveros J. enny. An interactive tool for comparing lists with Venn’s diagrams. [Internet].

60. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009; 4:44–57.

61. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009; 37:1–13.

62. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, Kuhn M, Bork P, Jensen LJ, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015; 43:D447–52.

63. von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, Jouffre N, Huynen MA, Bork P. STRING: known and predicted protein-protein associations, integrated and transferred across organisms. Nucleic Acids Res. 2005; 33:D433–7.

64. Mi H, Poudel S, Muruganujan A, Casagrande JT, Thomas PD. PANTHER version 10: expanded protein families and functions, and analysis tools. Nucleic Acids Res. 2015; 44:D336–42.

65. Mi H, Muruganujan A, Casagrande JT, Thomas PD. Large-scale gene function analysis with the PANTHER classification system. Nat Protoc. 2013; 8:1551–66.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 12390