DNA Methylation Profiles of Ovarian Epithelial Carcinoma Tumors and Cell Lines

Abstract

Background

Epithelial ovarian carcinoma is a significant cause of cancer mortality in women worldwide and in the United States. Epithelial ovarian cancer comprises several histological subtypes, each with distinct clinical and molecular characteristics. The natural history of this heterogeneous disease, including the cell types of origin, is poorly understood. This study applied recently developed methods for high-throughput DNA methylation profiling to characterize ovarian cancer cell lines and tumors, including representatives of three major histologies.

Methodology/Principal Findings

We obtained DNA methylation profiles of 1,505 CpG sites (808 genes) in 27 primary epithelial ovarian tumors and 15 ovarian cancer cell lines. We found that the DNA methylation profiles of ovarian cancer cell lines were markedly different from those of primary ovarian tumors. Aggregate DNA methylation levels of the assayed CpG sites tended to be higher in ovarian cancer cell lines relative to ovarian tumors. Within the primary tumors, those of the same histological type were more alike in their methylation profiles than those of different subtypes. Supervised analyses identified 90 CpG sites (68 genes) that exhibited ‘subtype-specific’ DNA methylation patterns (FDR<1%) among the tumors. In ovarian cancer cell lines, we estimated that for at least 27% of analyzed autosomal CpG sites, increases in methylation were accompanied by decreases in transcription of the associated gene.

Significance

The significant difference in DNA methylation profiles between ovarian cancer cell lines and tumors underscores the need to be cautious in using cell lines as tumor models for molecular studies of ovarian cancer and other cancers. Similarly, the distinct methylation profiles of the different histological types of ovarian tumors reinforces the need to treat the different histologies of ovarian cancer as different diseases, both clinically and in biomarker studies. These data provide a useful resource for future studies, including those of potential tumor progenitor cells, which may help illuminate the etiology and natural history of these cancers.

About the Authors:

Sahar Houshdaran

Contributed equally to this work with: Sahar Houshdaran, Sarah Hawley

Affiliation: Department of Biochemistry and Molecular Biology, University of Southern California, Los Angeles, California, United States of America

Sarah Hawley

Contributed equally to this work with: Sahar Houshdaran, Sarah Hawley

Affiliation: Canary Foundation, Palo Alto, California, United States of America

Chana Palmer

Affiliation: Canary Foundation, Palo Alto, California, United States of America

Mihaela Campan

Affiliations Department of Biochemistry and Molecular Biology, University of Southern California, Los Angeles, California, United States of America, Department of Surgery, University of Southern California, Los Angeles, California, United States of America

Mari N. Olsen

Affiliation: Department of Biochemistry, Stanford University, Stanford, California, United States of America

Aviva P. Ventura

Affiliation: Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America

Beatrice S. Knudsen

Affiliation: Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America

Charles W. Drescher

Affiliation: Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America

Nicole D. Urban

Affiliation: Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America

Patrick O. Brown

Affiliation: Department of Biochemistry, Stanford University, Stanford, California, United States of America

Peter W. Laird

* E-mail: [email protected]

Affiliations Department of Biochemistry and Molecular Biology, University of Southern California, Los Angeles, California, United States of America, Department of Surgery, University of Southern California, Los Angeles, California, United States of America, University of Southern California Epigenome Center, University of Southern California, Los Angeles, California, United States of America

Introduction

Ovarian cancer is the leading cause of death among all gynecological cancers in the United States [1], and is the sixth leading cause of all cancer deaths among women. There are four major histological types of ovarian cancer (serous, endometrioid, mucinous, and clear cell), each with distinct histopathological, clinical and molecular characteristics. The natural history of ovarian tumors, including their cell type of origin and steps of carcinogenesis, is poorly understood and has been the subject of much debate and discussion. Many studies of ovarian cancer are based on ovarian cancer derived cell lines, but it is unclear to what extent these cell lines accurately model the disease with respect to various molecular characteristics, including DNA methylation status.

DNA methylation is one of the epigenetic mechanisms that plays a role in many important biological processes including X-inactivation [2]–[4], silencing parasitic DNA elements [5], genomic imprinting [6], aging [7], male infertility [8], and cancer. Previous studies have shown CpG island DNA hypermethylation in various cancers, including ovarian tumors, as well as reduced levels of global DNA methylation associated with cancer [9]–[15]. Furthermore, the DNA methylation profile of a tumor cell is a reflection of its somatic lineage, environmental exposure, genetic predisposition, and cell-type specific chromatin structure. Recent technological advances have enabled quantitative assessment of the DNA methylation status of thousands of loci at once, providing an unprecedented opportunity to investigate the epigenetic signature of a cell. This rich source of information could be utilized to better understand the cell of origin of the subtypes and the steps of carcinogenesis, to identify appropriate model systems for future studies, and to discover candidate biomarkers for disease detection, classification, and monitoring.

Results

DNA Methylation Profiles of Ovarian Surface Epithelial Primary Tumors and Cell Lines

We assessed the DNA methylation profile of 1,505 CpG sites (associated with 808 genes) in 15 ovarian cell lines and 27 primary tumors (15 serous, 9 endometrioid, and 3 clear cell), using the Illumina GoldenGate Cancer Panel I [16] (see Supplemental File S1 for methylation array details). This panel has been used in several previous studies by others and us [8], [16]–[22]. Several investigators have validated the DNA methylation results of this panel by pyrosequencing [17], [19], Methylation Sensitive PCR (MSP) [16], [21], and bisulfite genomic sequencing [16]. The reproducibility of this assay has also been reported previously [16], [19]. The clinical and histopathological characteristics of the tumors and cell lines are detailed in Tables 1 and 2 respectively.

[Figure omitted. See PDF.]

Table 1. Histology and clinical characteristics of primary ovarian tumors.

https://doi.org/10.1371/journal.pone.0009359.t001

[Figure omitted. See PDF.]

Table 2. Characteristics of ovarian cancer cell lines.

https://doi.org/10.1371/journal.pone.0009359.t002

Of the 1,505 CpG sites (808 genes) surveyed, 1,184 sites (686 genes) exhibited sufficient variation across the 42 specimens for further analysis (see Methods). Unsupervised hierarchical clustering of both specimens and CpG sites using these 1,184 probes is shown in Figure 1. The primary division in the clustering of the specimens was between cell lines and tumors, each of which formed a distinct cluster. Among the tumors, the primary serous and endometrioid tumors each formed a distinct cluster (with one exception), whereas the three clear cell tumors were distributed between the two main clusters. There was no clear pattern of sub-clustering among the cell lines, despite the diversity of histological types nominally represented. The most notable cluster of CpG sites was a large group of CpG sites on the X chromosome with highly correlated methylation levels across all specimens. Because of the special role of DNA methylation on the X chromosome in dosage compensation, we excluded X-linked genes from further analyses. DNA methylation measurements from multiple probes representing different CpG sites associated with the same gene tended to be highly correlated (mean r value: 0.64 for related pairs of probes and 0.04 for unrelated pairs; see Figures S1 and S2).

[Figure omitted. See PDF.]

Figure 1. DNA methylation profiles of ovarian cancer cell lines and tumors.

DNA methylation profiles for 27 ovarian primary tumors and 15 ovarian cell lines are shown. Only the 1,184 CpG sites (686 genes) with variable DNA methylation levels across these specimens are shown (see Methods). Methylation values (β-values) were mean centered for each gene (across all specimens) then CpG sites and specimens were hierarchically clustered. Red indicates high methylation relative to the site-specific mean, green indicates low methylation relative to the mean. Rows of data corresponding to selected CpG sites (those with the greatest magnitude of difference between tumors and cell lines, or between pairs of tumor histologies) are indicated to the right of the panel. CL = Cell line, S = Serous tumor, E = Endometrioid tumor, CC = Clear cell tumor.

https://doi.org/10.1371/journal.pone.0009359.g001

