Gene expression reveals two distinct groups of anal carcinomas with clinical implications

Human papillomavirus (HPV) is a major aetiological agent in anal carcinomas. We here present a study of global gene expression using microarray hybridisation in a collection of anal carcinoma biopsies. Quantitative PCR was used to verify expression of selected genes. All biopsies contained integrated DNA of human papillomavirus subtype 16 (HPV16) and expressed HPV16 E7 mRNA. No other subspecies of HPV were detected in these 13 biopsies as assessed by PCR amplification and DNA sequencing. Unsupervised cluster analysis, based on global mRNA expression, divided the tumour biopsies into two distinct groups. Cluster analysis based on a number of high-risk HPV and/or E2F-regulated genes reproduced this biopsy grouping, suggesting that integrated HPV16 substantially influenced global gene expression in approximately half the biopsies studied. The levels of HPV16 E7 mRNA were significantly different between the two groups, but with considerable overlap. Thus, influence on global gene expression could not be absolutely ascribed to the expression level of HPV16. To investigate whether this distinction in gene expression had prognostic impact, we studied protein expression in an independent cohort of 55 anal carcinomas not included in the microarray study of two differentially expressed candidate genes, minichromosome maintenance complex component 7 (MCM7) and cyclin-dependent kinase inhibitor 2A (CDKN2A or p16). HPV status was assessed by in situ hybridisation. There was a significant association between in situ staining for HPV E7 mRNA and immunostaining for CDKN2A (p16) and MCM7 protein. CDKN2A (p16) mRNA was found significantly differentially expressed between the two tumour groups. However, cluster analysis on genes directly regulated by CDKN2A (p16) could not reproduce this split of biopsies into two groups, suggesting that the transcriptional regulatory activity of CDKN2A in these biopsies is inhibited. Furthermore, protein expression of CDKN2A (p16) could not be associated with survival. MCM7 is directly regulated by E2F and induced by HPV, and its mRNA was found differentially expressed between the two tumour groups. High level of MCM7 protein was found to be associated with both improved relapse-free survival (RFS, P=0.02) and cancer-specific survival (CSS, P=0.03) in anal cancer patients treated with radiation with or without additional chemotherapy.

. The concept that loss of E2 repressor function may be critical for malignant progression is supported by experiments showing that re-expression of E2 in cervical cancer cell lines causes growth suppression (Thierry and Yaniv, 1987). Other factors, including sequence variants of HPV subtypes, seem to be important for the oncogenic potential of the virus (Fujinaga et al, 1994;Ellis et al, 1995;Hecht et al, 1995).
The retinoblastoma protein (pRB) pathway plays a central role in regulating the progression through the G1 phase of the mammalian cell cycle (Sherr, 1996). The core members of this pathway include, in addition to pRB (and its family members p107 and p130), the D-type cyclins that, in association with CDK4 and CDK6, promote proliferation of the cell cycle and the INK4 family of cyclin-dependent kinase inhibitors (including cyclin-dependent kinase inhibitor 2A (CDKN2A) or p16) that specifically binds and inhibits the activity of CDK4 and CDK6. Also, the most studied and the best-understood targets for pRB are members of the E2F transcription factor family (Dyson, 1998;Muller et al, 2001). The transcription factor E2F is a key transcriptional regulator belonging to a family of transcription factors responsible for transcriptional control of genes encoding proteins involved in several biological processes (Cobrinik, 2005). This includes DNA replication, enzymes involved in deoxynucleotide biosynthesis, proteins that assemble to form functional origin of replication complexes, and kinases involved in the activation and initiation of cell-cycle progression (Ishida et al, 2001) as well as transcription of negative regulators of cell-cycle progression such as CDKN2C (p18) (Blais et al, 2002). The high frequency of alterations that have been identified in the pRB pathway in human cancer, a pathway central for the regulation of cell proliferation, suggest that the deregulation of the pRB pathway is an obligatory event in human cancer (Hanahan and Weinberg, 2000).
Given the many HPV subtypes, and the transforming function being dependent on sequence variants and integration events, establishing the oncogenic potential in a given HPV-infected cancer is a complex task. Rather than identifying the exact viral species or variant, we propose that it will be advantageous to investigate the effects of the integrated virus on global gene expression, thus possibly establishing potential markers that can be associated with disease outcome. Although HPV is strongly associated with anal cancer (Daling et al, 2004), it is not known how the presence of high-risk HPV affects global gene expression in anal cancer in vivo. To elucidate potential prognostic or predictive markers for HPV-associated anal cancer, we analysed the mRNA expression profile in 13 HPV16-positive anal cancer biopsies and 4 HPV16-negative normal anal mucosa samples. Furthermore, we tested the associations between protein expression of two isolated target genes, MCM7 and CDKN2A (p16) with HPV E7 mRNA expression and with the clinical outcome in 55 anal cancer patients.

