Skip to main content
Download PDF

An IL6-correlated signature in serous epithelial ovarian cancer associates with growth factor response

BMC Genomics volume 14, Article number: 508 (2013) Cite this article

Abstract

Background

Epithelial ovarian cancer (EOC) is one of the most lethal gynecological cancers; the majority of EOC is the serous histotype and diagnosed at advanced stage. IL6 is the cytokine that has been found most frequently associated with carcinogenesis and progression of serous EOCs. IL6 is a growth-promoting and anti-apoptotic factor, and high plasma levels of IL6 in advanced stage EOCs correlate with poor prognosis. The objective of the present study was to identify IL6 co-regulated genes and gene network/s in EOCs.

Results

We applied bioinformatics tools on 7 publicly available data sets containing the gene expression profiles of 1262 EOC samples. By Pearson's correlation analysis we identified, in EOCs, an IL6-correlated gene signature containing 40 genes mainly associated with proliferation. 33 of 40 genes were also significantly correlated in low malignant potential (LMP) EOCs, while 7 genes, named C5AR1, FPR1, G0S2, IL8, KLF2, MMP19, and THBD were IL6-correlated only in advanced stage EOCs. Among the 40-gene signature EGFR ligand HBEGF, genes of the EGR family members and genes encoding for negative feedback regulators of growth factor signaling were included. The results obtained by Gene Set Enrichment and Ingenuity Pathway Analyses enabled the identification, respectively, of gene sets associated with ‘early growth factor response’ for the 40-gene signature, and a biological network related to ‘thrombosis and cardiovascular disease’ for the 7-gene signature. In agreement with these results, selected genes from the identified signatures were validated in vitro by real time RT-PCR in serous EOC cell lines upon stimulation with EGF.

Conclusions

Serous EOCs, independently of their aggressiveness, co-regulate IL6 expression together with that of genes associated to growth factor signaling, arguing for the hypothesis that common mechanism/s driven by EGFR ligands characterize both advanced-stage and LMP EOCs. Only advanced-stage EOCs appeared to be characterized by a scenario that involves genes which are so far associated with thrombosis and cardiovascular disease, thus suggesting that this pathway is implicated in the growth and/or spread of more aggressive tumors. We have discovered novel activated signaling pathways that drive the expression of IL6 and of co-regulated genes and are possibly involved in the pathobiology of EOCs.

Background

Epithelial ovarian cancer (EOC) is the second most common and the most deadly malignancy of the female reproductive tract. Serous, endometrioid, clear-cell, and mucinous ovarian cancers are the four most common histotypes [1]. The majority of EOCs are diagnosed at stage III and IV when the tumor cells are spread in the peritoneum along with the presence of malignant ascites. The serous histotype accounts for about 80% of EOCs, and the majority show an inactivating mutation of the tumor suppressor gene TP53. Low malignant potential (LMP) serous EOCs are thought to arise by the transformation of tumors of borderline malignancy, and activating mutations in members of the RAS pathway (KRAS, BRAF, and ErbB2) are found in the majority of these tumors [2]. LMP EOCs show a relatively high growth capacity, are usually not invasive but resistant to conventional chemotherapy [1].

A number of studies suggest that factors related to the inflammation of the ovarian surface epithelium (OSE) such as ovulation, endometriosis, and pelvic inflammatory diseases are associated with an increased risk for EOC [3]. The most important hypothesis regarding EOC carcinogenesis is the ovulation theory, which relates the risk of ovarian cancer to incessant ovulation. Recently, it has been hypothesized that high grade serous, endometrioid and clear cell ovarian cancers arise from the fallopian tube epithelium and share a common pathogenic mechanism, i.e., iron-induced oxidative stress derived from retrograde menstruation [4]. Both the incessant ovulation and oxido-reductive fallopian tube epithelial damage hypotheses have provided evidence that inflammatory responses induced under physiological conditions may foster the development of EOC. In accordance with these hypotheses of ovarian tumorigenesis, a number of cyto/chemokines has been found at detectable levels in ascites from EOC patients [5]. Among those molecules, IL6 is the cytokine that has been most frequently associated with EOC carcinogenesis and progression [6]. Preclinical evidence has shown that IL6 enhances tumor cell survival and increases resistance to chemotherapy via JAK/STAT signaling in tumor cells [7] and IL6 receptor alpha trans-signaling on tumor endothelial cells [8, 9]. In addition, IL6 has pro-angiogenic properties [7], regulating immune cell infiltration, a stromal reaction, and the tumor-promoting actions of Th17 lymphocytes [10]. In patients with advanced disease, high plasma levels of IL6 correlate with poor prognosis [11] and elevated levels are also present in malignant ascites [12]. Treatment of EOC cells with the anti-IL6 antibody (Ab) siltuximab has been shown to reduce constitutive cyto/chemokine production and inhibit IL6 signaling, tumor growth, the tumor-associated macrophage infiltrate, and angiogenesis in IL6–producing intraperitoneal ovarian cancer xenografts [13]. IL6 stimulates inflammatory cytokine production, tumor angiogenesis and the tumor macrophage infiltrate in ovarian cancer and these actions can also be inhibited by a neutralizing anti-IL6 Ab in clinical studies [14]. However, further knowledge on IL6-expressing EOCs is needed to select patients who are possibly responsive to IL6-dependent therapies.

We have recently found that IL6 can be co-expressed together with plasminogen activator inhibitor (PAI)-1, encoded by SERPINE1, in a subset of advanced stage serous EOCs due to the activation of the ligand-dependent EGFR/NFkB signaling cascade [15]. Ex vivo, using 23 EOCs from advanced-stage patients with malignant ascites at surgery, we observed co-expression of EGFR, IL6, and PAI-1 in 57% of primary tumors and concomitant expression of both IL6 and PAI-1 in the corresponding ascites. Computational analysis on four publicly available data sets of EOC gene expression showed a correlation between the expression of the IL6 and SERPINE1 genes in advanced stage EOC patients, which in one case was associated with shorter progression-free survival [15]. These results further highlight the involvement of IL6 in the progression of EOC.