We further explored differences in DNA methylation between cell lines and tumors, using a two-sample t-test to compare DNA methylation levels of selected CpG sites. Only the 1,110 CpG sites that met a minimal level of variability in methylation levels across these specimens (see Methods) were included in this analysis. We thus identified 489 CpG sites (associated with 337 genes) that varied significantly in DNA methylation level between the ovarian cancer derived cell lines and the primary ovarian cancer tumors at a false discovery rate (FDR) less than 1% [23] (see Figure 2 and Supplemental File S2). Among CpG sites whose mean methylation levels differed between cell lines and tumors, the vast majority were more highly methylated in the cell lines (445 of 489 CpG sites, 299 of 337 genes). The proposed cell of origin of ovarian carcinoma (Ovarian Surface Epithelium) is currently under debate [24]. Whole ovary, which is comprised of various cell types, including granulosa and follicular cells, is not an appropriate control. On the other hand, the normal cells that contaminate tumor samples include blood cells. Therefore, we compared DNA methylation levels in leukocytes (obtained from buffy coat) from two healthy females (age>60). The inclusion of leukocytes could also serve as a first step to assess the feasibility of blood-based biomarker development. For the majority of probes, DNA methylation levels in the leukocytes were similar to those observed in the tumor samples (see Figure 2). The elevated levels of DNA methylation in the leukocytes among the CpG sites that had elevated methylation in tumors relative to cell lines was particularly striking. The 39 CpG sites (associated with 35 genes) whose methylation levels differed most significantly between cell lines from tumors (FDR<0.1%) and had the largest differences in mean methylation levels between cell lines and tumors are shown in Table 3.

[Figure omitted. See PDF.]

Figure 2. Differential DNA methylation between ovarian cell lines and tumors.

DNA methylation levels in tumors, cell lines and buffy coat (leukocytes) are shown for the 489 autosomal CpG sites (337 genes) with significant differences in DNA methylation between cell lines and tumors. Methylation values (β-values) were mean centered for each gene (across tumor and cell line specimens). Data were hierarchically clustered only for the CpG sites. The displayed CpG sites were selected based on the adjusted p-values associated with differential expression between tumors and cell lines. Red indicates high methylation relative to the site-specific mean, green indicates low methylation relative to the mean. Rows of data corresponding to sites with the greatest magnitude of difference in methylation levels (β-values) between tumors and cell lines are labeled (top 5 sites in each direction). CL =  Cell line, S = Serous tumor, E = Endometrioid tumor, CC = Clear cell tumor, BC = Buffy coat.

https://doi.org/10.1371/journal.pone.0009359.g002

[Figure omitted. See PDF.]

Table 3. Loci with differential methylation between ovarian tumors and cell lines.

https://doi.org/10.1371/journal.pone.0009359.t003

To explore differences in DNA methylation profiles between tumor histologies, we compared the DNA methylation levels of 818 autosomal CpG sites selected for variability in DNA methylation levels across the tumors (see Methods). We applied a two-sample t-test to compare pairs of tumor histologies, and adjusted for multiple comparisons by controlling the false discovery rate [23] (see Methods for details). In this way, we identified 90 CpG sites (68 genes) that showed significant differences in methylation levels between at least one pair of tumor histologies at an FDR<1% (see Figure 3, Supplemental File S2). Note that a small fraction of the CpG sites with ‘subtype-specific’ DNA methylation were identified in more than one pairwise comparison (11 of 90 CpG sites; 12 of 68 genes). Serous and endometrioid tumors were distinguished by the largest number of CpG sites (70 CpG sites, 52 genes), of which 49 CpG sites (36 genes) were more highly methylated on average in endometrioid tumors. Nineteen CpG sites (18 genes) and 12 CpG sites (10 genes) differed significantly in mean methylation between clear cell and serous tumors, or clear cell and endometrioid tumors, respectively. Nearly all of the 26 CpG sites (23 genes) with differential methylation between clear cell tumors and serous or endometrioid tumors (or both), were more highly methylated on average in clear cell tumors (24 of 26 CpG sites, 21 of 23 genes). For each pairwise comparison, the CpG sites with the most significant differences in methylation level between the two tumor histologies are shown in Tables 4, 5 and 6 (see Supplemental Figure S3 for extent of overlap between CpG sites and genes across lists, and Supplemental File S3, for CpG island information and distance to transcription start site for the probes shown in Tables 4, 5 and 6).

[Figure omitted. See PDF.]

Figure 3. Differential DNA methylation among ovarian tumor histologies.

DNA methylation levels are shown for the 90 autosomal CpG sites (68 genes) with a significant difference in methylation between one or more pairs of tumor histologies. Methylation values (β-values) were mean centered for each gene (across tumors). Data were hierarchically clustered only for the CpG sites. The displayed CpG sites were selected based on the adjusted p-values associated with differential expression between tumors and cell lines. Red indicates high methylation relative to the site-specific mean, green indicates low methylation relative to the mean. Rows of data corresponding to sites with the greatest magnitude of difference in methylation levels (β-values) between pairs of tumor histologies are labeled (up to 3 sites in each direction for each pairwise comparison). S = Serous tumor, E = Endometrioid tumor, CC = Clear cell tumor.

https://doi.org/10.1371/journal.pone.0009359.g003

[Figure omitted. See PDF.]

Table 4. Loci with differential methylation between endometrioid versus serous tumors (FDR<0.1%).

https://doi.org/10.1371/journal.pone.0009359.t004

[Figure omitted. See PDF.]

Table 5. Loci with differential methylation between clear cell versus serous tumors (FDR<1%).

https://doi.org/10.1371/journal.pone.0009359.t005

[Figure omitted. See PDF.]

Table 6. Loci with differential methylation between clear cell versus endometrioid tumors (FDR<5%).

https://doi.org/10.1371/journal.pone.0009359.t006

Age-associated DNA hypermethylation is a well-established phenomenon [7]. Therefore, we examined age-associated methylation changes in the tumor specimens. After adjusting for histological type and controlling the false discovery rate [23], no CpG sites showed a significant association with age.

Relationship between Gene Expression and DNA Methylation Status

To investigate the relationship between DNA methylation status and gene expression, we measured levels of transcripts from the same set of genes in the same cell lines and primary tumors, using HEEBO, a 70mer oligonucleotide DNA microarray [25]. We analyzed the cell line and tumor data sets separately, and in each case restricted the analysis to autosomal genes selected for variation in methylation levels across the specimen set (see Methods). Promoter DNA methylation and levels of the corresponding gene's transcript have been found to be inversely correlated in some cases, but the relationship is not necessarily linear. Therefore, we used Spearman correlation, rather than Pearson correlation to examine this relationship for each locus in each of the two datasets (Figure 4). In cell lines, we found a clear trend towards a negative correlation between DNA methylation levels and gene expression for the 828 CpG sites whose methylation levels varied among the cell lines (Figure 4A, mean ρ = −0.22). We estimated that for at least 27% of these 828 CpG sites, increases in methylation were accompanied by decreases in transcription of the associated gene (Figure 4A). For the tumor tissue samples, the correlation between DNA methylation levels and the corresponding transcript abundance, for the 758 CpG sites with significant variation in methylation across tumors, was also negative, but much weaker (Figure 4B, mean ρ = −0.08). We estimated that for at least 16% of these 758 CpG sites, increases in methylation were accompanied by decreases in transcription of the associated gene (Figure 4B). As expected, CpG sites at which methylation was positively correlated with expression of the corresponding transcript were rare. Indeed, the proportion of CpG sites exhibiting positive correlation (ρ>0.5) between methylation values and expression values was significantly less than that observed in randomly paired CpG sites (in the data from ovarian cancer cell lines, approximately 42% fewer than would be expected if there were no interaction). Scatterplots of DNA methylation levels versus gene expression levels for the 8 CpG sites with the strongest negative correlations between these values are shown in Figure 5A and 5B, for cell lines and tumors, respectively. Notably, the association between methylation of a given CpG and transcription of the corresponding gene across the cell lines was only weakly correlated to the same association across the tumors, (Pearson r = 0.19) (Supplemental Figure S4).

[Figure omitted. See PDF.]

Figure 4. Correlations between DNA methylation and gene expression levels.

