A model for co-expression pattern analysis of genes implicated in angiogenesis and tumour cell invasion in cervical cancer

To date, numerous genes have been identified which are involved in both tumour neovascularisation (angiogenesis) and tumour cell invasion, and most of them are also expressed to some extent under normal physiological conditions. However, little is known about how these genes co-express in these settings. This study was undertaken to quantitate mRNA levels in normal and malignant cervical tissues of nine selected genes (VEGF121, VEGF165, VEGF189, VEGF-C, eIF-4E, b-FGF, TSP-2, MMP-2 and MMP-9) implicated in the above processes using real-time quantitative RT–PCR. In addition, the Spearman's rank correlation was used to determine their co-expression patterns. The transcript levels for the different VEGF-A splice variants (VEGF121, VEGF165, VEGF189) were at least 10-fold higher in the cancer cases, with the highest levels in the primary tumours demonstrating lympho-vascular space involvement. The lymphangiogenic factor VEGF-C and MMP-9 were upregulated 130- and 80-fold respectively in cervical cancers. The highest levels of VEGF-C mRNA were found in the lymph-node positive group. The transcript levels for b-FGF were similar in normal cervical tissue and early-stage cervical cancer, however, higher levels were found in the cervical cancers with advanced stage disease. Comparing gene transcript levels between recurrent and non-recurrent cervical cancer patients revealed significant differences (P=0.038) in transcript levels for the angiogenesis inhibitor TSP-2, with the highest levels in non-recurrent cases. Co-expression pattern analysis in normal cervical tissue revealed highly significant co-expressions (P<0.0001) between TSP-2 and most other genes analysed (VEGF121, VEGF165, VEGF-C, b-FGF and MMP-2). In cervical cancer, TSP-2 appears only to be highly co-expressed with MMP-2 (P<0.0001). In contrast to normal cervical tissue, we found a highly significant co-expression (P<0.0001) between MMP-9 and VEGF189 in cervical cancer. The combined application of real-time quantitative RT–PCR and Spearman's rank correlation identifies gene transcripts which are simultaneously co-expressed. Our results revealed a significant co-expression between the angiogenesis inhibitor TSP-2 and most other genes analysed in normal cervical tissue. In cervical cancer, we found a strong upregulation of VEGF-C and MMP-9 mRNA, with a highly significant co-expression between MMP-9 and VEGF189. British Journal of Cancer (2002) 87, 537–544. doi:10.1038/sj.bjc.6600471 www.bjcancer.com © 2002 Cancer Research UK