Patient and tissue samples
Microarray expression analysis The biopsies used for gene expression analyses were taken from 13 patients with anal carcinoma, prior to treatment, and are summarised in Table 1. One sample was obtained at the University Hospital of Tromsø (Tromsø, Norway). The rest of the samples were obtained at Haukeland University Hospital (Bergen, Norway). Eleven of these patients were women. The mean age was 66.3 years (range: 45 -88 years). The histology was SCC in all, except for one case with Paget's disease (ca-29), and one had a basaloid type of SCC (ca-36). The TNM stage for each patient is given in Table 1. All patients except ca-20 and ca-47 were treated with radiation therapy combined with concomitant chemotherapy using mitomycin C (10 mg m À2 ) and 5-fluorouracil (1000 mg m À2 , days 1 -4). Patient ca-20 was a 63-year-old woman presenting with T4N2 disease, and patient ca-47 was a 54-year-old man with T3N1 disease and a possible retroperitoneal lymph node metastasis (M1). These two had first two courses of chemotherapy with 5-fluorouracil and cisplatin (80 mg m À2 ) with 3 weeks interval, and then they were given concomitant chemotherapy with cisplatin (60 mg m À2 ) and 5-fluorouracil during the first week of radiation therapy. Primary surgery was not performed in the 13 patients. The mean follow-up was 40 months (range: 25 -58 months for patients alive).
Three patients had a relapse (as indicated in Table 1) and died from anal carcinoma. Patient ca-41 had a distant relapse 52 months from diagnosis and died 3 months later. Patient ca-24 had a local relapse at 11 months and died 1 month later, while patient ca-47 (with possible M1 disease initially) died from distant metastases 31 months after diagnosis. Among the 10 nonrelapsing patients, two have experienced breast cancer during follow-up (patients ca-22 and ca-39).
Fresh tissue samples were taken with a conchotome, directly snap-frozen in liquid nitrogen and stored at À801C. All biopsy specimens were examined by a histopathologist and classified according to WHO criteria (Hermanek, 1992). Clinical data were collected from the hospital files. Cytological smears obtained from the biopsies immediately prior to tissue lysis and RNA extraction were evaluated for the presence of tumour cells.
Tissue microarray Independent samples of formalin-fixed and paraffin-embedded (FFPE) tumour tissues from the time of diagnosis, prior to therapy, from 62 patients with anal carcinoma were selected from a previous patient cohort. A tissue microarray (TMA) was performed as described below. Seven of them were treated with surgery alone. The 55 patients included in the survival analyses were treated with radiation alone or combined with concomitant chemotherapy (5-fluorouracil (1000 mg m À2 days 1 -4) and mitomycin C (10 mg m À2 ) (n ¼ 30); UKCCCR ACTWP, 1996; Bartelink et al, 1997). Among these 55 patients, the median radiation dose was 50 Gy (range: 28 -60 Gy).
Of 46 patients treated with concomitant chemotherapy, 28 received one course and 18 patients received two courses of chemotherapy. Mean follow-up time was 7.2 years.

RNA extraction
Total RNA was extracted using the Qiagen RNeasy minikit (QIAGEN GmbH, Hilden, Germany), as recommended by the manufacturer. The first flow-through fraction was collected for subsequent DNA extraction. Three micrograms of total RNA was subjected to DNase I treatment to remove any contaminating genomic DNA, using TURBO DNaset (Ambion, cat. no. 2238). The RNA quality and quantity were evaluated by Agilent Bioanalyzer and OD260-280, respectively.

DNA extraction
Genomic DNA was extracted from the first flow-through following the RNA purification protocol. Briefly, 350 ml flow-through were mixed with 1 ml linear acrylamide (Ambion, cat. no. 9520) and 500 ml 100% EtOH, vortexed and centrifuged at 13 000 g in a tabletop centrifuge. The precipitate was air-dried, and subsequently purified by Qiagen DNA minikit (cat. no. 13323).