The distributions of Spearman correlations between DNA methylation levels (β-values) and gene expression levels [log2(sample/reference)]. Cell lines and tumors were analyzed separately and was restricted to autosomal CpG sites with variable DNA methylation levels across the relevant specimen set and a minimum number of non-missing values (see Methods). A) Cell lines: Density plots of observed and simulated Spearman correlations for all filtered loci. N = 828 CpG sites (502 genes); B) Tumors: Density plots of observed and simulated Spearman correlations for all filtered loci. N = 758 CpG sites (482 genes); C) Cell Lines: Density plot of observed and simulated Spearman correlations for loci meeting CpG island criteria. N = 585; D) Tumors: Density plot of observed and simulated Spearman correlations for loci meeting CpG island criteria. N = 492.

https://doi.org/10.1371/journal.pone.0009359.g004

[Figure omitted. See PDF.]

Figure 5. Top correlations between DNA methylation and gene expression.

The 8 CpG sites with the most negative correlations between DNA methylation and gene expression levels are shown for each of the two specimen sets: A) Cell lines; B) Tumors.

https://doi.org/10.1371/journal.pone.0009359.g005

Hypermethylation of CpG dinucleotides located in promoter 5′CpG islands has been shown to be associated with gene silencing [26]–[29]. We therefore examined the relationship between promoter DNA hypermethylation and gene expression in the context of CpG islands (as defined according to either the Takai and Jones criteria [30] or the Gardiner-Garden criteria [31]). In the cell lines, we found that the correlation between DNA methylation and gene expression levels tended to be more strongly negative among CpGs contained within CpG islands (mean β-values: −0.25 versus −0.12 for loci that met one or more CpG island criteria (n = 585) versus loci that met zero criteria respectively (n = 243), KS test p-value <0.001). We estimated that for at least 29% of the 585 loci meeting the CpG island criteria, increases in methylation were accompanied by decreases in transcription of the associated gene (Figure 4C). The difference in correlations in CpG island loci versus non-CpG island loci was not borne out in the tumor dataset (mean β-values: −0.08 versus −0.10 for loci that met one or more CpG islands criteria versus those that met zero, KS test p-value = 0.052). In tumors, we estimated that for at least 10% of the 492 loci meeting the CpG island criteria, increases in methylation were accompanied by decreases in transcription of the associated gene (Figure 4D). We found no consistent relationship between the distance of a methylation locus from the transcription start site and the strength of the correlation between its methylation levels and expression levels of the associated gene (not shown).

Discussion

We quantitatively assessed the DNA methylation status of 1,505 CpG sites, (associated with 808 genes), in 15 commonly used ovarian cancer cell lines and 27 primary ovarian tumors, 15 serous (55.5%), 9 endometrioid (33.3%) and 3 clear cell (11.1%). These proportions roughly reflect the relative frequencies of the different histologies of ovarian carcinoma diagnosed in the US (serous tumors 50%, endometrioid 25–30%, mucinous 10–15%, clear cell 5%, and others <5%) [32].

We found that the DNA methylation profiles of the ovarian cancer derived cell lines were distinctly different from those of the primary ovarian tumors. These results are consistent with reports describing vastly different gene expression patterns between cell lines and primary tumors in ovarian and other cancers [33]–[35]. While cell lines can be useful for molecular mechanistic studies, this study reinforces the need for skepticism in considering them as models for the cancers from which they originated, both in general and in the specific case of epigenetic profiling of ovarian cancer. One notable general difference between cell lines and tumors was the tendency for DNA methylation measurements to be higher in the cell lines than in the tumors. These findings are consistent with reports of an enhanced hypermethylation phenotype in cancer cell lines, compared to primary tumors for other types of cancer [36]. One prior study of ovarian cancer found fewer hypermethylated CpG islands in cell lines versus primary tumors [37]. Our results do not support this finding. Because the methylation sites represented in the arrays used for this study were biased towards sites known to be methylated in some cancers, this finding may be consequence of the fact that cell lines represent a ‘pure’ culture of cancerous cells whereas the primary tumors inevitably contained a mixture of tumor cells and other cells. This ‘normal contamination’ effect is supported by our finding that nearly all of the CpG sites that showed higher methylation in tumors than in cell lines also had elevated levels of DNA methylation in leukocytes from healthy women. Accumulation of epigenetic changes during prolonged cell culture may also contribute to the higher frequency of hypermethylation in the cell lines. This was observed for cancer cell lines [36], as well as hES cell lines and their derivatives [38], [39]. Furthermore, DNA methylation profiles of the hESC-derived neural precursor cells [38], [39] and hESC-derived liver cells [40] were shown to be different from their in vivo derived counterparts.

Although the numbers of some subtypes, such as clear cell tumors, were small, unsupervised and supervised analyses of the DNA methylation profiles of the primary tumors indicate clear differences between the histological types. We identified 90 CpG sites (associated with 68 genes) with significantly subtype-specific methylation patterns, the majority of which distinguished serous from endometrioid tumors. The distinctness of the DNA methylation profiles of the histological subtypes is suggestive of different tumorigenic mechanisms and/or cells of origin and underscores the need to view the different histological types of ovarian cancer as different diseases. While this notion is common knowledge among clinicians, and is reflected in studies of tumor gene expression [41] and protein expression [42], regrettably, it is often ignored in molecular diagnostic studies, and all cases of ovarian cancer are still treated similarly in clinic despite their acknowledged differences in pathology. Successful detection and treatment strategies for ovarian cancer are likely to be histology-specific, and will require a better understanding of the biological characteristics of each subtype.

As one of the largest studies to examine correlations between gene expression and DNA methylation, our results in this regard are noteworthy. We found for some genes that increased DNA methylation was correlated with decreased levels of transcription of the associated gene, but that there were also many genes for which such a relationship was not evident. We estimated that increases in methylation levels were associated with decreases in transcription of the corresponding gene for 27% of the 828 CpG sites with variable DNA methylation across 15 ovarian cancer cell lines (Figure 4A). Among loci in CpG islands, we estimated that increases in DNA methylation were accompanied by decreases in transcription of the associated gene for at least 29% of 585 CpG sites across cell lines (Figure 4C). Across all CpGs and CpG island loci, the association between DNA methylation and gene expression was much stronger among the cell lines than among the primary tumor tissues, most likely due to the confounding effect of cellular heterogeneity in the tumor specimens. Our findings are consistent with a recent study of large B-cell lymphoma [43], in which a wide range of correlations between gene expression and methylation were reported. Interestingly, all three of the genes evaluated in this study that were highlighted in the previous study as having a strong correlation between hypermethylation and expression showed a similar negative correlation in our cell line dataset (MGMT: −0.63, RBP1: −0.78, IGSF4: −0.69). Our results support the generality of the recognized relationship between promoter CpG island hypermethylation and gene silencing. However, the fact that this inverse relationship was not observed for the majority of loci, similar to another report [43], further highlights the complexity of this relationship.

The data presented here provide a first glimpse of the extent of variation in levels of DNA methylation among ovarian cancer specimens and as such provide a useful new lens through which to examine their defining characteristics. In principal, DNA methylation profiling could also be useful for discovery of biomarkers for detection of ovarian cancer, based on assays using sensitive PCR-based approaches for detection of methylated DNA [44]–[51]. However, even if we can identify a DNA-methylation marker that distinguishes tumors from normal cells with high sensitivity and specificity, additional criteria would need to be satisfied in order for this approach to be viable: adequate quantities of tumor-derived DNA would need to reach and persist in an assayable fluid (like blood or proximal fluids) in women with occult ovarian cancers still at a curable stage, and the abundance of the specifically methylated DNA sequences in fluids from these women would need to significantly exceed background levels in cancer-free individuals. These criteria seem difficult to satisfy in a blood-based early detection test for ovarian cancer, given the small size of the tumors that we need to detect [52], [53]. To this end, and as a primary step, we compared the DNA methylation profiles of normal leukocytes and ovarian tumors. But the potential of DNA methylation markers as early-detection markers using proximal fluids or for detection of recurrent cancer in a blood-based assay remains to be evaluated. Additional profiling of primary tumors and corresponding normal tissues (fallopian tube, ovarian surface epithelium, peritoneum or endometrium) as well as measurements of methylated DNA in candidate detection media (blood or proximal fluids) in both health and disease will help define the natural history of ovarian cancers and critically assess the potential utility of DNA methylation based biomarkers in combating ovarian cancer.