Herein, to give further insight in the biology of IL6-expressing serous EOC we utilized a bioinformatics approach, described in the flowchart of Figure 1, to identify IL6 co-regulated genes and signaling pathway/s in which they are involved. First, we identified a list of genes representing a molecular signature for both advanced-stage and LMP serous EOC which recapitulate the so-called ‘early growth factor response’. We also identified an IL6-correlated signature of seven genes involved in vascular thrombosis specific for advanced-stage serous EOCs.

Figure 1
figure 1

Flowchart describing the analysis workflow.

Full size image

Results

IL6 expression significantly correlates with a defined gene set in advanced stage serous EOCs

Pearson’s correlation analysis of seven data sets containing the expression profiles of 1262 samples from serous EOCs (Table 1) was performed to identify genes whose expression was significantly correlated with IL6 expression in each data set. Correlation scores of each gene pair were computed using the R program essentially as described [16–18]. For genes represented by multiple probes in the same array format, the probe with the highest correlation to IL6 in the data set with the highest number of patients was chosen and considered for the other data sets when present (Additional file 1a). This analysis allowed the identification of genes whose expression positively correlated with IL6 along the seven data sets with a p-value ≤0.05 and a Pearson’s correlation coefficient (r) exceeding 0.4 (Additional file 1b). A further analysis across the seven data sets yielded 40 concordant correlated genes in at least four data sets (Additional file 1b and Figure 2). Of note, 38 of 40 genes were correlated in data set I (204 samples) obtained with an Affymetrix platform, and in data set VII (110 samples) obtained with an Agilent platform. Among the identified IL6-correlated genes CXCL2, HBEGF, SERPINE1, DUSP1, ZFP36, and IER3 were common to all data sets. The correlation between IL6 and HBEGF, an EGFR ligand, and SERPINE1, encoding PAI-1, is in agreement with our previously published results on co-expression of IL6 and PAI-1 in high grade EOCs due to EGFR activation [15]. The majority of genes are associated to the biological process ‘proliferation’ (50%) (Table 2). Among the genes associated with proliferation, there were a number of growth factor early response genes (EGR1, EGR3, NR4A1, FOSB, IER3). The IL6-correlated signature also included genes associated with ‘inflammation’ (20%), and the remaining genes were associated with ‘cell cycle and apoptosis’, ‘metabolism’ and ‘migration and invasion’.

Table 1 List of EOC data sets of gene expression analyzed in the present study
Full size table
Figure 2
figure 2

Heatmap of IL6-correlated genes. The heatmap of Pearson’s correlation coefficient (r) of the genes with IL6 was drawn by using R programming language. The r scores are represented in grayscale as reported in the color key. IL6 self-correlation was artificially set to the maximum score. Correlation score below 0.4 were considered not significant (NS). Genes not spotted on the array were defined NA (not available). The number of data sets in which the gene resulted significantly correlated with IL6 is reported on the right.

Full size image
Table 2 Biological functions of the IL6-correlated genes
Full size table

Thus, we identified a gene signature of IL6 correlated genes in serous EOC containing mainly proliferation-associated genes.

Advanced stage EOC-specific IL6-correlated gene signature functionally associated with control of cell morphology and cardiovascular disease

Next, to determine whether the identified gene signature was specific for advanced stage EOC or could also be associated with LMP EOC, Pearson’s correlation analysis to IL6 was applied to gene expression data of LMP EOCs reported in data sets I, II, and IV (see Table 1). The density plot of IL6 intensities showed a similar trend of expression in advanced stage and LMP EOCs (Additional file 2). The data obtained comparing advanced stage EOCs and LMP EOCs were similar in the three data sets, and were more reliable (number of cases, genes identified, significance level) in data set I. Among the above identified advanced-stage EOCs IL6-correlated genes, 33 were also significantly correlated in LMP EOCs, while 7 genes (C5AR1, FPR1, G0S2, IL8, KLF2, MMP19, and THBD) were specific for advanced-stages only (Additional file 3). Among these genes, IL8 has already been associated with aggressiveness and progression of malignant EOC [27], while the others have not previously been associated with EOC biology and clinical outcome.

To provide insight into the possible biological significance of the 40-gene signature, functional analysis of positively correlated genes (41, including IL6) was carried out by Ingenuity Pathway Analysis software (IPA) [28]. The top two functions (N1 and N2), associated with the highest score network, were ‘Cell death, cellular function and maintenance, hematological system development and function’ and ‘Cell death, cellular development, cellular growth and proliferation’ (Figure 3 and Additional file 4). When IPA analysis was performed on the seven-gene signature specific for advanced stage EOCs, the top function, associated with the highest score network, was ‘Cell morphology, cell function, cardiovascular disease’. As shown in Figure 3 and listed in Additional file 4, all seven genes are included in this network (N3) together with genes already known to have a role in the progression of EOCs such as VEGF, the receptor tyrosine kinases EGFR and HER2, and the PI3K complex [1]. In addition to the input genes, it is noteworthy that IL6 is not present in the identified networks (Additional file 4), but when added manually to each network establishes a connection with some of the correlated genes (Figure 3).

Figure 3
figure 3

Graphical representation of the top score networks identified by IPA. Molecular interactions between IL6-correlated genes in at least four data sets are reported. The top two networks (N1 and N2) were identified by loading all IL6-correlated genes in Fig. 1. Network 3 (N3) was identified by loading the 7 IL6-correlated genes specific for advanced stage EOCs. IL6 (highlighted in blue) was manually added to each network. IL6-correlated genes are highlighted in red and the intensity indicates the number of data sets where the gene is correlated. The name of the network is reported below the graph.

Full size image

The NFkB complex was included in networks N1 and N3 (Figure 3), highlighting its possible pivotal role in EOC progression.