Detection of genomic human papillomavirus DNA
Quantification of HPV16 and 18 DNA was done as described (Lamarcq et al, 2002). All samples were run in triplicate. Estimation of the relative copy number of HPV16 or HPV18 DNA was performed by the comparative CT method relative to the peripheral myelin protein 22 gene (pmp22). In addition, subspecies of HPVs were assessed by PCR amplification with degenerate PCR primers (MY09 and MY11) (Manos et al, 1989) followed by DNA sequencing (data not shown).

Microarray
Gene expression analysis was performed on the ABI 1700 Expression array system (ABI, Foster City, CA, USA) with the Applied Biosystems Chemiluminescent RT-IVT Labelling Kit (cat. no. PN 4339629) and Human Genome microarray (ABI, cat. no. 4337467) using 2 mg high-quality total RNA. Fifteen micrograms of mRNAs were hybridised to the ABI microarray. The Applied Biosystems Human Genome Survey Microarray v. 1.0 provides 31 700 probes for interrogation of 27 868 genes. Each probe is a synthetic oligodeoxynucleotide of 60 nucleotides (nt). The oligonucleotide probes are modified at the 3 0 end by the addition of a C-6 spacer and are covalently attached to the support via a terminal amino group. In almost all cases, probes are placed within 1500 nt of the 3 0 end, where labelling is more robust. The Applied Biosystems 1700 Chemiluminescent Microarray Analyzer is a component of the integrated Applied Biosystems Expression Array System, which uses the power and sensitivity of chemiluminescence to identify and measure gene expression levels.
All procedures were performed according to the manufacturer's recommendations.
Analysis of microarray data J-Express Pro v. 2.7 (Dysvik and Jonassen, 2001) was used for all filtering, clustering and statistical analysis of the microarray experiments. The compiled set consisted of 22 385 individual probes. Euclidean hierarchical clustering with average linkage (weighted pair group method with arithmetic mean, WPGMA) was performed on this data set for the total unsupervised clustering. This unsupervised clustering divided the tumour biopsies into two distinct clusters, group 1 (Gr1) and group 2 (Gr2). Biopsies belonging to Gr1 were labelled in blue and in Gr2 labelled in red and are consistently illustrated so throughout this paper. Additional filtering left 8449 spots remaining. All further analysis was done on this set. To isolate the genes accounting for the observed clustering, we grouped the biopsies according to the result of the unsupervised cluster analysis and performed supervised T-test and significance of microarray (SAM) (Tusher et al, 2001) analysis of the two groups (Gr1 and Gr2), in both cases returning 500 probes with the lowest P-value. The resulting 500 probes from the T-test analysis were sorted according to the biological process classification by the Human protein reference database (HPRD) (http:// www.hprd.org/). Of the 500 probes, 344 returned a unique HPRD reference, the rest being either duplicate probes or hypothetical genes that had no reference in this particular database. The HPRD classification divided the 344 probes into eight distinct biological processes. Euclidean hierarchical clustering with average linkage was performed for each of the eight subgroups.

Tissue microarray
Sixty-two formalin-fixed, paraffin-embedded tumour samples of representative tumour regions from individual patients were used for the preparation of a tumour TMA block. The TMA block was constructed using a Manual Tissue Arrayer 1 (Beecher Instruments, Silver Spring, MD, USA) with a 1.0-mm diameter core biopsy needle. A pathologist reviewed the H&E-stained slides to control for tissue quality and confirm the diagnosis. Following the construction of the array block, 5 mm sections were cut with a microtome. Each slide used for immunohistochemistry (IHC) was de-paraffinised and rehydrated following standard protocols.
For the detection of CDKN2A (p16) protein, target retrieval was performed in 10 mM citrate-buffer pH 6 in a microwave oven. Immunohistochemistry was performed using a CDKN2A (p16)specific antibody (BD Biosciences, Sparks, MD, USA, P/N 554070) diluted 1 : 200 in Antigen Diluent Solution, using a Dako Autostainer and the EnVision DAB detection system (Dako, P/N K5007).
Immunohistochemical analysis of p16 and MCM7 expression was scored independently by two of the authors. The staining intensity of MCM7 was categorised as negative (0), weak (1), moderate (2) or strong (3). An index was made from the product of staining intensity (0 -3) and from the percentage of positive tumour cells (0 -100). This index then was categorised as low (index o140, n ¼ 27) or high (index X140, n ¼ 28). The cutoff value of 140 was the median value of the complete data set, dividing the tumour biopsies essentially into two equal numbered groups. The impact of the protein-staining index on clinical outcome then was assessed on the 55 patients treated with radiation with or without additional chemotherapy (that is the 7 patients treated with surgery alone were left out). The Kaplan -Meier method (Bland and Altman, 1998) was used to estimate relapse-free survival (RFS) and cancer-specific survival (CSS), and the log-rank method was used to compare survival between groups.

In situ hybridisation
The HPV in situ hybridisation (ISH) was performed using reagents and INFORM s HPV III Family 16 DNA Probe (B) (P/N 800-4295) provided by Ventana Medical Systems Inc. (Tucson, AZ, USA). Developed to stain FFPE tissue sections, this probe is designed to detect HPV genotypes 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58 and 66. This probe was used in conjunction with a Discovery automated slide stainer, the RiboMap kit (P/N 760-102) and BlueMap-detection kit (P/N 760-120), using a streptavidin alkaline phosphatase conjugate that gives a blue colour in the chromogenic reaction with NBP (nitroblue terazolium) and BCIP (5-bromo-4chloro-indolyl-phosphate). The 5-mm-thick sections were deparaffinised and rehydrated on the stainer. DNA unmasking was performed with RiboCC buffer (P/N 760-107) and protease 3 (P/N 760-2020) for 20 min. A mix of 200 ml HPV III probe, 75 ml RiboHybe (P/N 760-104) and 25 ml 2 Â SSC was applied on each slide for hybridisation. The probe and target were denatured at 951C for 8 min and incubated for 2 h at 521C for hybridisation. Stringent wash were performed in 2 Â SSC at 721C. A 20-min incubation with a rabbit anti-DNP antibody (P/N 780-4335) was included prior to the 20-min incubation with the Universal Secondary Antibody (P/N 760-4205). The slides were counterstained with ISH Red Counterstain (P/N 780-2186), dehydrated through graded alcohols and xylene, and mounted. HPV-positive cervix tumour was used as a positive control. A negative control, omitting the probe, was included.
For HPV-ISH, nucleolar E7 mRNA was scored as negative, weak, moderate, strong or intense by two independent investigators. For analysis of association between MCM7 protein staining and E7 mRNA, E7 mRNA was recoded as negative or positive, and analysed by w 2 statistics including all 62 patients.
The Regional Ethics Committee Region West approved this study.

Patient characteristics
Histopathological evaluation of the biopsies used for RNA extraction is described in Table 1. Biopsies on the tissue array included 62 patients with histologically confirmed primary anal SCC.

Identification of HPV virus and evaluation of HPV16 viral load
All cancer biopsies included in the microarray study tested positive for HPV16 mRNA by quantitative real-time PCR analysis (Figure 1), while the four normal mucosa biopsies tested negative. No other subspecies of HPV virus were detected in this selection of biopsies when performing PCR amplification with the degenerate primers followed by DNA sequencing. The relative amount of HPV16 E7 mRNA was higher in Gr2 compared with Gr1 (P ¼ 0.05). In two biopsies, ca-29 (Gr1) and ca-24 (Gr2), only small amounts of HPV E7 mRNA were detected. For two of the biopsies (ca-3 in Gr1 and ca-11 in Gr2), we were unable to isolate genomic DNA. For the rest of the biopsies, the relative amount of viral DNA was estimated by comparison to the single-copy gene, pmp22. The estimated copy number of HPV16 genomes compared with host cell genome varied from approximately six copies per cell to one copy per 20 cells (HPV16 genome: host cell genome) (Figure 1).

Microarray analysis
Unsupervised hierarchical clustering on 22 385 genes readily distinguished tumours from the normal anal mucosa, and further identified two main subclasses of anal carcinomas based on distinct patterns of gene expression (Figure 3Aa). Further filtering, eliminating genes not differentially regulated between the samples, reduced the total number of probes to 8449. When performing a SAM analysis (Tusher et al, 2001) Figure 1 Taqman quantitative PCR for Gr1 and Gr2 anal cancer biopsies of the relative amounts of HPV16 E7 DNA relative to the single-copy gene pmp22 (P ¼ 0.12) and mRNA relative to b-actin mRNA (P ¼ 0.05).  normal tissue group with the tumour group, 1296 probes showed a difference in expression level of more than twofold. Comparison of the two tumour groups resulted in 2420 probes with expression levels differing more than twofold. In all, 344 individual genes were identified by performing a T-test between the two groups, followed by annotation according to the HPRD and subsequent division into eight groups depending on biological functions (Figure 2). Euclidean hierarchical clustering was performed for the 344 genes as a whole (Figure 2A), and subsequently for each of the eight subgroups ( Figure 2B). The gene symbol and location in the cluster of a subset of the 344 genes is illustrated in Figure 2B. Most notable, genes encoding proteins involved in DNA replication, DNA repair, DNA binding and key regulators of the cell cycle were upregulated in Gr2, while genes involved in cell -cell adhesion, desmosome formation and apoptosis as well as several protein kinases and phosphatases, members of the Rho-Ras GTPase families and structural proteins were downregulated. Among the 344 genes identified, several were regulated by the transcription factor E2F. Therefore, we determined expression profiles on a selection of 114 genes specifically regulated by E2F (Ishida et al, 2001) and subsequently performed hierarchical clustering. The resulting hierarchical clustering (Figure 3Ad) was highly similar to both the hierarchical clustering of our unsupervised analysis ( Figure 3Aa) and to the hierarchical clustering of genes selected by performing two independent statistical analyses between the two tumour groups (Figure 3Ab and Ac). Similarly, to investigate the possible transcriptional influence of high-risk HPV, we performed hierarchical cluster analysis on genes reported differentially  expressed in several systems related to high-risk HPV infection. This included gene expression profiling of a mouse SCC model (Loercher et al, 2004), effect on overall gene expression associated with integration of HPV16 in cervix SCC (Alazawi et al, 2002), gene expression profiling of primary HPV16-and HPV18-infected early stage cervical cancer (Santin et al, 2005), effect of HPV18 E6 and E7 on global gene expression in an organotypic keratinocytes culture system (Garner-Hamrick et al, 2004), effect of HPV16 E6 and/or E7 on global gene expression in differentiating cervical keratinocytes (Nees et al, 2001) and, finally, effect of HPV31 on global gene expression in normal human keratinocytes (Chang and Laimins, 2000). Again, the resulting hierarchical cluster (Figure 3Ae-Aj) was practically identical compared to our unsupervised analysis ( Figure  3Aa). Figure 3B illustrates the number of genes reported by others and used in this study as supporting information, and the distribution of these genes within our selection.

MCM7
P=0.0081   To test whether this division into two groups was particular to our collection of biopsies, we tested the mRNA expression of a number of genes listed in Figure 2B in a new series of anal carcinoma biopsies from six individual patients. The transcription of these genes in these biopsies also divided them into two groups with three biopsies in each group (data not shown).
Finally, we investigated expression profiling of genes selectively regulated by either CDKN2A (p16) or by E2F as reported by Vernell et al (2003). We found that hierarchical clustering based on CDKN2A (p16)-regulated genes could not reproduce the unsupervised clustering (Figure 3Ca), while genes regulated directly by E2F reproduced the unsupervised cluster analysis (Figure 3Cb).
Quantitative real-time PCR analyses confirmed the results from the microarray analysis for a selection of these genes (Figure 4).

Immunohistochemistry and tissue microarray
Following the results from the microarray analysis, CDKN2A (p16) and MCM7 were selected as a typical gene overexpressed in Gr2 to be investigated in an independent cohort of 55 anal cancer biopsies. CDKN2A (p16) protein staining was not significantly related to survival in Kaplan -Meier analyses (log-rank test, P ¼ 0.70) (data not shown). However, we found that high expression of MCM7 in this cohort of anal carcinomas correlated with favourable clinical outcome (log-rank test, P ¼ 0.017 for probability of RFS and P ¼ 0.011 for probability of CSS), as illustrated in Figure 5Bb.

In situ hybridisation and HPV detection
Out of 62 cases with available representative tumour tissue, E7 mRNA was detected in 34 anal carcinoma biopsies (55%). The intensity of the nucleolar staining was recorded as negative in 28 (45%), weak in 12 (19%), moderate in 8 (13%), strong in 8 (13%) and intense in 6 (10%). There was a significant association between E7 mRNA signal detected by ISH and MCM7 protein staining detected by IHC. Out of 62 anal carcinoma patients, 33.3% of those with a low MCM7 staining intensity (product o140) were HPV positive, as compared with 68.4% of those with a high MCM7 product (P ¼ 0.007, w 2 statistics). Also, there was a significant trend for increasing MCM7 immunostaining with increasing E7 mRNA staining intensity (recorded 0 -4) (P ¼ 0.007 w 2 test for trend). In addition, there was a significant association between presence of HPV mRNA and high CDKN2A (p16) immunostaining (P ¼ 0.001). While a high MCM7 protein staining was significantly associated with an improved CSS in 55 anal carcinoma patients treated with radiation with or without chemotherapy ( Figure 5), presence of HPV mRNA assessed by ISH was not significantly related to survival in Kaplan -Meier analyses (log-rank test P ¼ 0.92).

DISCUSSION
To our knowledge, this is the first description of how high-risk HPV infection influences in vivo global gene expression patterns in anal cancer. However, the effect of high-risk HPV on global gene expression has been studied and characterised in several model systems on different microarray platforms (Chang and Laimins, 2000;Nees et al, 2001;Alazawi et al, 2002;Vernell et al, 2003;Garner-Hamrick et al, 2004;Loercher et al, 2004;Santin et al, 2005). Several biomarkers predicting both HPV status and response to treatment for several cancers have been proposed. MCM7 and CDKN2A (p16) were proposed as biomarkers for HPVpositive head and neck cancer (Strati et al, 2006), while CDKN2A (p16) alone was suggested as a marker in HPV-infected oropharyngeal cancers with favourable prognosis (Weinberger et al, 2006). MCM7 was also proposed as an informative biomarker in cervical cancer (Brake et al, 2003). Both MCM7 and CDKN2A (p16) are among a number of genes regulated by the transcription factor E2F, and should therefore be transcriptionally induced in cells infected by high-risk HPV.
By investigating gene expression profiles in anal cancer biopsies, we identified two major subsets of anal cancers with unique and distinct molecular fingerprints. All cancer biopsies included in the microarray series tested positive for the presence of HPV16 mRNA by quantitative PCR. No other species of HPV was detected by PCR amplification followed by DNA sequencing in these particular biopsies. This does not reflect the distribution of HPV in Norwegian tumour samples reported by others, although HPV16 is reported to be the most common high-risk HPV in Norwegian specimens, accounting for more than 50% of all cases (Kraus et al, 2004;Molden et al, 2005). Owing to the technical limitation of the PCR/sequencing assay, we neither can exclude the presence of other subspecies of HPV within our samples, nor can we exclude the presence of viral variants, sequence disruptions or very low HPV levels that escaped detection by the methods used.
Cyclin-dependent kinase inhibitor 2A (p16) and MCM7 were significantly upregulated in one of the tumour subsets, indicating that this subset (Gr2) could be characterised by E2F-regulated transcription. This distinction could not be ascribed solely to the presence of HPV16, as there were extensive variations in both the relative amount of genomic viral DNA present and the level of viral mRNA expressed. However, the distinction could be explained by a major impact in penetrance of transcriptional regulation specifically controlled by E6 and E7, two of the genes encoded for by the HPV16 genome. Furthermore, one of the two subsets (Gr2) could be characterised by transcriptional regulation of many genes involved in biological processes similar to those controlled by E2F. Among the genes highly upregulated in Gr2, we identified several proteins involved in the initiation of DNA replication, as MCM2, MCM6, MCM7 and MCM8 (Figure 2), and three of the four cyclindependent kinase inhibitors, namely, CDKN2A (p16), CDKN2C (p18) and CDKN2D (p19) ( Figure 4B).
To further characterise differences in the transcriptional response in the two tumour groups, we isolated genes transcriptionally controlled by CDKN2A (p16). However, cluster analysis based on these genes could not reproduce the unsupervised clustering ( Figure 3C). This indicates that, although CDKN2A were highly upregulated in Gr 2, the effect of CDKN2A protein on gene expression does not follow the mRNA expression pattern of the cyclin-dependent kinase inhibitors.
Analysis of the association between MCM7 and CDKN2A (p16) protein expression, as accessed by IHC, and prognosis were performed on a separate cohort of anal carcinomas consisting of biopsies from 55 patients arranged in a tissue array. Human papillomavirus statuses of these biopsies were assessed by ISH. We found a strong correlation between high-risk HPV E7 staining and both CDKN2A (p16) and MCM7 expression, but no association between HPV expression and survival. This could, however, be due to sensitivity limitations of the HPV assay. On the other hand, elevated protein levels of MCM7 were significantly associated with favourable prognosis in patients treated with concomitant radiation and chemotherapy, both regarding RFS and CSS.
Our results suggest that transcriptional control exercised by the cyclin-dependent kinase inhibitor CDKN2A is inhibited in anal carcinomas. We also find that there is no association between CDKN2A protein expression and probability for survival. Furthermore, we find no correlation between HPV expression and survival, probably due to the limitations of the HPV assay, both regarding sensitivity and the fact that it cannot detect variants of HPV. We did, however, find a strong correlation between the MCM7 gene product (E2F regulated, HPV induced and absolute required for DNA replication) and both CSS and RFS. These findings may have clinical implications.