Materials and Methods

Sample Collection and Processing

Ovarian tumors.

Tumor tissues were obtained from patients with informed consent, which was approved by the Institutional Review Board of the Swedish Hospital in Seattle. Tissue acquisition, processing and storage were conducted by the Pacific Ovarian Cancer Research Consortium (POCRC) at the Fred Hutchinson Cancer Research Center. The quality of tumor tissues, the percentage of cancer cells and the histological analysis of the tumors were determined via a centralized pathology review at the POCRC repository. All samples were distributed in a de-identified manner and cannot be traced back to patients. Tumor and patient characteristics are summarized in Table 1.

Cell lines.

The sources and characteristics of the 15 ovarian cancer-derived cell lines analyzed are summarized in Table 2. All cell lines were cultured according to the specifications outlined by ATCC. Cell lines were grown to 80% confluence, then serum-deprived in 0.5% FBS for 24 hours before harvesting. Cells were washed twice with Hank's Balanced Salt Solution (Invitrogen, Carlsbad, California) and drained. 1 mL of SDS extraction buffer (0.1M NaCl, 20 mM Trizma base, 25 mM EDTA, 0.5% w/v SDS) was added directly to the plate to lyse the cells. Cellular lysate was scraped off the plates, pipetted into a cryovial, and stored at −80°C until DNA extraction. Lysates were treated with 200 ug/ml ProteinaseK at 50°C overnight. DNA was precipitated with one volume of Isopropanol and dissolved in TE−4 buffer.

Buffy coat.

Whole blood from two females (age>60) was purchased from HemaCare Corporation. Blood was collected in anticoagulant containing blood collection bags and kept on ice. Samples were processed immediately upon receipt.

DNA and RNA extraction from ovarian tumors.

The tumor tissues were disrupted in a Biospec Tissue-Tearor Model 985370-395 in the presence of TRIzol (Invitrogen) with a ratio of 10 ml TRIzol per 50 mg of tissue. Following an extraction step of 0.2 ml of Chloroform per 1 ml of TRIzol, the upper aqueous phase was removed for RNA isolation. The remaining inter-phase and the phenol phase was briefly stored in 4°C for DNA isolation. The upper aqueous phase was processed for total RNA according to Invitrogen protocol and further purified using RNeasy Mini kit (Qiagen). Genomic DNA was isolated by back-extraction from the TRIzol inter-phase and the phenol phase as described on Ambion website (page 4) [54]. Detailed amplification and labeling protocols are available at the Brown lab website [55].

DNA extraction from buffy coat.

Blood samples were centrifuged for 10 min at 300 g in 50 ml conical tubes, followed by another 10 min centrifugation at 1,600 g. Buffy coat was isolated and DNA was extracted using Qiagen QIAamp DNA Blood Kit and stored in −30 C until sodium bisulfite conversion.

RNA extraction from ovarian cell lines.

Cell lines were lysed in Phenol followed by a Phenol/Chloroform extraction. The aqueous layer was added to 70% ethanol and further purified by loading the RNA/Ethanol mix into an RNeasy Mini kit (Qiagen) column. Purification was continued and finalized as described in the RNeasy kit (Qiagen) Manual. Detailed amplification and labeling protocols are available at the Brown lab website [55].

DNA Methylation Analysis

DNA methylation assay platform.

DNA methylation analysis was performed using Illumina's Golden Gate Cancer Panel 1 bead array [16]. These arrays contain DNA probes representing 1505 CpG islands near 808 genes. Probes were selected for cancer-related genes based on data from Illumina's scientific collaborators, literature, and public databases.

Experimental methods.

Briefly, genomic DNA from 27 primary ovarian tumors and 15 ovarian cell lines was subject to sodium bisulfite conversion, labeled with fluorescent dyes, and hybridized to Illumina Golden Gate Cancer Panel 1 bead arrays at the University of Southern California Epigenome Center. We utilized a number of DNA methylation control reactions to assess the extent of bisulfile conversion completion as explained in detail elsewhere [56]. Methods and reagents for sodium bisulfite conversion and array hybridizations were as previously described [8], [16], [57].

Data processing.

DNA methylation values from the bead array hybridizations were scored as β-values [16]. β-values for each probe in each experiment were calculated as the ratio of the methylated signal over the total fluorescent signal on a scale between 0 and 1 [16]. β-values thus indicate a locus' observed methylation relative to the maximum potential methylation of that site. For quality control, we obtained non-background corrected signal intensities of the methylated (M) and unmethylated (U) probes and the mean negative control cy5 (red) and cy3 (green) signal intensities. We also obtained the Detection p-values, which provide an indication of DNA methylation measurement quality for each locus and are calculated based on the difference in signal intensity of each probe compared to the set of negative control probes. Measurements for which the detection p-value is less than 0.05 are considered to have a signal intensity significantly above background. Data points with a detection p-value>0.05 are masked as “NA”, and represent beta values with non-significant detection of DNA methylation compared to background.

Gene Expression Analysis

Assay platform.

We measured transcript levels in the 27 primary ovarian tumors and 15 cell lines using the HEEBO [16] oligonucleotide microarrays, containing 44,544 70-mer probes and printed at the Stanford Functional Genomics Facility. Detailed information on the HEEBO arrays is available at the Stanford Functional Genomics Facility website [25].

Experimental methods.

Microarray experiments were performed as described at the Brown lab website [58]. Briefly, 500 µg of total RNA from each of the ovarian tumors were amplified using Amino Allyl Message Amp™ II aRNA Kit (Ambion, Austin, TX, USA). The aRNAs were labeled with Cy5 and co-hybridized with Cy3 labeled Stratagene reference aRNA. For some samples, the mRNA was amplified and hybridized in duplicate or triplicate. The samples were then hybridized to HEEBO microarrays. The arrays were scanned in a low-ozone environment using a GenePix 4000A microarray scanner and images were analyzed with Genepix 5.0 (Axon instruments, Union City, CA). The raw data were deposited into Stanford Microarray Database [59] and can be viewed at the SMD website [60].

Data processing.

1) By Spot: A spot quality filter was applied to every spot on each array. Only spots with a Cy5/Cy3 ratio of intensity/background >1.5 in either channel were included. Spots that were flagged manually for poor quality by visual inspection of spot uniformity were excluded. The log (base2) of Cy5/Cy3 normalized ratio (mean) was calculated for each spot. Ratios were normalized for each array such that the distribution of log ratios for the array had a mean of 0.0. Empty and control spots were dropped after normalization. Spots with the same Entrez Gene ID were averaged. 2) By Array: Replicate hybridizations of the same amplified mRNA were averaged by Entrez Gene ID. Hybridizations of independent amplifications of the same sample mRNA were averaged by Entrez Gene ID.

Statistical Methods

DNA methylation clustering and visualization.

To visualize the DNA methylation patterns across all ovarian cancer cell lines and tumors, we first filtered for genes with some variation across the sample set. Specifically, we restricted our analysis to CpG sites with minimum β-value greater than 0.2, ratio of maximum β-value to minimum β-value greater than or equal to 2 and 1 or fewer missing β-values. We then mean-centered the selected 1184 β-values by gene and applied two-way average-linkage unsupervised hierarchical clustering on the resulting data. Note that while this data are filtered, no comparisons were made and data were clustered for both samples and CpG sites. We visualized the clusters in Java Treeview [61], [62].

Identification of loci with biopecimen-specific DNA methylation levels (cell line versus tumor).

To identify loci that showed significantly different levels of DNA methylation in cell lines versus tumors, we first restricted our analysis to the 1,110 autosomal CpG sites whose β-values were variable across the set of 15 cell lines and 27 tumors. We defined variable CpG sites as those CpG sites with minimum β-value greater than 0.2, ratio of maximum β-value to minimum β-value greater than or equal to 2 and 1 or fewer missing β-values. We then performed a two-sample t-test comparing the 15 cell lines to the 27 tumors. P-values from the t-test were adjusted for multiple comparisons by controlling the false discovery rate [23]. All statistical analyses were performed using the R software package [63].