It thus appears that two signatures are related to IL6 in EOCs: a 33-gene signature common to advanced stage and LMP EOCs and associated to control of cell growth and death, while the 7-gene signature, associated only to advanced-stage EOCs likely presenting NFkB transcriptional activation, might be a determinant of tumor aggressiveness, and may be associated with a pathway regulating vascular thrombosis.

An IL6-correlated gene set recapitulates the early growth factor response

To give further insight in the biology of EOCs expressing IL6, a GSEA [29] analysis was performed for each data set listed in Table 1. By the “use a gene (IL6) as phenotype” analysis, GSEA first ranks the genes according to their correlation to IL6. It then determines whether a priori defined set of genes, in this instance those belonging to the C2 curated catalogue of functional gene sets, are randomly distributed throughout the gene list or primarily found at the top or bottom. Common significant gene sets obtained from GSEA analysis of the two largest data sets (I and VI) were selected and analyzed in the other datasets. This yielded 20 significantly enriched gene sets for all datasets. Normalized enrichment scores and FDR values in the different datasets are listed in Table 3. A literature search was conducted to identify signaling pathways previously implicated in the progression of EOCs and/or in epithelial–derived malignancies. Among the most significant gene sets, the BILD_KRAS_ONCOGENIC_SIGNATURE [21] which includes genes whose expression is induced by the activation of H-RAS oncogene, was originally derived from the herein named data set III and can be considered a positive control. Three additional gene sets, AMIT_EGF_RESPONSE_60_HELA, AMIT_EGF_RESPONSE_120_HELA [30] and NAGASHIMA_NRG1_SIGNALING_UP [31], were considered possible candidates of signaling pathways associated with EOC, and are associated with ‘growth factor response’. These gene sets comprise early response genes, i.e. the EGR family members, and the negative feedback regulators of the growth factor signaling, i.e. ZFP36 and KLF2. The fifth selected gene set, named KIM_WT1_TARGET_UP in some ways also recapitulates the growth factor response, since among WT1 target genes the EGF family ligands EREG, AREG and HBEGF are included [32].

Table 3 Significant IL6 correlated gene sets identified by GSEA analysis.
Full size table

Furthermore, among the WT1 target genes SERPINE1 was also identified in the same study. Enrichment plots related to the above described gene-sets in data set I are shown in Figure 4. It is of note that IL6 is not included in the selected gene sets (Additional file 5) as well as other genes that are included in network 3 identified by analysis using IPA. Based on the results obtained by the above-described computational analysis and on our recent demonstration that IL6 is up-modulated in EOC cells upon EGF stimulation in time-dependent manner [15], in vitro validation of 12 genes selected from the IL6-correlated gene sets was performed with real time RT-PCR using total RNA from EGF-stimulated serous EOC cell lines (Figure 5). The IL6 was up-modulated in all EOC cells analyzed upon EGF stimulation. Concordantly, 75%, 58%, and 75% of the gene transcripts were up-modulated in IGROV1, OAW42, and SKOV3 cells, respectively (Figure 5). Among the correlated genes common to 7 data sets (see Figure 2), CXCL2, HBEGF, SERPINE1 and DUSP1 were increased in all three EOC cell lines analyzed. Additionally, NR4A1, a correlated gene in 6 data sets, was up-modulated upon EGF stimulation in all EOC cells. THBD and KLF2 transcripts, associated with ’Cardiovascular disease’ by IPA analysis, were up-modulated in 2 of 3 EGF-stimulated EOC cells. In contrast, the MMP19 transcript, whose relevant protein is associated with invasion and tumor progression [33], was not up-modulated in EGF-stimulated EOC cells. Interestingly, in non-transformed ovary cells, named IOSE- HTERT64 [34], although IL6 was slightly up-modulated by EGF stimulation, only 25% of the transcripts analyzed were up-modulated.

Figure 4
figure 4

GSEA enrichment plots for the five gene sets enriched in EOCs. On the top of each plot, the name of the gene set is reported. For each gene set, the enrichment plot was extracted from the GSEA output results and each gene set showed significant enrichment in IL6 expressing advanced stage EOC (FDR Q value = 0.0; Fig. 3). Genes with higher expression in IL6-positive tumors have higher enrichment scores, and are therefore plotted on the left side of the graph, whereas those with lower expression in IL6-positive tumors have lower enrichment scores and are plotted on the right side of the graph. The bottom portion of the plot shows the value of the ranking metric moving down the list of ranked genes. A positive ranking metric indicates that a gene is correlated with the IL6 positive phenotype. The results from dataset 1 are reported.

Full size image
Figure 5
figure 5

In vitro validation of selected IL6 correlated genes. Real time RT-PCR on selected IL6-correlated genes was performed using total RNA of starved EOC cell lines untreated (white bars) or treated (grey bars) for 4 hr (IGROV1, OAW42 and IOSE 64 hTERT) or 8 hr (SKOV3) with EGF (20 ng/ml). The number of data sets in which the gene resulted significantly correlated with IL6 is reported on the bottom. Data are mean values (± SD) presented as relative expression normalized for GAPDH mRNA levels. Asterisks indicate significant positive variations (Student’s t test).

Full size image

These data indicate that ligand-dependent EGFR activation in serous EOC cells induces the transcription of genes correlated with IL6 expression.

Discussion

Microarray technology has developed very rapidly, and it has become relatively easy to analyze the expression levels of thousands of genes within cancer cells. However, genes do not act in isolation, but each acts in complexes and builds networks and activated pathways that ultimately give rise to a specific cell phenotype. Thus, the search of co-regulated genes applying bioinformatics approaches may spread light on the biology of a tumor and its development. Previously, by applying this kind of ‘in silico’ approach on gene expression profiles of ovarian and thyroid carcinomas [16, 17] and melanomas [18], we have been able to identify novel signaling pathways activated in those tumors. The present study, by applying similar bioinformatics tools, highlights possible novel signaling pathways activated in IL6-expressing EOCs. Among those, growth factor-dependent signaling was also experimentally validated in vitro in selected cellular EOC models.

