Skip to main content
  • Research article
  • Open access
  • Published:

Identification of differentially expressed genes using an annealing control primer system in stage III serous ovarian carcinoma

Abstract

Background

Most patients with ovarian cancer are diagnosed with advanced stage disease (i.e., stage III-IV), which is associated with a poor prognosis. Differentially expressed genes (DEGs) in stage III serous ovarian carcinoma compared to normal tissue were screened by a new differential display method, the annealing control primer (ACP) system. The potential targets for markers that could be used for diagnosis and prognosis, for stage III serous ovarian cancer, were found by cluster and survival analysis.

Methods

The ACP-based reverse transcriptase polymerase chain reaction (RT PCR) technique was used to identify DEGs in patients with stage III serous ovarian carcinoma. The DEGs identified by the ACP system were confirmed by quantitative real-time PCR. Cluster analysis was performed on the basis of the expression profile produced by quantitative real-time PCR and survival analysis was carried out by the Kaplan-Meier method and Cox proportional hazards multivariate model; the results of gene expression were compared between chemo-resistant and chemo-sensitive groups.

Results

A total of 114 DEGs were identified by the ACP-based RT PCR technique among patients with stage III serous ovarian carcinoma. The DEGs associated with an apoptosis inhibitory process tended to be up-regulated clones while the DEGs associated with immune response tended to be down-regulated clones. Cluster analysis of the gene expression profile obtained by quantitative real-time PCR revealed two contrasting groups of DEGs. That is, a group of genes including: SSBP1, IFI6 DDT, IFI27, C11orf92, NFKBIA, TNXB, NEAT1 and TFG were up-regulated while another group of genes consisting of: LAMB2, XRCC6, MEF2C, RBM5, FOXP1, NUDCP2, LGALS3, TMEM185A, and C1S were down-regulated in most patients. Survival analysis revealed that the up-regulated genes such as DDAH2, RNase K and TCEAL2 might be associated with a poor prognosis. Furthermore, the prognosis of patients with chemo-resistance was predicted to be very poor when genes such as RNase K, FOXP1, LAMB2 and MRVI1 were up-regulated.

Conclusion

The DEGs in patients with stage III serous ovarian cancer were successfully and reliably identified by the ACP-based RT PCR technique. The DEGs identified in this study might help predict the prognosis of patients with stage III serous ovarian cancer as well as suggest targets for the development of new treatment regimens.

Peer Review reports

Background

Ovarian cancer is a complex disease, characterized by successive accumulation of multiple molecular alterations in both the cells undergoing neoplastic transformation and host cells [1]. These anomalies disturb the expression of genes that control critical cell processes, leading to the initiation of tumorigenesis and development. At the time of diagnosis most patients with ovarian cancer have advanced stage disease (i.e., stage III-IV) where surgery and chemotherapy results in an approximately 25% overall 5-year survival rate. Consequently, ovarian cancer is the leading cause of death from a gynecological malignancy. Epithelial ovarian cancer (EOC) accounts for 90% of all ovarian cancers; there is significant heterogeneity within the EOC group. For example, histologically defined subtypes such as serous, endometrioid, mucinous, and low- and high-grade malignancies all have variable clinical manifestations and underlying molecular signatures [2].

Gene expression has been extensively applied to screening for the prognostic factors associated with ovarian cancer; identification of such factors would help to determine patient prognosis. Studies have focused on differential gene expression between tumor and normal tissues [3], distinguishing between histological subtypes [4] and identifying differences between invasive tumors and those with low malignant potential [5]. However, to date, the use of differentially expressed genes (DEG) have not been implemented in ovarian cancer therapies; this is mainly because their reliability and validity have not yet been well established. Microarray technology permits large scale analysis of expression surveys to identify the genes that have altered expression as a result of disease. However, microarray data is notorious for its unreliable reproducibility of DEGs across platforms and laboratories, as well as validation problems associated with prognostic signatures [6]. In addition, identification of a gene responsible for a specialized function during a certain biological stage can be difficult to determine because the gene might be expressed at low levels, whereas the bulk of mRNA transcripts within a cell are abundant [7].

To screen DEGs in low concentrations, while minimizing false positive results, the polymerase chain reaction (PCR) based technique has been used. One screening method, differential display, requires PCR using short arbitrary primers. This method is simple, rapid and only requires small amounts of total RNA. However, many investigators have reported significantly high false-positive rates [8] and poor reproducibility of the results [9] because of nonspecific annealing by the short arbitrary primers. Recently, the annealing control primer (ACP) system has been developed; this technique provides a primer with annealing specificity to the template and allows only genuine products to be amplified [10]. The structure of the ACP includes (i) a 3" end region with a target core nucleotide sequence that substantially complements the template nucleic acid of hybridization; (ii) a 5" end region with a non-target universal nucleotide sequence; and (iii) a polydeoxyinosine [poly(dI)] linker bridging the 3" and 5" end sequences. Because of the high annealing specificity during PCR using the ACP system, the application of the ACP to DEG identification generates reproducible, accurate, and long (100 bp to 2 kb) PCR products that are detectable on agarose gels.