Identification of loci with histological subtype-specific DNA methylation levels.

To identify ovarian cancer subtype-specific loci, we first restricted our analysis to the 818 autosomal CpG sites whose β-values were variable among the 27 primary tumors. We defined variable CpG sites as those CpG sites with minimum β-value greater than 0.2, ratio of maximum β-value to minimum β-value greater than or equal to 2 and 1 or fewer missing β-values. We performed pair-wise two-sample t-tests comparing each subtype independently to the other subtypes (serous vs. endometrioid, serous vs. clear cell, endometrioid vs. clear cell). P-values from each pair-wise analysis were adjusted for multiple comparisons by controlling the false discovery rate [23]. All statistical analyses were performed using the R software package [63].

Association of DNA methylation and gene expression.

We paired loci used in methylation analysis with gene expression data from the appropriate loci according to Entrez Gene ID. For gene expression data, we used the normalized log ratios of Cy5/Cy3 (see Microarray Methods). For DNA methylation, normalized β-values were used. We excluded X-linked loci from this analysis because of the X-chromosome associated dosage compensation mechanisms could cause a decoupling of methylation levels and gene expression levels. We further restricted our analysis to loci with at least 5 (of 15) non-missing gene expression measurements in the cell line analysis and 9 (of 27) non-missing gene expression measurements in the tumor analysis. For the 758 CpG sites that were sufficiently variable across all tumor specimens (see filtering above), we calculated the Spearman correlation between DNA methylation levels and gene expression levels. Similarly, we calculated Spearman correlations for the 828 CpG sites that were sufficiently variable across all cell line samples (see filtering above). For each dataset, the false discovery rate was estimated by calculating the median number of loci with correlations below a given threshold over 500 permutations of the gene labels. All statistical analyses were performed using the R software package [63].

Supporting Information

[Figure omitted. See PDF.]

File S1.

Methylation array probe annotations and associated methylation and gene expression data. List of probes and associated annotations (locus, gene, distance from transcription start, etc) for the Illumina bead array used to profile specimens. Observed methylation values (β-values) for each probe and specimen are recorded as are the gene expression levels (log base2 ratio) of the associated transcripts.

https://doi.org/10.1371/journal.pone.0009359.s001

(2.20 MB XLS)

File S2.

Loci with differential DNA methylation. List of CpG sites with significant differences in DNA methylation between cell lines and tumors or between tumor histologies. Associated probe data are also included.

https://doi.org/10.1371/journal.pone.0009359.s002

(0.40 MB XLS)

File S3.

CpG island information and distance from the transcription start site of probes with differential methylation between tumor histologies.

https://doi.org/10.1371/journal.pone.0009359.s003

(0.06 MB XLS)

Figure S1.

Correlation in DNA methylation values between pairs of probes. The distribution of Pearson correlations is shown for i) pairs of related probes (2 CpG sites, 1 gene) ii) pairs of unrelated probes (2 CpG sites, 2 genes). For related probes, only the 1184 probes (686 genes) which exhibited sufficient variation across the 42 specimens (see Methods) were included.

https://doi.org/10.1371/journal.pone.0009359.s004

(0.36 MB TIF)

Figure S2.

Correlation matrix for DNA methylation values across loci. Pearson correlations in methylation values for pairs of CpG sites are shown colorimetrically. CpG sites are sorted by their location in the genome. The 1,184 CpG sites (686 genes) selected for variation in methylation across all specimens were used (Figure 1). Correlation calculations were based on 42 specimens. Red indicates positive correlation, green indicates negative correlation.

https://doi.org/10.1371/journal.pone.0009359.s005

(5.37 MB PDF)

Figure S3.

Overlap between tumor histology lists. A) Venn diagram showing number of CpG sites in each histology-specific list and overlap between lists; B) Venn diagram showing number of genes in each histology-specific list and overlap between lists. S = Serous tumor, E = Endometrioid tumor, CC = Clear cell tumor.

https://doi.org/10.1371/journal.pone.0009359.s006

(0.13 MB TIF)

Figure S4.

Methylation-gene expression correlation across cell lines versus methylation-gene expression correlation across tumors.

https://doi.org/10.1371/journal.pone.0009359.s007

(0.17 MB TIF)

Appendix S1.

Pacific Ovarian Cancer Research Consortium Menopausal Determination

https://doi.org/10.1371/journal.pone.0009359.s008

(0.02 MB DOC)

Acknowledgments

We thank Kathy O'Briant for her generous help with the collection and distribution of the samples. We also thank Tiffany I. Long, Dr. Dan Weisenberger and Dr. Fei Pan for their assistance in conducting and analyzing the Illumina GoldenGate assays. POB is an investigator of the Howard Hughes Medical Institute.

Author Contributions

Conceived and designed the experiments: POB PWL. Performed the experiments: MC MNO APV BSK. Analyzed the data: SH SH CP. Contributed reagents/materials/analysis tools: CWD NDU POB PWL. Wrote the paper: SH.

Citation: Houshdaran S, Hawley S, Palmer C, Campan M, Olsen MN, Ventura AP, et al. (2010) DNA Methylation Profiles of Ovarian Epithelial Carcinoma Tumors and Cell Lines. PLoS ONE 5(2): e9359. https://doi.org/10.1371/journal.pone.0009359

References

1. Jemal A, Siegel R, Ward E, Hao Y, Xu J, et al. (2008) Cancer statistics, 2008. CA Cancer J Clin 58: 71–96.A. JemalR. SiegelE. WardY. HaoJ. Xu2008Cancer statistics, 2008.CA Cancer J Clin587196

2. Boumil RM, Ogawa Y, Sun BK, Huynh KD, Lee JT (2006) Differential methylation of Xite and CTCF sites in Tsix mirrors the pattern of X-inactivation choice in mice. Mol Cell Biol 26: 2109–2117.RM BoumilY. OgawaBK SunKD HuynhJT Lee2006Differential methylation of Xite and CTCF sites in Tsix mirrors the pattern of X-inactivation choice in mice.Mol Cell Biol2621092117

3. Zuccotti M, Monk M (1995) Methylation of the mouse Xist gene in sperm and eggs correlates with imprinted Xist expression and paternal X-inactivation. Nat Genet 9: 316–320.M. ZuccottiM. Monk1995Methylation of the mouse Xist gene in sperm and eggs correlates with imprinted Xist expression and paternal X-inactivation.Nat Genet9316320

4. Barr H, Hermann A, Berger J, Tsai HH, Adie K, et al. (2007) Mbd2 contributes to DNA methylation-directed repression of the Xist gene. Mol Cell Biol 27: 3750–3757.H. BarrA. HermannJ. BergerHH TsaiK. Adie2007Mbd2 contributes to DNA methylation-directed repression of the Xist gene.Mol Cell Biol2737503757

5. Walsh CP, Bestor TH (1999) Cytosine methylation and mammalian development. Genes Dev 13: 26–34.CP WalshTH Bestor1999Cytosine methylation and mammalian development.Genes Dev132634

6. Li E, Beard C, Jaenisch R (1993) Role for DNA methylation in genomic imprinting. Nature 366: 362–365.E. LiC. BeardR. Jaenisch1993Role for DNA methylation in genomic imprinting.Nature366362365

7. Issa JP, Ottaviano YL, Celano P, Hamilton SR, Davidson NE, et al. (1994) Methylation of the oestrogen receptor CpG island links ageing and neoplasia in human colon. Nat Genet 7: 536–540.JP IssaYL OttavianoP. CelanoSR HamiltonNE Davidson1994Methylation of the oestrogen receptor CpG island links ageing and neoplasia in human colon.Nat Genet7536540

8. Houshdaran S, Cortessis VK, Siegmund K, Yang A, Laird PW, et al. (2007) Widespread epigenetic abnormalities suggest a broad DNA methylation erasure defect in abnormal human sperm. PLoS ONE 2: e1289.S. HoushdaranVK CortessisK. SiegmundA. YangPW Laird2007Widespread epigenetic abnormalities suggest a broad DNA methylation erasure defect in abnormal human sperm.PLoS ONE2e1289

