Variation in gene expression patterns in effusions and primary tumors from serous ovarian cancer patients
© Schaner et al. 2005
Received: 22 April 2005
Accepted: 21 July 2005
Published: 21 July 2005
While numerous studies have characterized primary ovarian tumors, little information is available regarding expression patterns of metastatic sites of this cancer. To define sets of genes that distinguish primary and metastatic ovarian tumors, we used cDNA microarrays to characterize global gene expression patterns in 38 effusions (28 peritoneal, 10 pleural) and 8 corresponding primary ovarian tumors, and searched for associations between expression patterns and clinical parameters.
We observed multidimensional variation in expression patterns among the cancers. Coordinate variation in expression of genes from two chromosomal regions, 8q and 19q, was seen in subsets of the cancers indicating possible amplifications in these regions. A set of 112 unique genes of known function was differentially expressed between primary tumors and effusions using supervised analysis. Relatively few differences were seen between effusions isolated from the pleural and peritoneal cavities or between effusions from patients diagnosed with stage III and stage IV cancers. A set of 84 unique genes was identified that distinguished high from lower grade ovarian cancers. The results were corroborated using immunocytochemistry, mRNA in situ hybridization, and immunoblotting.
The extensive variation in expression patterns observed underscores the molecular heterogeneity of ovarian cancer, but suggests a similar molecular profile for ovarian carcinoma cells in serosal cavities.
Epithelial ovarian carcinoma claims more lives than any other gynecologic malignancy, largely because it frequently escapes detection after it has metastasized . Ovarian carcinoma initially metastasizes primarily to the serosal surface of the peritoneal cavity and abdominal organs. The pleural space is often involved as well, either at diagnosis or, more commonly, at later stages of clinical progression. Pleural effusion is the most common presentation of stage IV disease . A number of metastasis-associated molecules have been reported to be differentially expressed between primary ovarian tumors and tumor cells in effusions [3–12], but little is known regarding the mechanism of metastases.
Molecular characterization of ovarian carcinoma using DNA microarrays has so far focused on primary tumors [13–22]. The paucity of data regarding the biological characteristics of ovarian carcinoma cells in effusions at both the phenotypic and genotypic level limits our understanding of tumor progression in this disease. Specifically, we do not know how ovarian carcinoma cells in ascites and pleural effusions differ from those in the corresponding solid primary tumors, or whether and how carcinoma cells in peritoneal and pleural effusions differ. Moreover, molecular analysis of malignant effusions might contribute to better predictions of survival and treatment response.
To identify genes whose expression may be associated with this metastatic behavior, we analyzed global gene expression patterns of ovarian cancer cells obtained from 3 distinctive anatomic sites: 28 peritoneal, 10 pleural and 8 primary tumors (see supplementary Table S1.xls). A valuable feature of this dataset is that it includes 8 paired samples of primary tumors and malignant effusions from the same patients. We were able to define a number of genes that differentiate primary tumors from effusions.
Overview of global gene expression patterns among ovarian cancers
Genes involved in cell cycle progression and cell proliferation were variably expressed among the cancers presumably reflecting variation in proliferation rates among the tumors in a coordinated manner (Figure 1C, panel c) [23, 24]. Co-expression of previously identified markers for ovarian cancer including kallikreins 6 and 10 as well as the S100 calcium-binding proteins S100A1 and S100A13 (Figure 1C, panel d) was also seen. Co-clustering of CEACAM5 and CEACAM7 with syndecan-1 (Figure 1C, panel e) was observed. These genes were part of a larger cluster with expression of Keratins 5, 6B, 7, 15 and 19, and S100A10, caveolin-2 and SLPI (Figure 1C, panel e). In another cluster (Figure 1C, panel f) the sprouty homologs: 1, 2, and 4, were co-expressed with the dual specificity phosphatase 6 (DUSP 6) and two ets variant genes 4 and 5 (ETV4 and ETV5), indicating possible involvement of the MAP kinase pathway in some tumors. Also, co-expression of caveolin 1 and EGFR was seen, in agreement with data from experimental models [25, 26] (Figure 1C, panel a). The complete gene list and clustered file is found in supplementary Figure S1 CDT.cdt.
Expression differences between primary tumors, peritoneal effusions and pleural effusions
Selected genes identified by SAM that are differentially expressed in primary tumors vs. all effusions. The complete list is available in Table S2.
Genes more highly expressed in Effusions
CLDN7 : claudin 7 : Hs.278562
KRT7 : keratin 7 : Hs.23881
CRIP1 : cysteine-rich protein 1 (intestinal) : Hs.423190
L1CAM : L1 cell adhesion molecule MASA: Hs.1757
ADM : adrenomedullin : Hs.394
CRYAB : crystallin, alpha B : Hs.391270
IL18RAP : interleukin 18 receptor accessory protein : Hs.158315
KRT19 : keratin 19 : Hs.182265
Genes more highly expressed in Primary Tumors
COL6A3 : collagen, type VI, alpha 3 : Hs.80988
COL1A2 : collagen, type I, alpha 2 : Hs.179573
MYB : v-myb myeloblastosis viral oncogene homolog (avian) : Hs.1334
BGN : biglycan : Hs.821
IGFBP7 : insulin-like growth factor binding protein 7 : Hs.119206
SPARC : secreted protein, acidic, cysteine-rich (osteonectin) : Hs.111779
AEBP1 : AE binding protein 1 : Hs.118397
APOD : apolipoprotein D : Hs.75736
CDH11 : cadherin 11, type 2, OB-cadherin (osteoblast) : Hs.75929
Genes identified by SAM that are differentially expressed in Pleural vs. Peritoneal Effusions
Positive Significant Genes: Higher in pleural vs. ascites
NR4A1: nuclear receptor subfamily 4, group A, member 1: Hs.1119
CYR61: cysteine-rich, angiogenic inducer, 61: Hs.8867
CXCL2: chemokine (C-X-C motif) ligand 2: Hs.75765
RAB21: member RAS oncogene family: Hs.184627
CTGF: connective tissue growth factor: Hs.75511
CXCL3: chemokine (C-X-C motif) ligand 3: Hs.89690
TCEB3: transcription elongation factor B (SIII): Hs.155202
IGLL1: immunoglobulin lambda-like polypeptide 1: Hs.348935
CTGF: connective tissue growth factor: Hs.75511
IGHG3: immunoglobulin heavy constant gamma 3: Hs.413826
C1QB: complement component 1, q subcomponent beta: Hs.8986
CYR61: cysteine-rich, angiogenic inducer, 61: Hs.8867
C1QG: complement component 1, q subcomponent gamma:Hs.94953
TAGLN: transgelin: Hs.433399
GSTA4: glutathione S-transferase A4: Hs.169907
Negative Significant Genes: Lower in pleural vs. ascites
N33: Putative prostate cancer tumor suppressor: Hs.71119 100791
PLEC1: plectin 1, intermediate filament binding protein 500 kDa: Hs.79706
SLC40A1: solute carrier family 40, member 1: Hs.5944
We then applied a supervised statistical method, PAM  to see whether it was possible to find a set of genes that could classify the primary tumors from peritoneal and pleural effusions. By comparing the three groups, a set of genes (615; 436 unique genes) was identified that correctly classified the primary tumors and most of the peritoneal (error rate 10%) and pleural effusions (error rate 20%). Cross-validation was not as successful, suggesting again that the difference between the 2 groups of effusions is not as clear as the distinction between all effusions and the primary tumors, which is in accordance to the results of the SAM analyses. Some differences were however detected. Clustering of the 615 genes (Figure 4 and supplementary Tables S3a, S3b.xls) illustrates that the expression of a number of chemokine ligands, including CCL2, CCL8 and CCL3L1 was more frequent in pleural effusions possibly reflecting the larger number of leukocytes in these specimens (Figure 4C, panel b). Higher expression of a number of genes on Chromosome 8q24 was observed among most of the ascites (Figure 4C, panel a highlighted in red). In addition, the expression of genes whose proteins have previously been shown to be produced by both ovarian carcinoma and stromal cells [28–30], including TIMP-2, vimentin and basic fibroblast growth factor separated the primary tumors from the effusions, suggesting that the cancer-stroma crosstalk is associated with different biological pathway activation than that observed in effusions (Figure 4C, panel c and supplementary Table S2.xls).
Hierarchical clustering of all 38 effusions was carried out as described previously. Results were very similar to those obtained when clustering both the primary tumors and effusions (Web Supplement, Figure S2.pdf) with the same 11 samples residing on the left branch as in Fig. 1 and with a strong overexpression of the immune response cluster of genes (Figure 1h). Interestingly, the majority of the pleural effusions (6/10) clustered to this branch (p = 0.02).
Identification of Expression Profiles Based on 'Intrinsic' Genelist
representative list is shown below and the complete list is available in Table S5.
DERP6: S-phase 2 protein: Hs.417029
VARS2L: valyl-tRNA synthetase 2-like: Hs.102910
SC5DL: sterol-C5-desaturase (ERG3 delta-5-desaturase homolog, fungal)-like: Hs.287749
MYLIP: myosin regulatory light chain interacting protein: Hs.443793
HMGN3: high mobility group nucleosomal binding domain 3: Hs.77558
MAGEA6: melanoma antigen, family A, 6: Hs.441113
ATP5L2: ATP synthase, H+ transporting, mitochondrial F0 complex, subunit g, isoform 2
GPR48: G protein-coupled receptor 48: Hs.160271
MASK: multiple ankyrin repeats, single KH-domain (MASK) homolog: Hs.528646
PRO1853: hypothetical protein PRO1853: Hs.433466
YKT6: SNARE protein Ykt6: Hs.296244
ZNF-kaiso: kaiso: Hs.143604
NCK1: NCK adaptor protein 1: Hs.54589
We examined global gene expression patterns in primary ovarian cancers and malignant effusions and found differences that may be related to tumor progression. Our results are consistent with many reports that document the heterogeneity of gene expression patterns noted in ovarian cancer [14, 16–18, 20–22]. These data underscore the heterogeneity of this disease and the profound molecular differences within tumor sub-groups with comparable morphology.
Several clusters of variably expressed genes may have relevance to the biology of ovarian cancer. A large cluster of genes from chromosome 8q21-24 was more highly expressed among a subset of the cancers, suggesting an amplification of this region on chromosome 8. This was seen in both the primary and the effusions from the corresponding samples with these abnormalities. Transcripts over-expressed in this cluster include YWHAZ, that encodes the zeta isoform of 14-3-3 protein (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein) family. 14-3-3 proteins are expressed in a number of cancers and are involved in the cell cycle and also in prolonging cell survival. TPD52 (hD52) that has been shown to be in the peak of the 8q21 amplicon in breast cancer cell lines was also over-expressed in a number of cancers in this cluster. A recent report suggests that TPD52 is a candidate target gene and putative oncogene at 8q21 . In addition, co-clustering of a number of genes in this region was also observed in the cluster resulting from the intrinsic genelist (Figure 5C, panel a).
Involvement of the MAP kinase pathway may be indicated in a large number of the cancers in our study based on the co-expression of the 3 sprouty transcripts (1, 2, and 4) (Figure 1C panel f). Expression of all three sprouty homologs has been observed before in mouse development , but not in human cancers. Human sprouty 2 has been shown to inhibit the mitogen-activated protein (MAP) kinase pathway , and sprouty proteins have been implicated in the negative regulation of the receptor tyrosine kinase-induced MAP kinase pathway. DUSP, another gene in our sprouty cluster, belongs to a family of 9 cytoplasmic and/or nuclear enzymes, that function by dephosphorylating threonine and tyrosine residues of p38, JNK and ERK [36, 37]. The transcriptional regulation of Ets transcription factors, including that of PEA3 (ETV4), on the other hand, involves signals initiated by growth factor signaling through tyrosine kinase receptors that may be mediated by MAPK . We have previously shown that PEA3 and DUSP member PAC-1 predict poor outcome in ovarian carcinoma, whereas expression of all three MAPK predicts improved survival [39–41]. The coordinated expression of 3 sprouty family members, DUSP6, and ETV4 and 5 in a subset of ovarian cancers therefore raises the possibility that the activity of a specific, MAPK-related signaling pathway may have a role in ovarian cancer.
We have previously shown that clinical and molecular markers that are of established prognostic role in primary ovarian cancer have little or no significance in effusions . Besides the obvious fact that cells at this site represent tumor progression, they also appear to be biologically different than tumor cells in both primary tumors and solid metastases. Using PAM to analyze the three different groups, primary tumors, and pleural and peritoneal effusions, we identified a subset of genes that is more highly expressed in the primary tumors, but fewer differences between pleural and peritoneal effusions. These genes are mainly characteristic of previously defined 'stromal' signature , but that are produced by both tumor and stromal cells including collagens, TIMP2, bFGF, vimentin and SPARC [28–30]. To study tumor progression we evaluated the expression levels in relation to clinical parameters using SAM, focusing on FIGO stage and grade. The most notable distinctions were based on grade rather than stage.
A large cluster of genes on chromosome segment 19q13.1 showed covarying expression among the cancers in this study, consistent with previous reports that the 19q13.1 region is amplified in some ovarian carcinomas [44, 45]. The AKT2 oncogene contained in this region of chromosome 19, has been shown to be associated with the progression of ovarian cancer [44, 45]. Interestingly, we found that some of the genes in the 'chromosome 19 cluster' including PDCD5 and PEN2, were more highly expressed in the Stage IV cancers than earlier-stage supporting the notion that amplification of the 19q13.1 region and concomitant elevated expression of these genes may play a role in the progression of the disease . Furthermore in a larger cohort of advanced ovarian carcinomas amplification of this region was found to be associated with poor survival (Wang et. al., manuscript to be submitted).
In summary, we have examined the relationship between primary tumors in ovarian cancer and their corresponding effusions. The most notable difference was observed in expression patterns between effusions and the primary tumors. There is significant molecular variation among the cancers. There are some hints of differences related to tumor grade and effusion versus primary tumor that will need further investigation to see whether these are significant. Finally, some consistent features of expression patterns in subsets of the cancers may suggest possible molecular alterations involved in the biology of the tumors.
Materials and methods
Material consisted of 38 fresh non-fixed peritoneal and pleural effusions submitted from the Department of Gynecological Oncology to the Division of Cytology, Department of Pathology, The Norwegian Radium Hospital, during the period of April 1998-May 2002. Specimens were obtained pre-operatively, intra-operatively, or at disease recurrence, from 36 patients diagnosed with ovarian carcinoma (35 of the serous type, one clear cell type) and one patient diagnosed with primary peritoneal carcinoma (PPC). From one patient, two ascites specimens were obtained one month apart. Effusion specimens consisted of 28 peritoneal and 10 pleural effusions. Patient age ranged from 35 to 73 years (mean = 58 years). Twenty patients were diagnosed with FIGO stage III disease and 16 with stage IV disease. The remaining patient had a stage I tumor. Tumor grade for 36 serous carcinomas was as follows: 5 grade 1, 10 grade 2 and 21 grade 3 carcinomas. All relevant clinical data were obtained from the Department of Gynecologic Oncology, and are presented in detail in web supplementary Table S1.xls. In order to preserve physiological activity, specimens submitted to our laboratory arrived within minutes after collection and were processed immediately. Cells were suspended and frozen in RPMI+DMSO at -70°C. Smears and cellblock sections from all specimens underwent morphological evaluation by three experienced cytopathologists, and were further characterized using immunocytochemistry with broad antibody panels against epithelial and mesothelial epitopes, as previously detailed [5, 6, 47]. In all specimens included in this study, cancer cells comprised 50% or more of the entire cell population based on cytology smears. All patients in this study were treated with platinum based chemotherapy according to current guidelines, and the samples were collected under an IRB approved protocol (S-01127; June 22, 2001).
We obtained samples of primary tumors from 8 of the patients whose effusions were analyzed. These were snap-frozen in liquid nitrogen upon removal and stored at -70°C. Frozen sections were obtained from all biopsies in order to evaluate the percentage of tumor cells and tissue viability. The former exceeded 50% of cells in all cases.
RNA isolation, Labeling and Hybridization
Total RNA was isolated from effusions and solid tumors using the TRIZOL Reagent (Gibco BRL, Life Technologies). mRNA isolation from total RNA was undertaken using d(T) coated Dynabeads (Dynal, Oslo, Norway). One to two μg of experimental sample mRNA was used for labeling with Cy5-dUTP. mRNA was reverse transcribed with Superscript II (Invitrogen Life Technologies) using an oligo dT primer (Operon Technologies, Alameda, CA). Each sample was comparatively hybridized to cDNA microarrays together with a common reference labeled with Cy3-dUTP (Stratagene). Fluorescent dyes were purchased from Amersham Pharmacia Biotech, Piscataway, NJ. Hybridizations were carried out using the standard protocol described previously. [23, 24]. Complete experimental details may be found at: http://brownlab.stanford.edu/protocols.html.
DNA Microarrays: All arrays were printed at the Stanford Functional Genomics Facility. DNA clones on the microarrays are based primarily the sequence verified IMAGE clones from the Research Genetics Corporation (Huntsville, AL) http://www.invitrogen.com/content/sfs/manuals/sequenceverifiedclones_man.pdf and the CGAP clone set http://cgap.nci.nih.gov/Genes/PurchaseReagents as well as a small percentage of custom spots and control spots. Complete details regarding the clones on the arrays may be found at http://www.microarray.org/sfgf/jsp/servicesFrame.jsp#productionArrays. These microarrays were comprised of 41,805 elements (42 K), representing an estimated 25,695 genes as judged by the number of unique Unigene symbols.
Data Analysis and Clustering
Data Selection: Data were analyzed by using either the GenePix 3.0 or GenePix 4.0 software (Axon Instruments). Spots with aberrant measurements due to obvious array artifacts or poor technical quality were manually flagged and removed from further analysis. A filter was applied to omit measurements where fluorescent signal from the DNA spot was less than 50% above the measured background fluorescence surrounding the printed DNA spot in both the Cy3 and Cy5 Channels. Genes that did not meet these criteria for at least 80% of the measurements across the cases were excluded from further analysis. Data were retrieved as log2(Cy5/Cy3). The (Cy5/Cy3) ratio is defined in Stanford Microarray Database (SMD) as the normalized ratio of the background-corrected intensities . Genes whose expression level differed by more than 2.5-fold from their mean expression level in the sample set in at least 3 samples, were selected for further analysis.
Significance Analysis of Microarrays (SAM)
SAM is a statistical approach to identify genes whose expression patterns are significantly associated with specific characteristics of sample sets  SAM analysis was applied to the ovarian dataset to examine differences between primary, pleural and peritoneal effusions, and to examine different clinical parameters, including stage and grade. A two-way, unpaired test was carried out comparing the two groups of interest. A 10-nearest neighbor imputation engine was applied to estimate missing data , and 500 permutations were carried out to compute expected values and to calibrate false positive calls.
Prediction Analysis for Microarrays (PAM)
Prediction Analysis for Microarrays (PAM) http://www-stat.stanford.edu/~tibs/PAM/ was carried out using the Excel version of the program. This method is applied to gene expression data to provide sample classification by 'shrunken centroids'. Data were filtered as described earlier for clustering analyses. The three groups used for analysis were defined as primary tumors, peritoneal effusions (ascites) and pleural effusions .
An 'intrinsic' genelist, comprising genes whose expression varied significantly more between samples from different patients than between replicate samples for the same patient, was selected based on the 8 primary tumor/effusion pairs and the paired effusion cancers using methods as described previously . A score was constructed which was the average within-pair variation relative to the between-pair variation. This was the ratio of the variance of the differences, relative to the variance of the averages.
Immunocytochemical analysis (ICC)
ICC was performed using antibodies directed against Syndecan-1 (CD138; SDC1), the integrin β3 subunit (platelet glycoprotein IIIa, antigen CD61), and CD44 antigen, all part of identified gene clusters showing variable degrees of expression among the samples. Negative controls consisted of sections that underwent a similar staining procedure, with the exclusion of primary antibody application, or that were stained with mouse myeloma protein of the same isotype as the primary antibody used. Biopsies in which immunoreactivity for the studied antigens had previously been demonstrated were used as positive controls.
Frozen effusion specimens were thawed and washed twice in phosphate buffered saline (PBS). Samples were subsequently lysed in 1% NP-40, 20 mM Tris HCl (pH 7.5), 137 mM NaCl, 10% glycerol, 1 mM Phenylmethylsulfonyl-fluoride, and 1 mM Sodium Orthovanadate, with 0.02 mg/ml aprotinin, leupeptin and pepstatin and 10 μg/ml phosphatase inhibitor cocktail I. All inhibitors were from Sigma-Aldrich (Saint Louis, MO). After centrifugation, the supernatant was collected and protein content was evaluated by the Bradford assay. Twenty-five μg total protein lysate/lane was resolved by SDS polyacrylamide gel electrophoresis (7.5 or 12 % gels) and transferred on to PVDF immobilon membranes (Millipore, Bedford, MA). Successful transfer was evaluated by staining of membranes withNaphtol-blue-black (Sigma-Aldrich). Filters were blocked in TBST containing 5 % dried skimmed milk and 0.05% Tween-20 for 1 hr at room temperature. Thereafter, the filters were incubated over night at 4°C with primary antibodies diluted in TBST containing 5 % skimmed milk and 0.05% Tween-20. Primary antibodies directed against the laminin γ-2 chain (Santa Cruz Biotechnology, Santa Cruz, CA) and Muc 5B (Santa Cruz) were chosen as examples since they were part of novel clusters displaying variation in expression among the cancers. Filters were washed 3 times 10 minutes each with TBST (0.05% Tween-20). They were subsequently incubated with HRP-conjugated secondary antibody diluted 1:5000 in TBST containing 5% dried skimmed milk and 0.05% Tween-20 for 45 minutes at room temperature. Immunoreactivity was detected using the ECL-plus western blotting system (Amersham-Pharmacia)
mRNA In Situ hybridization (ISH)
Three genes were chosen for ISH as representative examples for the evaluation of mRNA expression levels in the tumor cells themselves: MMP-14, MMP-9, and PEA3 (ETV4). The following antisense oligonucleotide probes were obtained from Research Genetics (Huntsville, AL) [1–28, 51, 52]:
MMP-9: 5' CCGGTCCACCTCGCTGGCGCTCCGGU 3', PEA3: 5' TGA ATT ATG AGA AGC TGA GCC G 3', MMP-14: 5' TCC ATC ACT TGG TTA TTC CTC ACC CGC C 3'.
A poly d(T)20 oligonucleotide (Research Genetics) was used to verify the integrity and lack of degradation of mRNA in each sample. The DNA probes were hyperbiotinylated. Stock dilution was prepared with a resulting equal concentration for both probes. The stock dilution was diluted with probe diluent (Research Genetics) immediately before use. Specific sense oligonucleotides were used for the evaluation of non-specific activity for each probe.
Cellblock sections (4 micron-thick) of formalin-fixed, paraffin-embedded specimens were mounted on ProbeOn Plus slides (Fisher Scientific, Pittsburgh, PA). Sectioning was performed in RNase-free water. Hybridization using the probes was carried out as previously described and by using the microprobe manual staining system (Fisher Scientific). Known positive controls were used in each hybridization reaction. These consisted of 2 ovarian carcinomas for which positive hybridization was reproducible in a previous study. Controls for endogenous alkaline phosphatase for all probes included treatment of the sample in the absence of the probe and use of chromogen alone.
This work was supported by grants from the Norwegian Cancer Society (D 99061), The Research Council of Norway (155218/300), NIH Grant 2HFZ542 and the Marsha Rivkin Scholar Award (M.E.S.). We wish to thank Mike Fero, and the staff of the Stanford Functional Genomics Facility, members of the Brown lab for helpful discussions, and the Stanford Microarray Database, with special thanks to Jeremy Gollub and Gavin Sherlock. Trevor Hastie provided advice and constructed the intrinsic genelist. We also would like to thank Pat Brown for helpful discussions, comments and support.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.