During tumour progression, two phases can be recognised with regard to neovascularization: a prevascular and a vascular phase. The transition from one phase to the other is referred to as the 'angiogenic switch' (Folkman, 1971(Folkman, , 2000Carmeliet and Jain, 2000;Hanahan and Folkman, 1996). The current working hypothesis is that the 'switch' involves either the induction of a positive regulator and/or the loss of a negative regulator. In the healthy adult organism, endothelial cell turnover is very low. This is thought to be due to the dominance of endogenous negative regulators, since positive regulators are frequently detected in adult tissues in which there is apparently no angiogenesis. Endogenous negative regulators identified to date include members of the thrombospondin (TSP) family, angiostatin and endostatin (Volpert et al, 1995;O'Reilly et al, 1994O'Reilly et al, , 1997. Among the positive endogenous regulators of tumour angiogenesis are the vascular endothelial growth factor (VEGF) and fibroblast growth factor (FGF) families of cytokines.
Molecular cloning has revealed at least five different isoforms of VEGF which are generated from a single mRNA by alternative splicing, and which have significantly different biochemical features and biological effects Houck et al, 1991;Veik-kola and Alitalo, 1999;Yancopoulos et al, 2000). There is recent evidence in breast cancer patients that levels of VEGF-A mRNA correlate with mRNA levels of the translation initiation factor eIF-4E, implying co-regulation of these genes (Scott et al, 1998). In several cancer types (breast, head and neck and prostate cancer) eIF-4E causes increased translational efficiency of several oncogene transcripts (e.g. cyclin D1 and c-Myc) and growth factors (e.g. basic-FGF and VEGF), leading to overexpression of their products (De Benedetti and Harris, 1999). Recent studies have shown that the growth factor VEGF-C plays a crucial role in the growth of lymphatic vessels (lymphangiogenesis), an effect which is mediated in part through VEGFR-3 (Jeltsch et al, 1997;Lymboussaki et al, 1998).
Many biological processes, including angiogenesis and tumour cell invasion, require extracellular proteolysis for degradation of the extracellular matrix. This allows endothelial and tumour cells to migrate through tissue stroma and subsequently to enter and spread via the blood and/or lymphatic systems. The matrix metalloproteinases (MMPs), including MMP-2 and MMP-9 (gelatinase-A and -B), are amongst the best studied proteases in this respect. Enhanced mRNA and protein levels of both MMP-2 and MMP-9 have been detected in breast, colon and pancreatic cancer (Campo et al, 1992). A recent study in transgenic mice suggests that MMP-9 triggers the angiogenic switch during carcinogenesis (Bergers et al, 2000).
Until recently, studies on angiogenesis in several cancers have largely been limited to analysis of microvascular density (MVD). However, little information is available regarding the co-expression of genes implicated in angiogenesis in different tumour types. Both angiogenesis and lymphangiogenesis are likely to be key processes in disease progression in cervical malignancy, given the well documented pattern of direct spread through parametrial tissues, lymphatic spread to pelvic lymph nodes, and haematogenous dissemination to distant organs. We have recently shown that approximately 50% of early-stage cervical cancers shed tumour cells to the pelvic lymph nodes (Van Trappen et al, 2001). The lymphangiogenic factor VEGF-C may be a key element in this process.
Recent technical advances have provided us with opportunities to assess levels of gene transcription and to identify patterns of co-expression. Although microarray technology can measure the expression of thousands of genes simultaneously, it is a semi-quantitative technique (Schena et al, 1995;Eisen et al, 1998;Duggan et al, 1999;Lander, 1999;Liotta and Petricoin, 2000). As the genes involved in tumour angiogenesis are also expressed in normal tissues, a real understanding of gene expression in carcinogenesis requires an accurate and fully quantitative technique. Real-time quantitative RT -PCR fulfils this requirement and has recently been used to accurately quantitate transcription levels of a variety of genes (Heid et al, 1996).
The aim in our study was to identify which genes are coexpressed in tumours in which angiogenesis and/or lymphangiogenesis are prominent, and to determine whether a different coexpression pattern exists between normal and malignant cervical tissue. Transcripts of the following nine genes were analysed: VEGF 121 , VEGF 165 , VEGF 189 , VEGF-C, eIF-4E, basic-FGF (b-FGF), TSP-2, MMP-2 and MMP-9.

Patients and samples
The patient population consisted of 52 women of whom 34 had malignant and 18 normal cervical tissue. Ethical approval was obtained from the Local Research Ethics Committee of the East London and City Health Authority. During the period between January 1998 to March 2000, radical hysterectomy and pelvic lymphadenectomy were carried out in 29 cervical cancer patients with early-stage disease confined to the uterus. Five patients were treated with primary radiotherapy for advanced stage disease (stage IIB -III) during the same time period, and cervical biopsies were taken before radiotherapy. None of the cancer patients had undergone any previous treatment to the cervix. Tissue specimens were taken from the periphery of the tumour, as it is known that neovascular hot spots occur predominantly at tumour margins (Tokumo et al, 1998). One portion was snap frozen for RT -PCR analysis while a neighbouring portion was fixed in 10% formaldehyde for histological examination to determine whether the biopsy specimen contained 480% tumour tissue. The histological cell types of the 34 cervical cancers were assigned according to the World Health Organisation (WHO) classification: 23 were classified as squamous cell carcinoma, 8 as adenocarcinoma, 2 as clear cell carcinoma and 1 as carcinoid. Clinical staging was determined using the International Federation of Gynaecology and Obstetrics (FIGO) classification: 6 were stage IA2, 15 were staged as stage IB1, 6 were stage IB2, 2 were stage IIA, 3 were stage IIB, 2 were stage III (Table 1: patient characteristics).
As controls, 18 biopsy specimens were obtained from normal cervices from patients who underwent a simple hysterectomy for non-malignant lesions in the corpus of the uterus. None of these patients had conditions, such as endometriosis, in which angiogenesis may be an important component. Light microscopic examination confirmed the normal histology.

RNA extraction
Total RNA was extracted from all samples using the acid-guanidinium isothiocyanate-phenol-chloroform method as described previously (Chomczynski and Sacchi, 1987). Disintegration of the tissue was achieved by vortexing and by using a Dounce homogeniser. The RNA samples were stored in 50 ml diethyl pyrocarbonate (DEPC) treated distilled water at 7808C. RNAs were further purified using an RNeasy Mini Protocol for RNA Clean up (Qiagen Ltd., Crawley, UK) and quantitated in a Genequant spectrophotometer (Pharmacia, St Albans, UK).

Primers and probes
Primers and probes were designed using Primer Express software (PE-ABI). The probes were designed specifically to span an intron in order to avoid potential amplification of small amounts of contaminating DNA in the analysed samples. The probes were labelled at the 5' end with 6-carboxy fluorescein (FAM) and at the 3' end with 6-carboxy-tetramethyl rhodamine (TAMRA). Primers and probes were stored at 7208C until use. Table 2

Single-step real-time quantitative RT -PCR
Real-time quantitative RT -PCR analysis was performed using the Taqman PCR Core Reagent Kit (PE Biosystems). Prior to reverse transcription, the RNA template was heated for 2 min at 508C in the presence of 0.01 U l 71 uracil N-glycosylase. To prevent carry-over of contaminating DNA, the reaction was carried out in the presence of dUTP. Total RNA (100 ng) was reverse transcribed in a 25ml reaction mixture at 608C for 30 min. After 5 min denaturation at 928C, PCR was carried out for 40 cycles with denaturation at 928C for 20 s and extension at 628C for 60 s in the presence of the fluorescent oligonucleotides indicated above. RT -PCR analyses were carried out at least twice for each gene transcript to determine consistency of results, and GAPDH mRNA was assessed as internal control for RNA quality. Reactions were recorded and analysed using the ABI 7700 Prism Sequence detection system (Perkin-Elmer Applied Biosystems, Warrington, UK). Accurate quantitation was achieved through the generation of standard curves by serially diluting a known amount of cDNA from corresponding synthetic RNA and performing Taqman PCR on the series alongside patient samples. All gene transcript levels are expressed per mg of RNA.

Preparation of cDNA
Two mg of DNase treated RNA was mixed with 2 ml (1 mg) of random hexamers in DEPC water, in a total volume of 14 ml. The mixture was heated at 758C for 15 min and cooled on ice for 5 min. Subsequently, 5 ml of M-MLV RT reaction buffer, 5 ml of nucleotide pool and 1 ml of M-MLV RT enzyme was added. This mixture was left at room temperature for 10 min and incubated at 408C for 50 min. Seventy five ml of DEPC water was added to make up a total volume of 100 ml.

Real-time PCR for 18S rRNA
The housekeeping gene 18S rRNA was evaluated in all samples as internal control. In a total volume of 25 ml the following was added as reaction mixture: 12.5 ml of TaqMan Universal Master Mix (Applied Biosystems), 2.5 ul of human 18S rRNA Assay Reagents (Pre-developed; Applied Biosystems), 1 ul (100 ng) of cDNA and 9 ml of DEPC water. After 2 min at 508C, denaturation was performed at 958C for 10 min; PCR was carried out for 40 cycles with denaturation at 958C for 15 s and annealing at 608C for 1 min. The threshold cycle (Ct) for each sample was determined.

Statistical Analysis
The Mann -Whitney U-test was used to compare the means of the gene expressions between normal and malignant tissues (SPSS software, Version 9.0 for Windows, SPSS, Inc, Chicago, IL, USA). We have used the Spearman's rank correlation to compare the coexpression patterns between the normal and the cancer group. A correlation with a P-value of 50.05 was considered to be significant. In addition, the z-test was performed to compare the Spearman correlation matrices between normal and malignant (MedCalc, Version 6.0 for Windows, 1998). The two compared correlation coefficients (r) are transformed into two z values by this formula: z ¼ 0:5 Â ½lnð1 þ rÞ À lnð1 À rÞ which almost follows a normal distribution of mean m, and variance 1/(n-3), with n the number of individuals in the sample. Table 2 Primers/probe sequences used in the real-time quantitative RT -PCR assays to analyse transcription levels of nine genes involved in angiogenesis. The housekeeping gene GAPDH was also included. Then, for the statistical value D/SD D a P-value is given:

Transcript levels in normal and malignant cervical tissue
Real-time quantitative RT -PCR analysis was performed on mRNA transcripts of nine genes (VEGF 121 , VEGF 165 , VEGF 189 , VEGF-C, eIF-4E, b-FGF, TSP-2, MMP-2 and MMP-9) involved in angiogenesis and/or lymphangiogenesis as well as tumour cell invasion. Experiments were repeated at least twice for each gene. GAPDH was initially considered as an internal control. However, considerable variation in mRNA levels were found for this housekeeping gene. In contrast, all samples showed a threshold cycle between 17 and 18 for 18S rRNA, demonstrating its usefulness as internal control for quality of RNA. Figures 1A, B, C show the transcript levels in normal and malignant cervical tissues for the nine genes analysed. All tumour and normal samples expressed the VEGF-A isoforms in order of decreasing abundance, 16541214189 ( Figure 1A). In normal cervical tissue, a greater distribution (95% confidence interval (CI)) was observed for the mRNA copy numbers of the different VEGF-A splice variants compared to cervical cancer. This may be due in part to the influence of the hormonal cycle. The levels of the different VEGF-A splice variants were significantly higher in malignant compared to normal cervical tissues. There was an 11fold increase in the mean copy number of VEGF 165 , a 23-fold increase in VEGF 121 and a 13-fold increase in VEGF 189 in malignant compared to normal tissues. The highest levels of VEGF 121 , VEGF 165 and VEGF 189 mRNA were found in cervical tumours with lympho-vascular space involvement in the primary tumour, which is known to be a poor prognostic factor (lympho-vascular space includes all small vessels, capillaries as well as lymphatic vessels, present in the primary tumour and the surrounding tissue). All cervical cancers expressed high levels of VEGF-C, and the levels were more than 130-fold higher than in normal cervical tissues ( Figure 1B). The highest levels of VEGF-C mRNA were found in the lymph node-positive group. Of the 18 normal cervical tissues analysed, only 10 expressed 4100 copies of VEGF-C mRNA/mg total RNA, suggesting that VEGF-C is not implicated in normal physiological processes in cervical tissue. Levels of b-FGF mRNA were similar in normal cervical tissue and early-stage cervical cancers. However, high levels of b-FGF mRNA were found in the five cervical cancers with advanced stage disease. The levels of TSP-2 and of eIF-4E mRNA were increased four-and seven-fold respectively in cervical cancer. No differences in the level of MMP-2 mRNA were found between normal and malignant cervical tissue. However, an approximately 80-fold increase in MMP-9 mRNA was observed in cervical cancers ( Figure 1C).
Using the Mann -Whitney U-test the differences in the means between the two groups (normal vs cancer) are highly significant (P-values 40.006), except for the genes b-FGF and MMP-2 (respective P-values: 0.336, 0.584) ( Figure 1A, B, C).

Gene co-expression patterns in normal and malignant cervical tissue
Tables 3A and B show the Spearman correlation matrices of gene co-expressions (with rho and P-values) in normal and malignant cervical tissue. In normal cervical tissue, a strong co-expression is observed between the angiogenesis inhibitor TSP-2 and the other genes analysed. Furthermore, MMP-2 is also strongly co-expressed with the other genes analysed in normal cervical tissue, except for MMP-9. In cervical cancer, TSP-2 seems to be strongly co-expressed only with MMP-2. In both normal and malignant cervical tissue, significant co-expressions were found between the different VEGF-A splice variants. In normal cervical tissue, VEGF-C is strongly co-expressed with all other genes analysed, except for MMP-9 and eIF4E. However, in cervical cancer, VEGF-C is strongly co-expressed with MMP-9, VEGF 189 and b-FGF. In addition, MMP-9 is strongly co-expressed with VEGF 189 in cervical cancer but not in normal cervical tissue.

Molecular and Cellular Pathology
To compare the co-expressions between the two groups (normal vs malignant), the z-statistic of Fischer was used on the Spearman correlation matrices (Table 4). Of the 36 correlation coefficients, 18 (50%) showed significant differences in favour of a higher correlation in the normal group, except for MMP-9/VEGF 189 .
Although the sample size was smaller for the normal group and the distribution (95% confidence interval) of quantitative gene expression levels was broader, we observed a highly consistent co-expression pattern in the normal cases. This could reflect the housekeeping mode in normal cases.

Differences in gene (co-)expressions between recurrent and non-recurrent cervical cancers
Six patients (18%) with primary cervical cancer developed recurrent disease during a mean follow-up of 27 months. Three patients were initially diagnosed with early-stage disease, the other three patients had advanced disease at initial presentation. No significant differences in co-expression patterns were found between recurrent and non-recurrent cervical cancers. However, a significant difference was found in the mean of expression levels for TSP-2, between the recurrent and the non-recurrent group (P-value=0.038; Mann -Whitney U-test). The mean value of TSP-2 in non-recurrent cases was 4.9610 5 copies mg 71 RNA (s.d.: +/77.7610 5 ), and in recurrent cases 1.8610 5 copies/mg RNA (SD: +/73610 5 ). This suggests that TSP-2 might play a crucial role in progression of disease in cervical cancer. In contrast to normal cervical tissue, TSP-2 seems to be only strongly co-expressed with MMP-2 in malignant cervical tissue. The main difference between normal and malignant cervical tissue is the strong co-expression between MMP-9 and VEGF 189 in cervical cancer. P-values in bold represent highly significant co-expressions (P50.0001).

Molecular and Cellular Pathology
Co-expression patterns of genes in cervical cancer PO Van Trappen et al

DISCUSSION
We show here that combined application of real-time quantitative RT -PCR with correlation matrices of gene expression data allows us to identify differences in gene co-expression patterns between normal and malignant cervical tissue. In vivo and in vitro studies using semi-quantitative techniques have recently shown that tumours exhibiting high expression levels of vascular endothelial growth factor (VEGF) have a high metastatic potential (Potgens et al, 1996;Ohta et al, 1997). In vivo studies in cervical cancer have demonstrated a significant correlation between MVD and VEGF mRNA expression, with the highest levels in stage I and stage IV disease (Fujimoto et al, 1999;Kodama et al, 1999). VEGF 121 and VEGF 165 were found to be the two dominant VEGF-A subtypes with similar expression levels in normal cervix as well as in cervical cancer. VEGF 189 was barely detectable. In our study, using a fully quantitative RT -PCR technique, mRNA levels of the different VEGF-A splice variants were significantly higher in malignant compared to normal cervical tissues, with a more than 10-fold increase in VEGF 121 , VEGF 165 , and VEGF 189 . Interestingly, the highest levels of VEGF 121 , VEGF 165 and VEGF 189 mRNA were found in cervical tumours histologically characterised by lympho-vascular space involvement, which is known to be a poor prognostic factor. In normal cervical tissue, a greater distribution was observed for the mRNA copy numbers of the different VEGF-A splice variants compared to cervical cancer. This may be due in part to the influence of the hormonal cycle, however, the co-expression patterns were stronger and more consistent than in cervical cancer.
We have shown that all tumour and normal samples expressed the VEGF-A isoforms, 16541214189, in order of decreasing abundance. This order differs from normal and malignant breast tissue in which the VEGF 121 isoform is most abundant (Scott et al, 1998). There is evidence that b-FGF mRNA levels are higher in advanced cervical cancers, regardless of histological type, suggesting that it may predominantly play a role in the later molecular events which regulate angiogenesis in this particular tumour (Fujimoto et al, 1997). Furthermore, it has been demonstrated that VEGF is essential for the initial but not continued growth of human breast carcinoma cells in vivo, and that other angiogenic factors, such as b-FGF, can substitute for VEGF during disease progression (Yoshiji et al, 1997). This is consistent with our findings that early-stage cervical cancers had no significant increase in b-FGF mRNA, whereas a significant increase was observed in advanced stage cancers. Elevated levels of eIF-4E have been associated with increased growth rates and cell transformation in HeLa cells (De Benedetti and Rhoads, 1990). It has been suggested that eIF-4E may regulate expression of genes including angiogenic factors such as VEGF and b-FGF as well as subsequent tumour growth (Kevil et al, 1996). In our study, there was a seven-fold increase in transcript level of eIF-4E in cervical cancer.
Clinical and pathological observations have revealed that lymph node involvement is one of the earliest features of metastatic disease, and this certainly holds true for cervical cancer. Bloodborne tumour cells seed out in the lymphatic hilum, and with an increasing number of tumour cells in the blood, the frequency of nodal involvement rises. This occurs particularly in patients with large volume disease and increased numbers of microvessels in the primary tumour. A second route is via lymphatic vessels, and this pathway may be of particular importance in early microdissemination to lymph nodes (Ohta et al, 2000;reviewed by Van Trappen and Pepper, 2002). It has recently been demonstrated that primary solid tumours expressing vascular endothelial growth factors -C and -D (VEGF-C and -D), induce lymphangiogenesis de novo, thereby providing a direct conduit for tumour cell dissemination to lymph nodes (Mandriota et al, 2001;Skobe et al, 2001;Stacker et al, 2001). In prostatic carcinoma, the expression of VEGF-C mRNA in the primary tumour is significantly higher in lymph node-positive patients compared to lymph node-negative patients (Tsurusaki et al, 1999).
Our results have revealed that VEGF-C transcript levels are more than 130-fold higher in cervical cancer compared to normal cervical tissues. Using real-time quantitative RT -PCR we found the highest VEGF-C mRNA levels in the lymph node-positive cervical cancers.
Although VEGF-C is the principal lymphangiogenic growth factor described to date, it also has angiogenic properties under certain circumstances, including tumorigenesis (Olofsson et al, 1999). This implies that VEGF-C may serve to increase both tumour vascularisation and lymphangiogenesis, thereby contributing to the increase in tumour growth and the formation of metastases.
Proteolytically-mediated extracellular matrix (ECM) degradation is crucial for angiogenesis as well as for the detachment of malig- Table 4 Comparison of the co-expression patterns between the two groups (normal vs malignant) using the z-statistic of Fischer on the Spearman correlation matrices. Of the 36 correlation coefficients, 18 (50%) showed significant differences in favour of a higher correlation in the normal group. the P-values of the z-statistic with one star show highly significant differences between the two groups, with a significant level of 0.1% (P-value 50.001). These correlation coefficients are positive and only highly significant in the normal group. The two stars represent P-value with strong significant differences between the two correlations, with a significant level of 1% (P-value50.01). In all cases, correlation coefficients are positive but higher in the normal group. The P-values of the z-statistic with three stars show significant differences between the two groups with a significant level of 5% (P-value 50.05). Almost all correlations are positive and higher in the normal group, except for MMP9-VEGF189.

Molecular and Cellular Pathology
Co-expression patterns of genes in cervical cancer PO Van Trappen et al nant cells from the primary tumour and migration through the surrounding stroma (Pepper, 2001). Matrix metalloproteinases (MMPs) have been heavily implicated in these processes, and are indispensable for the degradation of native triple-helical collagen, a major component of ECMs. Two of the best studied MMP's in angiogenesis and cancer are MMP-2 and MMP-9 (gelatinase-A and -B), which amongst other molecules degrade collagen IV, one of the major components of the basement membrane. MMP-2 and MMP-9 play an important role in triggering the angiogenic switch in null mice (reviewed by Pepper, 2001). In addition, enhanced mRNA and protein expression of both enzymes has been reported in breast, colon and pancreatic cancer (Campo et al, 1992). In our study, no differences in MMP-2 were found between normal and malignant cervical tissue at the mRNA level. However, an approximately 80-fold increase was observed for MMP-9 mRNA in cervical cancers. The fact that we observed no increase in MMP-2 and a marked increase in MMP-9, provides no information about the status of MMP activation, which from a mechanistic point of view may be more relevant. Nonetheless, our findings have revealed a clear relationship between increased MMP-9 mRNA and cervical cancer, increasing the probability that MMP-9 plays a key role in the progression of early-stage cervical cancer.
In normal cervical tissue, gene co-expression pattern analysis revealed a strong co-expression between TSP-2, MMP-2 and the other genes analysed, except for MMP-9. In cervical cancer, TSP-2 appears only to be co-expressed with MMP-2. The lymphangiogenic factor VEGF-C is strongly co-expressed with MMP-9, VEGF 189 and b-FGF in cervical cancer. In addition, MMP-9 is strongly co-expressed with VEGF 189 in cervical cancer but not in normal cervical tissue. Comparison of the co-expression patterns between the two groups revealed that almost all correlations are positive and stronger in the normal group, except for the coexpression MMP-9/VEGF 189 . Recently, in vivo studies have shown that colon cancer metastasis results from an alteration in the balance between the angiogenesis inhibitor thrombospondin-2 and the angiogenic growth factor VEGF 189 . We found a significant difference in the mean of expression levels for thrombospondin-2 between recurrent and non-recurrent cancer cases, with the lowest levels in the recurrent cases. VEGF 189 has been shown to correlate with poor prognosis in colon and lung cancer Tokunaga et al, 1998), and VEGF-C may facilitate lymphatic spread (Mandriota et al, 2001;Skobe et al, 2001;Stacker et al, 2001). In our study, both VEGF 189 and VEGF-C appear to be strongly co-expressed with MMP-9 in cervical cancer but not in normal cervical tissue.
In this report, we described the combined application of fully quantitative gene expression data with correlation matrices. This approach revealed different co-expression patterns between normal and malignant cervical tissue. This combined molecular-statistical model could be applied to several other cancers or other sets of genes.