9. Jones PA, Laird PW (1999) Cancer epigenetics comes of age. Nat Genet 21: 163–167.PA JonesPW Laird1999Cancer epigenetics comes of age.Nat Genet21163167

10. Cheng P, Schmutte C, Cofer KF, Felix JC, Yu MC, et al. (1997) Alterations in DNA methylation are early, but not initial, events in ovarian tumorigenesis. Br J Cancer 75: 396–402.P. ChengC. SchmutteKF CoferJC FelixMC Yu1997Alterations in DNA methylation are early, but not initial, events in ovarian tumorigenesis.Br J Cancer75396402

11. Strathdee G, MacKean MJ, Illand M, Brown R (1999) A role for methylation of the hMLH1 promoter in loss of hMLH1 expression and drug resistance in ovarian cancer. Oncogene 18: 2335–2341.G. StrathdeeMJ MacKeanM. IllandR. Brown1999A role for methylation of the hMLH1 promoter in loss of hMLH1 expression and drug resistance in ovarian cancer.Oncogene1823352341

12. Yu Y, Xu F, Peng H, Fang X, Zhao S, et al. (1999) NOEY2 (ARHI), an imprinted putative tumor suppressor gene in ovarian and breast carcinomas. Proc Natl Acad Sci U S A 96: 214–219.Y. YuF. XuH. PengX. FangS. Zhao1999NOEY2 (ARHI), an imprinted putative tumor suppressor gene in ovarian and breast carcinomas.Proc Natl Acad Sci U S A96214219

13. Catteau A, Harris WH, Xu CF, Solomon E (1999) Methylation of the BRCA1 promoter region in sporadic breast and ovarian cancer: correlation with disease characteristics. Oncogene 18: 1957–1965.A. CatteauWH HarrisCF XuE. Solomon1999Methylation of the BRCA1 promoter region in sporadic breast and ovarian cancer: correlation with disease characteristics.Oncogene1819571965

14. Baldwin RL, Nemeth E, Tran H, Shvartsman H, Cass I, et al. (2000) BRCA1 promoter region hypermethylation in ovarian carcinoma: a population-based study. Cancer Res 60: 5329–5333.RL BaldwinE. NemethH. TranH. ShvartsmanI. Cass2000BRCA1 promoter region hypermethylation in ovarian carcinoma: a population-based study.Cancer Res6053295333

15. Barton CA, Hacker NF, Clark SJ, O'Brien PM (2008) DNA methylation changes in ovarian cancer: implications for early diagnosis, prognosis and treatment. Gynecol Oncol 109: 129–139.CA BartonNF HackerSJ ClarkPM O'Brien2008DNA methylation changes in ovarian cancer: implications for early diagnosis, prognosis and treatment.Gynecol Oncol109129139

16. Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, et al. (2006) High-throughput DNA methylation profiling using universal bead arrays. Genome Res 16: 383–393.M. BibikovaZ. LinL. ZhouE. ChudinEW Garcia2006High-throughput DNA methylation profiling using universal bead arrays.Genome Res16383393

17. Christensen BC, Houseman EA, Marsit CJ, Zheng S, Wrensch MR, et al. (2009) Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context. PLoS Genet 5: e1000602.BC ChristensenEA HousemanCJ MarsitS. ZhengMR Wrensch2009Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context.PLoS Genet5e1000602

18. Katari S, Turan N, Bibikova M, Erinle O, Chalian R, et al. (2009) DNA methylation and gene expression differences in children conceived in vitro or in vivo. Hum Mol Genet 18: 3769–3778.S. KatariN. TuranM. BibikovaO. ErinleR. Chalian2009DNA methylation and gene expression differences in children conceived in vitro or in vivo.Hum Mol Genet1837693778

19. Ladd-Acosta C, Pevsner J, Sabunciyan S, Yolken RH, Webster MJ, et al. (2007) DNA methylation signatures within the human brain. Am J Hum Genet 81: 1304–1315.C. Ladd-AcostaJ. PevsnerS. SabunciyanRH YolkenMJ Webster2007DNA methylation signatures within the human brain.Am J Hum Genet8113041315

20. Martinez R, Martin-Subero JI, Rohde V, Kirsch M, Alaminos M, et al. (2009) A microarray-based DNA methylation study of glioblastoma multiforme. Epigenetics 4: 255–264.R. MartinezJI Martin-SuberoV. RohdeM. KirschM. Alaminos2009A microarray-based DNA methylation study of glioblastoma multiforme.Epigenetics4255264

21. Bibikova M, Chudin E, Wu B, Zhou L, Garcia EW, et al. (2006) Human embryonic stem cells have a unique epigenetic signature. Genome Res 16: 1075–1083.M. BibikovaE. ChudinB. WuL. ZhouEW Garcia2006Human embryonic stem cells have a unique epigenetic signature.Genome Res1610751083

22. The Cancer Genome Atlas Research Network (2008) Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 455: 1061–1068.The Cancer Genome Atlas Research Network 2008Comprehensive genomic characterization defines human glioblastoma genes and core pathways.Nature45510611068

23. Klipper-Aurbach Y, Wasserman M, Braunspiegel-Weintrob N, Borstein D, Peleg S, et al. (1995) Mathematical formulae for the prediction of the residual beta cell function during the first two years of disease in children and adolescents with insulin-dependent diabetes mellitus. Med Hypotheses 45: 486–490.Y. Klipper-AurbachM. WassermanN. Braunspiegel-WeintrobD. BorsteinS. Peleg1995Mathematical formulae for the prediction of the residual beta cell function during the first two years of disease in children and adolescents with insulin-dependent diabetes mellitus.Med Hypotheses45486490

24. Dubeau L (2008) The cell of origin of ovarian epithelial tumours. Lancet Oncol 9: 1191–1197.L. Dubeau2008The cell of origin of ovarian epithelial tumours.Lancet Oncol911911197

25. Stanford Functional Genomics Facility website (2009) Stanford Functional Genomics Facility website2009http://microarray.org/sfgf/heebo.do. http://microarray.org/sfgf/heebo.do.

26. Felsenfeld G, Groudine M (2003) Controlling the double helix. Nature 421: 448–453.G. FelsenfeldM. Groudine2003Controlling the double helix.Nature421448453

27. Chiang JW, Karlan BY, Cass L, Baldwin RL (2006) BRCA1 promoter methylation predicts adverse ovarian cancer prognosis. Gynecol Oncol 101: 403–410.JW ChiangBY KarlanL. CassRL Baldwin2006BRCA1 promoter methylation predicts adverse ovarian cancer prognosis.Gynecol Oncol101403410

28. Roman-Gomez J, Jimenez-Velasco A, Castillejo JA, Agirre X, Barrios M, et al. (2004) Promoter hypermethylation of cancer-related genes: a strong independent prognostic factor in acute lymphoblastic leukemia. Blood 104: 2492–2498.J. Roman-GomezA. Jimenez-VelascoJA CastillejoX. AgirreM. Barrios2004Promoter hypermethylation of cancer-related genes: a strong independent prognostic factor in acute lymphoblastic leukemia.Blood10424922498

29. Brock MV, Gou M, Akiyama Y, Muller A, Wu TT, et al. (2003) Prognostic importance of promoter hypermethylation of multiple genes in esophageal adenocarcinoma. Clin Cancer Res 9: 2912–2919.MV BrockM. GouY. AkiyamaA. MullerTT Wu2003Prognostic importance of promoter hypermethylation of multiple genes in esophageal adenocarcinoma.Clin Cancer Res929122919

30. Takai D, Jones PA (2002) Comprehensive analysis of CpG islands in human chromosomes 21 and 22. Proc Natl Acad Sci U S A 99: 3740–3745.D. TakaiPA Jones2002Comprehensive analysis of CpG islands in human chromosomes 21 and 22.Proc Natl Acad Sci U S A9937403745

