| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology and Genetics |
Ontario Cancer Institute [J. B., B. B., D. N., J. K., B. Z., J. A. S.] and Departments of Medical Oncology and Hematology [J. D. B., B. Z.] and Gynecological Oncology [B. R., J. M., S. L.], Princess Margaret Hospital, University Health Network; Departments of Medical Biophysics [B. Z., J. A. S.] and Laboratory Medicine and Pathobiology [J. B., B. B., J. A. S.], University of Toronto; and Microarray Centre [P. F. M., M. A.], Clinical Genomics Center, University Health Network, Toronto M5G 2M9, Canada
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Cytogenetic analyses of ovarian carcinomas show both simple numerical and structural changes and complex aberrant changes (4, 5, 6, 7, 8, 9, 10, 11) . A cytogenetic survey of 244 ovarian carcinomas identified a clustering of chromosomal translocation breakpoints occurring in the regions 1p1*, 1q1*, 1p2*, 1q2*, 1p3*, 1q3, 3p1*, 1q4*, 6q1*, 6p2, 6q2, 7p1, 7p2*, 11p1*, 11q2*, 12p1, 12q2*, 13p1, and 19q1 [where the asterisk (*) denotes the most commonly involved regions of chromosomal rearrangement] (12, 13, 14) . The presence of translocations in the regions shown above was associated with reduced patient survival. Chromosomal alterations in the regions 1p1 and 3p1 were found to confer an independent deleterious effect. One disadvantage of this study was that individual breakpoints were assigned to very large chromosomal regions (containing several hundred genes), making specific recurrent structural aberration identification difficult.
In recent years, molecular cytogenetic studies such as CGH6 (15) and SKY (16) have demonstrated their power in identifying recurrent chromosomal aberrations. To date, CGH has been widely used (17, 18, 19, 20, 21, 22, 23, 24, 25) in studying ovarian carcinomas and has identified increased copy numbers at 1q32, 3q26, 8q24.1-q24.2, 20q13.2-qter and frequent chromosomal losses identified at 5q21, 9q, 17p, 17q12-q21, 4q26-q31, 16q, and 22q. Sites of amplification have been identified at 7q36, 8q24.1-q24.2, 3q26.3, 17q25, 19q13.1-q13.2, and 20q13.2-qter by CGH.
Microarray analysis is a relatively new technique (26 , 27) , which allows the simultaneous expression analysis of thousands of genes from a single sample. Profiling and clustering of data from solid tumors have suggested new molecular classifications in lymphoma and breast cancer and generated hypotheses for metastatic markers in melanoma. Similar approaches have identified deregulated genes in ovarian cancer (13 , 28, 29, 30, 31, 32) ; however, the relationship between these candidates and prognosis or classification is not yet clear. Microarray studies on nonovarian malignancies show that cancers with identical histological phenotype have marked differences in global expression patterns. Given this variability, we hypothesized that correlating expression analysis with chromosomal dosage change or rearrangement would allow better identification of key genetic changes in ovarian cancer.
In this study, we have performed parallel analysis using SKY, CGH, and microarray methods on a small cohort of ovarian carcinomas. Our objective was to identify a correlation between deviations in gene expression and chromosomal abnormalities identified by SKY and CGH, in combination with a pilot microarray expression profiling of ovarian carcinoma.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Preparation of Metaphase Spreads
Tumor specimens were finely minced and collagenase treated (Life Technologies, Inc., Gaithersburg, MD) for 24 h at 37°C in a CO2 incubator. The following day, the disaggregated tissue was processed for short-term culture in
-MEM supplemented with FCS (20%; Sigma), penicillin/streptomycin (1%; Life Technologies, Inc.), and L-glutamine (1%; Life Technologies, Inc.). The cultures were harvested within 4 days using 0.1 µg/ml colcemid treatment (Life Technologies, Inc.) for 24 h, hypotonic treatment with 0.075 M KCl, and fixed in 3:1 methanol:acetic acid. The slides were prepared and aged for several days at room temperature.
SKY Method and Nomenclature
SKY analysis was carried out on metaphase slides aged for <1 month. The assay was carried out using the SKY Paint according to the manufacturers instructions (ASI, Carlsbad, CA; Ref. 36
). The full nomenclature describing the karyotypes is presented as supplementary information,2 and is also summarized in the National Center for Biotechnology Information SKY and CGH Database (2001). These operations represent only the clonal aberrations.7
A gain of a chromosome was described when identified in at least two metaphase spreads, a loss when identified in three or more cells, and a chromosomal rearrangement when identified in two or more cells. To assess the karyotype complexity, the number of cytogenetic abnormalities detected by SKY was determined by counting the aberrations described in the karyotype in each case according to the following criteria: Multiple copies of structural or numerical aberration were counted only once; in derivative chromosomes, each structural abnormality was counted once; and each breakpoint in complex rearrangements was counted once. In cases 13A/B, 15A/B, and 21A/B, where there are paired primary and metastatic samples, the stemline karyotype was counted once, with only the addition of the unique clonal changes in each of the primary and metastatic sample added to Fig. 2
. Thus, the relative comparisons of overall structural aberrations and imbalances shown in Fig. 2
and summarized in Table 2
refer to the total number of SKY breakpoints and CGH copy number changes detected along the length of a given chromosome. In determining the mean number of breakpoints present between endometrioid and serous subtypes, the full karyotypes (i.e., stemline and unique breakpoints) were averaged giving a total of 443 breakpoints for tumors OCA2, OCA3A, OCA8, OCA13A/B, OCA15A/B, and OCA21A/B.
|
|
Microarray Analysis
Microarrays used in this study were purchased from the Microarray Centre, University Health Network (Toronto, Ontario, Canada). These arrays contain 1718 human cDNAs selected from the Swissprot database, with 74 positive control and 74 negative control features. All features are spotted in duplicate for a total of 3840 spots. A complete list of the cDNA collection used for these arrays ("Human 1.7K2" and "Human 1.7K3") and protocols used for array construction can be found at the University Health Network Microarray Centre web site.8
Human 1.7K2 and Human 1.7K3 contain the same set of human ESTs, and the difference is the cDNA spot locations for some ESTs. All reverse transcriptions (direct and indirect labeling) and hybridizations conditions are available.9
Two series of microarray experiments were carried out as described below.
Microarray Analysis with a Nine Cell Line Reference RNA.
For the microarray analysis of 13 tumor samples against RNA derived from a pool of nine cell lines, the RNA samples were labeled using amino-allyl indirect labeling, using a 2:1 ratio of amino-allyl dUTP:dTTP. Tumor RNAs were labeled with Cy5, and the pooled reference RNA was labeled with Cy3, followed by column purification (High Pure PCR Product Purification kit; Roche). The eluates were lyophilized and stored at -70°C. Duplicate arrays were hybridized for each tumor sample. For a subset of 4 tumor samples, there was sufficient RNA for comparisons to two normal ovarian reference RNAs. In these experiments, all RNAs were labeled using direct labeling. In one set of arrays, all 4 tumor RNAs and the Stratagene normal ovary RNA were labeled with Cy3, and the reference normal ovary RNA (Ambion) was labeled with Cy5.
Microarray Analysis Using a Normal Ovarian RNA Reference.
In a second set of arrays, all tumor RNAs and the Stratagene normal ovary RNA were labeled with Cy3, and the reference RNA (Ambion) was labeled with Cy5 (dye switches) to account for possible labeling bias. Slides were scanned using either an Axon GenePix 4000A (Axon, Foster City, CA) or ScanArray 4000 scanner (Packard BioScience) dual laser scanner and quantified with GenePix Pro 3.0 software (Axon). Features that were not quantified because they were flagged bad or absent by GenePix Pro were manually reviewed for accuracy. Filtering, normalization of the raw data, and complete data analysis of the data sets were carried out with an algorithm developed in house that will be published elsewhere.10
Briefly, each of the 16 subgrids on each array was independently normalized by equalizing the Cy3 intensities, with respect to the Cy5 intensities, while excluding features flagged bad or absent by GenePix Pro. Other features excluded by the algorithm included saturated spots and spots with foreground:background intensity ratio < 2. The Cy5:Cy3 normalized intensity ratio was determined for each spot, and these values from the duplicate spots within each array were averaged. Subsequently, Cy5:Cy3 ratios for the same spots between replicate experiments of each sample were averaged together as a single project file. Finally, the project files representing each sample and its replicates were combined into a single text file for processing by the Cluster software package for hierarchical clustering, and the Treeview software was used for generating the graphical visualization of the clustering (39)
.11
Clustering was performed on genes showing expression values present in
80% of samples. For the first series of microarray experiments, we used median centering of the genes and tumors to emphasize differences between tumors rather than changes with respect to the standard comparator RNA (39
, 40)
. The agglomerative hierarchical clustering algorithm used in the Cluster software successively joins gene expression profiles to form a dendrogram based on their pairwise similarities. The same procedure is followed when clustering by experiment, i.e., the similarity between each experiment is determined over the total number of genes as an average, and experiments with similar averages are grouped together (41)
. Two-dimensional hierarchical clustering first reorders genes and then reorders tumors based on similarities of gene expression between samples. The datasets used for hierarchical clustering generated by this study are available as tab-delimited format text files.12
Custom software was also developed to automate retrieval of chromosome localizations of microarray cDNAs from the UniGene database (at present, build 144). This allowed for arrangement of cDNAs into sequential order in megabases along each chromosome. The mapping distribution of the 1718 human cDNAs (1.7K2 and 1.7K3) used in this study, as well as the software, to carry out retrieval of chromosomal localizations are available.13
Confirmation of Microarray Results by Semi-quantitative RT-PCR.
cDNA was prepared from 20 ng of total RNA extracted from samples OCA17, OCA19, OCA21A, and OCA21B by reverse transcriptase-PCR using specific primers for the following genes using the Sigma One-Step RT-PCR kit (Sigma-Aldrich, St. Louis, MO): OPN (sense: 5'ACAGCATCGTCGGGACCAGA3', antisense: 5'CATCCATCCACCTTCATCCACC3'); decorin (sense: 5'CAGTGTTCTGATTTGGGTCTGGAC3', antisense: 5'CGTAAGGGAAGGAGGAAGACC3'); and ribosomal 60S (sense: CCGCCATAACAAGGACCGAAAGT3', antisense: 5'ACCCGGCCGTCTTGTTTTCC3'). RT-PCR reactions were also run in parallel with PCR primer pairs for the 18S rRNA (Ambion) as an internal standard. The following reaction/cycling conditions were used: 45 min at 45°C; 3 min at 94°C; 30 cycles of 45 s at 94°C, 30 s at 55°C, and 1 min at 72°C; and 5 min at 72°C. An aliquot of each reaction was run onto a 1% agarose gel, and the gel was stained with ethidium bromide. The intensities of each band were quantified using a Gel Documentation AlphaImager System (Alpha Innotech, San Leandro, CA) and were normalized using the intensities obtained for the 18S internal standard. The results are displayed in Fig. A.2
, 14
| RESULTS |
|---|
|
|
|---|
|
40 clonal aberrations detected by SKY in the right and left ovary were concordant; however, the right ovary (OCA15A) possessed additional aberrations, including deletions at 1q22, 3p21, and 8p10 and complex rearrangements affecting chromosomes 2, 5, and 11 that were not detected in the left ovary (OCA15B). Similarly, the left ovary exhibited unique aberrations, including deletions at 2p11 and 11q11 and a translocation der(11)
t(11;15)(p11;q12) not evident in the tumor present at the right ovary. OCA21A and OCA21B were primary and metastatic tumors, respectively. OCA21A was found to be tetraploid, whereas OCA21B was found to be diploid, but the majority of structural and numerical aberrations were concordant and duplicated in the primary tumor (Fig. 1, B and C)
q24::17q22. However, several unique clonal rearrangements were detected in the metastatic tumor OCA21B such as the three translocations +der (1)
t(1;18)(q41;q22), +der(8)
t(7;8)(p12;p12?), der(13)
t(9;13)(q21p12), and der (13)
t(9;13)(p11;p11). Similarly for the tumor pair OCA13A and OCA13B, the majority of the
68 aberrations found were common to both samples. Aberrations that were unique to the metastatic tumor OCA13B included del(X)(q11) and two translocations der(4)
t(4;5)(q24;?p11), and der(18)
t(18;20)(q11.2;?). In summary, it was possible to demonstrate by SKY that three of the four paired samples (OCA13A/B, OCA15A/B, and OCA21A/B) were closely related in terms of their pattern of cytogenetic aberration determined by SKY, and that this concordance indicates an origin from a common progenitor.
Chromosomal Aberrations Identified by both CGH and SKY.
In Fig. 1A
, a representative comparison between the CGH findings and SKY results using tumor OCA13A is shown. It can be seen that extensive areas of CGH gain on chromosomes 1, 2, 3, and X arise as a result of large unbalanced translocations detected by SKY. Smaller regions of CGH loss and gain are also detected in other chromosomal regions that appear to be generated by much smaller rearrangements. In cases where an unbalanced chromosomal rearrangement was determined by SKY, but the specific sub region could not be recognized because of a lack of identifiable banding information, CGH was used to infer the probable region of genomic imbalance. These data are illustrated for OCA13A in Fig. 1A
, where a net gain was detected at 8q21-q22. SKY analysis detected aberrations on chromosome 8: der(8)
t(8;10)(p11;p11); and chromosome 19: der(19)
t(8;19)(?;p11),der (19)
t(8;19) (?p11;q13),+der(19)
(9pter
p13::19p13
q13::8?p11
p?). On the basis of the CGH findings, it is likely that the deleted region of 8p arises as a result of a complex translocation involving 19 and 13 (inverted DAPI banding suggests that this segment derives from 8p). The remainder of chromosome 8 was involved in at least 4 breakpoints on 8 from the 8q21 region because this was detected as a regional gain by CGH. Imbalances identified by CGH analysis (Fig. 2)
detected frequent chromosomal gains at 3q (58.3%), 8q (50%), 7q (50%), 20 (41.6%), 1q (33.3%), and 6p (25%). The most frequent losses were detected in chromosomes 4 (58.3%), 3p (41.6%), 6q (33.3%), and 18q (33.3%). The minimal common regions of chromosomal gains were identified at 1q21-q44, 3q13.3-q29, 6p21.2-p23, 7q32-q36, 8q23-q24.3, and 20p, whereas the minimal region of common loss was identified as 3p14-p26, 4q, 6q21-q27, and 18q12-q22.
To determine an integrated overall chromosomal distribution of genomic rearrangements as identified by both SKY and CGH, the total number of aberrations per chromosome shown in Fig. 2
was enumerated, and 396 rearrangements were identified in the study group. When the distribution was normalized with respect to relative genomic length, chromosomes 3, 8, 11, 17, and 21 had the highest frequencies of structural and numerical aberrations identified by SKY and CGH. Analysis of these combined CGH and SKY data demonstrates that regions of chromosomal imbalance (as detected by CGH) were concordant with the locations of chromosomal breakpoints (as identified by SKY), particularly for chromosomes 3 and 8. Consistent losses distal to 3p21 and gains of 3q (detailed in Fig. 1D
) in tumors OCA2, OCA3A, OCA15A/B, and OCA19 were observed, whereas tumors OCA3, OCA5, OCA19, and OCA21A/B all had both SKY and CGH aberrations of the 8q23-q24 regions (Fig. 1E)
.
Microarray Expression Analysis of Study Group.
To determine whether patterns of gene expression could be associated with the chromosomal features associated with the study group and to explore expression profiling of the tumor samples, two series of microarray experiments were carried out. At the time we began this study, normal ovary RNA was not available to us (commercially or otherwise). Therefore, in a first series of experiments, 13 ovarian tumors were analyzed by microarray using a pool of RNA derived from nine cell lines as an internal standard. When used in conjunction with two-dimensional average linkage unsupervised clustering, this approach allows the determination of relative gene expression profiles in multiple samples and is commonly used in microarray studies (40
, 43, 44, 45)
. The results of the two-dimensional hierarchical clustering for this first series of experiments is presented in Fig. 3A
. Profiles from primary tumor and metastasis isolated from the same patient, irrespective of biopsy time and anatomical site, showed the greatest similarity and always clustered close together as shown previously (27
, 40
, 46, 47, 48, 49, 50)
, indicating the robustness of the clustering algorithm. Consistent grouping of poorly differentiated tumors within this tumor set was seen independent of histological subtypes. Surprisingly, OCA15A/B, also a poorly differentiated tumor, clustered with the well-differentiated tumors. In addition, the three tumor pairs of karyotypically related tumors (OCA13A/B, OCA15A/B, and OCA21A/B) exhibited similar profiles and each tumor pair clustered together.
|
|
Among the genes that showed up-regulation in the tumor samples (cluster I) was the gene encoding the OPN precursor, which showed a 53-fold increase in expression in sample OCA17 and a 2-fold increase in the other tumor samples. Five genes encoding metallothionein showed up-regulation (average ratios between tumor samples between 7.01 for MT-II and 2.87 for MT-1H), as well as gene encoding for the VEGF receptor 2 precursor (average ratio, 3.17), caspase-8 precursor (average ratio, 3.13), MYB-related protein (average ratio, 2.57), and six HLA class II histocompatibility antigens (average ratios between 3.66 and 2.23). Among the genes that showed down-regulation in the tumor samples (cluster II), there were 19 transcripts encoding 60S and 40S ribosomal proteins (average ratios between 0.49 and 0.23), IGFBP-3 (average ratio, 0.34), IGFBP-4 precursor (average ratio of 0.28), tumor necrosis factor
-induced protein 1 (average ratio, 0.26), G2/mitotic specific cyclin 2 (average ratio, 0.25), and decorin (average ratio, 0.15). The expression level of a subset of the differentially expressed genes was analyzed by semi-quantitative RT-PCR. For OPN, decorin, and the ribosomal 60S L19 gene, the levels of overexpression or underexpression relative to a normal ovary were consistent with the levels determined by microarray expression analysis.2
Microarray Expression Analysis of cDNAs Mapping to Chromosomes 3 and 8.
Because both the literature and our parallel SKY and CGH analysis have identified elevated frequencies of rearrangement and imbalance on chromosomes 3 and 8, we performed a detailed analysis of the results of ovarian cancer expression levels relative to a normal ovary for serous-derived tumors OCA17, OCA19, and OCA21A. A total of 94 and 51 targets represented genes with specific Blast and UniGene chromosomal band assignments on chromosomes 3 and 8, respectively, on the 1.7K3 cDNA array (Table 3)
. There was a contiguous distribution of reduced expression of seven genes mapping to 3p25.53p21.3 interval and increased expression of four genes from 3q13.333q28. These expression data were consistent with the observed pattern of CGH loss in distal 3p and gain in 3q observed both in this study and in the literature. Similarly, CGH data have consistently showed frequent gain and amplification in the 8q2324 region in ovarian cancer, and the Exostosin Type I gene mapping to 8q23 exhibited increased expression in all three tumors.
|
| DISCUSSION |
|---|
|
|
|---|
Previous conventional cytogenetic analysis has identified frequent breakage at 1p1, 1q1, 1p2, 1q2, 1p3, 1q3, 3p1, 1q4, 6q1, 6p2, 6q2, 7p1, 7p2, 11p1, 11q2, 12p1, 12q2, 13p1, and 19q (12) . Similarly, loss of heterozygosity studies by Launonen et al. (51) have identified an association between loss of sequences at 3p14.2, 11p15.5, 11q23.2-q24, 16q24.3, and 17p13.1 and adverse outcome. The resulting loss of heterozygosity in these and other regions can result from simple deletions or deletions as a result of translocations. Because SKY can identify cryptic translocations and poorly resolved "marker chromosomes" with greater certainty, a much more accurate picture of overall cytogenetic change is possible. Many chromosomal translocations were centromeric or pericentromeric, suggesting the involvement of repetitive centromeric sequences or Alu repeats that may facilitate rearrangements leading to chromosomal imbalance (52) . Chromosomal breakage at 1p1* and 3q* occur in ovarian cancer (5 , 7 , 11) and confers a poor outcome (13) .
Ovarian cancer often presents in more than one site or has metastasized at the time of surgery. We studied the chromosomal and gene expression patterns in tumors arising in two sites to determine their common origin. In three instances, two separate tumors were available for simultaneous analysis of patterns of gene expression and karyotypes. In all three pairs, there was a marked similarity in the chromosomal constitution of the tumors and each tumor pair clustered together on the dendrogram obtained with microarray analysis. One tumor pair presented an interesting opportunity to determine cytogenetic changes during metastasis. OCA21 possessed a tetraploid primary tumor and a diploid karyotype in the metastatic tumor, suggesting that either a reduction in chromosomal ploidy occurred in the primary tumor before metastasis or that the reduction occurred at the metastatic site. The cytogenetic description of the metastatic biopsy is similar to components of structural alterations observed in the primary tumor, suggesting that a selection of a dominant clone with specific structural abnormalities arose in the metastatic tumor.
The minimal regions of chromosomal imbalance identified by CGH were gains at 1q21-q44, 3q13.3-q29, 6p21.2-p23, 7q32-q36, and 8q23-q24.3 and losses at 3p14-p26, 4q, 6q21-q27, and 18q12-q22. Previous CGH studies (17, 18, 19, 20, 21, 22, 23, 24) have identified a similar pattern of genomic change, and attempts have been made to correlate these data with clinical endpoints (53) . Studies by Yonescu et al. (54) and Jenkins et al. (11) show that structural and numerical changes predominantly occur on chromosome 11 endometrioid tumors. The high percentage of chromosome 4 loss (58.3%) is consistent with literature reports of chromosome 4 involvement in advanced staged tumors (18) .
The expression clustering of genes/pathways may elucidate those critical genes that promote both a more aggressive disease course and worse outcome and may provide a marker for treatment resistance. The first series of microarray experiments described here consisted of analyzing 13 tumor samples using a pool of nine cell lines as a common reference. Two-dimensional hierarchical clustering identified a grouping of predominantly poorly differentiated subtypes (OCA13A/B, OCA19, and OCA21A/B), which were independent of histological grouping. The only discrepancy was the tumor OCA15 that failed to cluster with the other poorly differentiated tumors in the study group. It is possible that histological heterogeneity or the lower quantity of RNA obtained for this particular tumor sample resulted in this grouping obtained by cluster analysis.
The second series of microarray experiments carried out on a subset of tumor samples and using normal ovarian tissue as a common reference identified 194 genes differentially expressed at a level of at least 2 (up or down) between tumor samples and the reference normal sample. Two-dimensional hierarchical clustering and statistical analysis revealed a very similar expression profile for all 4 tumor samples, with only 13% of genes showing a SD > 1 across all tumor samples, consistent with their similar serous histology (OCA17, OCA19, and OCA21A/B). As expected, the normal ovary sample used as a negative control displayed very few significantly differentially expressed genes when compared with the normal reference (Fig. 4)
.
One of the objectives of this study was to identify correlations between gene expression deviations observed by microarray analysis and chromosomal imbalances identified by CGH. Because chromosomes 3 and 8 have been previously shown to be subject to rearrangement and copy number change (12, 13, 14) and were found to be commonly involved in rearrangements by SKY and/or CGH, we performed a more detailed expression evaluation of cDNA mapping to these chromosomes. We identified a contiguous distribution of reduced expression of seven genes mapping to 3p25.53p21.3 interval and increased expression of four genes from 3q13.333q28 that were consistent with loss of 3p and gain of 3q in this and other CGH studies of ovarian cancer (17, 18, 19, 20, 21, 22, 23, 24, 25) . Within the 3q13.333q28 interval, it is noteworthy that there was elevated expression of the calcium-sensing receptor (55) that has been previously shown to induce human ovarian surface epithelial cells to modulate extracellular calcium and induce proliferation.
Two-dimensional hierarchical clustering revealed two large clusters of expression profiles (I and II, Fig. 3B
). The first cluster consists of up-regulated genes and most striking is the very high value (ratio of 53) obtained for the OPN precursor in sample OCA17. This gene is also overexpressed in the other tumor samples but to a much lesser extent. OPN is a noncollagenous bone-related protein that has been detected by immunohistochemical and in situ hybridization in calcified psammona bodies, which frequently occur in ovarian serous papillary tumors (56)
. This protein may be causally related with the calcification of the psammona bodies of the ovarian serous papillary tumors such as OCA17. Five genes encoding metallothionein are also up-regulated in all 4 tumor samples. Metallothionein is a key component of platinum resistance in epithelial ovarian cancer and is thought to be of prognostic significance (57)
.
Cluster I also contains five genes encoding HLA class II histocompatibility antigens. All five genes are located in the 6p region, which was involved in genomic imbalance and structural rearrangement, so it was of interest that MHC class II genes of the leukocyte antigen (HLA) complex emerges as a strong cluster of expression. This observation is in agreement with a recent high-density filter array analysis of ovarian cancer, which also reported differential expression of immune response mediators (58) . Previously, an association between the T-cell response and alterations to MHC expression has been observed in ovarian cancer (59 , 60) . However, it is conceivable that infiltrating lymphocyte may also contribute to overexpression of immune response genes.
The VEGF receptor 2 precursor is up-regulated in all tumor samples, and this is consistent with a previous cDNA microarray study (58) on four poorly differentiated serous papillary tumors that showed an overall increase in angiogenesis-related markers, including VEGF.
Among the down-regulated genes, we identified 19 genes encoding ribosomal proteins. Ribosomal proteins have been found to be up-regulated in some tumor tissues (61) . However, a recent study (62) using a data mining tool called Digital Differential Display has shown that distinct ribosomal proteins were found to be up-regulated and down-regulated in a tumor type-dependent manner. In that study, some ribosomal proteins (S15A, S17, and L35) were found to be up-regulated in breast and prostate carcinoma-derived libraries, whereas others such as L9, L23, and L37A were down-regulated in ovarian cancer-derived libraries. Our present findings are consistent with these results. IGFBP-3 precursor showed down-regulation in all tumor samples as observed by others (63) . Finally, the gene encoding a decorin precursor showed the highest level of down-regulation. Decorin is a small proteoglycan protein, which is part of the cellular matrix and has been shown to inhibit the growth of two ovarian cell lines in vitro (64) . Decorin inhibits transforming growth factor ß and induces p21, resulting in inhibition of proliferation (65) . A marked decrease in the endogenous level of decorin precursor might therefore result in an increase in cell proliferation in tumor cells.
A more complete cytogenetic and microarray analysis of a larger collection of ovarian tumor and normal samples and using cDNA microarrays spotted with 19,200 genes and ESTs is under way in our laboratories that will allow us to uncover additional differentially expressed genes and to further explore the ones identified in the present study,16 as well as additional regions of chromosomal instability associated with ovarian cancer. In this study, we performed a comparative analysis using SKY, CGH, and microarrays on a small cohort of ovarian carcinomas. We showed that chromosomes 3, 8, 11, 17, and 21 had the highest overall frequencies of structural and numerical aberrations and losses distal to 3p21, and consistent gains of 3q were associated with a pattern of reduced expression of genes mapping to 3p25.53p21.3 and increased expression of genes from 3q13.333q28. Thus, although there is inherently limited resolution of all cytogenetically based techniques, we have shown that a parallel strategy of genomic and expression analysis can facilitate the identification of smaller subsets of genes pertinent to ovarian cancer.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Supported in part by the Ovarian Fashion Show Committee. ![]()
2 Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org). ![]()
3 Present address: Department of Oncology (Box 193), School of Clinical Medicine, University of Cambridge, Addenbrookes Hospital, Cambridge CB2 2QQ, United Kingdom. ![]()
4 Present address: The Cross Cancer Institute, 11560 University Avenue, Edmonton, Alberta, T6G 1Z2 Canada. ![]()
5 To whom requests for reprints should be addressed, at Ontario Cancer Institute, Division of Cellular and Molecular Biology, 610 University Avenue, Room 9-721, Toronto, Ontario, M5G 2M9 Canada. Phone: (416) 946-4509; Fax: (416) 920-5413; E-mail: jeremy.squire{at}utoronto.ca ![]()
6 The abbreviations used are: CGH, comparative genomic hybridization; SKY, spectral karyotyping; EST, expressed sequence tag; OPN, osteopontin; DAPI, 4',6-diamidino-2-phenylindole; RT-PCR, reverse transcription-PCR; IGFBP, insulin-like growth factor binding protein; VEGF, vascular endothelial growth factor. ![]()
7 Internet address: http://www.ncbi.nlm.nih.gov/sky/skyweb.cgi. ![]()
8 Internet address: http://www.uhnres.utoronto.ca/services/microarray. ![]()
9 Internet address: http://www.uhnres.utoronto.ca/services/microarray/protocols/index.html. ![]()
10 B. Beheshti, I. Braude, P. Marrano, P. T. Thorner, M. Zielenska, and J. A. Squire, manuscript in preparation. ![]()
11 Internet address: http://rana.lbl.gov. ![]()
12 Internet address: http://www.utoronto.ca/cancyto/OVCA2001CR/. ![]()
13 Internet address: http://www.utoronto.ca/cancyto/. ![]()
14 Fig. A to be available as supplementary information on the web edition. RT-PCR analysis of OPN, decorin, and ribosomal 60S using ovarian cancer samples OCA17, OCA19, and OCA21A/B. ![]()
15 Table A to be available as supplementary on the web edition. Full description of SKY and CGH results. ![]()
16 D. Grisaru, M. Albert, J. Goncalves, J. Woodgett, J. A. Squire, K. J. Murphy, B. Rosen, A. Covens, R. Osborne, H. Begley, P. Shaw, and P. F. Macgregor. Prognostic prediction in advanced serous epithelial ovarian cancer using molecular profiling, manuscript in preparation. ![]()
Received 6/14/01. Accepted 4/30/02.
| REFERENCES |
|---|
|
|
|---|
and TNF-
. Immunol. Lett., 64: 153-160, 1998.[Medline]
This article has been cited by other articles:
![]() |
A. D. PANANI, C. ARAVIDIS, Z. KOSMAIDOU, A. RODOLAKIS, and A. ANTSAKLIS Simple Structural Chromosomal Abnormalities in Advanced Stage of Ovarian Cancer In Vivo, May 1, 2009; 23(3): 425 - 428. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. K. Zorn, T. Bonome, L. Gangi, G. V.R. Chandramouli, C. S. Awtrey, G. J. Gardner, J. C. Barrett, J. Boyd, and M. J. Birrer Gene Expression Profiles of Serous, Endometrioid, and Clear Cell Subtypes of Ovarian and Endometrial Cancer Clin. Cancer Res., September 15, 2005; 11(18): 6422 - 6430. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. M. Levin, D. Ghosh, K. R. Cho, and S. L. R. Kardia A model-based scan statistic for identifying extreme chromosomal regions of gene expression in human tumors Bioinformatics, June 15, 2005; 21(12): 2867 - 2874. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Hibbs, K. M. Skubitz, S. E. Pambuccian, R. C. Casey, K. M. Burleson, T. R. Oegema Jr, J. J. Thiele, S. M. Grindle, R. L. Bliss, and A. P.N. Skubitz Differential Gene Expression in Ovarian Carcinoma: Identification of Potential Biomarkers Am. J. Pathol., August 1, 2004; 165(2): 397 - 414. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. A. Heinzelmann-Schwarz, M. Gardiner-Garden, S. M. Henshall, J. Scurry, R. A. Scolyer, M. J. Davies, M. Heinzelmann, L. H. Kalish, A. Bali, J. G. Kench, et al. Overexpression of the Cell Adhesion Molecules DDR1, Claudin 3, and Ep-CAM in Metaplastic Ovarian Epithelium and Ovarian Cancer Clin. Cancer Res., July 1, 2004; 10(13): 4427 - 4436. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Varma, T. Rollason, J. K Gupta, and E. R Maher Endometriosis and the neoplastic process Reproduction, March 1, 2004; 127(3): 293 - 304. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. D. Pack, O. M. Alper, K. Stromberg, M. Augustus, M. Ozdemirli, A. M. Miermont, G. Klus, M. Rusin, R. Slack, N. F. Hacker, et al. Simultaneous Suppression of Epidermal Growth Factor Receptor and c-erbB-2 Reverses Aneuploidy and Malignant Phenotype of a Human Ovarian Carcinoma Cell Line Cancer Res., February 1, 2004; 64(3): 789 - 794. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. K. Gruvberger-Saal, P. Eden, M. Ringner, B. Baldetorp, G. Chebil, A. Borg, M. Ferno, C. Peterson, and P. S. Meltzer Predicting continuous values of prognostic markers in breast cancer from microarray gene expression profiles Mol. Cancer Ther., February 1, 2004; 3(2): 161 - 168. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Dobbin, J. H. Shih, and R. Simon Questions and Answers on Design of Dual-Label Microarrays for Identifying Differentially Expressed Genes J Natl Cancer Inst, September 17, 2003; 95(18): 1362 - 1369. [Full Text] [PDF] |
||||
![]() |
L. A. Loeb, K. R. Loeb, and J. P. Anderson Multiple mutations and cancer PNAS, February 4, 2003; 100(3): 776 - 781. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. F. Macgregor and J. A. Squire Application of Microarrays to the Analysis of Gene Expression in Cancer Clin. Chem., August 1, 2002; 48(8): 1170 - 1177. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |