Histotype-specific copy-number alterations in ovarian cancer
BMC Medical Genomics volume 5, Article number: 47 (2012)
Epithelial ovarian cancer is characterized by multiple genomic alterations; most are passenger alterations which do not confer tumor growth. Like many cancers, it is a heterogeneous disease and can be broadly categorized into 4 main histotypes of clear cell, endometrioid, mucinous, and serous. To date, histotype-specific copy number alterations have been difficult to elucidate. The difficulty lies in having sufficient sample size in each histotype for statistical analyses.
To dissect the heterogeneity of ovarian cancer and identify histotype-specific alterations, we used an in silico hypothesis-driven approach on multiple datasets of epithelial ovarian cancer.
In concordance with previous studies on global copy number alterations landscape, the study showed similar alterations. However, when the landscape was de-convoluted into histotypes, distinct alterations were observed. We report here significant histotype-specific copy number alterations in ovarian cancer and showed that there is genomic diversity amongst the histotypes. 76 cancer genes were found to be significantly altered with several as potential copy number drivers, including ERBB2 in mucinous, and TPM3 in endometrioid histotypes. ERBB2 was found to have preferential alterations, where it was amplified in mucinous (28.6%) but deleted in serous tumors (15.1%). Validation of ERBB2 expression showed significant correlation with microarray data (p=0.007). There also appeared to be reciprocal relationship between KRAS mutation and copy number alterations. In mucinous tumors where KRAS mutation is common, the gene was not significantly altered. However, KRAS was significantly amplified in serous tumors where mutations are rare in high grade tumors.
The study demonstrates that the copy number landscape is specific to the histotypes and identification of these alterations can pave the way for targeted drug therapy specific to the histotypes.
Ovarian cancer is often dubbed a ‘silent’ killer because of its non-specific symptoms and late clinical onset which contribute to overall poor prognosis. There has been a steady increase in incidence over the last three decades with 204,000 new cases diagnosed each year globally . It ranked fifth in mortality among cancers in women and has the highest case-fatality rate in gynecological cancers . The 5-year survival rate for women with advanced disease remains at 29% with estimated 125,000 deaths annually [1, 3].
About 90% ovarian cancers are epithelial ovarian cancer (EOC) , histologically subtyped as serous, mucinous, endometrioid, or clear cell. Further subtyping include the borderline cases, such as mucinous or serous borderline, often presented as stable diseases with more favorable outcome compared to the non-borderline subtypes. It is now recognized that epithelial ovarian cancer is a spectrum of diseases with varied genetic mutations among histotypes . Genetically, mutations differ between the grades of the disease. Low-grade serous carcinoma have high frequency of KRAS and BRAF mutations but few p53 mutations while high grade serous carcinoma shows the inverse in frequency of these mutations [6, 7]. Mucinous histotype has KRAS mutations  and endometrioid has PTEN mutations . Despite the molecular heterogeneity, the treatment standard remains as taxane/platinum-based chemotherapy for all histotypes.
Genomic alterations such as copy number alterations (CNA) have been known to harbor drivers in carcinogenesis. Driver genes are genes that confer growth advantage on the cancer cells . Several known copy number alteration drivers in cancers include receptor tyrosine kinases such as EGFR, FGFR, and ERBB2, which are targets for drugs therapy . Successful incorporation of genomics alterations studies in cancer treatment has been evident in breast, leukemia, and lung cancers, where targeted therapies are part of the standard treatment protocols [12–14]. For example, Trastuzumab, a targeted therapy that can significantly reduce risk of disease recurrence and improve overall survival, is now standard of care for early-stage patients with Her2-positive breast cancer . To date, targeted drug therapy has not been successfully incorporated in EOC.
One of the challenges in elucidating copy number alterations in EOC is the disproportional prevalence of histotype. Serous is the most prevalent (70-85%), followed by endometrioid (5-10%), clear cell (5-10%), and mucinous (~5%). The lower prevalent diseases tend to suffer from small sample size for statistical analyses to identify copy number alterations. To elucidate CNA in each histotype, we combined data from multiple studies with similar platforms to identify copy number alterations that are specific to the histotypes. The motivation is to identify high confidence histotype-specific alterations that may otherwise be obscured due to disproportionate prevalence of the histotypes.
Two major results are presented here; the histotype-specific copy number alterations for EOC and the identification of potential driver genes. Three datasets with corresponding gene expression and copy number profiling on similar platforms were used for this study. A hypothesis-driven approach using stringent false discovery rate (FDR) filtering was used to identify potential copy number driver genes (Methods). The genes were compared with other studies as well as cancer genes reported in literature. In addition, we also validated the expression of ERBB2, a driver gene, via quantitative real-time PCR (qPCR).
Distinct copy number alterations in EOC histotypes
Figure 1a shows the global frequency of copy number aberrations for EOC across the genome for the merged dataset. In concordance with other reports, the four commonly reported chromosomes of 3q, 8, 17p, and 20q showed broad copy number alterations in EOC [15–17]. When this was de-convoluted into histotypes (Figure 1b, shown in 4 tracks for clear cell, endometrioid, mucinous, and serous), distinct differences were observed. It is evident that the general frequency for EOC was mirrored in serous tumors, not surprising since it has the largest sample size. In the four most commonly reported altered regions, 3q, 8, and 20q amplifications were observed in serous, clear cell, endometrioid but not mucinous tumors. Endometrioid and serous tend to harbor more copy number alterations, with more broad regions of alterations involving the p- or q-arm. The genomics landscape for clear cell and mucinous tumors appeared different from the other histotypes, with lesser broad regions of alterations and in lower frequency.
To assess the significance of copy number altered regions, we used a 2-pronged approach using merged and individual datasets (Additional file 1: Figure S1). Figure 2 shows significant copy number altered regions (see Methods) in the histotypes. Broad 3p amplification and 8p deletion were observed in serous tumors; 8q amplifications in clear cell and serous tumors; 17p deletions in mucinous and serous tumors; and chr20 amplifications in serous tumors. The nature of the alterations also differs, e.g. focal versus broad alterations. For example in 3q, serous tumors showed broader amplifications in the region of 3q13.31-29 while focal amplification was observed for clear cell tumors at 3q26.2-26.32. This is interesting as it has been reported that overlapping broad and focal aberrations can have distinct functional consequences . In other chromosomes, alterations were specific to histotypes as well; such as 9p21 focal deletions in mucinous histotype, reportedly harboring homozygous deletions in EOC [19–22]. There are regions that displayed opposite trend in alterations between histotypes. One particular region, 8p23.1, showed amplification in clear cell but deletion in serous tumors. Another region which showed opposite trend in alteration was 17q12 which harbor the oncogene ERBB2; the gene was significantly amplified in mucinous but deleted in serous tumors. Excluding borderline cases, 28.6% (4/14) mucinous samples had amplifications and 15.2% (15/99) of serous tumors had deletions of ERBB2. 1/5 mucinous borderline also showed ERBB2 amplification. It should be of note that although ERBB2 was found significantly deleted in serous tumors, 5/99 (5.1%) of serous samples harbored the amplification, close to the 3% reported by TCGA for high grade serous . No ERBB2 deletion was observed in the mucinous samples. As serous tumors have a comparatively larger sample size than the other histotypes, we would expect more significant regions for this histotype. Nevertheless, using stringent criteria, we were able to identify some significant CNA for the lower prevalent tumors.
To quantify genes that were altered in each histotype, we mapped genes to regions that were identified in each histotype. 6375 unique genes were found to be altered: 2682 amplifications and 3712 deletions (Additional file 2: Table S1). 91% of genes were amplified and 97.1% deleted in serous tumors, 19.1% amplified and 1.5% deleted in clear cell, 14.3% amplified in endometrioid, and 0.5% amplified and 11.5% deleted in mucinous. A total of 5360 genes were specific to each histotype (amplification=2014, deletion=3346), 5011 in serous tumors, 193 in endometrioid tumors, 79 in clear cell tumors, and 77 in mucinous tumors. Within each histotype, the type of alterations varied. Clear cell tumors had more amplified genes than deleted genes, while mucinous and serous tumors had more deleted genes than amplified genes. Only amplified genes were found in endometrioid tumors. Figure 3 shows the Venn diagram of overlapping genes in the lower prevalent histotypes. Due to the differences in sample size of each histotype, comparisons between overlapped genes were limited to genes found only in the non-serous tumors where the sample sizes are more comparable. A small number of overlapped amplified and deleted genes between clear cell and mucinous tumors were observed; none with endometrioid. This suggests that most of the CNA are specific to histotypes, adding to the mounting evidence of molecular differences between histotypes. Despite larger sample size in serous histotype, comparison of CNA found in the lower prevalent histotypes with serous is also of interest (Additional file 3: Figure S2). Clear cell tumors had the highest number of common altered genes with serous tumors (7.6%) while endometrioid tumors had the lowest number of common altered genes; 3.0% common altered genes with serous tumors only. There is also more overlap on amplifications than deletions.
Consistency with other studies
Three studies have reported copy number changes in EOC [16, 17, 24] (Additional file 4: Table S2). These studies focused on global changes in EOC rather than histotype-specific alterations. Nevertheless, histotype-specific alterations from our study should overlap to some degree with these reported genes. We compiled a list of 551 significant genes reported in the 3 studies of which 545 genes were reported by either Haverty et al. or Gorringe et al. Eleven genes (2%) were commonly identified in at least 2 of the studies. In comparison, 39.2% (216/551) of genes were found in our study, indicating that our approach can identify copy number altered genes. We also compared our findings in serous tumors with TCGA high grade serous study (n=489)  where they reported 63 regions of gains and 50 regions of deletion. Based on overlapping of genes (if available) or genomic regions, 29/63 (46%) amplified and 27/50 (54%) deleted regions were also found in our study (n=101) (Additional file 5: Table S3).
Copy number alterations in known cancer genes
300 cancer genes were previously reported , of which 76 cancer genes were found within the altered regions (Table 1). We found that copy number alterations of these cancer genes were specific to histotypes as well; e.g. TPM3 amplification in endometrioid tumors; JAK2 deletion in mucinous; RB1 deletion in both clear cell and serous tumors; TP53 and MAP2K4 deletion in mucinous and serous tumors. ERBB2, a gene implicated in breast and EOC showed significant focal amplification in mucinous tumors but deletions in serous tumors. Evaluation of ERBB2 expression between mucinous and serous tumors in the 3 datasets showed the trend of over expression of ERBB2 in mucinous compared to serous (Additional file 6: Table S4). The focal amplification of ERBB2 has been observed in various studies [26–28], supporting our findings.
Identification of candidate driver genes in EOC histotypes
To identify candidate driver genes that might contribute to carcinogenesis of EOC, we looked for genes that showed association between copy number and gene expression (Methods, Additional file 1: Figure S1). Table 1 summarized the alterations and potential cancer driver genes (highlighted in bold) based on cytoband. Pathway analysis of potential driver genes in Table 1 showed top molecular functions of these genes to be involved with cell cycle, cellular development, growth and proliferation.
Candidate drivers in known cancer genes
Among 76 identified cancer genes listed in Table 1, several genes were potential drivers (highlighted in bold) in EOC in the association analysis. These include: FH, GMPS, PIK3CA, EIF4A2, ZNF384, and SS18L1 (amplifications in serous tumors); TET2, FGFR10P, NF1, ERBB2, and SH3GL1 (deletions in serous tumors) and ERBB2 (amplifications in mucinous tumors). Figure 4 shows the association between copy number and gene expression of ERBB2. The correlation was significant for all 3 datasets (Dataset1: R=0.80, p=3.47E-09, Dataset2: R=0.74, p=9.44E-6, Dataset3: R=0.79, p=1.14E-6, meta-p=3.90E-17), suggesting a driver mechanism. Amplification was observed in mucinous and deletion in serous tumors. Another interesting observation was MYC, TP53, KRAS, and BRCA1, genes reportedly to be commonly mutated in cancers but did not show significant association between copy number and gene expression. Similarly, two other genes reported to be mutated in EOC (PIK3R1 and STK11) also did not show potential driver mechanism.
Validation of ERBB2 expression
It is not within the scope of this study to validate the candidate copy number driver genes. As ERBB2 has potential for targeted therapy, we validated the expressions of ERBB2 via qPCR of 7 samples in Dataset1 that were found to be amplified or deleted. It should be noted that no more samples were available in Dataset2 for validation. Figure 5 shows the scatter plot between gene expression from microarray and qPCR (measured as fold-change). Significant correlation was observed (p=0.007) between microarray and qPCR. Four ERBB2 amplified samples (1 serous, 1 mucinous, 1 mucinous borderline, and 1 clear cell) were expressed correspondingly higher than the 3 serous samples with deleted ERBB2 (p=0.06, Wilcoxon). All these data support ERBB2 as a copy number driver gene in EOC.
The study reveals genomics diversity in EOC. It is conceivable that some of these alterations are involved in the tumorigenesis of EOC but the pathogenesis is likely regulated by aberrations of histotype specific alterations. By stratifying based on histotypes, we were able to identify alterations in the lower prevalent clear cell, endometrioid, and mucinous samples. An example is ERBB2. Several groups have investigated alterations of ERBB2 in EOC with mixed results [29, 30]. However, when stratified into histotypes, the high prevalence of ERBB2 amplification in mucinous is clearly evident in our results and other studies [26–28]. In addition, our results support ERBB2 as a potential copy number driver gene in EOC. The differences of ERBB2 copy number alterations amongst the histotypes could be due to the origin of histotypes in EOC. Our study demonstrates the importance of histotype-specific analyses where the differing copy number landscape amongst the histotypes adds to the mounting evidence that EOC should not be treated as one disease.
76 cancer genes listed in Table 1 were found to be copy number altered in EOC. They include ERBB2, TPM3, BRCA1, BRAF, KRAS, and PIK3CA; some of which are potential copy number drivers e.g. PIK3CA and BRAF in serous histotypes. Another interesting observation was KRAS, a gene reported to be mutated in mucinous tumors. In our study, KRAS was not found significantly altered in mucinous tumors where mutations are common. However, KRAS was significantly amplified in serous tumors where mutations are rare in high grade tumors. The reciprocal relationship between KRAS mutation and copy number alterations is also observed in gastric cancer . The 8q24.21 region harboring MYC, on the other hand, was altered in most histotypes other than mucinous. TCGA study also indicates MYC is highly amplified in high grade serous tumors. This suggests that MYC inhibitors may be applicable for these histotypes. For cancer genes which have been reported to have somatic or germline mutations in EOC, 3 were found to harbor copy deletions in serous histotype: PIK3R1, BRCA1, and STK11. For BRCA1, our finding was concordant with previous report that BRCA1 locus could be lost via either deletion or epigenetic silencing other than mutation in sporadic EOC .
A number of the candidate drivers in Table 1 are also implicated in translocation aberrations, e.g. TPM3, BCL9, GMPS, ZNF384 and SS18L1. It’s interesting that these genes were amplified in endometrioid and/or serous tumors. TPM3 and BCL9 reside in 1q21, a frequent site for chromosomal rearrangements. TPM3 was specifically amplified in endometrioid tumors and the gene has been shown to constitute a fusion gene with NTRK1 which belongs to the group of TRK oncogenes reported for papillary thyroid carcinoma . Interestingly, NTRK1 is also significantly amplified only in endometrioid tumors and further investigation is required to ascertain if this is due to gene fusion. BCL9 is a novel oncogene in Wnt signaling pathway, playing a critical role in epithelial-mesenchymal transition in colon epithelium and adenocarcinomas [33, 34]. Translocation of BCL9 has been reported with 14q32  and the gene was amplified in both endometrioid and serous tumors. Translocations for GMPS, ZNF384, and SS18L1 were also found in leukemia and synovial sarcomas [36–40] and all were amplified in serous tumors.
There are several drugs targeting the genes, e.g. for ERBB2, inhibitors include Trastuzumab, Lapatinib, and Pertuzumab. Lately, a clinical trial on combination of Pertuzumab, Trastuzumb and Docetaxel improved outcome of patients with HER2 positive metastatic breast cancer . BRAF mutations are more common in low grade serous while BRAF amplification is more common in high grade serous. Our data showed that it is a potential copy number driver and hence may be targetable by BRAF-inhibitors in serous tumors. Most BRAF inhibitors target various mutations and its efficacy on amplified BRAF is not yet well understood. A study has shown that BRAF amplified colorectal cancer cells acquired resistance to the MEK1/2 inhibitors selumetinib . PIK3CA is significantly amplified in serous histotypes and could be a potential target for PI3K inhibitors. In a study of PI3K inhibitor on breast and gynecologic malignancies harboring PIK3CA mutations, patients with the mutations treated with the inhibitor showed higher response rate than patients without the mutations .
In combining the 3 datasets, there was concern with regards to the genetic diversity amongst the Chinese Japanese, and Caucasian samples. The Hapmap  and Human Genome Diversity projects  have showed that these ethnic groups are different, though Chinese (CHB) and Japanese (JPT) tend to have high similarity in population structure. As genetic differences can be evaluated via principal component analysis (PCA) , we used PCA to assess the copy number data of the 3 cohorts. No distinct clustering between the groups (Additional file 7: Figure S3) was observed, suggesting that in this particular copy number landscape, the genetic effect is not evident and therefore has minimal effect in the analyses. We also used ANOVA test to assess ~200 housekeeping genes between the 3 datasets; none of the genes showed significance (Additional file 8: Table S5). Note that ERBB2 also did not show any significance. Nevertheless, genetic differences were taken into consideration in the preprocessing protocol. Individual dataset was normalized with respect to the relevant ethnic group from Hapmap data, i.e. Dataset1 with JPT, and Dataset2 with CHB. Ethnic-specific common structural polymorphism was also filtered out (see Data analysis) to ensure the copy number alterations identified in this study are de novo alterations in tumors.
We recognize that the regions identified could still be limited by the individual sample size of the histotypes. The larger number of copy number altered genes in serous tumors could be attributed to the larger sample size in this collection. We performed sub-sampling analyses to ascertain the effects and in addition, to ensure robustness of results, we used stringent criteria to filter the regions as well as criteria to consider CNA genes if it were supported by at least 2 dataset (Additional file 1: Figure S1). The flip side of this filtering was that true regions of alterations could be filtered out (as shown in the green area in Figure 1), leading to probably more false negatives. Nevertheless, we observed that despite the filtering and limited sample size of some histotypes, significant regions were still observed in the less prevalent histotypes; e.g. the 1p36.33, 2p11.1, 19q13.31, and 20q13.33 amplification and 9q32 deletion in clear cell tumors (n=29); 1q21.2-3 amplification in endometrioid tumors (n=20); 17q12 amplification and p24.1 deletion in mucinous tumors (n=19). Note that endometrioid tumors were not available in Dataset 2 although the total number of tumors was comparable with clear cell and mucinous. The concordance criteria of agreement on 2 datasets in the analytical workflow would thus bias the identification of regions for this histotype. Despite this, significant alterations were still observed for endometroid (e.g. TPM3) and given the stringent criteria; these are likely high confidence alterations. It should be noted that the samples were stratified according to the 4 main histotypes, including some borderline cases. Although borderline cases are presented clinically as a different subtype, they were included to simplify the stratification of histotypes and analyses. The significance of this approach can be seen in ERBB2, where both mucinous and mucinous borderline cases harbor amplification and corresponding up regulation of expression as well. This was similarly observed in other studies [26–28]. To assess if copy number alterations differ between borderline and non-borderline tumors and would thus cause bias in our analyses, we evaluated PCA of these samples (Additional file 9: Figure S4). No distinct clustering was observed between the borderline and non-borderline groups.
In summary, our study showed genomic diversity in EOC and highlighted distinct copy number alterations in histotypes that may have potential for drug targeted therapy. ERBB2 is significantly amplified in mucinous tumors and is a candidate copy number driver gene. By merging multiple datasets of similar platform, we demonstrated that CNA in the lower prevalent histotypes could be elucidated, even with limited sample size.
56 archived frozen tumor samples from the Department of Gynecology & Obstetrics, Kyoto University Graduate School of Medicine, Japan were profiled on microarray. It contained 12 clear cell carcinoma, 6 endometrioid adenocarcinoma, 2 mucinous adenocarcinoma, 5 mucinous-borderline tumors, 26 serous adenocarcinoma, and 5 serous-borderline tumors.
46 archived frozen tumor samples collected from Department of Obstetrics and Gynecology, Tri-Service General Hospital, Taiwan, containing 9 clear cell, 6 mucinous, and 31 serous.
GSE19539 consisting of 8 clear cell, 14 endometrioid, 6 mucinous, and 39 serous . Blood normal available in the dataset was used for normalization in concordance with the paper.
Human ovarian carcinoma samples
Two collections of archived flash frozen ovarian carcinoma samples were obtained from the Department of Gynecology & Obstetrics, Kyoto University Graduate School of Medicine, Japan (Dataset1, N= 56) and Department of Obstetrics and Gynecology, Tri-Service General Hospital, Taiwan (Dataset2, N=46). All samples were collected with the donor’s written informed consent. Ethical clearance has been approved by the Institutional Review Board for both institutes. All samples were reviewed by at least one pathologist from the respective institutes on histopathological typing and purity of samples. Tumor genomic DNAs were extracted by using phenol-chloroform extraction method. Tumor RNAs were extracted by using Qiazol followed by column clean-up using miRNeasy kit (Qiagen).
Copy number profiling
Affymetrix Genome-Wide Human SNP Arrays 6.0 (Affymetrix, Santa Clara, California) were used for copy number analysis according to the cytogenetics protocol from the manufacturer. Data was pre-processed and normalized with Hapmap JPT or CHB for Japan and Taiwan samples respectively using the Affymetric Genotyping Console. Copy number segments were obtained from the circular binary segmentation (CBS) algorithm  implemented in R package DNAcopy using default settings.
Gene expression profiling
Affymetrix GeneChip Human Gene 1.0ST Array was used for gene expression analysis according to the protocols from the manufacturer. Data was pre-processed and RMA normalized  using Affymetrix Gene Expression Console. Expressions for genes were mean-aggregated for each gene based on Affymetrix probes annotation. Note: 3 mucinous samples in Dataset2 and 2 serous samples in Dataset3 do not have corresponding gene expression data.
The gene expression and copy number datasets are MIAME compliant and have been submitted to National Centre for Biotechnology Information’s (NCBI) Gene Expression Omnibus (GEO) website, series accession number GSE30311.
Quantitative real time PCR
Total and miRNA was isolated from the ovarian carcinoma tissues using the miRNeasy kit (QIAGEN), of which 500ng were used to generate cDNA using the RT2 first strand kit (QIAGEN). For the qPCR run, 200ng of first strand cDNA was used per gene analysis. To determine the expression profile, ERBB2 transcript expression levels were normalized against the averaged expression levels of 5 housekeeping genes (ACTB, B2M, GADPH, HPRT and RPL13A).
Delta-Ct (ΔCt) and fold-change determination
Ct was determined using the SDS software (version 2.3, Applied Biosystems). Briefly, Ct values were determined by setting the baseline between cycle 2 of the run (total run: 40 cycles) and 2 cycles before the start of the first log-phase amplification. The threshold was set by positioning the limit to the lower third of the earliest amplification. ΔCt was calculated by the formulae:
whereby: Ct (GOI): Ct value of the respective gene of interest (GOI), Ct (HKG): average Ct values of the 5 housekeeping genes (HKG) used in the assay, Fold-change of the transcript is determined by the following formula:
To identify significant copy number altered regions, we used a 2-pronged workflow employing the GISTIC algorithm . GISTIC identifies copy number alterations based on the frequency as well as the log relative ratio (LRR) signals to compute the q value (false discovery rate). Default settings were used in the GISTIC analysis, and amplification and deletion thresholds were set at 0.2 and −0.2 respectively. Additional file 1: Figure S1 shows the 2-pronged workflow involving merged and individual copy number datasets to identify copy number alterations. Alterations were considered significant if it passed the following filtering criteria: (i) q < 0.25 (individual dataset), (ii) q < 0.05 (merged dataset), and (iii) concordance in 2 or more datasets. The significant regions were than mapped to genes (hg18 Refseq) by averaging the segments within each gene. ANOVA test was used to identify histotype-specific alterations. The analyses resulted in a list of significant gain and loss genes for each histotype, summarized in Figure 2 and Table 1 (known cancer genes).
To identify potential driver genes, non-parametric Spearman correlation was used to assess association between gene expression and copy number alterations of individual gene for each dataset (Additional file 1: Figure S1). Fisher’s combined probability test (meta-p)  was then used to combine the correlation statistics from each dataset to identify potential driver genes. This hypothesis driven association approach has been used to identify potential cancer driver genes [50, 51]. Potential driver genes of known cancer genes are listed in Table 1 (in bold).
PCA plots were generated using Partek Genomics Suite (Partek, Missouri, USA). Pathway analyses were performed using Ingenuity Pathway Analysis software (Ingenuity, California, USA). The frequency plot for copy number altered regions in Figure 1 was generated using the threshold of LRR ≥ |0.2|. All statistical analyses and plots were done using the R programming package (http://www.r-project.org).
Sub-sampling analyses to ascertain effects of sample size
To assess the effects of disparate sample size of the histotypes in the merged copy number data, multiple sub-sampling (with replacement) on the merged serous tumors was performed to ascertain the false positive and negative. The results showed that ≥97% of genes identified in sample size of 20–30 were also found in sample size of 101. However, 43-57% of genes found in the larger sample size were not identified in the smaller sample size datasets. In view of this, we have mainly confined comparison of genes found in the non-serous tumors.
Epithelial ovarian cancer
Polymerase chain reaction
Quantitative real-time polymerase chain reaction
False discovery rate
The cancer genome atlas
Principal component analysis
Circular binary segmentation
National Centre for Biotechnology Information
Gene expression omnibus
Gene of interest
Genomic identification of significant targets in cancer
Analysis of variance.
Sankaranarayanan R, Ferlay J: Worldwide burden of gynaecological cancer: the size of the problem. Best Pract Res Clin Obstet Gynaecol. 2006, 20 (2): 207-225. 10.1016/j.bpobgyn.2005.10.007.
Bonome T, Levine DA, Shih J, Randonovich M, Pise-Masison CA, Bogomolniy F, Ozbun L, Brady J, Barrett JC, Boyd J, et al: A gene signature predicting for survival in suboptimally debulked patients with ovarian cancer. Cancer Res. 2008, 68 (13): 5478-5486. 10.1158/0008-5472.CAN-07-6595.
Jemal A, Siegel R, Ward E, Murray T, Xu J, Smigal C, Thun MJ: Cancer statistics, 2006. CA Cancer J Clin. 2006, 56 (2): 106-130. 10.3322/canjclin.56.2.106.
Chuaqui RF, Cole KA, Emmert-Buck MR, Merino MJ: Histopathology and molecular biology of ovarian epithelial tumors. Ann Diagn Pathol. 1998, 2 (3): 195-207. 10.1016/S1092-9134(98)80007-5.
Khalil I, Brewer MA, Neyarapally T, Runowicz CD: The potential of biologic network models in understanding the etiopathogenesis of ovarian cancer. Gynecol Oncol. 116 (2): 282-285.
Sieben NL, Macropoulos P, Roemen GM, Kolkman-Uljee SM, Jan Fleuren G, Houmadi R, Diss T, Warren B, Al Adnani M, De Goeij AP, et al: In ovarian neoplasms, BRAF, but not KRAS, mutations are restricted to low-grade serous tumours. J Pathol. 2004, 202 (3): 336-340. 10.1002/path.1521.
Ahmed AA, Etemadmoghadam D, Temple J, Lynch AG, Riad M, Sharma R, Stewart C, Fereday S, Caldas C, Defazio A, et al: Driver mutations in TP53 are ubiquitous in high grade serous carcinoma of the ovary. J Pathol. 2010, 221 (1): 49-56. 10.1002/path.2696.
Enomoto T, Weghorst CM, Inoue M, Tanizawa O, Rice JM: K-ras activation occurs frequently in mucinous adenocarcinomas and rarely in other common epithelial tumors of the human ovary. Am J Pathol. 1991, 139 (4): 777-785.
Obata K, Morland SJ, Watson RH, Hitchcock A, Chenevix-Trench G, Thomas EJ, Campbell IG: Frequent PTEN/MMAC mutations in endometrioid but not serous or mucinous epithelial ovarian tumors. Cancer Res. 1998, 58 (10): 2095-2097.
Stratton MR, Campbell PJ, Futreal PA: The cancer genome. Nature. 2009, 458 (7239): 719-724. 10.1038/nature07943.
Deng N, Goh LK, Wang H, Das K, Tao J, Tan IB, Zhang S, Lee M, Wu J, Lim KH, et al: A comprehensive survey of genomic alterations in gastric cancer reveals systematic patterns of molecular exclusivity and co-occurrence among distinct therapeutic targets. Gut. 2012, 61 (5): 673-684. 10.1136/gutjnl-2011-301839.
Janku F, Stewart DJ, Kurzrock R: Targeted therapy in non-small-cell lung cancer–is it becoming a reality?. Nat Rev Clin Oncol. 2010, 7 (7): 401-414. 10.1038/nrclinonc.2010.64.
Barnes DJ, Melo JV: Management of chronic myeloid leukemia: targets for molecular therapy. Semin Hematol. 2003, 40 (1): 34-49. 10.1016/S0037-1963(03)70041-9.
Verma S, Lavasani S, Mackey J, Pritchard K, Clemons M, Dent S, Latreille J, Lemieux J, Provencher L, Verma S, et al: Optimizing the management of her2-positive early breast cancer: the clinical reality. Curr Oncol. 2010, 17 (4): 20-33.
Ramakrishna M, Williams LH, Boyle SE, Bearfoot JL, Sridhar A, Speed TP, Gorringe KL, Campbell IG: Identification of candidate growth promoting genes in ovarian cancer through integrated copy number and expression analysis. PLoS One. 2010, 5 (4): e9983-10.1371/journal.pone.0009983.
Gorringe KL, Ramakrishna M, Williams LH, Sridhar A, Boyle SE, Bearfoot JL, Li J, Anglesio MS, Campbell IG: Are there any more ovarian tumor suppressor genes? A new perspective using ultra high-resolution copy number and loss of heterozygosity analysis. Genes Chromosomes Cancer. 2009, 48 (10): 931-942. 10.1002/gcc.20694.
Haverty PM, Hon LS, Kaminker JS, Chant J, Zhang Z: High-resolution analysis of copy number alterations and associated expression changes in ovarian tumors. BMC Med Genomics. 2009, 2: 21-10.1186/1755-8794-2-21.
Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, Linhart D, Vivanco I, Lee JC, Huang JH, Alexander S, et al: Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. Proc Natl Acad Sci USA. 2007, 104 (50): 20007-20012. 10.1073/pnas.0710052104.
Rodabaugh KJ, Biggs RB, Qureshi JA, Barrett AJ, Welch WR, Bell DA, Berkowitz RS, Mok SC: Detailed deletion mapping of chromosome 9p and p16 gene alterations in human borderline and invasive epithelial ovarian tumors. Oncogene. 1995, 11 (7): 1249-1254.
Feltmate CM, Lee KR, Johnson M, Schorge JO, Wong KK, Hao K, Welch WR, Bell DA, Berkowitz RS, Mok SC: Whole-genome allelotyping identified distinct loss-of-heterozygosity patterns in mucinous ovarian and appendiceal carcinomas. Clin Cancer Res. 2005, 11 (21): 7651-7657. 10.1158/1078-0432.CCR-05-1008.
Campbell IG, Foulkes WD, Beynon G, Davis M, Englefield P: LOH and mutation analysis of CDKN2 in primary human ovarian cancers. Int J Cancer. 1995, 63 (2): 222-225. 10.1002/ijc.2910630213.
Devlin J, Elder PA, Gabra H, Steel CM, Knowles MA: High frequency of chromosome 9 deletion in ovarian cancer: evidence for three tumour-suppressor loci. Br J Cancer. 1996, 73 (4): 420-423. 10.1038/bjc.1996.75.
Cancer Genome Atlas Research Network: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474 (7353): 609-615. 10.1038/nature10166.
Gorringe KL, Jacobs S, Thompson ER, Sridhar A, Qiu W, Choong DY, Campbell IG: High-resolution single nucleotide polymorphism array analysis of epithelial ovarian cancer reveals numerous microdeletions and amplifications. Clin Cancer Res. 2007, 13 (16): 4731-4739. 10.1158/1078-0432.CCR-07-0502.
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer. 2004, 4 (3): 177-183. 10.1038/nrc1299.
Yan B, Choo SN, Mulyadi P, Srivastava S, Ong CW, Yong KJ, Putti T, Salto-Tellez M, Lim GS: Dual-colour HER2/chromosome 17 chromogenic in situ hybridisation enables accurate assessment of HER2 genomic status in ovarian tumours. J Clin Pathol. 2011, 64 (12): 1097-1101. 10.1136/jclinpath-2011-200082.
McAlpine JN, Wiegand KC, Vang R, Ronnett BM, Adamiak A, Kobel M, Kalloger SE, Swenerton KD, Huntsman DG, Gilks CB, et al: HER2 overexpression and amplification is present in a subset of ovarian mucinous carcinomas and can be targeted with trastuzumab therapy. BMC Cancer. 2009, 9: 433-10.1186/1471-2407-9-433.
Han CP, Hsu JD, Yao CC, Lee MY, Ruan A, Tyan YS, Yang SF, Chiang H: HER2 gene amplification in primary mucinous ovarian cancer: a potential therapeutic target. Histopathology. 2010, 57 (5): 763-764. 10.1111/j.1365-2559.2010.03689.x.
Nakayama K, Nakayama N, Jinawath N, Salani R, Kurman RJ, Shih Ie M, Wang TL: Amplicon profiles in ovarian serous carcinomas. Int J Cancer. 2007, 120 (12): 2613-2617. 10.1002/ijc.22609.
Pastor T, Popovic B, Gvozdenovic A, Boro A, Petrovic B, Novakovic I, Puzovic D, Lukovic L, Milasin J: Alterations of c-Myc and c-erbB-2 genes in ovarian tumours. Srp Arh Celok Lek. 2009, 137 (1–2): 47-51.
Senturk E, Cohen S, Dottino PR, Martignetti JA: A critical re-appraisal of BRCA1 methylation studies in ovarian cancer. Gynecol Oncol. 2010, 119 (2): 376-383. 10.1016/j.ygyno.2010.07.026.
Greco A, Miranda C, Pierotti MA: Rearrangements of NTRK1 gene in papillary thyroid carcinoma. Mol Cell Endocrinol. 2010, 321 (1): 44-49. 10.1016/j.mce.2009.10.009.
Deka J, Wiedemann N, Anderle P, Murphy-Seiler F, Bultinck J, Eyckerman S, Stehle JC, Andre S, Vilain N, Zilian O, et al: Bcl9/Bcl9l are critical for Wnt-mediated regulation of stem cell traits in colon epithelium and adenocarcinomas. Cancer Res. 2010, 70 (16): 6619-6628. 10.1158/0008-5472.CAN-10-0148.
Mani M, Carrasco DE, Zhang Y, Takada K, Gatt ME, Dutta-Simmons J, Ikeda H, Diaz-Griffero F, Pena-Cruz V, Bertagnolli M, et al: BCL9 promotes tumor progression by conferring enhanced proliferative, metastatic, and angiogenic properties to cancer cells. Cancer Res. 2009, 69 (19): 7577-7586. 10.1158/0008-5472.CAN-09-0773.
Huret JL, Senon S, Bernheim A, Dessen P: An Atlas on genes and chromosomes in oncology and haematology. Cell Mol Biol (Noisy-le-Grand). 2004, 50 (7): 805-807.
Storlazzi CT, Mertens F, Mandahl N, Gisselsson D, Isaksson M, Gustafson P, Domanski HA, Panagopoulos I: A novel fusion gene, SS18L1/SSX1, in synovial sarcoma. Genes Chromosomes Cancer. 2003, 37 (2): 195-200. 10.1002/gcc.10210.
Martini A, La Starza R, Janssen H, Bilhou-Nabera C, Corveleyn A, Somers R, Aventin A, Foa R, Hagemeijer A, Mecucci C, et al: Recurrent rearrangement of the Ewing's sarcoma gene, EWSR1, or its homologue, TAF15, with the transcription factor CIZ/NMP4 in acute leukemia. Cancer Res. 2002, 62 (19): 5408-5412.
La Starza R, Aventin A, Crescenzi B, Gorello P, Specchia G, Cuneo A, Angioni A, Bilhou-Nabera C, Boque C, Foa R, et al: CIZ gene rearrangements in acute leukemia: report of a diagnostic FISH assay and clinical features of nine patients. Leukemia. 2005, 19 (9): 1696-1699. 10.1038/sj.leu.2403842.
Zhong CH, Prima V, Liang X, Frye C, McGavran L, Meltesen L, Wei Q, Boomer T, Varella-Garcia M, Gump J, et al: E2A-ZNF384 and NOL1-E2A fusion created by a cryptic t(12;19)(p13.3; p13.3) in acute leukemia. Leukemia. 2008, 22 (4): 723-729. 10.1038/sj.leu.2405084.
Pegram LD, Megonigal MD, Lange BJ, Nowell PC, Rowley JD, Rappaport EF, Felix CA: t(3;11) translocation in treatment-related acute myeloid leukemia fuses MLL with the GMPS (GUANOSINE 5' MONOPHOSPHATE SYNTHETASE) gene. Blood. 2000, 96 (13): 4360-4362.
Baselga J, Cortes J, Kim SB, Im SA, Hegg R, Im YH, Roman L, Pedrini JL, Pienkowski T, Knott A, et al: Pertuzumab plus trastuzumab plus docetaxel for metastatic breast cancer. N Engl J Med. 2012, 366 (2): 109-119. 10.1056/NEJMoa1113216.
Little AS, Balmanno K, Sale MJ, Newman S, Dry JR, Hampson M, Edwards PA, Smith PD, Cook SJ: Amplification of the driving oncogene, KRAS or BRAF, underpins acquired resistance to MEK1/2 inhibitors in colorectal cancer cells. Sci Signal. 2011, 4 (166): ra17-10.1126/scisignal.2001752.
Janku F, Tsimberidou AM, Garrido-Laguna I, Wang X, Luthra R, Hong DS, Naing A, Falchook GS, Moroney JW, Piha-Paul SA, et al: PIK3CA mutations in patients with advanced cancers treated with PI3K/AKT/mTOR axis inhibitors. Mol Cancer Ther. 2011, 10 (3): 558-565. 10.1158/1535-7163.MCT-10-0994.
International HapMap Consortium: The International HapMap Project. Nature. 2003, 426 (6968): 789-796. 10.1038/nature02168.
Jakobsson M, Scholz SW, Scheet P, Gibbs JR, VanLiere JM, Fung HC, Szpiech ZA, Degnan JH, Wang K, Guerreiro R, et al: Genotype, haplotype and copy-number variation in worldwide human populations. Nature. 2008, 451 (7181): 998-1003. 10.1038/nature06742.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38 (8): 904-909. 10.1038/ng1847.
Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5 (4): 557-572. 10.1093/biostatistics/kxh008.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.
Fisher RA: Statistical methods for research workers. 1925, Edinburgh: Oliver and Boyd
Woo HG, Park ES, Lee JS, Lee YH, Ishikawa T, Kim YJ, Thorgeirsson SS: Identification of potential driver genes in human liver carcinoma by genomewide screening. Cancer Res. 2009, 69 (9): 4059-4066. 10.1158/0008-5472.CAN-09-0164.
Lando M, Holden M, Bergersen LC, Svendsrud DH, Stokke T, Sundfor K, Glad IK, Kristensen GB, Lyng H: Gene dosage, expression, and ontology analysis identifies driver genes in the carcinogenesis and chemoradioresistance of cervical cancer. PLoS Genet. 2009, 5 (11): e1000719-10.1371/journal.pgen.1000719.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/5/47/prepub
We thank Duke-NUS Genome Biology Facility for their help in profiling the samples in Dataset1. This study is supported by Khoo Discovery Award KDP/2008/0002 and KDP/2009/0006, and Duke-NUS core grant (LG), and National Research Foundation (RYH, JPT). We like to thank the reviewers for their constructive comments which have greatly improved the manuscript.
The authors declare no competing interests.
RYJH, JPT, and LG drafted the manuscript. GBC and LG performed the analyses. NM, HCL, IK collected the tissues. RYJH, SM, MKW extracted the RNA and DNA for profiling. RYJH and MKW carried out the PCR. LG conceived of the study. All authors read and approved the final manuscript.
Ruby YunJu Huang, Geng Bo Chen and Liang Goh contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Data Analysis Workflow. Two pronged approach for individual and merged datasets through selective threshold of GISTIC q-value and concordance in copy number analysis. As some histotypes have lower prevalence, filtering thresholds for individual and merged dataset were set at q<0.25 and q<0.05 respectively to overcome differences in sample size. In addition, any genomic alterations are supported by at least 2 datasets (i.e. concordance criteria). Specifically, the filtering criteria for histotype-specific regions were: (i) q < 0.25 (individual dataset), (ii) q < 0.05 (merged dataset), and (iii) concordance in 2 or more datasets. This resulted in a list of significant gains and loss regions. To identify copy number driver genes that are specific to histotype, copy number segments were mapped to genes and ANOVA was used to identify the differentially altered genes. This resulted in a list of histotype-specific altered genes. Spearman correlation between gene expression and copy number was then used to assess potential driver genes in each individual dataset. (JPEG 123 KB)
Additional file 3: Figure S2: Venn diagram of copy number altered genes between the 4 histotypes. Left: amplified genes; Right: deleted genes. Clear cell tumors had the highest number of common altered genes with serous tumors while endometrioid tumors had the lowest number of common altered genes. (JPEG 743 KB)
Additional file 7: Figure S3: Principal component analysis of copy number altered gene from the merged datasets. The plot shows that there is minimal copy number alterations difference between the 3 datasets. (JPEG 338 KB)
Additional file 9: Figure S4: Principal component analysis of copy number altered genes from the merged datasets showing borderline and non-borderline tumors. Borderline tumors were available only in serous and mucinous histotypes. No distinct clustering was observed between borderline and non-borderline tumors for (a) mucinous and (b) serous. (JPEG 472 KB)
About this article
Cite this article
Huang, R.Y., Chen, G.B., Matsumura, N. et al. Histotype-specific copy-number alterations in ovarian cancer. BMC Med Genomics 5, 47 (2012). https://doi.org/10.1186/1755-8794-5-47