31. Gardiner-Garden M, Frommer M (1987) CpG islands in vertebrate genomes. J Mol Biol 196: 261–282.M. Gardiner-GardenM. Frommer1987CpG islands in vertebrate genomes.J Mol Biol196261282

32. NCI SEER website (2009) NCI SEER website2009http://seer.cancer.gov/statfacts/html/ovary.html. http://seer.cancer.gov/statfacts/html/ovary.html.

33. Hough CD, Sherman-Baust CA, Pizer ES, Montz FJ, Im DD, et al. (2000) Large-scale serial analysis of gene expression reveals genes differentially expressed in ovarian cancer. Cancer Res 60: 6281–6287.CD HoughCA Sherman-BaustES PizerFJ MontzDD Im2000Large-scale serial analysis of gene expression reveals genes differentially expressed in ovarian cancer.Cancer Res6062816287

34. Furlong MT, Hough CD, Sherman-Baust CA, Pizer ES, Morin PJ (1999) Evidence for the colonic origin of ovarian cancer cell line SW626. J Natl Cancer Inst 91: 1327–1328.MT FurlongCD HoughCA Sherman-BaustES PizerPJ Morin1999Evidence for the colonic origin of ovarian cancer cell line SW626.J Natl Cancer Inst9113271328

35. Ross DT, Scherf U, Eisen MB, Perou CM, Rees C, et al. (2000) Systematic variation in gene expression patterns in human cancer cell lines. Nat Genet 24: 227–235.DT RossU. ScherfMB EisenCM PerouC. Rees2000Systematic variation in gene expression patterns in human cancer cell lines.Nat Genet24227235

36. Smiraglia DJ, Rush LJ, Fruhwald MC, Dai Z, Held WA, et al. (2001) Excessive CpG island hypermethylation in cancer cell lines versus primary human malignancies. Hum Mol Genet 10: 1413–1419.DJ SmiragliaLJ RushMC FruhwaldZ. DaiWA Held2001Excessive CpG island hypermethylation in cancer cell lines versus primary human malignancies.Hum Mol Genet1014131419

37. Ahluwalia A, Yan P, Hurteau JA, Bigsby RM, Jung SH, et al. (2001) DNA methylation and ovarian cancer. I. Analysis of CpG island hypermethylation in human ovarian cancer using differential methylation hybridization. Gynecol Oncol 82: 261–268.A. AhluwaliaP. YanJA HurteauRM BigsbySH Jung2001DNA methylation and ovarian cancer. I. Analysis of CpG island hypermethylation in human ovarian cancer using differential methylation hybridization.Gynecol Oncol82261268

38. Shen Y, Chow J, Wang Z, Fan G (2006) Abnormal CpG island methylation occurs during in vitro differentiation of human embryonic stem cells. Hum Mol Genet 15: 2623–2635.Y. ShenJ. ChowZ. WangG. Fan2006Abnormal CpG island methylation occurs during in vitro differentiation of human embryonic stem cells.Hum Mol Genet1526232635

39. Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, et al. (2008) Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature 454: 766–770.A. MeissnerTS MikkelsenH. GuM. WernigJ. Hanna2008Genome-scale DNA methylation maps of pluripotent and differentiated cells.Nature454766770

40. Brunner AL, Johnson DS, Kim SW, Valouev A, Reddy TE, et al. (2009) Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver. Genome Res 19: 1044–1056.AL BrunnerDS JohnsonSW KimA. ValouevTE Reddy2009Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver.Genome Res1910441056

41. Zorn KK, Bonome T, Gangi L, Chandramouli GV, Awtrey CS, et al. (2005) Gene expression profiles of serous, endometrioid, and clear cell subtypes of ovarian and endometrial cancer. Clin Cancer Res 11: 6422–6430.KK ZornT. BonomeL. GangiGV ChandramouliCS Awtrey2005Gene expression profiles of serous, endometrioid, and clear cell subtypes of ovarian and endometrial cancer.Clin Cancer Res1164226430

42. Kobel M, Kalloger SE, Boyd N, McKinney S, Mehl E, et al. (2008) Ovarian carcinoma subtypes are different diseases: implications for biomarker studies. PLoS Med 5: e232.M. KobelSE KallogerN. BoydS. McKinneyE. Mehl2008Ovarian carcinoma subtypes are different diseases: implications for biomarker studies.PLoS Med5e232

43. Pike BL, Greiner TC, Wang X, Weisenburger DD, Hsu YH, et al. (2008) DNA methylation profiles in diffuse large B-cell lymphoma and their relationship to gene expression status. Leukemia 22: 1035–1043.BL PikeTC GreinerX. WangDD WeisenburgerYH Hsu2008DNA methylation profiles in diffuse large B-cell lymphoma and their relationship to gene expression status.Leukemia2210351043

44. Kawakami K, Brabender J, Lord RV, Groshen S, Greenwald BD, et al. (2000) Hypermethylated APC DNA in plasma and prognosis of patients with esophageal adenocarcinoma. J Natl Cancer Inst 92: 1805–1811.K. KawakamiJ. BrabenderRV LordS. GroshenBD Greenwald2000Hypermethylated APC DNA in plasma and prognosis of patients with esophageal adenocarcinoma.J Natl Cancer Inst9218051811

45. Esteller M, Sanchez-Cespedes M, Rosell R, Sidransky D, Baylin SB, et al. (1999) Detection of aberrant promoter hypermethylation of tumor suppressor genes in serum DNA from non-small cell lung cancer patients. Cancer Res 59: 67–70.M. EstellerM. Sanchez-CespedesR. RosellD. SidranskySB Baylin1999Detection of aberrant promoter hypermethylation of tumor suppressor genes in serum DNA from non-small cell lung cancer patients.Cancer Res596770

46. Wong IH, Lo YM, Zhang J, Liew CT, Ng MH, et al. (1999) Detection of aberrant p16 methylation in the plasma and serum of liver cancer patients. Cancer Res 59: 71–73.IH WongYM LoJ. ZhangCT LiewMH Ng1999Detection of aberrant p16 methylation in the plasma and serum of liver cancer patients.Cancer Res597173

47. Wong IH, Lo YM, Yeo W, Lau WY, Johnson PJ (2000) Frequent p15 promoter methylation in tumor and peripheral blood from hepatocellular carcinoma patients. Clin Cancer Res 6: 3516–3521.IH WongYM LoW. YeoWY LauPJ Johnson2000Frequent p15 promoter methylation in tumor and peripheral blood from hepatocellular carcinoma patients.Clin Cancer Res635163521

48. Fleischhacker M (1999) Meeting report - First international symposium on circulating nucleic acids in plasma/serum - Implications in cancer diagnostics, prognosis or follow up and in prenatal diagnosis - Menthon France, August 18-20, 1999. Eur J Med Res 4: 488–490.M. Fleischhacker1999Meeting report - First international symposium on circulating nucleic acids in plasma/serum - Implications in cancer diagnostics, prognosis or follow up and in prenatal diagnosis - Menthon France, August 18-20, 1999.Eur J Med Res4488490

49. Leon SA, Shapiro B, Sklaroff DM, Yaros MJ (1977) Free DNA in the serum of cancer patients and the effect of therapy. Cancer Res 37: 646–650.SA LeonB. ShapiroDM SklaroffMJ Yaros1977Free DNA in the serum of cancer patients and the effect of therapy.Cancer Res37646650

50. Goessl C, Krause H, Muller M, Heicappell R, Schrader M, et al. (2000) Fluorescent methylation-specific polymerase chain reaction for DNA-based detection of prostate cancer in bodily fluids. Cancer Res 60: 5941–5945.C. GoesslH. KrauseM. MullerR. HeicappellM. Schrader2000Fluorescent methylation-specific polymerase chain reaction for DNA-based detection of prostate cancer in bodily fluids.Cancer Res6059415945

51. Weisenberger DJ, Trinh BN, Campan M, Sharma S, Long TI, et al. (2008) DNA methylation analysis by digital bisulfite genomic sequencing and digital MethyLight. Nucleic Acids Res 36: 4689–4698.DJ WeisenbergerBN TrinhM. CampanS. SharmaTI Long2008DNA methylation analysis by digital bisulfite genomic sequencing and digital MethyLight.Nucleic Acids Res3646894698