First of all, Pearson’s correlation analysis allowed the identification of genes co-regulated with IL6 in aggressive EOC providing evidence that co-regulated genes can encode proteins involved in common signaling pathways. To identify IL6-coregulated genes we adopted thresholds which allowed to obtain a good balance among the statistical significance, the strength of the correlation and the biological reproducibility. Furthermore, we performed the analysis on 7 different data sets, containing the gene expression profiles of more than 1200 EOC samples, obtained on different array platforms, to increase the robustness on the bioinformatics results. We found a gene signature common to both advanced stage and LMP serous EOCs, and another 7-gene signature specific for advanced stage EOCs. The integration of the results obtained by IPA and GSEA, allowed us to determine that all EOCs, independently of their aggressiveness, co-regulate IL6 together with genes associated with cell growth and early growth factor response, arguing for the hypothesis of common mechanism/s of transformation. Only advanced-stage EOCs appeared to be characterized by a scenario that involves genes such as FPR1, KLF2 and THBD, to date associated with thrombosis and cardiovascular disease, thus suggesting that this pathway contributes to the growth and/or the spread of this type of tumor.

Our data indicate the existence of a biological interaction between IL6 expression and that of the co-regulated genes as resulted upon IPA and GSEA analyses. On the other hand, since knockdown of IL-6 by specific siRNA did not affect the amount of the transcripts analyzed (data not shown), the regulation of the expression of the identified genes appeared not directly dependent to that of IL6. Although the IL6 gene was not associated with the networks identified by IPA or with the gene sets selected by GSEA, these results are in agreement with our previous observations in a subset of advanced stage EOC where ligand-dependent EGFR activation induced NFkB-dependent transcription of IL6 together with PAI-1, encoded by the SERPINE1 gene [15]. NFkB also emerged to be a possible transcriptional regulator of 13 out of 40 genes according to the reported informations [35], data which might further indicate that a growth factor-dependent NFkB signaling is activated in a subset of EOC. It is noteworthy that IL6 and 19 of the 40 correlated genes were found up-modulated upon 2 hr serum stimulation of quiescent keratinocytes [36]. We can therefore argue that the activation of growth factor activated signaling can either directly or indirectly induce the expression of IL6 and genes which likely play a role in the growth of EOCs. Furthermore, this growth factor-induced signaling pathway induces positive regulators of cellular function that are in turn regulated by negative feedback regulators such as ZFP36 and KLF2 [30]. The HBEGF gene, encoding for an EGFR ligand, was also highly significantly correlated in all seven data sets analyzed, indicating the prevalence of ligand-dependent EGFR activation. The regulation of growth factor signaling pathways by negative feedback is a universal mechanism for limiting the duration and intensity of signaling output. While negative feedback is a key component of normal cellular signaling, its role in cancer cells is more complex. Indeed, the loss of some negative feedback regulators might contribute to tumor progression, but might also be expressed at considerably higher levels in oncogene-mutant tumors as observed in BRAF-mutated melanomas [37]. Interestingly, the presence of a feedback negative mechanism has also been associated with greater efficacy of growth factor receptor-targeted therapy [38]. The fact that in EOC cells, with active EGFR/NFkB/IL6 signaling, EGFR-targeted therapy was more effective might be due to the up-regulation of feedback negative regulators of growth factor signaling [15]. Taken together, these data suggest that the IL6-associated signature might have a translational impact helping to select EOC patients who are likely responsive to EGFR-targeted therapy. Experiments are now ongoing to verify this hypothesis.

Nuclear expression of the Wilm’s tumor suppressor is found in OSE cells and in the majority of serous EOCs [39]. However, the corresponding gene, named WT1, has been never associated with IL6 gene expression. WT1 is required for kidney development, and the report in which the relevant gene set has been derived particularly emphasized the finding that the genes encoding the EGF family ligands EREG, AREG, and HBEGF may be transcriptionally regulated by WT1, orchestrating a fine-tuning of the EGF signaling pathway [32]. Altogether these observations support that the EGF signaling pathway is pivotal in the biology of EOC.

The gene signature common to advanced stage serous and LMP EOCs is not unexpected if one considers the theory that LMP EOCs derive from serous low grade EOC with a borderline morphologic phenotype [40]. However, if this is the case, advanced-stage and LMP tumors might share common genetic alterations that induce aberrant growth. In addition, in vitro validation experiments performed on gene transcripts of non-transformed surface ovary cells argue for the notion that the signature associated with a growth factor response is not expressed and/or EGF-dependent in normal ovary cells.

The bioinformatics approach also produced hypothesis-generating results. The association of the 7-gene signature with advance stage EOCs is novel. At present only angiogenesis-related genes and proteins, such as VEGF and its receptor, have a well documented role in EOC biology and are already well-exploited targets in the therapy of more aggressive EOCs [39]. Our findings open new questions on the role of genes associated with thrombosis and cardiovascular disease in the progression of EOCs. It has been recently hypothesized that low-dose aspirin as antithrombotic therapy may inhibit progression rather than the induction of EOC [41]. Indeed, aspirin and selective COX inhibitors could reduce progression not only by inhibiting prostaglandin production, thus reducing inflammation, but also by negatively modulating thrombosis-associated genes. Therefore, the inhibition of both pathways synergistically might be an interesting approach to block the growth and dissemination of advanced stage EOCs.

MMP19, a gene of the 7-gene signature specific for malignant EOCs and encoding the metallo-protease (MMP) 19 was present in network 3 of IPA analysis, but was not in any of the gene sets selected by GSEA analysis. MMPs are key molecules of tumor cell invasion, including EOCs [42] and, since the majority of samples were advanced stage EOCs, MMP19 could be a new player in the dissemination of these tumors and experiments are now ongoing to test its presence and role in advanced stage EOCs.

Conclusions

By applying a bioinformatics approach we identified genes co-regulated with IL6 expression in clinically-relevant subtypes of EOC, their interactions in networks and pathways as well as their functional association to growth factor response. IL6 gene expression together with that of the correlated gene signature could help identifying EOC patient’s subgroups in which the identified signaling pathways might be biologically relevant during the progression of the disease and, in the long term, might represent new pharmacological targets.

Methods

Computational analysis

Seven EOC data sets, six arrayed on Affymetrix platforms and one on an Agilent platform, were analyzed (Table 1). Raw data of data sets I, II, and III [19–21] were downloaded from the NCBI Gene Expression Omnibus (GEO) repository (IDs GSE9891, GSE12172 and GSE3149, respectively) and those of data set VI were downloaded from the proprietary repository [43]. Data sets IV and V [23, 24] were downloaded from the Duke Institute website as suggested in the original publications. The raw data from Affymetrix were normalized through the RMA method using the Expression Console software developed by Affymetrix. Upon quality control, probes were annotated with the current annotation files (version 32) for the proper array format. Normalized data of data set VII [25], obtained on Agilent platform, were downloaded from GEO (ID GSE17260).

For each data set, the expression data from serous histotype cases were selected. Since in all but one (IV), data sets the percent of cases at early stage (I-II) ranged from 0 to 10%, no stage selection was applied; in the case of data set IV, in which stage I and II represented 50% of case material, to avoid difficulty in comparison with the others, only advanced stages (III-IV) were selected. According to these selection criteria, we considered our overall case material to be composed of advanced stage EOC. Each data set was analyzed separately and the gene expression intensity of IL6, represented by a single probe in all the analyzed array formats, was correlated to the remaining probes across all EOCs samples in the array. The Pearson’s correlation coefficients (r), p and FDR values were calculated using cor, cor.test and p.adjust (using the Benjamini & Hochberg method) functions, respectively, from the Stats package in R programming language (version 2.12.0). For genes represented by multiple probes in the same array format, the probe with the highest correlation to IL6 in the data set with the highest number of patients was chosen and considered for the other data sets when present. Only genes exhibiting a p value ≤ 0.05 and r ≥ 0.4 in at least 4 of the 7 data sets were considered significant (Additional file 1). In three studies (I, II and IV), serous LMP EOCs were also profiled and their expression data analyzed as described above. Correlation values to IL6 corresponding only to the list of genes significant in advanced stage EOC were further considered (Additional file 1). IPA (Ingenuity Systems, 2012 release), a software leveraging a manually reviewed repository of biological interactions and functional annotations was used to analyze the signalling pathways, cellular location, function and, network connections of the identified genes [28].

Gene Set Enrichment Analysis (GSEA) [29], was used to find whether a set of genes defined based on prior biological knowledge (e.g., those in a common signaling pathway) shows statistically significant correlations with IL6. Briefly, for each of the seven EOCs datasets, through the “use a gene as phenotype” option, GSEA ranks the genes according to their correlation with IL6. This ranked lists is then interrogated against gene sets contained within the C2 curated gene sets (c2.all.v3.0.symbols.gmt), a collection of 2516 gene sets that are part of the Molecular Signatures Database (MSigDB) v3.0 (12, 13). The primary GSEA result is the enrichment score (ES), which reflects the degree to which a gene set is overrepresented at either the top or bottom of the ranked list of genes. To estimate the statistical significance of the ES, a nominal p value is calculated by permuting the genes 1,000 times. The ES score is normalized to account for the gene set sizes (NES). Gene sets associated to a false positive rate (FDR) of less than 0.25 were considered significant.

Reagents

Recombinant human EGF was from Peprotech. Taqman® Gene Expression Assays were from Applied Biosystems (Foster City, CA, USA).

Ovarian cancer cell lines

SKOV3, IGROV1 (serous histotype) cell lines were obtained from ATCC and maintained in RPMI 1640 medium (Sigma Aldrich) with 10% fetal calf serum (FCS) (Hyclone, Logan, UT) and 2 mmol/L glutamine, in a 5% CO2 humidified atmosphere at 37°C. OAW42 (serous histotype, kindly provided by Dr. A. Ullrich, Max Planck Institute of Biochemistry, Martinsried, Germany) cells were cultured in MEM (Sigma Aldrich) and supplemented as above. IOSE-64 hTERT cells were maintained and prepared as described [34]. All cell lines used in this study were subjected to short tandem repeat (STR) analysis and the profiles were compared to publically available databases to verify their authenticity. For the in vitro validation, a time course (up to 24 hr) with EGF stimulation was performed and IL6 expression was monitored by real timer RT-PCR in order to assess the shorter time necessary to detect IL6 up-modulation. Based on this method IGROV1, OAW42 and IOSE 64 hTERT were EGF stimulated for 4 hr and and SKOV3 for 8 hr.

RNA Extraction and real time RT-PCR

Real time RT-PCR on selected IL6-correlated genes was performed on total RNA extracted from EOC cell lines stimulated for 4 hr (IGROV1 and OAW42 cells) and 8 hr (SKOV3 cells) with EGF (20 ng/ml). Total RNA from cell lines was extracted using a commercial kit (Amersham Bioscience-GE Healthcare). RT-PCR analysis was performed as described [17]. Human GAPD (GAPDH) Endogenous Control (VIC/MGB Probe) (RefSeq NM_002046.3) was used as housekeeping gene for normalization among samples. The Taqman Assays used for amplification were: Hs00174131_m1 for IL6; Hs00236966_m1 for CXCL2; Hs00181813_m1 for HBEGF; Hs01126604_m1 for SERPINE1; Hs00610256_g1 for DUSP1; Hs00185658_m1 for ZFP36; Hs00174674_m1 for IER3; Hs00171851_m1 for FOSB; Hs00374230_m1 for NR4A1; Hs00166165_m1 for EGR2; Hs00275699_ for MMP19; Hs00264920_s1 for THBD; Hs003604396_g1 for KLF2 (Applied Biosystems). Data analysis was performed by the Sequence Detection System (SDS) 2.2.2 software (Applied Biosystems).