In this study, the ACP-based PCR method was used to identify the DEGs of patients with stage III serous ovarian cancer, the findings were compared to normal ovarian tissue. A total of 60 arbitrary ACPs were used and 114 DEGs were identified by sequencing differentially expressed bands. For the confirmation of differential expression of the DEGs, quantitative real-time PCR was performed on 38 selected DEGs; the results showed good agreement with the ACP findings. These results could be used as preliminary data for further study of the molecular mechanism underlying stage III serous ovarian cancer.

Methods

Patient information

After obtaining written informed consent from all patients included in the study, samples of primary epithelial ovarian cancer were snap frozen in liquid nitrogen and stored at -80°C. Analysis of tissues from patients was approved by the Institutional Review Board of The Catholic University of Korea (Seoul, Korea). The histopathological diagnoses were determined using the WHO criteria, and the tumor histotype was serous adenocarcinoma in all patients. Classification of cancer stage and grade was performed according to the International Federation of Gynecology and Obstetrics (FIGO). A total of 16 patients with serous ovarian carcinoma were enrolled in this study, all patients were diagnosed as stage IIIC with high-grade cancer (grade 3).

ACP-based GeneFishingâ„¢ reverse transcription polymerase chain reaction

Total RNAs from the ovarian tissues of the serous carcinoma were isolated by gentle homogenization using Trizol®. The normal human ovary total RNA was purchased from Stratagene (Total RNA Human Ovary, #540071). The RNA was used for the synthesis of first-strand cDNAs by reverse transcriptase. Reverse transcription was performed 1.5 hours at 42°C in a final reaction volume of 20 μl containing 3 μg of the purified total RNA, 4 μl of 5" reaction buffer (Promega, Madison, WI, USA), 5 μl of dNTPs (2 mmol each), 2 μl of 10 μM dT-ACP1 (5"-CTGTGAATGCTGCGACTACGATIIIII(T)18)-3", where "I" represents deoxyinosine), 0.5 μL of RNasin® RNase Inhibitor (40 U/μl, Promega), and 1 μl of Moloney murine leukemia virus reverse transcriptase (200 U/μl, Promega). First-strand cDNAs were diluted by the addition of 80 μL of RNase-free water for the GeneFishing PCR and stored -20°C until use.

DEGs were screened by ACP-based PCR method using the GeneFishing™ DEG kits (Seegene, Seoul, South Korea). Briefly, second-strand cDNA synthesis was conducted at 50°C (low stringency) during one cycle of first-stage PCR in a final reaction volume of 49.5 μl containing 3-5 μl (about 50 ng) of diluted first-strand DNA cDNA, 5 μl of 10x PCR buffer plus Mg (Roche Applied Science, Mannheim, Germany), 5 μl of dNTP (each 2 mM), 1 μl of 10 μM dT-ACP2 (5'-CTGTGAATGCTGCGACTACGATIIIII(T)15-3"), and 1 μl of 10 μM arbitrary ACP. The tube containing the reaction mixture was kept at 94°C while 0.5 μl of Taq DNA polymerase was added to the reaction mixture (5 U/μl, Roche Applied Science). Sixty PCR reactions for each sample were carried out with 60 arbitrary ACPs, respectively. The PCR protocol for second-strand synthesis was one cycle at 94°C for 1 min, followed by 50°C for 3 min, and 72°C for 1 min. After completion of second-strand DNA synthesis, 40 cycles were performed. Each cycle involved denaturation at 94°C for 40 sec, annealing at 65°C for 40 sec, extension at 72°C for 40 sec, and a final extension at 72°C to complete the reaction. The amplified PCR products were separated in 2% agarose gel and stained with ethidium bromide. The overall scheme of the experiment is shown in Figure 1.

Figure 1
figure 1

A schematic diagram of the experimental procedure.

Cloning and sequencing

The differentially expressed bands were extracted from the gel using the GENCLEAN® II Kit (Q-BIO gene, Carlsbad, CA.,USA), and directly cloned into a TOPO TA® cloning vector (Invitrogen, Karlsruhe, Germany) according to the manufacturer's instructions. The cloned plasmids were sequenced with an ABI PRISM® 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA., USA). Complete sequences were analyzed by searching for similarities using the Basic Local Alignment Search Tool (BLAST) search program at the Genbank database of the National Center for Biotechnology Information (NCBI).

Quantitative real-time polymerase chain reaction and statistical analysis

For the confirmation of the differential expression of DEGs, quantitative real-time PCR was carried out for 38 DEGs selected from 114 DEGs. The concentrations of the reagents were adjusted to reach a final volume of 20 μL containing 5 ng of cDNA template, 10 μl of SYBR® Premis Ex Taq™ II (Takara Bio,Otsu, Japan), 0.4 μl of ROX™ reference Dye II, 0.4 μl of 10 μM forward and reverse primers of DEGs with β-actin as an internal control (Table 1). The cDNA templates were constructed with the total RNA extracted from 16 ovarian cancer tissues and normal human ovary total RNA. The PCR amplification protocol was 50°C for 2 min and 95°C for 10 min followed by 40 cycles of 95°C for 30 sec, 60°C for 30 sec, and 72°C for 30 sec. The real-time PCR analysis was performed on an Applied Biosystems Prism 7900 Sequence Detection System (Applied Biosystems). Relative quantification with the data obtained was performed according to the user's manual. The fold change for gene expression, between the cancer and normal samples, was calculated by using the threshold cycle (CT): fold change = 2-ΔΔ C T , ΔΔC T = [(CT of gene of interest - CT of β-actin)cancer sample- (CT gene of interest - CT of β-actin)normal sample)]. The fold change was log2 transformed for the cluster and survival analysis. The R packages mclust and survival http://www.r-project.org were used for the cluster and survival analysis, respectively.

Table 1 Primer sequences of 38 DEGs and β-actin used for the quantitative real-time PCR.

Results

Differentially expressed genes (DEGs) in stage III serous ovarian cancers

The patients included in this study ranged in age between 38 and 69 (mean age 53.3 ± 7.5 years). All 16 patients had a diagnosis of FIGO stage III papillary serous ovarian carcinoma. To identify genes that showed a predominant change of expression in patients with stage III serous ovarian cancer, the total RNAs from 16 serous ovarian tissues of stage III and the normal human ovary total RNA (Stratagene, #540071) were individually subjected to ACP-based RT PCR analysis using a combination of 60 arbitrary primers and two anchored oligo (dT) primers (dT-ACP1 and dT-ACP2). All PCR amplicons were compared on agarose gels (Figure 2). When the bands generated by the normal sample showed a clear difference compared to the bands generated on the cancer sample, the band was defined as a differentially expressed band. After all poor appearing bands were excluded, the differentially expressed bands were extracted, amplified using the TOPO TA Cloning Kit (Invitrogen, Cat. #K4500-1) and sequenced. The sequences of 114 DEGs were obtained. The DNA sequence of each DEG was analyzed by searching for similarities using the BLASTX program at the Genbank database (NIH, MD, USA). Table 2 shows the 114 DEGs assessed by Genbank and the best homologues.

Table 2 Annotation of the 114 DEGs by the BLAST search.
Figure 2
figure 2

An example of GeneFishingâ„¢ using an arbitrary ACP in combination with an oligo (dT) ACP as indicated in the Methods section. M and N represents a 100-bp size marker generated by Forever 100-bp Ladder Personalizer (Seegene, Seoul, South Korea) and a normal sample, respectively. The lanes 1-16 include each of the cancer samples from the 16 patients. The bands showing a clear difference between the normal and cancer samples, here marked with a red arrow, were excised from the gel for further cloning and sequencing.

The DEGs identified from the clones that were up-regulated included: AP000659 (U4, ARHGEF12), NM_013974.1 (U8, DDAH2), NM_022873.2 (U16, IFI6), NM_020529 (U21, NFKBIA), BC023599 (U30, TFG), and NG_011683 (U39, SEPT9). The up-regulation of SEPT9 mRNA was reported in a bank of ovarian tumors, which included benign, borderline and malignant tumors [11]. The genes including DDAH2, IFI6 and NFKBIA are known to be involved in the apoptosis inhibitory process while ARHGEF12 and TGF have been implicated in signaling pathways. The DEGs identified from down-regulated clones included: AK291099 (D7, BST2), NM_001025159.1 (D10, CD74), BC004186 (D18, GNB1), NG_011554 (D34, HTRA1), BC095490 (D36, IGKC), and BC001693 (D41, LGALS1). The down-regulation of HTRA1 was associated with ovarian cancer metastasis [12]. The genes including BST2, CD74 and IGKC were associated with the immune system; while GNB1 and LGALS1 were related to the modulation of cell-cell interaction and the G protein coupled receptor protein signaling pathway, respectively.

Confirmation of ACP observation by quantitative real-time PCR and cluster analysis

To confirm the efficacy of the ACP system, confirmation of the differential expression of DEGs was performed with quantitative real-time PCR for 38 DEGs selected from the total 114 DEGs using a specific primer pair for each gene (Table 1 and Additional file 1). The expression ratio of the cancer to normal sample was calculated by using CT and then was log2 transformed (see Methods section for detail). The DEGs were considered differentially expressed if the log2 ratio was > 1.0 or < -1.0. Differential expression was clearly observed in all 38 DEGs, which indicates a high reliability of the ACP system (Figure 3).

Figure 3
figure 3

Clustering of the log 2 expression ratios of the cancer to normal samples measured by the quantitative real-time PCR for the 38 DEGs from the 16 patients and its representation as a heat map. Patients are ordered along the X-axis and genes along the Y-axis.

For the detection of more conserved expression patterns in patients with stage III serous ovarian cancer, cluster analysis was performed using the R package mclust http://www.r-project.org. Thirty eight DEGs were divided into four groups according to their expression profiles with assignment of each gene to a group (Figure 3). A clear contrast in expression patterns was noted between groups 1 and 4. That is, the overall up- and down-regulation in group 1 was 86.8% and 5.6%, respectively; the overall up- and down-regulation in group 4 was 9.3% and 60.3%, respectively, when up-regulation corresponds to log2 ratio > 1 and down-regulation to log2 ratio < -1 (Table 3). This regulation pattern in each group was well maintained at various thresholds used for definition of differential expression. These findings suggest that the genes in groups 1 and 4 might be used as potential markers for prognosis in patients with stage III serous ovarian cancer. The group 1 consisted of: NM_003143.1 (SSBP1), NM_022873.2 (IFI6), NM_001355.3 (DDT), NM_005532 (IFI27), NM_207429.2 (C11orf92), NM_020529 (NFKBIA), NM_032470 (TNXB), EF177379 (NEAT1) and BC02359 (TFG); group 4 was composed of: NM_002292.3 (LAMB2), AK025219, AK293439 (XRCC6), NM_001131005 (MEF2C), AK302766 (RBM5), NM_032682 (FOXP1), NG_001229 (NUDCP2), BC001120 (LGALS3), NM_032508.1 (TMEM185A), and NM_201442 (C1S) (Table 3). All genes except for NEAT1 in group 1 were identified from up-regulated clones while all genes except for NUDCP2 in group 4 were identified from down-regulated clones (Tables 2 and 3). These results were in good agreement with the ACP findings.

Table 3 Clustering of 38 DEGs according to the expression profiles of 16 patients with serous ovarian cancer; Up-regulation corresponds to log2 ratio > threshold while down-regulation to log2 ratio < minus value of threshold.

Survival analysis

The Kaplan-Meier method was performed using the 38 DEGs that were up- and down-regulated. The up and down regulated genes were considered according to the log2 expression ratio > 0 and < 0, respectively. A significant difference, in the overall survival (p values < 0.05) between the up- and down-regulated group, was observed for three DEGs including NM_013974.1 (dimethylarginine dimethylaminohydrolase 2, DDAH2), NM_001004333.3 (ribonuclease K, RNase K) and NM_080390 (transcription elongation factor A (SII)-like 2, TCEAL2) (Figure 4). The overall survival decreased with up-regulation of these genes. DDAH2 predominates in the vascular endothelium, which is the site of endothelial nitric oxide synthase (eNOS) expression [13, 14]. TCEAL2 is a nuclear phosphoprotein that modulates transcription in a promoter context-dependent manner and has been recognized as an important nuclear target for intracellular signal transduction. The function of RNase K is unknown but might be related to an enhanced degradation of the tumor suppressor gene mRNAs, which leads to the development of cancer. The up-regulation of these three genes was observed in more than 60% of the total number of patients.

Figure 4
figure 4

Kaplan-Meier estimates of overall survival stratified by up- and down-regulation for three genes including DDAH2 , RNase K and TCEAL2.

The survival analysis also was performed for chemo-resistance. Following debulking surgery, all patients received platinum-based chemotherapy, considered the standard of care for patients with advanced ovarian cancer. The patients that had either progression during chemotherapy or relapse within six months of treatment were considered chemo-resistant. Among 16 patients, eight were classified as chemo-resistant and the others were categorized as chemo-sensitive. The difference in overall survival between these two groups was significant (p value < 0.05, Figure 5). The shorter survival time of patients with chemo-resistance is consistent with a prior report [15]. To consider chemo-resistance and gene expression simultaneously, in the prediction of overall survival, the Cox multivariate analysis was carried out with the expression information of 38 DEGs and chemo-resistance information from 16 patients. Multivariate analysis demonstrated a significant difference in overall survival between the chemo-resistant and sensitive groups for four DEGs including: NM_001004333.3 (RNase K), NM_032682 (forkhead box transcription factor family, FOXP1), NM_002292.3 (a family of extracellular matrix glycoproteins, LAMB2), and NM_001100167 (murine retrovirus integration site 1 homolog, MRVI1) (p values < 0.05, Figure 6). The RNase K showed significance in both univariate (gene expression) and multivariate (gene expression and chemo-resistance) analysis while DDAH2 and TCEAL2 showed significance only in univariate analysis. This was mainly due to the similar proportion of up- and down-regulation of DDAH2 and TCEAL2 between the chemo-resistant and chemo-sensitive groups. The overall survival of patients with chemo-resistance was significantly decreased with up-regulation of these four genes; while the chemo-sensitive patients had down-regulation of these genes and a good prognosis. FOXP1 has a diverse repertoire of functions ranging from the regulation of B-cell development and monocyte differentiation to the facilitation of cardiac valve and lung development [16, 17]. LAMB2 might be involved in the cell adhesion or motility of prostate cancer cells [18]. The down-regulation of FOXP1 and LAMB2 was observed in more than 60% of all patients while MRVI was up-regulated among 63%.

Figure 5
figure 5

Kaplan-Meier estimates of overall survival stratified by chemo-resistance.

Figure 6
figure 6

The overall survival estimated by the Cox proportional hazards multivariate model including gene expression and chemo-resistance. The up- and down-regulation is represented by red and blue, respectively.

Discussion

Ovarian cancer is the most common cause of death among all gynecological malignancies. The five-year survival rates in patients with ovarian cancer are about 80-90% for stage Ia-Ic, 70-80% for stage IIa-IIc, 30-50% for stages IIIa-IIIc and 13% for stage IV [19]. The high rate of death is due to the fact that most of patients (>60%) present with advanced stage disease (FIGO stages III/IV). Despite an initial response rate of 65%-80% to first-line chemotherapy, most ovarian carcinomas relapse. Acquired resistance to further chemotherapy is generally responsible for treatment failure. Several studies have sought to identify gene expression signatures that correlate with clinical outcome to identify those genes that are associated with survival and relapse and to use as predictive biomarkers for response to chemotherapy [20–22].

There are many types of ovarian cancer. EOC accounts for 85%-90%; half of such cases are serous EOC. As with many cancers thought to be of epithelial origin, it is important to establish an appropriate control for evaluating differential gene expression between "normal" and cancer. Expression profiling studies of ovarian cancer have relied on a variety of sources of normal cells or comparison with tumors, including whole ovary samples (WO), ovarian surface epithelium (OSE), exposed to short-term culture, and immortalized OSE cell lines (IOSE) [23]. Direct comparison of the gene expression profiles generated from OSE brushings, WO samples, short-term cultures of normal OSE (NOSE), and telomerase-immortalized OSE (TIOSE) cell lines revealed that these "normal" samples formed robust, but very distinct groups in hierarchical clustering [24]. These indicate that the selection of normal control to compare epithelial ovarian samples in microarray studies can strongly influence the genes that are identified as differentially expressed. In this study, the total RNA from WO (Stratagene, http://www.stratagene.com) was used as control. WO samples have potential to obscure epithelial pattern due to large amounts of stroma, but they offer the advantages of avoidance of exposure to culture conditions and identification of differential gene expression pattern between tumor and normal tissue [25].

A total of 114 DEGs from patients with serous ovarian cancer stage III were identified using the ACP-based GeneFishingâ„¢ PCR system, which uses primers that anneal specifically to the template and allows only genuine products to be amplified. As the GeneFishingâ„¢ system is based on PCR, it can overcome the difficulty in identifying the genes responsible for a specialized function during a certain biological stage; this is because the gene is expressed at low levels, whereas most mRNA transcripts within a cell are abundantly expressed. Among the 114 DEGs, 42 were identified as up-regulated clones while 72 were down-regulated clones. These DEGs were involved in a variety of biological processes including apoptosis, signal transduction and the immune response. Apoptosis inhibitory processes were associated with genes such as NFKBIA, DDAH2 and IFI6 identified from up-regulated clones; while the immune system associated genes such as IGKC, CD74 and BST2 were found in down-regulated clones. The differential expression of DEGs identified by the Genefishingâ„¢ system showed good agreement with the results of the quantitative real-time PCR.

Cluster analysis based on gene expression profiles identified two groups showing a contrast in the expression pattern. That is, one group including: SSBP1, IFI6 DDT, IFI27, C11orf92, NFKBIA, TNXB, NEAT1 and TFG, was up-regulated in most patients with stage III serous ovarian cancer (Figure 3 and Table 3). These genes might be utilized as potential targets in patients with stage III serous ovarian cancer. IFI6 is involved in apoptosis inhibitory activity while TGF is implicated in up regulation of the I-κ B kinase/NF-κ B cascade. TNXB encodes TNX, a protein of unknown function that is mainly expressed in the peripheral nervous system and muscles. The promotion of tumor invasion and metastasis has been reported in mice deficient in TNX through the activation of the matrix metalloproteinase 2 (MMP2) and MMP9 genes [26]. This indicates that the up-regulation of TNXB, in patients with advanced stage of ovarian cancer, might induce low expression of MMP2 and MMP9. The decrease of MMP2 has been reported in liver metastases in advanced colorectal cancers [27]. Another group consists of LAMB2, XRCC6, MEF2C, RBM5, FOXP1, NUDCP2, LGALS3, TMEM185A, and C1S, which was down-regulated in most patients with stage III serous ovarian cancer (Figure 4 and Table 3). The down-regulation of LAMB2 and MEF2C might be involved with cell adhesion or motility in invasive prostate cancer cells [18] and apoptosis via BCL2 transformation [28], respectively. FOXP1 is a potential therapeutic target in cancer and can be considered either an oncogene and/or a tumor suppressor gene [29]. That is, its over-expression confers a poor prognosis in a number of types of lymphomas while the loss of its expression in breast cancer is associated with a poor outcome. The function of FOXP1 in serous ovarian cancer remains unclear. LGALS3 encoding galectin-3 has been implicated in advanced stage disease [30]. The distribution of galectin has been associated with stage III-V with cellular changes such as dysplasia, cancer cells' nest formation, breakage of the basement membrane, and infiltration of cells into non-native tissue. However, 13 out of the total 16 patients with stage III serous ovarian cancer had down-regulation of LGALS3. This is consistent with the report by van den Brule et al. [31] that showed that galectin-3 expression was decreased in 67% of cases compared to the normal epithelial cells. However, it conflicts with the observations of Lurisci et al. [32] that galectin-3 serum levels in patients with ovarian cancer were significantly elevated. The expression of LGALS3 might be affected by the stage of ovarian cancer.

The overall survival estimated by the Kaplan-Meier method was significantly different between the up- and down-regulated patient cohort with regard to three genes including DDAH2, RNase K and TCEAL2. The up-regulation of these genes was associated with a shorter overall survival (Figure 4). In rapidly growing cells like tumor cells, the activity of RNases is decreased [33], as dictated by the requirement of significant amounts of RNA for protein synthesis. However, high RNase activity has been reported in chronic myeloid leukemia [34] and pancreatic carcinoma [35]. In this study, 70% of the patients with stage III serous ovarian cancer showed up-regulation of RNase K. The function of RNase K in advanced ovarian cancer remains to be clarified. In addition, the overall survival of patients with chemo-resistance was significantly decreased with up-regulation of the genes including: RNase K, FOXP1, LAMB2 and MRVI1 (Figure 6). This might implicate these genes in chemoresistance.

Conclusion

One hundred and fourteen DEGs were identified from 16 patients with stage III serous ovarian carcinoma using the ACP-based RT-PCR technique. Fifteen percent of the total DEGs were associated with apoptosis, the immune response, cell adhesion, and signal pathways. The genes related to apoptosis inhibitory processes tended to be up-regulated while the genes associated with the immune response tended to be down-regulated. The up- and down-regulated genes were identified in most of the patients and might be used as predictive markers in stage III serous ovarian cancer.

References

  1. Chang XH, Zhang L, Yang R, Feng J, Cheng YX, Cheng HY, YE X, Fu TY, Cui H: Screening for genes associated with ovarian cancer prognosis. Chin Med J. 2009, 122 (10): 1167-1172.

    CAS  PubMed  Google Scholar 

  2. Bonome T, Lee JY, Park DC, Radonovich M, Pise-Masison C, Brady J, Gardner GJ, Hao K, Wong WH, Barrett JC, Lu KH, Sood AK, Gershenson DM, Mok SC, Birrer MJ: Expression profiling of serous low malignant potential, low-grade, and high-grade tumors of the ovary. Cancer Res. 2005, 65: 10602-10612. 10.1158/0008-5472.CAN-05-2240.

    Article  CAS  PubMed  Google Scholar 

  3. Welsh JB, Zarrinkar PP, Sapinoso LM, Kern SG, Behling CA, Monk BJ, Lockhart DJ, Burger RA, Hampton GM: Analysis of gene expression profiles in normal and neoplastic ovarian tissue samples identifies candidate molecular markers of epithelial ovarian cancer. Proc Natl Acad Sci USA. 2001, 98: 1176-1181. 10.1073/pnas.98.3.1176.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Schwartz DR, Kardia SL, Shedden KA, Kuick R, Michailidis G, Taylor JM, Misek DE, Wu R, Zhai Y, Darrah DM, Reed H, Ellenson LH, Giordano TJ, Fearon ER, Hanash SM, Cho KR: Gene expression in ovarian cancer reflects both morphology and biological behavior, distinguishing clear cell from other poor-prognosis ovarian carcinomas. Cancer Res. 2002, 62: 4722-4729.

    CAS  PubMed  Google Scholar 

  5. Gilks CB, Vanderhyden BC, Zhu S, van de Rijn M, Longacre TA: Distinction between serous tumors of low malignant potential and serous carcinomas based on global mRNA expression profiling. Gynecol Oncol. 2005, 96: 684-694. 10.1016/j.ygyno.2004.11.039.

    Article  CAS  PubMed  Google Scholar 

  6. Klebanov L, Yakovlev A: How high is the level of technical noise in microarray data?. Biol Direct. 2007, 2: 9-10.1186/1745-6150-2-9.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Wan JS, Sharp SJ, Poirier GM, Wagaman PC, Chambers J, Pyati J, Hom YL, Galindo JE, Huvar A, Peterson PA, Jackson MR, Erlander MG: Cloning differentially expressed mRNAs. Nat Biotechnol. 1996, 14: 1685-1691. 10.1038/nbt1296-1685.

    Article  CAS  PubMed  Google Scholar 

  8. Debouck C: Differential display or differential dismay?. Curr Opin Immunol. 1995, 6: 597-599.

    CAS  Google Scholar 

  9. Liang P, Pardee AB: Recent advances in differential display. Curr Opin Immunol. 1995, 7: 274-280. 10.1016/0952-7915(95)80015-8.

    Article  CAS  PubMed  Google Scholar 

  10. Kim YJ, Kwak CI, Gu YY, Hwang IT, Chun JY: Annealing control primer system for identification of differentially expressed genes on agarose gels. Biotechniques. 2004, 36 (3): 424-434.

    CAS  PubMed  Google Scholar 

  11. Scott M, McCluggage WG, Hillan KJ, Hall PA, Russell SE: Altered patterns of transcription of the septin gene, SEPT9, in ovarian tumorigenesis. Int J Cancer. 2006, 118 (5): 1325-1329. 10.1002/ijc.21486.

    Article  CAS  PubMed  Google Scholar 

  12. He X, Ota T, Liu P, Su C, Chien J, Shridhar V: Downregulation of HtrA1 promotes resistance to anoikis and peritoneal dissemination of ovarian cancer cells. Cancer Res. 2010, 70 (8): 3109-3118. 10.1158/0008-5472.CAN-09-3557.

    Article  CAS  PubMed  Google Scholar 

  13. Achan V, Tran CT, Arrigoni F, Whitley GS, Leiper JM, Vallance P: All trans retinoic acid increases nitric oxide synthesis by endothelial cells. Circ Res. 2002, 90: 764-769. 10.1161/01.RES.0000014450.40853.2B.

    Article  CAS  PubMed  Google Scholar 

  14. Onozato ML, Tojo A, Leiper J, Fujita T, Palm F, Wilcox CS: Expression of DDAH and PRMT isoforms in the diabetic rat kidney; effects of angiotensin II receptor blocker. Dibetes. 2008, 57 (1): 172-180. 10.2337/db06-1772.

    Article  CAS  Google Scholar 

  15. Cheng WF, Huang CY, Chang MC, Hu YH, Chiang YC, Chen YL, Hsieh CY, Chen CA: High mesothelin correlates with chemoresistance and poor survival in epithelial ovarian carcinoma. Br J Cancer. 2009, 100 (7): 1144-1153. 10.1038/sj.bjc.6604964.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Hu H, Wang B, Borde M, Nardone J, Maika S, Allred L, Tucker PW, Rao A: Foxp1 is an essential transcriptional regulator of B cell development. Nat Immunol. 2006, 7 (8): 819-826. 10.1038/ni1358.

    Article  CAS  PubMed  Google Scholar 

  17. Shi C, Zhang X, Chen Z, Sulaiman K, Feinberg MW, Ballantyne CM, Jain MK, Simon DI: Integrin engagement regulates monocyte differentiation through the forkhead transcription factor Foxp1. J Clin Invest. 2004, 114 (3): 408-418.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Ashida S, Nakagawa H, Katagiri T, Furihata M, Iiizumi M, Anazawa Y, Tsunoda T, Takata R, Kasahara K, Miki T, Fujioka T, Shuin T, Nakamura Y: Molecular features of the transition from prostatic intraepithelial neoplasia (PIN) to prostate cancer: genome-wide gene-expression profiles of prostate cancers and PINs. Cancer Res. 2004, 64 (17): 5963-5972. 10.1158/0008-5472.CAN-04-0020.

    Article  CAS  PubMed  Google Scholar 

  19. Oliver RI, van Burden M, van't Veer LJ: The role of gene expression profiling in the clinical management of ovarian cancer. Eur J of Cancer. 2006, 42: 2930-2938. 10.1016/j.ejca.2006.04.013.

    Article  Google Scholar 

  20. Berchuck A, Iversen ES, Lancaster JM, Pittman J, Luo J, Lee P, Murphy S, Dressman HK, Febbo PG, West M, Nevins JR, Marks JR: Patterns of gene expression that characterize long-term survival in advanced stage serous ovarian cancers. Clin Cancer Res. 2005, 11: 3686-3696. 10.1158/1078-0432.CCR-04-2398.

    Article  CAS  PubMed  Google Scholar 

  21. Helleman J, Jansen MP, Span PN, van Staveren IL, Massuger LF, Meijer-van Gelder ME, Sweep FC, Ewing PC, van der Burg ME, Stoter G, Nooter K, Berns EM: Molecular profiling of platinum resistant ovarian cancer. Int J Cancer. 2006, 118: 1963-1971. 10.1002/ijc.21599.

    Article  CAS  PubMed  Google Scholar 

  22. Dressman HK, Berchuck A, Chan G, Zhai J, Bild A, Sayer R, Cragun J, Clarke J, Whitaker RS, Li L, Gray J, Marks J, Ginsburg GS, Potti A, West M, Nevins JR, Lancaster JM: An integrated genomic-based approach to individualized treatment of patients with advanced-stage ovarian cancer. J Clin Oncol. 2007, 25: 517-525. 10.1200/JCO.2006.06.3743.

    Article  CAS  PubMed  Google Scholar 

  23. Farley J, Ozbun L, Birrer M: Genomic analysis of epithelial ovarian cancer. Cell Res. 2008, 18: 538-548. 10.1038/cr.2008.52.

    Article  CAS  PubMed  Google Scholar 

  24. Zorn KK, Jazaeri AA, Awtrey CS, Gardner GJ, Mok SC, Boyd J, Birrer MJ: Choice of normal ovarian control influences determination of differentially expressed genes in ovarian cancer expression profiling studies. Clin Cancer Res. 2003, 9: 4811-4818.

    CAS  PubMed  Google Scholar 

  25. Welsh JB, Zarrinkar PP, Sapinoso LM, Kern SG, Behling CA, Monk BJ, Lockhart DJ, Burger RA, Hampton GM: Analysis of gene expression profiles in normal and neoplastic ovarian tissue samples identifies candidate molecular markers of epithelial ovarian cancer. Proc Natl Acad Sci USA. 2001, 98: 1176-1181. 10.1073/pnas.98.3.1176.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Lévy P, Ripoche H, Laurendeau I, Lazar V, Ortonne N, Parfait B, Leroy K, Wechsler J, Salmon I, Wolkenstein P, Dessen P, Vidaud M, Vidaud D, Bièche I: Microarray-based identification of Tenascin C and Tenascin XB, genes possibly involved in tumorigenesis associated with neurofibromatosis type 1. Clin Cancer Res. 2007, 13 (2): 398-407. 10.1158/1078-0432.CCR-06-0182.

    Article  PubMed  Google Scholar 

  27. Chan CC, Menges M, Orzechowski HD, Orendain N, Pistorius G, Feifel G, Zeitz M, Stallmach A: Increased matrix metalloproteinase 2 concentration and transcript expression in advanced colorectal carcinomas. Int J Colorectal Dis. 2001, 16: 133-140. 10.1007/s003840100287.

    Article  CAS  PubMed  Google Scholar 

  28. Nagel S, Meyer C, Quentmeier H, Kaufmann M, Drexler HG, MacLeod RA: MEF2C is activated by multiple mechanisms in a subset of T-acute lymphoblastic leukemia cell lines. Leukemia. 2008, 22 (3): 600-7. 10.1038/sj.leu.2405067.

    Article  CAS  PubMed  Google Scholar 

  29. Koon HB, Ippolito GC, Banham AH, Tucker PW: FOXP1: a potential therapeutic target in cancer. Expert Opin Ther Targets. 2007, 11 (7): 955-965. 10.1517/14728222.11.7.955.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Balasubramanian K, Vasudevamurthy R, Venkateshaiah SU, Thomas A, Vishweshwara A, Dharmesh SM: Galectin-3 in urine of cancer patients: stage and tissue specificity. J Cancer Res Clin Oncol. 2009, 135: 355-363. 10.1007/s00432-008-0481-4.

    Article  CAS  PubMed  Google Scholar 

  31. van den Brule FA, Berchuck A, Bast RC, Liu FT, Gillet C, Sobel ME, Castronovo V: Differential expression of the 67-kD laminin receptor and 31-kD human laminin-binding protein in human ovarian carcinomas. Eur J Cancer. 1994, 30: 1096-1099. 10.1016/0959-8049(94)90464-2.

    Article  Google Scholar 

  32. Iurisci I, Tinari N, Natoli C, Angelucci D, Cianchetti E, Iacobelli S: Concentrations of galectin-3 in the sera of normal controls and cancer patients. Clin Cancer Res. 2000, 6: 1389-1393.

    CAS  PubMed  Google Scholar 

  33. Peclet C, de Lamirande G, Daoust R: Sensitivity to RNase treatment of ribosomes and rRNA from normal rat liver and Novikoff hepatoma. Br J cancer. 1982, 45: 140-143.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Doran G, Allen-Mersh TG, Reynolds KW: Ribonuclease as a tumor marker for pancreatic cancer carcinoma. J Clin Pathol. 1980, 33: 1212-1213. 10.1136/jcp.33.12.1212.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Naskalski JW, Celinski A: Determination of actual activities of acid and alkaline ribonuclease in human serum and urine. Mater med Pol. 1991, 23: 107-110.

    CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This study was supported by a grant from the National R&D Program for Cancer Control, Ministry for Health, Welfare and Family affairs, Republic of Korea (0820330).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Woong Shick Ahn.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YSK, DHB and WSA designed the study and provided the clinical background. DHB and WSA performed sample annotation and gathered follow-up of the patients. SB performed the experiments. JHD and SB carried out data analysis and JHD wrote the manuscript. All authors contributed to the manuscript and approved it.

Electronic supplementary material

12885_2010_2375_MOESM1_ESM.XLS

Additional file 1: The log2 expression ratios, of cancer to normal samples measured by the quantitative real-time PCR, for the 38 DEGs from 16 patients. (XLS 33 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kim, YS., Hwan Do, J., Bae, S. et al. Identification of differentially expressed genes using an annealing control primer system in stage III serous ovarian carcinoma. BMC Cancer 10, 576 (2010). https://doi.org/10.1186/1471-2407-10-576

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-10-576

Keywords