52. Brown PO, Palmer C (2009) The preclinical natural history of serous ovarian cancer: defining the target for early detection. PLoS Med 6: e1000114.PO BrownC. Palmer2009The preclinical natural history of serous ovarian cancer: defining the target for early detection.PLoS Med6e1000114

53. Lutz AM, Willmann JK, Cochran FV, Ray P, Gambhir SS (2008) Cancer screening: a mathematical model relating secreted blood biomarker levels to tumor sizes. PLoS Med 5: e170.AM LutzJK WillmannFV CochranP. RaySS Gambhir2008Cancer screening: a mathematical model relating secreted blood biomarker levels to tumor sizes.PLoS Med5e170

54. Ambion website (2009) Ambion website2009http://www.ambion.com/techlib/append/supp/tri.pdf. http://www.ambion.com/techlib/append/supp/tri.pdf.

55. Brown Laboratory website (2009) Brown Laboratory website2009http://cmgm.stanford.edu/pbrown/protocols/Amplified_Protocol1.html. http://cmgm.stanford.edu/pbrown/protocols/Amplified_Protocol1.html.

56. Campan M, Weisenberger DJ, Trinh B, Laird PW (2009) MethyLight. Methods Mol Biol 507: 325–337.M. CampanDJ WeisenbergerB. TrinhPW Laird2009MethyLight.Methods Mol Biol507325337

57. illumina website (2009) illumina website2009www.illumina.com. www.illumina.com.

58. Brown Laboratory website (2009) Brown Laboratory website2009http://cmgm.stanford.edu/pbrown/protocols/Array_Hyb_Wash1.html. http://cmgm.stanford.edu/pbrown/protocols/Array_Hyb_Wash1.html.

59. Sherlock G, Hernandez-Boussard T, Kasarskis A, Binkley G, Matese JC, et al. (2001) The Stanford Microarray Database. Nucleic Acids Res 29: 152–155.G. SherlockT. Hernandez-BoussardA. KasarskisG. BinkleyJC Matese2001The Stanford Microarray Database.Nucleic Acids Res29152155

60. Stanford MicroArray website (2009) Stanford MicroArray website2009http://genome-www4.stanford.edu/MicroArray/SMD/index.html. http://genome-www4.stanford.edu/MicroArray/SMD/index.html.

61. Eisen Laboratory website (2009) Eisen Laboratory website2009http://rana.lbl.gov/EisenSoftware.htm. http://rana.lbl.gov/EisenSoftware.htm.

62. Saldanha AJ (2004) Java Treeview–extensible visualization of microarray data. Bioinformatics 20: 3246–3248.AJ Saldanha2004Java Treeview–extensible visualization of microarray data.Bioinformatics2032463248

63. (2008) R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2008R: A language and environment for statistical computing.R Foundation for Statistical Computing

64. Lowe KA, Shah C, Wallace E, Anderson G, Paley P, et al. (2008) Effects of personal characteristics on serum CA125, mesothelin, and HE4 levels in healthy postmenopausal women at high-risk for ovarian cancer. Cancer Epidemiol Biomarkers Prev 17: 2480–2487.KA LoweC. ShahE. WallaceG. AndersonP. Paley2008Effects of personal characteristics on serum CA125, mesothelin, and HE4 levels in healthy postmenopausal women at high-risk for ovarian cancer.Cancer Epidemiol Biomarkers Prev1724802487

65. DiSaia PJ, Sinkovics JG, Rutledge FN, Smith JP (1972) Cell-mediated immunity to human malignant cells. A brief review and further studies with two gynecologic tumors. Am J Obstet Gynecol 114: 979–989.PJ DiSaiaJG SinkovicsFN RutledgeJP Smith1972Cell-mediated immunity to human malignant cells. A brief review and further studies with two gynecologic tumors.Am J Obstet Gynecol114979989

66. Eva A, Robbins KC, Andersen PR, Srinivasan A, Tronick SR, et al. (1982) Cellular genes analogous to retroviral onc genes are transcribed in human tumour cells. Nature 295: 116–119.A. EvaKC RobbinsPR AndersenA. SrinivasanSR Tronick1982Cellular genes analogous to retroviral onc genes are transcribed in human tumour cells.Nature295116119

67. Behrens BC, Hamilton TC, Masuda H, Grotzinger KR, Whang-Peng J, et al. (1987) Characterization of a cis-diamminedichloroplatinum(II)-resistant human ovarian cancer cell line and its use in evaluation of platinum analogues. Cancer Res 47: 414–418.BC BehrensTC HamiltonH. MasudaKR GrotzingerJ. Whang-Peng1987Characterization of a cis-diamminedichloroplatinum(II)-resistant human ovarian cancer cell line and its use in evaluation of platinum analogues.Cancer Res47414418

68. Lau DH, Lewis AD, Ehsan MN, Sikic BI (1991) Multifactorial mechanisms associated with broad cross-resistance of ovarian carcinoma cells selected by cyanomorpholino doxorubicin. Cancer Res 51: 5181–5187.DH LauAD LewisMN EhsanBI Sikic1991Multifactorial mechanisms associated with broad cross-resistance of ovarian carcinoma cells selected by cyanomorpholino doxorubicin.Cancer Res5151815187

69. Benard J, Da Silva J, De Blois MC, Boyer P, Duvillard P, et al. (1985) Characterization of a human ovarian adenocarcinoma line, IGROV1, in tissue culture and in nude mice. Cancer Res 45: 4970–4979.J. BenardJ. Da SilvaMC De BloisP. BoyerP. Duvillard1985Characterization of a human ovarian adenocarcinoma line, IGROV1, in tissue culture and in nude mice.Cancer Res4549704979

70. Buick RN, Pullano R, Trent JM (1985) Comparative properties of five human ovarian adenocarcinoma cell lines. Cancer Res 45: 3668–3676.RN BuickR. PullanoJM Trent1985Comparative properties of five human ovarian adenocarcinoma cell lines.Cancer Res4536683676

71. Hamilton TC, Young RC, McKoy WM, Grotzinger KR, Green JA, et al. (1983) Characterization of a human ovarian carcinoma cell line (NIH:OVCAR-3) with androgen and estrogen receptors. Cancer Res 43: 5379–5389.TC HamiltonRC YoungWM McKoyKR GrotzingerJA Green1983Characterization of a human ovarian carcinoma cell line (NIH:OVCAR-3) with androgen and estrogen receptors.Cancer Res4353795389

72. Schilder RJ, Hall L, Monks A, Handel LM, Fornace AJ Jr, et al. (1990) Metallothionein gene expression and resistance to cisplatin in human ovarian cancer. Int J Cancer 45: 416–422.RJ SchilderL. HallA. MonksLM HandelAJ Fornace Jr1990Metallothionein gene expression and resistance to cisplatin in human ovarian cancer.Int J Cancer45416422

73. Provencher DM, Lounis H, Champoux L, Tetrault M, Manderson EN, et al. (2000) Characterization of four novel epithelial ovarian cancer cell lines. In Vitro Cell Dev Biol Anim 36: 357–361.DM ProvencherH. LounisL. ChampouxM. TetraultEN Manderson2000Characterization of four novel epithelial ovarian cancer cell lines.In Vitro Cell Dev Biol Anim36357361

74. Wolf CR, Hayward IP, Lawrie SS, Buckton K, McIntyre MA, et al. (1987) Cellular heterogeneity and drug resistance in two ovarian adenocarcinoma cell lines derived from a single patient. Int J Cancer 39: 695–702.CR WolfIP HaywardSS LawrieK. BucktonMA McIntyre1987Cellular heterogeneity and drug resistance in two ovarian adenocarcinoma cell lines derived from a single patient.Int J Cancer39695702

75. Fogh J, Trempe G (1975) New human tumor cell lines. In: Fogh J, editor. Human Tumor Cells in Vitro. New York: Plenum Publishing Corp. pp. 115–159.J. FoghG. Trempe1975New human tumor cell lines.J. FoghHuman Tumor Cells in VitroNew YorkPlenum Publishing Corp115159