Combined array-comparative genomic hybridization and single-nucleotide polymorphism-loss of heterozygosity analysis reveals complex genetic alterations in cervical cancer
BMC Genomics volume 8, Article number: 53 (2007)
Cervical carcinoma develops as a result of multiple genetic alterations. Different studies investigated genomic alterations in cervical cancer mainly by means of metaphase comparative genomic hybridization (mCGH) and microsatellite marker analysis for the detection of loss of heterozygosity (LOH). Currently, high throughput methods such as array comparative genomic hybridization (array CGH), single nucleotide polymorphism array (SNP array) and gene expression arrays are available to study genome-wide alterations. Integration of these 3 platforms allows detection of genomic alterations at high resolution and investigation of an association between copy number changes and expression.
Genome-wide copy number and genotype analysis of 10 cervical cancer cell lines by array CGH and SNP array showed highly complex large-scale alterations. A comparison between array CGH and SNP array revealed that the overall concordance in detection of the same areas with copy number alterations (CNA) was above 90%. The use of SNP arrays demonstrated that about 75% of LOH events would not have been found by methods which screen for copy number changes, such as array CGH, since these were LOH events without CNA. Regions frequently targeted by CNA, as determined by array CGH, such as amplification of 5p and 20q, and loss of 8p were confirmed by fluorescent in situ hybridization (FISH). Genome-wide, we did not find a correlation between copy-number and gene expression. At chromosome arm 5p however, 22% of the genes were significantly upregulated in cell lines with amplifications as compared to cell lines without amplifications, as measured by gene expression arrays. For 3 genes, SKP2, ANKH and TRIO, expression differences were confirmed by quantitative real-time PCR (qRT-PCR).
This study showed that copy number data retrieved from either array CGH or SNP array are comparable and that the integration of genome-wide LOH, copy number and gene expression is useful for the identification of gene specific targets that could be relevant for the development and progression in cervical cancer.
The development of cervical cancer is a multi-step process. This process starts with infection of the epithelial layer with high-risk human papilloma virus (HPV). Cervical intraepithelial neoplasia (CIN) grade I, II and III, representing reversible stages of disease, is believed to develop towards neoplasm. From CIN I onwards genomic alterations and aneuploidy can be detected [1, 2] and when disease is progressing towards CIN II and III, genetic alterations become more evident . HPV E6 and E7 can induce cellular transformation due to interaction with cellular proteins involved in cell cycle control, apoptosis and genomic instability . The induction of genomic instability seems to be particularly important for the establishment and development of an invasive tumour .
So far, different studies have investigated genomic alterations in cervical cancer by means of metaphase comparative genomic hybridization (mCGH) and microsatellite marker analysis for the detection of loss of heterozygosity (LOH). Chromosomal changes involving loss of 2q, 3p, 4p, 4q, 5q, 6q, 11q, 13q and 18q regions and gain of 1q, 3q, 5p and 8q regions at various stages of disease have been identified [1, 6–8]. Some of these imbalances have been associated with decreased survival such as LOH of 18q  or with the transition from severe dysplasia to invasive carcinoma such as gain of 3q . However, few genes, targeted by these genomic imbalances and possibly involved in carcinogenesis, have been discovered due to the lack of resolution when applying mCGH  whereas classical LOH studies are laborious. Array comparative genomic hybridization (array CGH) is an established high resolution method to study the whole genome for chromosomal amplifications and deletions . Although array CGH has relatively few printed elements on the arrays, the strength of this technique is the high spatial resolution . A limitation of array CGH is the lack of genotype information. Therefore, it does not provide information about regions of LOH without copy number alteration (CNA) such as mitotic recombination or gene conversion events. The under-reprepresentation of LOH in copy number methods may be very high with up to 2/3 of LOH events being undetected because of copy-neutral and copy-gain LOH events, as was found in pancreatic cancer . Recently, high-density oligonucleotide-based single polymorphism arrays (SNP arrays) have been used to identify copy number and LOH of chromosomal regions [15, 16]. Appropriate software tools for the detection of CNA from SNP data are becoming available. The advantage of a combined SNP-CGH approach is the identification of allele specific gain and loss by SNP array and the robust copy number detection by array CGH.
Different factors, including changes in genome copy number, can disrupt proper gene functioning. A few studies showed correlations between chromosome CNA and gene expression changes of the affected regions [17, 18]. In this study, genome-wide LOH and copy number analyses of 10 cervical cancer cell lines were investigated using SNP and CGH arrays and the sensitivity to detect CNA by the 2 platforms was compared. We aimed at identifying common regions of gains and losses and determining the influence of CNA on gene expression to find genes involved in cervical carcinogenesis.
Genome-wide LOH and copy number analyses of 10 cervical cancer cell lines were investigated using SNP and CGH arrays to identify common regions of amplification. Also the two platforms were compared in their sensitivity to detect CNA. Furthermore, data from gene expression arrays were combined with array CGH data to investigate an association between gene expression and copy number.
CGH and SNP-LOH analysis
Various genomic alterations affecting the majority of chromosomes were found in cervical cancer cell lines. The CGH findings for all amplifications, gains and losses detected in these cell lines are displayed in Figure 1A. Gains or amplifications of 5p, 5q, 8q, 9q, 17q and 20q occurred in most of the cell lines. Physical loss of chromosomal regions at 4p, 8p, 11q, 13q and 18q, occurred in at least half of the cell lines. Regions frequently targeted by gain or loss with small overlapping regions were present at 8p, 8q, 9p, 11q and 17q. The minimal region of overlap at 8p23 involved TUSC3, a tumor suppressor gene that showed loss in 5 cell lines. At 8q24, the minimal region of overlap included c-Myc which showed gain in 8 out of 10 cell lines. Amplification at 9q33-q34 included endoglin, a gene involved in tumor angiogenesis and predictive for metastasis in cervical cancer . Loss at 11q25 involved 2 members of the IgLON family of cell adhesion molecules, OPCML and HNT. The minimal region of overlap at 17q11-q21 included ERBB2 and TOP2A, genes that are amplified in breast cancer . As for SNP-LOH, most chromosome arms were affected by LOH in the cell lines (Figure 1B). More than 5 overlapping regions of LOH were evident for chromosome region 5q, 6, 8p, 10q, 11q, 14q, 18 and 20p.
Comparison between SNP and BAC microarray
Concordance between array CGH and SNP array in the detection of CNA was investigated. Overall, the BAC and SNP arrays identified exactly the same regions of gain in 91% of cases, when taking an average of all cell lines (Table 1). Similarly, the detection of the same regions with loss between the 2 techniques was on average 94% (Table 1). Subsequently, we investigated the overlap between copy number loss and LOH and copy number gain and LOH for both platforms. For LOH to occur, measured as a reduction in copy number, physical loss is required. Thus all cases with copy number reductions are expected to have undergone LOH. In case of physical loss, around 60% was observed to encompass regions of LOH (Table 2). Overall, copy number reductions as measured by array CGH showed a higher overlap with LOH than SNP array (Table 2), namely 66% for array CGH and 57% for SNP array, when taking an average of all cell lines. The occurrence of physical gain involves a change in copy number and thus a change in the proportion of genotypes which would imply a call for LOH. Around 28% of copy number gain was found to cover regions of LOH (Table 2) and frequencies were comparable between SNP array and array CGH. Also, the overall percentage of LOH in copy neutral, copy gain and copy loss events was investigated for each cell line. The average frequency of LOH from all cell lines was 21–26% for loss, 17% for gain and 60% for copy neutral events (Table 3). Although array CGH and SNP array are comparable in the level of detection of copy number gain and loss, alterations in regions of LOH were often not found by copy number analysis because these LOH changes did not result in CNA (Figure 1 and Table 3). For example, 6 of 10 cell lines showed LOH of 6p while array CGH detected physical loss of 6p in only one cell line. This suggests that the other 5 cell lines underwent mitotic recombination upon loss of one arm, thereby not affecting copy number. Examples of complex genomic alterations are shown for chromosome 8 and 20 of SiHa and chromosome 3 of CaSki (Figure 2). In these examples (Figure 2A and 2B), LOH is caused by loss of one chromosome, accompanied by duplication of the remaining chromosome and subsequent deletions and translocations. The copy number data from SNP array as deduced from the dChip data (in blue lines) and array CGH (red lines) clearly follow the same pattern in these examples, supporting their compliance.
Correlation between DNA copy number and gene expression
To correlate the chromosomal changes to alterations in gene expression profiles, we have searched for changes in the expression of genes located in the regions targeted by chromosomal gains and deletions. This was done by combining gene expression data from Affymetrix focus arrays with array CGH data. The average gene expression values of all cell lines were used as a reference for each cell line to determine increases or decreases in gene expression. Genome-wide, no correlation was observed between copy-number and expression. Since 5p, 8p23, 8q24, 9q33-q34, 11q25 and 20q are regions with many alterations, we focused on these areas. For 20q we compared cell lines harboring amplifications with cell lines containing gains since no cell line was devoid of 20q gain or amplification but we did not find significant differences in gene expression. Neither did we find significant changes in gene expression at 8p23, 8q24 and 11q25. At 9q33-q34, cell lines with amplification (CSCC7, CC11+, CC10A, CC10B, SiHa and HeLa) showed a significantly higher expression compared to cell lines with gains or no CNA correlation for 9 out of 86 genes (10%); POMT1, FLJ35348, LHX3, NDOR, SH3GLB2, ST6GALNAC4, CDK9, PTAN1 and TRAF2 (p < 0.05, T-test), of which the latter 3 have been implicated in tumorigenesis. Fourteen of 64 genes (22%) at 5p showed a significantly higher expression in cell lines with amplification (CaSki, SiHa and HeLa) compared to cell lines with gains or no CNA (Table 4).
Validation of genetic alterations by FISH and LOH
The DNA copy number retrieved from the array CGHs was validated by FISH. Based on the gained, amplified or lost regions of 5p, 20q and 8p, BACs were selected that covered potential oncogenes or tumour suppressor genes. FISH analysis on all cell lines was performed using the BAC clones encompassing SKP2, TRIO, hTERT, GDNF and ANKH genes on 5p, TUSC3 on 8p and ZNF-217, CYP24A1, MYBL2 and AIB1 on 20q (Table 5). Besides additional copies of the area covered by BACs at 5p, other complex genetic alterations such as translocations and inversions were detected (Figure 3A). TUSC3 is a putative tumour suppressor gene at 8p that is covered by 2 overlapping BAC clones, RP11-44L18 and RP11-184C1 BAC. Loss of the area at 8p was concordant between array CGH and FISH for CC10A, CC11- and SiHa (Figure 3B). Signals for BACs at chromosome 20q were equal to the number of centromere signals or showed additional copies (Table 5 and Figure 3C), indicating that the whole arm or the chromosome was gained, which is in agreement with the array CGH data (Figure 1A). Overall, FISH assays confirmed the accuracy in which DNA CNA were detected by array CGH for these loci and highlighted the gains and amplifications of the ZNF-217, CYP24A1, MYBL2 and AIB1 genes at 20q, SKP2, TRIO, hTERT, GDNF and ANKH genes at 5p and the loss of the TUSC3 gene at 8p. LOH of 6p detected by SNP array in cell line CC11-, CC11+ and CC10A was confirmed by LOH analysis of the tumors from which the cell lines originate, using microsatellite markers (Figure 4).
Validation of gene expression by qRT-PCR
To validate the microarray data, upregulated expression of 3 potential oncogenes (SKP2, ANKH and TRIO) were analyzed by qRT-PCR. These genes were located on 5p that was shown to be gained or amplified by CGH in half of the cell lines. Microarray data showed upregulation of ANKH in 2 of the 3 cell lines with amplification of 5p which did not reach statistical significance. Five normal cervical epithelial cell cultures (NPE) were included in the analysis to investigate the difference in expression level between normal and tumour cells, with and without amplification, from the cervix. Expression was significantly different between cell lines with 5p amplification, cell lines with 5p gains or no 5p alterations, and NPE (Figure 5). At chromosome arm 8p, deletion of the tumour suppressor gene TUSC3, as found by FISH analysis (Table 5 and Figure 3B), was evident in CaSki, SiHa, CC10A, CC11- and CC11+. No probe set on the expression array was available for TUSC3. Therefore, qRT-PCR was carried out on all cell lines. Only SiHa and CaSki showed a substantial drop in expression, compared to the expression levels of the other cell lines and NPE (data not shown).
The present study describes genome wide LOH and copy number analyses of 10 cervical cancer cell lines using SNP and CGH arrays to identify common regions of amplification and deletion. Also, the sensitivity to detect CNA by the 2 platforms was compared. Furthermore, integrated copy number and gene expression analyses were performed to investigate an association between CNA and gene expression.
CNA retrieved from array CGH and genotype information analyzed from SNP arrays revealed very complex large-scale changes. Chromosome arms that were found to be gained and/or amplified in most of the cell lines by array CGH in this study were 5p, 8p, 8q, 9p, 9q, 17q and 20q. Especially gains of 5p, 8q and 20q have been frequently reported in cervical cancer [6, 7, 10, 21–25]. Also the 1q and 3q arms, often found to be gained [7, 10, 21, 24–26], were overrepresented in 5 of the 10 cell lines. Gains of 17q and chromosome 9 are less consistently reported [6, 22]. Gain of 8q24, present in 8 cell lines, included c-Myc in the minimal region of overlap. Amplification of the c-Myc oncogene is a frequent event in cervical tumors, found in 25% to more than half of the tumours [27, 28]. Chromosome arms affected by loss in half, or more, of the cell lines were 4p, 8p, 11q and 13q of which especially 4p and 13q have been reported regularly in cervical cancer [7, 10, 22, 25, 26, 29]. Six cell lines showed loss of 11q25 that includes OPCML in the minimal region of overlap. OPCML is a tumour suppressor gene in epithelial ovarian cancer . Although the genotype data have to be interpreted with some caution because of the use of cell lines, many of the frequent changes that we observed were in concordance with previous studies.
The SNP-LOH array data showed that most chromosome arms were affected by LOH in the cell lines. More than 5 overlapping regions of LOH were evident for chromosome region 5q, 6p, 6q, 8p, 10q, 11q, 14q, 18p, 18q and 20p. LOH of the majority of these chromosomes, specifically 5q, 6p, 8p, 10q, 11q and 18, has been reported previously as frequent events in cervical cancer [9, 31, 32]. Whereas LOH of 8p resulted in copy number loss, LOH of 6p did not result in CNA for most cases, indicative of mitotic recombination events.
A comparison between CGH and SNP array revealed that the overall concordance in detection of the same areas with CNA was on average 91% for gain and 94% for loss (Table 1). Thus both platforms are comparable in the detection of particular regions with gain or loss. Additionally, the use of high-density SNP arrays showed that cervical cancer cell lines harbor multiple regions of LOH which are not detected by methods which screen for CNA such as array CGH. Accordingly, comparison of array CGH and SNP array in genome-wide analysis shows that detection of CNA by SNP array is reliable and has the added advantage of providing genotype information. However, LOH was not consistently detected by SNP array because many copy number changes did not result in true reduction to homozygosity since most of the cervical cancer cell lines are multiploid (Figure 2C and Table 2). As expected, the areas with gain and amplification showed less LOH than regions with deletion. The overlap between copy-reduction regions and LOH was 10% higher for array CGH than for SNP array (Table 2), which suggests that array CGH is slightly better in assigning the correct regions for CNA. Although some regions of LOH are missed, thereby particularly underestimating the amount of LOH in areas with gain and amplification, an estimation was made for the distribution of LOH in regions with loss, gain and copy-neutral regions. We found that copy-neutral and copy-gain regions accounted for around 75% of all LOH identified (17% represented copy-gain LOH) while LOH associated with copy-reduction occurred at a rate of about 25% (Table 3). These percentages are similar to a recently published paper from Calhoun et al. that described only 32% of LOH events to be associated with copy-reduction . Further improvement in the data from the SNP arrays may be possible using higher density SNP arrays, thereby increasing resolution and the ability to assess subtle copy number changes.
Genome-wide, we did not find a correlation between copy-number changes and expression. Lack of an association may be due to the different histological types from which the cell lines originate (Squamous carcinoma, Adenocarcinomas and Adenosquamous carcinomas) with complex genetic aberrations. Besides increased DNA dosage which may result in overexpression, transcription is influenced by other factors such as the presence and quantity of transcription and repressive factors; CpG methylation; histone methylation and acetylation status, and possibly microRNA's . Thus, amplification or deletion on its own does not have to be a prerequisite for a change in gene expression. Amplification of 5p, however, correlated with a significant higher gene expression for 22% of the genes, including SKP2, TRIO, OSMR, RPL37 and PDCD6 (Table 4). The latter genes were previously reported in cancer, either associated with increased expression or with growth stimulatory properties [34–38].
Array CGH findings of CNA at 5p, 8p and 20q were confirmed by FISH analysis. Gains and amplifications at 5p, as found by array CGH, were validated by BACs covering candidate oncogenes, SKP2, TRIO, ANKH, GDNF and hTERT. SKP2 is important for cell cycle progression and its inhibition decreases proliferation of tumour cells . Besides gain of the region covering the SKP2 gene, we found a significant increase in expression between normal cervical epithelial cells, tumour cells without amplification and tumour cells with amplification. Correlation between amplification of chromosome 5p and increased expression of SKP2 was previously shown in HPV-immortalized keratinocytes . Upregulated expression of SKP2, ANKH and TRIO, as determined by qRT-PCR in cell lines with 5p amplification, confirmed gene expression data from Affymetrix focus arrays. TRIO, which has a putative role in cell cycle regulation, was found to be amplified and highly expressed in bladder cancer  that associated with proliferation and invasive phenotype. ANKH, which is involved in tissue calcification , was found to be amplified in bladder cancer  and small cell lung cancer cell lines . Thus, upregulation of expression through amplification may be the case for TRIO, ANKH and SKP2 in some neoplasms, including cervical cancer, as indicated by our data. This may hold for hTERT and PRKAA1 as well since it was reported that amplification at 5p in cervical cancer correlated with increased gene and/or protein expression [43, 44]. The telomeric region at 8p includes a candidate tumour suppressor gene, TUSC3 [45, 46]. However, physical loss of this region did not correlate with a decrease in gene expression of TUSC3, as determined by qRT-PCR. All cell lines showed gains or amplifications of 20q, a chromosome arm which is often amplified in cancer, including cervical cancer [10, 22, 24, 47]. At 20q increased copy numbers, detected by array CGH, were confirmed by FISH. Candidate oncogenes, ZNF-217, CYP24A1, MYBL2 and AIB2, encompassing the region analyzed by FISH, were previously reported to be amplified in other types of cancer as well [48–50]. However, no significant difference in gene expression was found between tumours with gain and tumours with amplification of 20q. Also in oesophageal adenocarcinoma, no correlation was detected between amplification of ZNF-217/CYP24A1 at 20q and increased mRNA expression .
From this study we conclude that copy number data retrieved from either array CGH or SNP are comparable. The SNP array provided a reliable copy number as well as genotype information, which was suitable for the analysis of complex genetic alterations present in cervical cancer cell lines. Frequently occurring gains, amplifications and deletions on chromosome arms 5p, 20q and 8q were verified by FISH. Amplification at 5p associated with overexpression of candidate oncogenes. Accordingly, the possibility to integrate genome-wide LOH, copy number and gene expression is useful for the identification of gene specific targets that could be relevant for the development and progression in cervical cancer.
Cell lines, cultures and patient material
Human cervical cancer cell lines, HeLa, SiHa, CaSki (ATCC, Rockville, MD), CSCC1, CSCC7, CC8, CC10A, CC10B, CC11-, and CC11+  were maintained in RPMI 1640 (Life Technologies, Grand Island, NY), supplemented with 10% fetal calf serum, streptomycin and penicillin at 37°C in a humid atmosphere containing 5% CO2. CC10A, CC10B, CC11- and CC11+ are sub clones from the same cell line, CC10 and CC11 respectively. Five normal primary cervical epithelial cell cultures (NPE) were prepared as described . Two archival formalin-fixed, paraffin-embedded cervical carcinomas were obtained from the tissue bank of the Department of Pathology (LUMC, Leiden, The Netherlands).
DNA and RNA isolation
DNA from all the cultured cell lines was isolated using the salting out procedure according to Miller and Polesky . Total RNA was isolated using TRIzol reagent (Gibco BRL/Life Technologies, Breda, The Netherlands) from cells that were grown to approximately 60% confluence in 250 ml culture flasks (Greiner, Frickenhausen, Germany). Total RNA was cleaned with Rneasy cleanup system columns (Qiagen GmbH, Hilden, Germany).
Array CGH was performed as previously described . Slides containing ~3,500 BACs in triplicate were made at the Leiden University Medical Centre. The particular BAC set used (Welcome Trust Sanger Institute (UK)) comprises clones spaced at ~1 Mb density over the full genome, a set of subtelomeric sequences for each chromosome arm, and a few hundred probes selected for their involvement in oncogenesis. Information regarding the full set is available in the "Cytoview" window of the Sanger Center mapping database site, Ensembl. The test/reference ratios were normalized for the median of the ratios of all features. The BioConductor package Limma (Smyth and Speed 2003) version 2.6.0 was used to perform background correction with method 'Edwards' and within array normalization with method 'Printtiploess', between array normalization was not performed. After normalization the replicate spots were averaged. Probes with less than 2 valid replicate values, and probes that showed a standard deviation above 0.1 between the replicate values were excluded. Chromosomal regions demonstrating a ratio of between 1.11 and 1.41 were scored as 'gained' and a ratio of 1.41 or greater was scored as 'amplified'. Chromosomal regions showing a ratio of between 0.84 and 0.73 were scored as lost whereas a second threshold for loss was set for regions showing a ratio of less than 0.73.
SNP 10K array
Affymetrix 10K SNP array (Santa Clara, CA) containing 11,555 biallelic polymorphic sequences randomly distributed throughout the genome, excluding the Y chromosome, was used. The median distance between SNPs is ~105 kb. The average heterozygosity for these SNPs is 0.37, with an average minor allele frequency of 0.25. DNA samples, including two control DNA samples from Affymetrix, were assayed according to the manufacturer's manual protocol (GeneChip Mapping Assay manual). Genotyping was performed using GDAS (Affymetrix, Santa Clara, CA) and copy number values were calculated using dChip version 1.3. Stretches of homozygous SNPs equal or longer than 24 calls were defined as LOH.
Gene expression microarray
Affymetrix focus arrays (Santa Clara, CA), representing 8,793 human sequences from the NCBI Refseq database, were used for mRNA expression profiling according to the manufacturer's instructions. Data files of all microarrays were normalized using the Robust Multi-Chip Average (GCRMA) method . Quality control was performed on the arrays using Affy-Probe Level Model (AffyPLM) . Microarrays for each cell line were performed in duplicate. The signal intensities of the replicates were averaged after normalization. To identify differentially expressed genes between cell lines with and without gain, amplification or deletion of a chromosome arm or region, unpaired two-sample t-tests were performed.
To compare the data from the 3 high dimensional techniques; array CGH, SNP array, and expression array, the copy number and expression values were transformed to contain log 2, median centered values. Access to complete microarray data sets will be available at NCBI Gene Expression Omnibus. Normalization of array CGH data and expression arrays resulted in values that are equivalent to log2 intensity values. The data retrieved from dChip version 1.3 from the SNP arrays showed a copy number of 2 for unaffected probes on a linear scale. These data were divided by 2 and log2 transformed. A per gene median centering to zero was performed. Array CGH smoothing was performed using a smoother that retains sharp transitions between regions with genomic alterations. The correlation between gene expression and physical genomic alterations was determined by computing the Pearson correlation between the expression of each gene to the smoothed array CGH value of the nearest BAC on the chromosome. The resulting p-values were subjected to the Benjamini and Yekutieli false discovery rate procedure . The data retrieved from dChip version 1.3 from the SNP arrays showed a copy number of 2 for unaffected probes on a linear scale.
Fluorescent in situ hybridization (FISH)
Metaphase preparations of all cell lines were obtained using colcemid according to standard procedures . BAC probes for the genes selected on chromosome 5, chromosome 8 and chromosome 20 (Table 5) were ordered from the BACPAC Resource Centre at Children's Hospital Oakland Research Institute (Oakland, CA). All BAC probes were fluorescein-12-dUTP- or digoxigenin-12-dUTP-labelled (Roche, Basel, Switzerland). The following centromere probes were used: the pG-A16 probe for centromere 5; the D8Z2 probe for centromere 8; and the p3.4 probe for chromosome 20. All centromere probes were kindly provided by J. Wiegant (Department of Molecular and Cellular Biology, LUMC) and labelled with biotin-16-dUTP (Roche, Basel, Switzerland). Triple-colour FISH experiments were performed as previously described . Centromere and BAC signals of 20 metaphases were analyzed per probe combination per cell line.
Genomic DNA was isolated from flow-sorted tumour cell subpopulations derived from formalin-fixed, paraffin-embedded cervical carcinomas as previously described . Eight fluorescein-labelled primer pairs of microsatellite markers comprising MOGC, D6S265, C125, TY2A, TNFα, D6S273, D6S294 and D6S1666 , spanning altogether 26 MB from 6p22.1 to 6p12.1, were used. PCR, electrophoresis and analysis was performed as described previously . Comparing normal and tumour DNA, a reduction of more than 50% of one allele was assigned as complete LOH.
Quantitative real-time PCR (qRT-PCR)
qRT-PCR was performed as described previously . Expression of the genes of interest was normalized by geometric averaging of multiple internal control genes (CPSF6, HNRPM, EEF1A1, RPL11 and RPL13) using the Genorm program . Out of 5 normalization genes the best 3 were selected using this program; EEF1A1, RPL13 and RPL11. Primer sequences for RPL11 are 5'-ACTTCGCATCCGCAAACTCT (FWD) and 5'-AAGGTGTTGGAGCAGXTCACA (REV). Primer sequences for EEF1A1, RPL13 were reported previously . Commercially available primers were used for TUSC3, ANKH and TRIO (Superarray, Frederick, USA). Primer sequences for SKP2 are: 5'-CCTATCACTCAGTCGGTGCTATGA (FWD) and 5'-AATCGTGCCAGATGGTACCCT (REV). To determine differential gene expression between normal cervical epithelial cell cultures (NPE), cell lines with and without deletion, gain or amplification, one-way ANOVA was performed.
Narayan G, Pulido HA, Koul S, Lu XY, Harris CP, Yeh YA, Vargas H, Posso H, Terry MB, Gissmann L, Schneider A, Mansukhani M, Rao PH, Murty VV: Genetic analysis identifies putative tumor suppressor sites at 2q35-q36.1 and 2q36.3-q37.1 involved in cervical cancer progression. Oncogene. 2003, 22: 3489-3499. 10.1038/sj.onc.1206432.
Melsheimer P, Klaes R, Doeberitz MV, Bastert G: Prospective clinical study comparing DNA flow cytometry and HPV typing as predictive tests for persistence and progression of CIN I/II. Cytometry. 2001, 46: 166-171. 10.1002/cyto.1101.
Kersemaekers AM, van de Vijver MJ, Kenter GG, Fleuren GJ: Genetic alterations during the progression of squamous cell carcinomas of the uterine cervix. Genes Chromosomes Cancer. 1999, 26: 346-354. 10.1002/(SICI)1098-2264(199912)26:4<346::AID-GCC9>3.0.CO;2-D.
Fan X, Chen JJ: Regulation of cell cycle progression and apoptosis by the papillomavirus E6 oncogene. Crit Rev Eukaryot Gene Expr. 2004, 14: 183-202. 10.1615/CritRevEukaryotGeneExpr.v14.i3.30.
Nishimura M, Furumoto H, Kato T, Kamada M, Aono T: Microsatellite instability is a late event in the carcinogenesis of uterine cervical cancer. Gynecol Oncol. 2000, 79: 201-206. 10.1006/gyno.2000.5940.
Rao PH, Arias-Pulido H, Lu XY, Harris CP, Vargas H, Zhang FF, Narayan G, Schneider A, Terry MB, Murty VV: Chromosomal amplifications, 3q gain and deletions of 2q33-q37 are the frequent genetic changes in cervical carcinoma. BMC Cancer. 2004, 4: 5-10.1186/1471-2407-4-5.
Huang FY, Kwok YK, Lau ET, Tang MH, Ng TY, Ngan HY: Genetic abnormalities and HPV status in cervical and vulvar squamous cell carcinomas. Cancer Genet Cytogenet. 2005, 157: 42-48. 10.1016/j.cancergencyto.2004.06.002.
Kersemaekers AM, Hermans J, Fleuren GJ, van de Vijver MJ: Loss of heterozygosity for defined regions on chromosomes 3, 11 and 17 in carcinomas of the uterine cervix. Br J Cancer. 1998, 77: 192-200.
Kersemaekers AM, Kenter GG, Hermans J, Fleuren GJ, van de Vijver MJ: Allelic loss and prognosis in carcinoma of the uterine cervix. Int J Cancer. 1998, 79: 411-417. 10.1002/(SICI)1097-0215(19980821)79:4<411::AID-IJC17>3.0.CO;2-7.
Heselmeyer K, Macville M, Schrock E, Blegen H, Hellstrom AC, Shah K, Auer G, Ried T: Advanced-stage cervical carcinomas are defined by a recurrent pattern of chromosomal aberrations revealing high genetic instability and a consistent gain of chromosome arm 3q. Genes Chromosomes Cancer. 1997, 19: 233-240. 10.1002/(SICI)1098-2264(199708)19:4<233::AID-GCC5>3.0.CO;2-Y.
Lichter P, Joos S, Bentz M, Lampel S: Comparative genomic hybridization: uses and limitations. Semin Hematol. 2000, 37: 348-357. 10.1016/S0037-1963(00)90015-5.
Oostlander AE, Meijer GA, Ylstra B: Microarray-based comparative genomic hybridization and its applications in human genetics. Clin Genet. 2004, 66: 488-495. 10.1111/j.1399-0004.2004.00322.x.
Ylstra B, van , Carvalho B, Brakenhoff RH, Meijer GA: BAC to the future! or oligonucleotides: a perspective for micro array comparative genomic hybridization (array CGH). Nucleic Acids Res. 2006, 34: 445-450. 10.1093/nar/gkj456.
Calhoun ES, Gallmeier E, Cunningham SC, Eshleman JR, Hruban RH, Kern SE: Copy-number methods dramatically underestimate loss of heterozygosity in cancer. Genes Chromosomes Cancer. 2006, 45: 1070-1071. 10.1002/gcc.20365.
Zhao X, Li C, Paez JG, Chin K, Janne PA, Chen TH, Girard L, Minna J, Christiani D, Leo C, Gray JW, Sellers WR, Meyerson M: An integrated view of copy number and allelic alterations in the cancer genome using single nucleotide polymorphism arrays. Cancer Res. 2004, 64: 3060-3071. 10.1158/0008-5472.CAN-03-3308.
Huang J, Wei W, Zhang J, Liu G, Bignell GR, Stratton MR, Futreal PA, Wooster R, Jones KW, Shapero MH: Whole genome DNA copy number changes identified by high density oligonucleotide arrays. Hum Genomics. 2004, 1: 287-299.
Tsafrir D, Bacolod M, Selvanayagam Z, Tsafrir I, Shia J, Zeng Z, Liu H, Krier C, Stengel RF, Barany F, Gerald WL, Paty PB, Domany E, Notterman DA: Relationship of gene expression and chromosomal abnormalities in colorectal cancer. Cancer Res. 2006, 66: 2129-2137. 10.1158/0008-5472.CAN-05-2569.
Carter SL, Eklund AC, Kohane IS, Harris LN, Szallasi Z: A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers. Nat Genet. 2006, 38: 1043-1048. 10.1038/ng1861.
Brewer CA, Setterdahl JJ, Li MJ, Johnston JM, Mann JL, McAsey ME: Endoglin expression as a measure of microvessel density in cervical cancer. Obstet Gynecol. 2000, 96: 224-228. 10.1016/S0029-7844(00)00864-4.
Harris LN, Yang L, Liotcheva V, Pauli S, Iglehart JD, Colvin OM, Hsieh TS: Induction of topoisomerase II activity after ErbB2 activation is associated with a differential response to breast cancer chemotherapy. Clin Cancer Res. 2001, 7: 1497-1504.
Hidalgo A, Monroy A, Arana RM, Taja L, Vazquez G, Salcedo M: Chromosomal imbalances in four new uterine cervix carcinoma derived cell lines. BMC Cancer. 2003, 3: 8-10.1186/1471-2407-3-8.
Harris CP, Lu XY, Narayan G, Singh B, Murty VV, Rao PH: Comprehensive molecular cytogenetic characterization of cervical cancer cell lines. Genes Chromosomes Cancer. 2003, 36: 233-241. 10.1002/gcc.10158.
Thein A, Trkova M, Fox M, Parrington J: The application of comparative genomic hybridization to previously karyotyped cervical cancer cell lines. Cancer Genet Cytogenet. 2000, 116: 59-65. 10.1016/S0165-4608(99)00108-9.
Matthews CP, Shera KA, McDougall JK: Genomic changes and HPV type in cervical carcinoma. Proc Soc Exp Biol Med. 2000, 223: 316-321. 10.1046/j.1525-1373.2000.22345.x.
Kirchhoff M, Rose H, Petersen BL, Maahr J, Gerdes T, Lundsteen C, Bryndorf T, Kryger-Baggesen N, Christensen L, Engelholm SA, Philip J: Comparative genomic hybridization reveals a recurrent pattern of chromosomal aberrations in severe dysplasia/carcinoma in situ of the cervix and in advanced-stage cervical carcinoma. Genes Chromosomes Cancer. 1999, 24: 144-150. 10.1002/(SICI)1098-2264(199902)24:2<144::AID-GCC7>3.0.CO;2-9.
Umayahara K, Numa F, Suehiro Y, Sakata A, Nawata S, Ogata H, Suminami Y, Sakamoto M, Sasaki K, Kato H: Comparative genomic hybridization detects genetic alterations during early stages of cervical cancer progression. Genes Chromosomes Cancer. 2002, 33: 98-102. 10.1002/gcc.1215.
Ocadiz R, Sauceda R, Cruz M, Graef AM, Gariglio P: High correlation between molecular alterations of the c-myc oncogene and carcinoma of the uterine cervix. Cancer Res. 1987, 47: 4173-4177.
Munzel P, Marx D, Kochel H, Schauer A, Bock KW: Genomic alterations of the c-myc protooncogene in relation to the overexpression of c-erbB2 and Ki-67 in human breast and cervix carcinomas. J Cancer Res Clin Oncol. 1991, 117: 603-607. 10.1007/BF01613296.
Aubele M, Zitzelsberger H, Schenck U, Walch A, Hofler H, Werner M: Distinct cytogenetic alterations in squamous intraepithelial lesions of the cervix revealed by laser-assisted microdissection and comparative genomic hybridization. Cancer. 1998, 84: 375-379. 10.1002/(SICI)1097-0142(19981225)84:6<375::AID-CNCR10>3.0.CO;2-1.
Sellar GC, Watt KP, Rabiasz GJ, Stronach EA, Li L, Miller EP, Massie CE, Miller J, Contreras-Moreira B, Scott D, Brown I, Williams AR, Bates PA, Smyth JF, Gabra H: OPCML at 11q25 is epigenetically inactivated and has tumor-suppressor function in epithelial ovarian cancer. Nat Genet. 2003, 34: 337-343. 10.1038/ng1183.
Mullokandov MR, Kholodilov NG, Atkin NB, Burk RD, Johnson AB, Klinger HP: Genomic alterations in cervical carcinoma: losses of chromosome heterozygosity and human papilloma virus tumor status. Cancer Res. 1996, 56: 197-205.
Mitra AB, Murty VV, Li RG, Pratap M, Luthra UK, Chaganti RS: Allelotype analysis of cervical carcinoma. Cancer Res. 1994, 54: 4481-4487.
Tchurikov NA: Molecular mechanisms of epigenetics. Biochemistry (Mosc ). 2005, 70: 406-423. 10.1007/s10541-005-0131-2.
Yokoi S, Yasui K, Mori M, Iizasa T, Fujisawa T, Inazawa J: Amplification and overexpression of SKP2 are associated with metastasis of non-small-cell lung cancers to lymph nodes. Am J Pathol. 2004, 165: 175-180.
Zheng M, Simon R, Mirlacher M, Maurer R, Gasser T, Forster T, Diener PA, Mihatsch MJ, Sauter G, Schraml P: TRIO amplification and abundant mRNA expression is associated with invasive tumor growth and rapid tumor cell proliferation in urinary bladder cancer. Am J Pathol. 2004, 165: 63-69.
Boing I, Stross C, Radtke S, Lippok BE, Heinrich PC, Hermanns HM: Oncostatin M-induced activation of stress-activated MAP kinases depends on tyrosine 861 in the OSM receptor and requires Jak1 but not Src kinases. Cell Signal. 2006, 18: 50-61. 10.1016/j.cellsig.2005.03.015.
Loging WT, Reisman D: Elevated expression of ribosomal protein genes L37, RPP-1, and S2 in the presence of mutant p53. Cancer Epidemiol Biomarkers Prev. 1999, 8: 1011-1016.
la Cour JM, Mollerup J, Winding P, Tarabykina S, Sehested M, Berchtold MW: Up-regulation of ALG-2 in hepatomas and lung cancer tissue. Am J Pathol. 2003, 163: 81-89.
Jiang F, Caraway NP, Li R, Katz RL: RNA silencing of S-phase kinase-interacting protein 2 inhibits proliferation and centrosome amplification in lung cancer cells. Oncogene. 2005, 24: 3409-3418. 10.1038/sj.onc.1208459.
Dowen SE, Neutze DM, Pett MR, Cottage A, Stern P, Coleman N, Stanley MA: Amplification of chromosome 5p correlates with increased expression of Skp2 in HPV-immortalized keratinocytes. Oncogene. 2003, 22: 2531-2540. 10.1038/sj.onc.1206296.
Ho AM, Johnson MD, Kingsley DM: Role of the mouse ank gene in control of tissue calcification and arthritis. Science. 2000, 289: 265-270. 10.1126/science.289.5477.265.
Coe BP, Henderson LJ, Garnis C, Tsao MS, Gazdar AF, Minna J, Lam S, Macaulay C, Lam WL: High-resolution chromosome arm 5p array CGH analysis of small cell lung carcinoma cell lines. Genes Chromosomes Cancer. 2005, 42: 308-313. 10.1002/gcc.20137.
Zhang A, Zheng C, Hou M, Lindvall C, Wallin KL, Angstrom T, Yang X, Hellstrom AC, Blennow E, Bjorkholm M, Zetterberg A, Gruber A, Xu D: Amplification of the telomerase reverse transcriptase (hTERT) gene in cervical carcinomas. Genes Chromosomes Cancer. 2002, 34: 269-275. 10.1002/gcc.10071.
Huang FY, Chiu PM, Tam KF, Kwok YK, Lau ET, Tang MH, Ng TY, Liu VW, Cheung AN, Ngan HY: Semi-quantitative fluorescent PCR analysis identifies PRKAA1 on chromosome 5 as a potential candidate cancer gene of cervical cancer. Gynecol Oncol. 2006, 103: 219-225. 10.1016/j.ygyno.2006.02.028.
Levy A, Dang UC, Bookstein R: High-density screen of human tumor cell lines for homozygous deletions of loci on chromosome arm 8p. Genes Chromosomes Cancer. 1999, 24: 42-47. 10.1002/(SICI)1098-2264(199901)24:1<42::AID-GCC6>3.0.CO;2-F.
Lerebours F, Olschwang S, Thuille B, Schmitz A, Fouchet P, Buecher B, Martinet N, Galateau F, Thomas G: Fine deletion mapping of chromosome 8p in non-small-cell lung carcinoma. Int J Cancer. 1999, 81: 854-858. 10.1002/(SICI)1097-0215(19990611)81:6<854::AID-IJC3>3.0.CO;2-1.
Wilting SM, Snijders PJ, Meijer GA, Ylstra B, van den Ijssel PR, Snijders AM, Albertson DG, Coffa J, Schouten JP, van de Wiel MA, Meijer CJ, Steenbergen RD: Increased gene copy numbers at chromosome 20q are frequent in both squamous cell carcinomas and adenocarcinomas of the cervix. J Pathol. 2006, 209: 220-230. 10.1002/path.1966.
Tanner MM, Grenman S, Koul A, Johannsson O, Meltzer P, Pejovic T, Borg A, Isola JJ: Frequent amplification of chromosomal region 20q12-q13 in ovarian cancer. Clin Cancer Res. 2000, 6: 1833-1839.
Yang SH, Seo MY, Jeong HJ, Jeung HC, Shin J, Kim SC, Noh SH, Chung HC, Rha SY: Gene copy number change events at chromosome 20 and their association with recurrence in gastric cancer patients. Clin Cancer Res. 2005, 11: 612-620.
Wong LJ, Dai P, Lu JF, Lou MA, Clarke R, Nazarov V: AIB1 gene amplification and the instability of polyQ encoding sequence in breast cancer cell lines. BMC Cancer. 2006, 6: 111-10.1186/1471-2407-6-111.
van Dekken H, Vissers K, Tilanus HW, Kuo WL, Tanke HJ, Rosenberg C, Ijszenga M, Szuhai K: Genomic array and expression analysis of frequent high-level amplifications in adenocarcinomas of the gastro-esophageal junction. Cancer Genet Cytogenet. 2006, 166: 157-162. 10.1016/j.cancergencyto.2005.11.002.
Koopman LA, Szuhai K, van Eendenburg JD, Bezrookove V, Kenter GG, Schuuring E, Tanke H, Fleuren GJ: Recurrent integration of human papillomaviruses 16, 45, and 67 near translocation breakpoints in new cervical cancer cell lines. Cancer Res. 1999, 59: 5615-5624.
Hazelbag S, Fleuren GJ, Baelde JJ, Schuuring E, Kenter GG, Gorter A: Cytokine profile of cervical cancer cells. Gynecol Oncol. 2001, 83: 235-243. 10.1006/gyno.2001.6378.
Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16: 1215-10.1093/nar/16.3.1215.
Knijnenburg J, Szuhai K, Giltay J, Molenaar L, Sloos W, Poot M, Tanke HJ, Rosenberg C: Insights from genomic microarrays into structural chromosome rearrangements. Am J Med Genet A. 2005, 132: 36-40.
Ensembl. 2007, [http://www.ensembl.org]
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.
Gautier L, Cope L, Bolstad BM, Irizarry RA: affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.
Eilers PH, de Menezes RX: Quantile smoothing of array CGH data. Bioinformatics. 2005, 21: 1146-1153. 10.1093/bioinformatics/bti148.
Kleinfeld RG, Sisken JE: Morphological and kinetic aspects of mitotic arrest by and recovery from colcemid. J Cell Biol. 1966, 31: 369-379. 10.1083/jcb.31.3.369.
Wiegant J, Ried T, Nederlof PM, van der PM, Tanke HJ, Raap AK: In situ hybridization with fluoresceinated DNA. Nucleic Acids Res. 1991, 19: 3237-3241. 10.1093/nar/19.12.3237.
Corver WE, Ter Haar NT, Dreef EJ, Miranda NF, Prins FA, Jordanova ES, Cornelisse CJ, Fleuren GJ: High-resolution multi-parameter DNA flow cytometry enables detection of tumour and stromal cell subpopulations in paraffin-embedded tissues. J Pathol. 2005, 206: 233-241. 10.1002/path.1765.
The GDB Human Genome Database. 2007, [http://www.GDB.org/]
Haven CJ, van Puijenbroek M, Karperien M, Fleuren GJ, Morreau H: Differential expression of the calcium sensing receptor and combined loss of chromosomes 1q and 11q in parathyroid carcinoma. J Pathol. 2004, 202: 86-94. 10.1002/path.1489.
Kloth JN, Fleuren GJ, Oosting J, de Menezes RX, Eilers PH, Kenter GG, Gorter A: Substantial changes in gene expression of Wnt, MAPK and TNFalpha pathways induced by TGF-beta1 in cervical cancer cell lines. Carcinogenesis. 2005, 26: 1493-1502. 10.1093/carcin/bgi110.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-10.1186/gb-2002-3-7-research0034.
The study was supported by the Dutch Cancer Society, grant no. RUL2001-2465.
JK performed RNA isolation, preparation of RNA for microarray hybridization of gene expression arrays, qRT-PCR experiments, analysis of array data and preparation of the manuscript. JO performed (integrative) statistical analysis of arrays, interpretation of data and contributed to the preparation of figures. TW supervised hybridization of SNP arrays, contributed to the design of the study and critically reviewed the manuscript. KS supervised technical aspects of array CGH hybridization and critically reviewed the manuscript. JK performed array CGH hybridization. AG contributed to the design of the study and critically reviewed the manuscript. GK contributed to the design of the study and critically reviewed the manuscript. GF contributed to the design of the study and critically reviewed the manuscript. EJ performed DNA isolation, preparation of DNA and microarray hybridization of Affymetrix 10K SNP arrays, FISH experiments, analysis of array data and contributed to the preparation of the manuscript. All authors read and approved the manuscript.
About this article
Cite this article
Kloth, J.N., Oosting, J., van Wezel, T. et al. Combined array-comparative genomic hybridization and single-nucleotide polymorphism-loss of heterozygosity analysis reveals complex genetic alterations in cervical cancer. BMC Genomics 8, 53 (2007). https://doi.org/10.1186/1471-2164-8-53
- Cervical Cancer
- Human Papilloma Virus
- Cervical Intraepithelial Neoplasia
- Copy Number Alteration
- Copy Number Gain