Abbreviations

EOC:

Epithelial ovarian cancer

LMP:

Low malignant potential

OSE:

Ovarian surface epithelium

IPA:

Ingenuity pathway analysis

GSEA:

Gene set enrichment analysis.

References

  1. Bast RC, Hennessy B, Mills GB: The biology of ovarian cancer: new opportunities for translation. Nat Rev Cancer. 2009, 9: 415-428. 10.1038/nrc2644.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Berns EM, Bowtell DD: The changing view of high-grade serous ovarian cancer. Cancer Res. 2012, 72: 2701-2704. 10.1158/0008-5472.CAN-11-3911.

    Article  CAS  PubMed  Google Scholar 

  3. Fleming JS, Beaugie CR, Haviv I, Chenevix-Trench G, Tan OL: Incessant ovulation, inflammation and epithelial ovarian carcinogenesis: revisiting old hypotheses. Mol Cell Endocrinol. 2006, 247: 4-21. 10.1016/j.mce.2005.09.014.

    Article  CAS  PubMed  Google Scholar 

  4. Vaughan S, Coward JI, Bast RC, Berchuck A, Berek JS, Brenton JD, Coukos G, Crum CC, Drapkin R, Etemadmoghadam D, Friedlander M, Gabra H, Kaye SB, Lord CJ, Lengyel E, Levine DA, McNeish IA, Menon U, Mills GB, Nephew KP, Oza AM, Sood AK, Stronach EA, Walczak H, Bowtell DD, Balkwill FR: Rethinking ovarian cancer: recommendations for improving outcomes. Nat Rev Cancer. 2011, 11: 719-725.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Kulbe H, Thompson R, Wilson JL, Robinson S, Hagemann T, Fatah R, Gould D, Ayhan A, Balkwill F: The inflammatory cytokine tumor necrosis factor-alpha generates an autocrine tumor-promoting network in epithelial ovarian cancer cells. Cancer Res. 2007, 67: 585-592. 10.1158/0008-5472.CAN-06-2941.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Lo CW, Chen MW, Hsiao M, Wang S, Chen CA, Hsiao SM, Chang JS, Lai TC, Rose-John S, Kuo ML, Wei LH: IL-6 trans-signaling in formation and progression of malignant ascites in ovarian cancer. Cancer Res. 2011, 71: 424-434. 10.1158/0008-5472.CAN-10-1496.

    Article  CAS  PubMed  Google Scholar 

  7. Wang Y, Niu XL, Qu Y, Wu J, Zhu YQ, Sun WJ, Li LZ: Autocrine production of interleukin-6 confers cisplatin and paclitaxel resistance in ovarian cancer cells. Cancer Lett. 2010, 295: 110-123. 10.1016/j.canlet.2010.02.019.

    Article  CAS  PubMed  Google Scholar 

  8. Yadav A, Kumar B, Datta J, Teknos TN, Kumar P: IL-6 promotes head and neck tumor metastasis by inducing epithelial-mesenchymal transition via the JAK-STAT3-SNAIL signaling pathway. Mol Cancer Res. 2011, 9: 1658-1667. 10.1158/1541-7786.MCR-11-0271.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Nilsson MB, Langley RR, Fidler IJ: Interleukin-6, secreted by human ovarian carcinoma cells, is a potent proangiogenic cytokine. Cancer Res. 2005, 65: 10794-10800. 10.1158/0008-5472.CAN-05-0623.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Wang L, Yi T, Kortylewski M, Pardoll DM, Zeng D, Yu H: IL-17 can promote tumor growth through an IL-6-Stat3 signaling pathway. J Exp Med. 2009, 206: 1457-1464. 10.1084/jem.20090207.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Plante M, Rubin SC, Wong GY, Federici MG, Finstad CL, Gastl GA: Interleukin-6 level in serum and ascites as a prognostic factor in patients with epithelial ovarian cancer. Cancer. 1994, 73: 1882-1888. 10.1002/1097-0142(19940401)73:7<1882::AID-CNCR2820730718>3.0.CO;2-R.

    Article  CAS  PubMed  Google Scholar 

  12. Lane D, Matte I, Rancourt C, Piche A: Prognostic significance of IL-6 and IL-8 ascites levels in ovarian cancer patients. BMC Cancer. 2011, 11: 210-10.1186/1471-2407-11-210.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Guo Y, Nemeth J, O'Brien C, Susa M, Liu X, Zhang Z, Choy E, Mankin H, Hornicek F, Duan Z: Effects of siltuximab on the IL-6-induced signaling pathway in ovarian cancer. Clin Cancer Res. 2010, 16: 5759-5769. 10.1158/1078-0432.CCR-10-1095.

    Article  CAS  PubMed  Google Scholar 

  14. Dijkgraaf EM, Welters MJ, Nortier JW, Van der Burg SH, Kroep JR: Interleukin-6/interleukin-6 receptor pathway as a new therapy target in epithelial ovarian cancer. Curr Pharm Des. 2012, 18: 3816-3827. 10.2174/138161212802002797.

    Article  CAS  PubMed  Google Scholar 

  15. Alberti C, Pinciroli P, Valeri B, Ferri R, Ditto A, Umezawa K, Sensi M, Canevari S, Tomassetti A: Ligand-dependent EGFR activation induces the co-expression of IL-6 and PAI-1 via the NFkB pathway in advanced-stage epithelial ovarian cancer. Oncogene. 2012, 31: 4139-4149. 10.1038/onc.2011.572.

    Article  CAS  PubMed  Google Scholar 

  16. Castellano G, Reid JF, Alberti P, Carcangiu ML, Tomassetti A, Canevari S: New potential ligand-receptor signaling loops in ovarian cancer identified in multiple gene expression studies. Cancer Res. 2006, 66: 10709-10719. 10.1158/0008-5472.CAN-06-1327.

    Article  CAS  PubMed  Google Scholar 

  17. Degl'Innocenti D, Alberti C, Castellano G, Greco A, Miranda C, Pierotti MA, Seregni E, Borrello MG, Canevari S, Tomassetti A: Integrated ligand-receptor bioinformatic and in vitro functional analysis identifies active TGFA/EGFR signaling loop in papillary thyroid carcinomas. PLoS One. 2010, 5: e12701-10.1371/journal.pone.0012701.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Sensi M, Catani M, Castellano G, Nicolini G, Alciato F, Tragni G, De SG, Bersani I, Avanzi G, Tomassetti A, Canevari S, Anichini A: Human cutaneous melanomas lacking MITF and melanocyte differentiation antigens express a functional Axl receptor kinase. J Invest Dermatol. 2011, 131: 2448-2457. 10.1038/jid.2011.218.

    Article  CAS  PubMed  Google Scholar 

  19. Tothill RW, Tinker AV, George J, Brown R, Fox SB, Lade S, Johnson DS, Trivett MK, Etemadmoghadam D, Locandro B, Traficante N, Fereday S, Hung JA, Chiew YE, Haviv I, Gertig D, DeFazio A, Bowtell DD: Novel molecular subtypes of serous and endometrioid ovarian cancer linked to clinical outcome. Clin Cancer Res. 2008, 14: 5198-5208. 10.1158/1078-0432.CCR-08-0196.

    Article  CAS  PubMed  Google Scholar 

  20. Anglesio MS, Arnold JM, George J, Tinker AV, Tothill R, Waddell N, Simms L, Locandro B, Fereday S, Traficante N, Russell P, Sharma R, Birrer MJ, DeFazio A, Chenevix-Trench G, Bowtell DD: Mutation of ERBB2 provides a novel alternative mechanism for the ubiquitous activation of RAS-MAPK in ovarian serous low malignant potential tumors. Mol Cancer Res. 2008, 6: 1678-1690. 10.1158/1541-7786.MCR-08-0193.

    Article  CAS  PubMed  Google Scholar 

  21. Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, Joshi MB, Harpole D, Lancaster JM, Berchuck A, Olson JA, Marks JR, Dressman HK, West M, Nevins JR: Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006, 19;439: 353-357.

    Article  Google Scholar 

  22. Cancer Genome Atlas Research Network: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474: 609-615. 10.1038/nature10166.

    Article  Google Scholar 

  23. Berchuck A, Iversen ES, Luo J, Clarke JP, Horne H, Levine DA, Boyd J, Alonso MA, Secord AA, Bernardini MQ, Barnett JC, Boren T, Murphy SK, Dressman HK, Marks JR, Lancaster JM: Microarray analysis of early stage serous ovarian cancers shows profiles predictive of favorable outcome. Clin Cancer Res. 2009, 15: 2448-2455. 10.1158/1078-0432.CCR-08-2430.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Dressman HK, Berchuck A, Chan G, Zhai J, Bild A, Sayer R, Cragun J, Clarke J, Whitaker RS, Li L, Gray J, Marks J, Ginsburg GS, Potti A, West M, Nevins JR, Lancaster JM: An integrated genomic-based approach to individualized treatment of patients with advanced-stage ovarian cancer. J Clin Oncol. 2007, 25: 517-525. 10.1200/JCO.2006.06.3743.

    Article  CAS  PubMed  Google Scholar 

  25. Yoshihara K, Tajima A, Yahata T, Kodama S, Fujiwara H, Suzuki M, Onishi Y, Hatae M, Sueyoshi K, Fujiwara H, Kudo Y, Kotera K, Masuzaki H, Tashiro H, Katabuchi H, Inoue I, Tanaka K: Gene expression profile for predicting survival in advanced-stage serous ovarian cancer across two independent datasets. PLoS One. 2010, 5: e9615-10.1371/journal.pone.0009615.

    Article  PubMed Central  PubMed  Google Scholar 

  26. The Genecards database. http://www.genecards.org. 2012.Ref Type: Generic

  27. Shahzad MM, Arevalo JM, Rmaiz-Pena GN, Lu C, Stone RL, Moreno-Smith M, Nishimura M, Lee JW, Jennings NB, Bottsford-Miller J, Vivas-Mejia P, Lutgendorf SK, Lopez-Berestein G, Bar-Eli M, Cole SW, Sood AK: Stress effects on FosB- and interleukin-8 (IL8)-driven ovarian cancer growth and metastasis. J Biol Chem. 2010, 285: 35462-35470. 10.1074/jbc.M110.109579.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. The Ingenuity Pathway Analysis Software. 2012,http://www.ingenuity.com,

  29. The Gene Set Enrichment Analysis software. 2012,http://www.broad.mit.edu/gsea,

  30. Amit I, Citri A, Shay T, Lu Y, Katz M, Zhang F, Tarcic G, Siwak D, Lahad J, Jacob-Hirsch J, Amariglio N, Vaisman N, Segal E, Rechavi G, Alon U, Mills GB, Domany E, Yarden Y: A module of negative feedback regulators defines growth factor signaling. Nat Genet. 2007, 39: 503-512. 10.1038/ng1987.

    Article  CAS  PubMed  Google Scholar 

  31. Nagashima T, Shimodaira H, Ide K, Nakakuki T, Tani Y, Takahashi K, Yumoto N, Hatakeyama M: Quantitative transcriptional control of ErbB receptor signaling undergoes graded to biphasic response for cell differentiation. J Biol Chem. 2007, 282: 4045-4056.

    Article  CAS  PubMed  Google Scholar 

  32. Kim HS, Kim MS, Hancock AL, Harper JC, Park JY, Poy G, Perantoni AO, Cam M, Malik K, Lee SB: Identification of novel Wilms' tumor suppressor gene target genes implicated in kidney development. J Biol Chem. 2007, 282: 16278-16287. 10.1074/jbc.M700215200.

    Article  CAS  PubMed  Google Scholar 

  33. Delassus GS, Cho H, Hoang S, Eliceiri GL: Many new down- and up-regulatory signaling pathways, from known cancer progression suppressors to matrix metalloproteinases, differ widely in cells of various cancers. J Cell Physiol. 2010, 224: 549-558. 10.1002/jcp.22157.

    Article  CAS  PubMed  Google Scholar 

  34. De Cecco L, Marchionni L, Gariboldi M, Reid JF, Lagonigro MS, Caramuta S, Ferrario C, Bussani E, Mezzanzanica D, Turatti F, Delia D, Daidone MG, Oggioni M, Bertuletti N, Ditto A, Raspagliesi F, Pilotti S, Pierotti MA, Canevari S, Schneider C: Gene expression profiling of advanced ovarian cancer: characterization of molecular signature involving the fibroblast growth factor 2. Oncogene. 2004, 23: 8171-8183. 10.1038/sj.onc.1207979.

    Article  CAS  PubMed  Google Scholar 

  35. NF-kB Target Genes.http://www.bu.edu/nf-kb/gene-resources/target-genes,

  36. Qi L, Higgins SP, Lu Q, Samarakoon R, Wilkins-Port CE, Ye Q, Higgins CE, Staiano-Coico L, Higgins PJ: SERPINE1 (PAI-1) is a prominent member of the early G0 --> G1 transition "wound repair" transcriptome in p53 mutant human keratinocytes. J Invest Dermatol. 2008, 128: 749-753. 10.1038/sj.jid.5701068.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Pratilas CA, Taylor BS, Ye Q, Viale A, Sander C, Solit DB, Rosen N: (V600E)BRAF is associated with disabled feedback inhibition of RAF-MEK signaling and elevated transcriptional output of the pathway. Proc Natl Acad Sci USA. 2009, 106: 4519-4524. 10.1073/pnas.0900780106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Chandarlapaty S: Negative feedback and adaptive resistance to the targeted therapy of cancer. Cancer Discov. 2012, 2: 311-319. 10.1158/2159-8290.CD-12-0018.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Masoumi MS, Amini A, Morris DL, Pourgholami MH: Significance of vascular endothelial growth factor in growth and peritoneal dissemination of ovarian cancer. Cancer Metastasis Rev. 2012, 31: 143-162. 10.1007/s10555-011-9337-5.

    Article  Google Scholar 

  40. Kurman RJ, Shih I: The origin and pathogenesis of epithelial ovarian cancer: a proposed unifying theory. Am J Surg Pathol. 2010, 34: 433-443. 10.1097/PAS.0b013e3181cf3d79.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Lo-Ciganic WH, Zgibor JC, Bunker CH, Moysich KB, Edwards RP, Ness RB: Aspirin, nonaspirin nonsteroidal anti-inflammatory drugs, or acetaminophen and risk of ovarian cancer. Epidemiology. 2012, 23: 311-319. 10.1097/EDE.0b013e3182456ad3.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Kenny HA, Lengyel E: MMP-2 functions as an early response protein in ovarian cancer metastasis. Cell Cycle. 2009, 8: 683-688. 10.4161/cc.8.5.7703.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. The Cancer Genome Atlas database. [https://tcga-data.nci.nih.gov/] 2012.

Download references

Acknowledgements

We thank Dr. Patrick Moore for English editing of the manuscript. Financial Support: Italian Association for Cancer Research (IG10302 to SC and IG13055 to AT) and to Italian Ministry of Health ("Progetto Oncologico di Medicina Molecolare: i Tumori Femminili" to SC and “Progetto Integrato in Oncologia” to AT).

Author information

Authors and Affiliations

  1. Unit of Molecular Therapies, Department of Experimental Oncology and Molecular Medicine, Fondazione IRCCS Istituto Nazionale dei Tumori, Milan, Italy

    Patrizia Pinciroli, Chiara Alberti, Silvana Canevari & Antonella Tomassetti

  2. Unit of Human Tumor Immunobiology, Department of Experimental Oncology and Molecular Medicine, Fondazione IRCCS Istituto Nazionale dei Tumori, Milan, Italy

    Marialuisa Sensi

Authors
  1. Patrizia Pinciroli
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chiara Alberti
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Marialuisa Sensi
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Silvana Canevari
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Antonella Tomassetti
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Silvana Canevari or Antonella Tomassetti.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PP carried out bioinformatics and statistical analysis. CA carried out the in vitro biological validation. MS contributed to the design of the study. SC and AT, conceived the study, participated in its design and coordination, and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2012_5214_MOESM1_ESM.xlsx

Additional file 1:Table containing: a. Selected probe sets for each platform; b. IL6-correlated genes in serous high malignant EOCs.(XLSX 45 KB)

12864_2012_5214_MOESM2_ESM.pptx

Additional file 2:Figure reporting IL6 distribution (density plot) in the three data sets containing expression data of both advanced stage (204, 60 and 40 patients in data set I, II and IV, respectively) and LMP (18, 30 and 19 patients in data set I, II and IV, respectively) EOCs.(PPTX 640 KB)

12864_2012_5214_MOESM3_ESM.xlsx

Additional file 3:Table reporting IL6-correlated genes in serous advanced stage and LMP EOCs from data set I.(XLSX 13 KB)

Additional file 4:Table reporting the networks identified by IPA software.(XLSX 11 KB)

Additional file 5:Table reporting IL6-correlated genes included in each the gene sets selected by GSEA.(XLSX 12 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Pinciroli, P., Alberti, C., Sensi, M. et al. An IL6-correlated signature in serous epithelial ovarian cancer associates with growth factor response. BMC Genomics 14, 508 (2013). https://doi.org/10.1186/1471-2164-14-508

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-14-508

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us