Skip to main content

Estimation of groin recurrence risk in patients with squamous cell vulvar carcinoma by the assessment of marker gene expression in the lymph nodes



Regional lymph node (LN) status is a well-known prognostic factor for vulvar carcinoma (VC) patients. Although the reliable LN assessment in VC is crucial, it presents significant diagnostic problems. We aimed to identify specific mRNA markers of VC dissemination in the LN and to address the feasibility of predicting the risk of nodal recurrence by the patterns of gene expression.


Sentinel and inguinal LN samples from 20 patients who had undergone surgery for stage T1-3, N0-2, M0 primary vulvar squamous cell carcinoma were analyzed. Gene expression profiles were assessed in four metastatic [LN(+)] and four histologically negative [LN(−)] lymph node samples obtained from four VC patients, by the Affymetrix U133 Plus 2.0 gene expression microarrays. Of the set of genes of the highest expression in the metastatic LNs compared to LN(−), seven candidate marker genes were selected: PERP, S100A8, FABP5, SFN, CA12, JUP and CSTA, and the expression levels of these genes were further analyzed by the real-time reverse transcription polymerase chain reaction (qRT-PCR) in 71 LN samples.


All of the seven genes in question were significantly increased in LN(+) compared to LN(−) samples. In the initial validation of the seven putative markers of metastatic LN, the Cox proportional hazard model pointed to SFN, CA12 and JUP expression to significantly relate to the time to groin recurrence in VC patients.


Our findings first provided evidence that SFN, CA12 and JUP have a potential of marker genes for the prediction of the groin recurrence LN in VC patients.

Peer Review reports


Vulvar carcinoma (VC) is a rare genital malignancy with age-standardized incidence rates (ASR) ranging worldwide between 0.5 and 1.5 per 100,000 [1]. The current ASR in Poland is 1.0, and 436 new cases of VC were diagnosed and 232 VC-related deaths were recorded in 2009 (age-standardized mortality rate 0.5) [2]. Squamous cell carcinoma is the most common histological type of VC, representing approximately 90% of lesions. VC spreads primarily by the local expansion and via the lymphatic system.

Lymph node (LN) status is the most important prognostic factor in VC patients [3]. Up to 24% of patients with LNs clinically considered to be normal have metastases, while more than 20% of patients with enlarged nodes in clinical examination are found metastases-free in the histopathologic assessment [4]. Therefore, in many centres the VC staging comprises an ultrasound-guided fine-needle biopsy of the suspicious groin LNs and cytologic evaluation of the regional LN status. In early stage VC patients, the probability of a positive inguinal LN finding is only 11-25% [5]. This means that if all women with early VC were assigned to lymphadenectomy, approximately 80% of them would undergo overtreatment. At present, the concept of sentinel LN biopsy in early stage VC is becoming increasingly accepted while the problem of optimal treatment of patients with metastatic regional LNs is still being discussed. Another issue is the assessment of the risk of metastases in non-sentinel LNs after removal of all positive sentinel LNs.

Finding a method of an accurate and sensitive pre- or perioperative determination of pathologic status of groin nodes would be of a great value to aid therapeutic decisions in early stages of VC. Time-consuming, routine pathological evaluation of hematoxylin- and eosin-stained (H&E) LN sections lacks sensitivity. Therefore, a real-time reverse-transcription PCR (qRT-PCR) is being used in various tumours to determine the presence of metastatic cells in LNs [6, 7]. However, before the qRT-PCR-based molecular staging of VC can be applied, clinically relevant markers [8] need to be identified. In this study, to identify specific mRNA markers of VC dissemination into LNs, expression microarrays were employed.


The material was obtained from 20 patients treated for vulvar carcinoma in the Maria Sklodowska-Curie Memorial Cancer Centre and Institute of Oncology in Warsaw between March 2003 and September 2006. Patients with microscopically confirmed vulvar squamous cell carcinoma in clinical stage T1-3, N0-2, M0 and with no prior treatment for this or any other malignancy were enrolled. The selected characteristics, including the TNM categories according to the AJCC (American Joint Committee on Cancer TNM staging system), of VC patients are presented in Table 1. Surgery and LN identification were performed as described previously [9]. The median LN count was 5 nodes in the right groin, range 2–13, while in the left groin the median LN count was 5.5 nodes, range 0–15. The study was approved by the Independent Ethics Committee of the Maria Sklodowska-Curie Memorial Cancer Centre and Institute of Oncology in Warsaw and all patients gave their informed consent. LN specimens were divided into two parts: one part was examined histologically as described previously [10], the other was frozen in liquid nitrogen immediately after collection and stored at −70°C until RNA isolation. Median follow-up time, determined from the date of surgery to the date of death or the date of the last interview, was 1.2 years (range 0.43 – 6.29).

Table 1 Individual patient data

Microarray analysis of gene expression profiles was performed in four pairs of LN samples from four VC patients. One metastatic LN [LN(+)] and one histologically negative [LN(−)], according to a routine pathological examination with H&E, was examined per each of the four patients. These patients had human papillomavirus (HPV)-negative tumours, as previously determined using the Linear Array HPV Detection Kit and Linear Array HPV Genotyping Test (Roche Molecular Systems, Inc) [10].

A total of 71 right and left LNs (sentinel and inguinal specimens), from 20 VC patients were included in the real-time RT-PCR analysis of expression of selected genes. In a routine pathological examination, 22 LNs were evaluated as LN(+) and the remaining 49 were histologically negative [LN(−)].


RNA isolation

Total RNA was isolated from approximately 200 mg of pulverised (with the Microdismembrator II, B Braun Biotech International) LN samples using Nucleospin RNA L kit (Macherey-Nagel), according to the manufacturer’s protocol. RNA quality was assessed using the Agilent 2100 Bioanalyzer and RNA 6000 Nano Chip Kit (Agilent Technologies).

Microarray procedure

cDNA synthesis was carried out from 5 μg of RNA with the One-Cycle cDNA Synthesis Kit (Affymetrix). After purification (GeneChip Ample Cleanup Module), 7 μl of double-stranded cDNA were used for biotinylated cRNA synthesis with IVT Labeling Kit (Affymetrix). The quantity and quality of the obtained cRNA were assessed using the Agilent 2100 Bioanalyser and RNA 6000 Nano Chip Kit (Agilent Technologies). Labeled cRNA was purified using a GeneChip Sample Cleanup Module, fragmented and hybridized with the Affymetrix GeneChip Human Genome U133 Plus 2.0 arrays. Washing, staining with streptavidin-phycoerythrin conjugate and scanning of the arrays in the Affymetrix GeneChip 3000 scanner were performed as recommended by the Affymetrix Gene Expression Analysis Technical Manual.

Microarray data analysis

All the arrays were normalized by the GCRMA algorithm using a Bioconductor [11] (version 2.8.1) package gcrma version 2.14.1. The expression levels were log2 transformed. We filtered out all the probe sets with signal level below 7.5 in at least 8 samples. Genes differentially expressed in involved [LN(+)] and uninvolved [LN(−)] lymph nodes were searched by a random-variance t-test, with the statistical significance threshold set at P<0.01. To reduce false-differential gene expression, a Benjamini-Hochberg [12] multiple testing correction was applied. The full dataset has been deposited in the Gene Expression Omnibus repository (accession no GSE28442).

The probe sets with differential expression in two samples - one LN(+) and one LN(−) - obtained from patient No 15 were chosen to be analyzed first with pumaDE function from PUMA package version 1.8.1 [13]. Differentially expressed probe sets were annotated with Gene Ontology (GO) terms (GO.db version 2.2.5) using the Bioconductor [11] packages GOstats (version 2.8.0) and package annotate (version 1.20.1). The significance of differential representation of GO terms between the specified lists of probe sets was determined by the hypergeometric test implemented in GOstats (version 2.8.0). P values returned by GOstats were corrected for testing of multiple hypotheses with the Benjamini–Hochberg method implemented in an R environment (version 2.8.1, The R Foundation for Statistical Computing; Adjusted P values of less than 0.05 were considered significant. The same procedure was applied to the three additional pairs of LN(+) and LN(−) samples from three patients (No. 21, 46 and 61).

This analysis has provided a set of candidate marker genes, which were further investigated to define their tumour specificity using Genevestigator V3, a web-based microarray database and analysis system [14].

Real time RT-PCR

RNA was extracted from 71 LN samples and 1 μg from each sample was reverse-transcribed by using the RT2 First Strand Kit (C-03) from SA Biosciences. Custom PCR arrays (SA Biosciences) were used to simultaneously examine the mRNA levels of seven genes of interest according to the manufacturer's protocol.

The arrays also included primers for two housekeeping genes and three internal controls. Quantitative real-time PCR analysis was performed with the RT2 Real-Time PCR Master Mix (SA Biosciences) in the 7500 Fast Real-Time PCR System (Applied Biosystems), according to the manufacturer's instructions (SA Biosciences). The collected data were analyzed using threshold-cycle (Ct) values for the genes with the SDS 2.1 software (Applied Biosystems). Normalization was performed based on the mean values of two housekeeping genes, HMBS and HSP90AB1, and the relative amounts of RNA for each gene were calculated by the 2−ΔCT method using DataAssist™ Software (Applied Biosystems). Expression levels for each gene were visualized using GraphPadPrism (La Jolla, CA, USA).

Statistical analysis

The significance of difference between the selected gene expression level in each of LN(+) and LN(−) samples was assessed by using Wilcoxon test. P < 0.01 was considered significant.

The Cox proportional-hazard model with Bonferroni correction for multiple hypotheses testing was applied to estimate the effect of gene expression levels in LN samples on the time to groin recurrence (TTR). The highest expression value of a given gene of all obtained in different LN samples of the same patient was chosen to be included in the statistical analysis. TTR was calculated from the date of primary surgery to the date of groin recurrence. MedCalc (Mariakerke, Belgium) software was used to generate Kaplan-Meier curves and to compare TTR using two sided log-rank analysis. The median of maximum gene expression levels (again, in one of the nodes examined per patient) in LNs from all the VC patients enrolled in the study was arbitrarily chosen as a cut-off value for patient stratification. P < 0.05 was considered significant.


Identification of expressed genes associated with LN metastasis

The microarray data from two LN samples, one LN(+) and one LN(−), from one patient (case No. 15) with a rapid disease progression, were analyzed using pumaDE function from the PUMA package 1.8.1 [13]. The filtered dataset comprised 907 probe-sets with P value below 0.05 (five percent Benjamini-Hochberg false discovery rate) differentially expressed in these two LN samples. This list included only up-regulated genes in LN(+) samples as no genes were found to be reduced in LNs(+) compared to LNs(−).

To evaluate which biological processes are represented within the obtained probe set list, the differentially expressed genes were annotated with the Gene Ontology (GO) terms. The most significantly over-represented ontology classes for the genes differentiating between LN(+) and LN(−) tissue are given in Table 2. The most up-regulated genes include those coding for epidermal molecules as well as molecules associated with cell adhesion, locomotory behaviour and inflammatory response.

Table 2 Significant Gene Ontology categories overrepresented among differentially expressed genes

Gene selection for validation

Based on the results of microarray analysis by the pumaDE function of PUMA package 1.8.1, we aimed to select genes with the highest fold-change of expression between LN(+) and LN(−) obtained from the patient No. 15, and lowest significant adjusted P values, to be assessed by the quantitative PCR. Genes with over two-fold differences in expression levels were selected from the list and further categorized into different functions and cellular processes. The genes significantly differentially expressed (at P value below 0.01) in the analyzed LN obtained from patient No. 15 were represented by 408 probe sets (Additional file: 1 Table S1). The genes selected were those with the highest fold-change values, represented by specific probe sets (named with the suffix ‘_at’), and those belonging to different GO categories. The group of genes of interest was further reduced to those which also proved to differentiate LNs(+) from LNs(−) in the PUMA analysis of microarray data of the three additional pairs of nodes from three consecutive patients (No. 21, 46 and 61). A list of the genes differentially expressed in LN(+) and LN(−) samples obtained from these four VC patients is presented in the Additional file: 2 Table S2.

Next, only the genes up-regulated in the four LNs(+) obtained from the four VC patients were analyzed in Genevestigator V3 system. Genevestigator allowed to validate their increased expression in vulva, as this repository contains high quality microarray expression data from experiments performed on 27 non-cancer vulvar tissue samples (including nine specimens of vulvar intraepithelial neoplasia, VIN). The genes that showed high signal intensities in these 27 microarray data sets were selected. The probe sets were chosen on the basis of the highest % Present” value generated in the Genevestigator. These values represent a fraction of arrays in which the signal for a given probe set is above the background. Finally, seven genes were selected from the microarray studies: TP53 apoptosis effector - PERP, S100 calcium binding protein A8 - S100A8, fatty acid binding protein 5 (psoriasis-associated) - FABP5, stratifin - SFN, carbonic anhydrase 12 - CA12, junction plakoglobin - JUP and cystatin A - CSTA. Five out of the seven selected genes, namely, S100A8 FABP5 SFN JUP and CSTA, have been described by Hsiao et al. [15] as “vulva-selective” .

Validation of gene expression levels by real-time RT-PCR

To investigate the reliability of the cDNA microarray results and to check if our gene selection was correct, the expression levels of seven genes, PERP, S100A8, FABP5, SFN, CA12, JUP and CSTA, were measured by the qRT-PCR in 22 LN(+) and 49 LN(−) samples.

We confirmed the significant differences in the expression levels between uninvolved and involved (as assessed by routine histopathological analysis) of the all the seven genes in question (see Figure 1). As shown in Additional file: 3 Table S3, the expression levels of all of the seven genes tested by the qRT-PCR, in agreement with the microarray data, was lower in LN(−) than in LN(+).

Figure 1
figure 1

The PERP , S100A8 , FABP5 , SFN , CA12 , JUP and CSTA expression levels in lymph node samples, measured by qRT-PCR. Abbreviations: LN(+) - involved lymph node, LN(−) - uninvolved lymph node sample, PP value computed by the Wilcoxon test.

Correlations between gene expression levels and time to groin recurrence in VC patients

Subsequently, we examined a correlation between the expression levels of the seven selected genes and time to groin recurrence in VC patients using Cox's proportional hazard analysis. The highest expression value of a given gene of all obtained in different LN samples of the same patient was taken into account. Three genes, SFN, CA12 and JUP, significantly correlated with TTR. High expression levels of the set of the four remaining genes in LNs were also associated with shorter TTR, but the trend did not reach statistical significance (Table 3).

Table 3 Impact of gene expression levels in LN samples on TTR and its significance

Finally, patients were stratified by the arbitrarily chosen cut-off values, i.e. the medians of maximum gene expression levels in LNs from all the VC patients enrolled in the study (again, one node with the highest gene expression was taken into account per patient). Kaplan-Meier analysis was used for comparison of TTR between the groups of patients stratified according to gene expression status, i.e. high (above the cut-off value) vs low gene expression status (below the cut-off value) and irrespective of the histopathological evaluation of their LNs status. Initially, Kaplan-Meier estimates of TTR were used for the evaluation of the prognostic value of the SFN, CA12 and JUP gene expression, the genes whose expression was found to correlate significantly with TTR in above-mentioned COX-model analysis. Further, an analogous Kaplan-Meier analysis of the remaining four genes revealed the prognostic value of the expression level of CSTA gene. The Kaplan-Meier plots show the association of TTR with the LN expression levels of SFN, CA12, JUP and CSTA in VC patients, i.e. the expression levels of these genes predicted TTR (see Figure 2). The results of comparison of the Kaplan-Meier estimates for TTR depending on the expression level of SFN, CA12, JUP and CSTA genes, together with the cut-off values used for this analysis, are presented in Table 4. The differences in TTR between the groups of VC patients stratified by PERP, S100A8 and FABP gene expression levels were not significant (Kaplan-Meier curves not shown).

Figure 2
figure 2

Kaplan-Meier estimates of TTR according to SFN , CA12 JUP and CSTA expression status (high vs low) in regional LNs. Abbreviations: PP value calculated using the Mantel-Cox log-rank test; high - expression level, i.e. above the cut-off value; low - expression level, i.e. below the cut-off value (i.e. the median of maximum gene expression levels - as measured in one node per VC patient.

Table 4 The results of the comparison of Kaplan-Meier estimates for TTR for VC patients according to their LN gene expression levels


Lymph node status is the most important prognostic factor in patients with VC [16]. Routine pathological staging is limited to the microscopic evaluation of H&E-stained LN sections. Unfortunately, small tumour foci may occasionally be missed by the pathologist in this examination [17]. As stated by Regauer [18], sentinel LNs of VC patients with even single tumour cells should be regarded as positive. Thus, to improve detection of LN metastases, pathological ultrastaging should ideally be performed [19, 20]. This histological work-up of LNs involves serial sectioning and immunohistochemical analysis; the procedure is labour-intensive and time-consuming. Real-time RT-PCR technology is more sensitive, as it allows to detect a very small number of cells in larger volumes of previously pulverized (or homogenized) tissues. However, this method carries some disadvantages, for example it makes morphologic evaluation of the analyzed samples impossible. Therefore, a careful marker selection for RT-PCR analysis is a fundamental issue.

Since there are no molecular markers for diagnosing VC dissemination, we aimed to identify the potential mRNA markers of VC dissemination into LNs, and to verify their prognostic value by correlating their levels with time to groin recurrence. We used oligonucleotide microarrays which measure the expression level of over 47,000 transcripts and variants of human genes.

At first, the PUMA package, enabling a comparison between small-sample gene expression data [21] was employed, to analyze two samples obtained from one patient, who underwent a rapid progression of the disease during the follow-up time, with local recurrence, and groin recurrence in the same groin where the LN(+), and died due to cancer progression. The PUMA analysis produced a list of genes differentially expressed in LN(+) and LN(−). The most significantly over-represented gene ontology (GO) classes in the list included epidermal molecules. Thus, some of the differentially expressed genes are those specifically expressed in the primary tumour site. This finding testifies the methodology used for the comparison of the single microarray experiments. Similar results were obtained for the three additional pairs of LNs obtained from three patients.

Interestingly, gene expression signature distinguishing between the four pairs of histologically positive and negative LNs comprised no genes down-regulated in the metastatic LNs of the same patients. A possible explanation of this observation is that only transcripts that were highly up-regulated in the metastatic cells could be distinguished from the background of the numerous non-cancer cells as an “added value” while down-regulated transcripts, i.e. transcripts’ deprivation, would not be noticed in this background environment.

The consecutive steps of narrowing-down of the list of the differentially expressed in the four LNs obtained from the four VC patients genes led us to the choice of seven genes, namely PERP S100A8 FABP5 SFN CA12 JUP and CSTA, to be validated by the real time RT-PCR. Gene expression data of 27 non-cancer vulvar tissue samples deposited in the Genevestigator V3 system [22, 23] have shown that these seven genes are regularly up-regulated in vulvar tissue. In addition, Hsiao et al. [15], who created a compendium of tissue-selective genes based on microarray data, described five of these genes as “vulva-selective.” This strongly supports tissue specificity of the selected genes for the primary tumour site.

To our knowledge, this is the first time that these genes have been linked to VC dissemination. PERP, as a component of intercellular desmosome junctions, plays a role in epithelial integrity and cell-cell adhesion [24] and constitutes a proapoptotic transcriptional target of TP53 [25]. PERP-deficiency promotes cancer by enhancing cell survival, desmosome loss, and inflammation [26]. S100 calcium binding protein A8, S100A8, belongs to S100 proteins, a family of EF-hand signalling proteins [27]. The complex of S100A8 and S100A9 called calprotectin induces a proinflammatory and thrombogenic response, and is involved in danger signalling; what is more, its stimulation results in a loss of cell–cell contacts [28, 29]. S100A8/S100A9 expression is increased in patients with various tumours, being involved in invasion and migration processes [30, 31]. S100A8/9 expression is minimal in normal epidermis and is elevated in skin diseases [32]. However, Dell'oste et al. [33] showed that in HPV-immortalized keratinocytes S100A8/9 expression was downregulated. Fatty acid binding protein 5 (psoriasis-associated, FABP5), found in epidermal cells, was first identified to be up-regulated in psoriasis keratinocytes [34]. In cancer cells, FABP5 increases cell proliferation and invasiveness, as recently demonstrated in oral squamous cell carcinoma [35] where its expression may be HPV-related [36]. In non-cancer keratinocytes FABP5 elevation may be necessary for the activation of cell motility during epidermal wound healing [37]. Stratifin, SFN (14-3-3 sigma), plays a role in various cellular processes being the most cancer-associated 14-3-3 isoform [38]. 14-3-3 sigma expression is lost in numerous tumours [38], including VC [39], while its increased expression could be associated with a loss of TP53 function [40, 41]. However, Wang et al. [42] found high levels of SFN detected immunohistochemically in over 70% of analyzed VC tumors, significantly correlating to large tumor diameter and deep invasion. Importantly, 14-3-3 sigma protein is present in various normal epithelia and absent in LNs [43, 44], and lymphocytes may express its mRNA at relatively low levels [45]. Carbonic anhydrase XII is one of the tumour-associated carbonic anhydrases. The enzyme is overexpressed under hypoxic conditions and constitutes a possible target for anticancer therapy [46, 47]. Junction plakoglobin, JUP, a member of the catenin family, is present both in desmosomes and in intermediate junctions [48]. The majority of studies suggest that plakoglobin has a tumour suppressor role [49]. Cystatin A (stefin A), CSTA, functions as a cysteine protease inhibitor and plays a role in epidermal development and maintenance. Serum level of CSTA has been proposed as a tool predicting nodal stage and poor prognosis in nasopharyngeal carcinoma [50].

The comparison of LNs(+) with LNs(−) confirmed that LNs(+) express significantly more transcripts of PERP, S100A8, FABP5, SFN, CA12, JUP and CSTA. The results of PCR analysis showed particularly increased expression of S100A8 and SFN in the metastatic LNs, where there was over 10-fold change of the mean expression level. Higher expression of PERP, FABP5, CA12, JUP and CSTA were also found, but the differences were lower.

The significance of gene expression levels in LN samples for time to groin recurrence (TTR) was statistically analyzed using Cox's proportional hazard model. A coefficient computed by this method for each variable expression of “predictor” genes indicates the direction and degree of flexing that the predictor has on the TTR curve. A positive value of coefficients indicated that larger values of the expression were associated with greater groin recurrence rates. Two genes, namely SFN and CA12, were found significantly predictive of such recurrence. High expression levels of the other genes of the set of five in LNs were also associated with shorter TTRs, but the trend did not reach statistical significance.

The log-rank test used for testing significance of the TTR functions in Kaplan-Meier analysis yielded significant P values for separating the two groups of VC patients with different TTR according to their LN gene expression levels of SFN, CA12, JUP and CSTA. Patients with low expression of these genes had a superior TTR values of not reached versus 6.9 to 8.7 months.

It is worth noting that our study had some important limitations. Firstly, only seven of the most differentially expressed genes in the microarray experiment were analyzed by the qRT-PCR. Notwithstanding, the remaining genes (see Additional file: 2 Table S2) are candidates for validation in further studies. Secondly, LN assessment was based on routine pathological evaluation with H&E staining, while the state-of-the-art would be to perform ultrastaging on the excised LNs. Lastly, Kaplan-Meier analysis of TTR that was applied to SFN, CA12, JUP and CSTA expression, should be considered as just an illustration of potential clinical use of the identified marker genes. Future studies should include more patients to assure a better assessment of the novel markers’ performance. Such studies should also provide better cut-off values for patients’ stratification than the arbitrarily chosen median expression levels in LNs, and enable the ROC curve analysis. Importantly, our microarray data were obtained for HPV-negative patients in order to exclude the interference of the infection with the expression results and to focus on metastasis-associated genes’ selection. However, patients with HPV infection should also be enrolled in future studies to assess the influence of HPV on the markers’ performance.

Groin recurrences from VC are often fatal. Unfortunately, “the optimum mode of treatment and predictive factors for patients with groin recurrences are unknown” [51]. The size of LN metastases correlates with survival in patients with early [52] and advanced [53] stage VC. On the other hand, “inguinofemoral lymphadenectomy can be avoided when the sentinel node is negative for disease” [54]. However, although Oonk and colleagues [52] reported that the risk of non-sentinel LN involvement increases with size of sentinel LN metastases, they were unable to determine the cut-off size below which the risk of non-sentinel LN metastases would be inconsiderable. Therefore, the authors concluded that all patients with sentinel LN metastases - regardless of their size - should undergo additional groin treatment. The risk of additional metastases when only isolated tumour cells are present in the sentinel LN is 4% [55]. Still, approximately 12% of early-stage VC patients with a negative sentinel LN develop local recurrence” [54, 56].

For the reasons discussed above, novel means to better stratify early-stage VC patients in order to optimize treatment modality, i.e. a tool to decide who should not undergo lymphadenectomy vs who should receive additional treatment to the groin would be highly beneficial. Further validation of markers and development of molecular tests is indispensable for the reliable up- or down-staging of LNs. Our results may ideally be used to develop a test that could be used during surgery to decide whether to remove inguinal LNs, similarly to the RT-PCR based, FDA-approved GeneSearch™ Breast Lymph Node (BLN) Assay (Veridex, LLC, Warren, NJ), used for the rapid intra-operative detection of sentinel LN metastases in breast cancer. Such tests may reduce the need for a second surgery for the axillary LN dissection [5759]. As sentinel node procedure was confirmed to be safe in the early-stage VC patients [54], tests similar to GeneSearch™ should also be of clinical utility for VC patients.


To conclude, the performance of qRT-PCR assays employing the LN(+) marker genes that we have identified, PERP, S100A8, FABP5, SFN, CA12, JUP and CSTA, may provide a promising tool for intra-operative sentinel LN evaluation in VC patients. Moreover, all the above marker genes have a potential of prognostic biomarkers, however, before their incorporation into clinical setting further studies are necessary to confirm their prognostic value in the qRT-PCR assays.

Conflicts of interest

The authors report no potential conflicts of interest.


  1. Cancer incidence in five continents vol. VIII. Edited by: Parkin DM, Ferlay J, Teppo L, Thomas DB. 2002, IARC Scientific Publications No. 155 IARC,

  2. National Cancer Registry. available at:

  3. Raspagliesi F, Hanozet F, Ditto A, Solima E, Zanaboni F, Vecchione F, Kusamura S: Clinical and pathological prognostic factors in squamous cell carcinoma of the vulva. Gynecol Oncol. 2006, 102: 333-337. 10.1016/j.ygyno.2005.12.027.

    Article  PubMed  Google Scholar 

  4. Selman TJ, Luesley DM, Acheson N, Khan KS, Mann CH: A systematic review of the accuracy of diagnostic tests for inguinal lymph node status in vulvar cancer. Gynecol Oncol. 2005, 99: 206-214. 10.1016/j.ygyno.2005.05.029.

    Article  CAS  PubMed  Google Scholar 

  5. Hacker NF, Berek JS, Lagasse LD, Leuchter RS, Moore JG: Management of regional lymph nodes and their prognostic influence in vulvar cancer. Obstet Gynecol. 1983, 61: 408-412.

    CAS  PubMed  Google Scholar 

  6. Lugo TG, Braun S, Cote RJ, Pantel K, Rusch V: Detection and measurement of occult disease for the prognosis of solid tumors. J Clin Oncol. 2003, 21: 2609-2615. 10.1200/JCO.2003.01.153.

    Article  PubMed  Google Scholar 

  7. Rutgers EJ: Sentinel node biopsy: interpretation and management of patients with immunohistochemistry-positive sentinel nodes and those with micrometastases. J Clin Oncol. 2008, 26: 698-702. 10.1200/JCO.2007.14.4667.

    Article  PubMed  Google Scholar 

  8. Knopp S, Trope C, Nesland JM, Holm R: A review of molecular pathological markers in vulvar carcinoma: lack of application in clinical practice. J Clin Pathol. 2009, 62: 212-218. 10.1136/jcp.2008.057240.

    Article  CAS  PubMed  Google Scholar 

  9. Kowalewska M, Radziszewski J, Kulik J, Barathova M, Nasierowska-Guttmajer A, Bidzinski M, Pastorek J, Pastorekova S, Siedlecki JA: Detection of carbonic anhydrase 9-expressing tumor cells in the lymph nodes of vulvar carcinoma patients by RT-PCR. Int J Cancer. 2005, 116: 957-962. 10.1002/ijc.21106.

    Article  CAS  PubMed  Google Scholar 

  10. Kowalewska M, Szkoda MT, Radziszewski J, Ptaszynski K, Bidzinski M, Siedlecki JA: The frequency of human papillomavirus infection in polish patients with vulvar squamous cell carcinoma. Int J Gynecol Cancer. 2010, 20: 434-437. 10.1111/IGC.0b013e3181d320f1.

    Article  PubMed  Google Scholar 

  11. Ihaka R, Gentleman R: R: a language for data analysis and graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314.

    Google Scholar 

  12. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B. 1995, 57: 289-300.

    Google Scholar 

  13. Liu X, Milo M, Lawrence ND, Rattray M: A tractable probabilistic model for Affymetrix probe-level analysis across multiple chips. Bioinformatics. 2005, 21: 3637-3644. 10.1093/bioinformatics/bti583.

    Article  CAS  PubMed  Google Scholar 

  14. Hruz T, Laule O, Szabo G, Wessendorp F, Bleuler S, Oertle L, Widmayer P, Gruissem W, Zimmermann P: Genevestigator v3: a reference expression database for the meta-analysis of transcriptomes. Adv Bioinformatics. 2008, 2008: 420747-

    Article  PubMed  PubMed Central  Google Scholar 

  15. Hsiao LL, Dangond F, Yoshida T, Hong R, Jensen RV, Misra J, Dillon W, Lee KF, Clark KE, Haverty P, et al: A compendium of gene expression in normal human tissues. Physiol Genomics. 2001, 7: 97-104.

    Article  CAS  PubMed  Google Scholar 

  16. Hacker NF: Revised FIGO staging for carcinoma of the vulva. Int J Gynaecol Obstet. 2009, 105: 105-106. 10.1016/j.ijgo.2009.02.011.

    Article  PubMed  Google Scholar 

  17. Knopp S, Holm R, Trope C, Nesland JM: Occult lymph node metastases in early stage vulvar carcinoma patients. Gynecol Oncol. 2005, 99: 383-387. 10.1016/j.ygyno.2005.06.023.

    Article  PubMed  Google Scholar 

  18. Regauer S: Histopathological work-up and interpretation of sentinel lymph nodes removed for vulvar squamous cell carcinoma. Histopathology. 2009, 55: 174-181. 10.1111/j.1365-2559.2009.03350.x.

    Article  PubMed  Google Scholar 

  19. Gietema HA, Vuylsteke RJ, de Jonge IA, van Leeuwen PA, Molenkamp BG, van der Sijp JR, Meijer S, van Diest PJ: Sentinel lymph node investigation in melanoma: detailed analysis of the yield from step sectioning and immunohistochemistry. J Clin Pathol. 2004, 57: 618-620. 10.1136/jcp.2003.011742.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Weaver DL, Le UP, Dupuis SL, Weaver KA, Harlow SP, Ashikaga T, Krag DN: Metastasis detection in sentinel lymph nodes: comparison of a limited widely spaced (NSABP protocol B-32) and a comprehensive narrowly spaced paraffin block sectioning strategy. Am J Surg Pathol. 2009, 33: 1583-1589. 10.1097/PAS.0b013e3181b274e7.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Pearson RD, Liu X, Sanguinetti G, Milo M, Lawrence ND, Rattray M: puma: a Bioconductor package for propagating uncertainty in microarray analysis. BMC Bioinformatics. 2009, 10: 211-10.1186/1471-2105-10-211.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Santegoets LA, Seters M, Helmerhorst TJ, Heijmans-Antonissen C, Hanifi-Moghaddam P, Ewing PC, van Ijcken WF, van der Spek PJ, van der Meijden WI, Blok LJ: HPV related VIN: highly proliferative and diminished responsiveness to extracellular signals. Int J Cancer. 2007, 121: 759-766. 10.1002/ijc.22769.

    Article  CAS  PubMed  Google Scholar 


  24. Ihrie RA, Attardi LD: A new Perp in the lineup: linking p63 and desmosomal adhesion. Cell Cycle. 2005, 4: 873-876. 10.4161/cc.4.7.1836.

    Article  CAS  PubMed  Google Scholar 

  25. Reczek EE, Flores ER, Tsay AS, Attardi LD, Jacks T: Multiple response elements and differential p53 binding control Perp expression during apoptosis. Mol Cancer Res. 2003, 1: 1048-1057.

    CAS  PubMed  Google Scholar 

  26. Beaudry VG, Jiang D, Dusek RL, Park EJ, Knezevich S, Ridd K, Vogel H, Bastian BC, Attardi LD: Loss of the p53/p63 regulated desmosomal protein Perp promotes tumorigenesis. PLoS Genet. 2010, 6: e1001168-10.1371/journal.pgen.1001168.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Santamaria-Kisiel L, Rintala-Dempsey AC, Shaw GS: Calcium-dependent and -independent interactions of the S100 protein family. Biochem J. 2006, 396: 201-214. 10.1042/BJ20060195.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Ehrchen JM, Sunderkotter C, Foell D, Vogl T, Roth J: The endogenous Toll-like receptor 4 agonist S100A8/S100A9 (calprotectin) as innate amplifier of infection, autoimmunity, and cancer. J Leukoc Biol. 2009, 86: 557-566. 10.1189/jlb.1008647.

    Article  CAS  PubMed  Google Scholar 

  29. Viemann D, Strey A, Janning A, Jurk K, Klimmek K, Vogl T, Hirono K, Ichida F, Foell D, Kehrel B, et al: Myeloid-related proteins 8 and 14 induce a specific inflammatory response in human microvascular endothelial cells. Blood. 2005, 105: 2955-2962. 10.1182/blood-2004-07-2520.

    Article  CAS  PubMed  Google Scholar 

  30. Gebhardt C, Nemeth J, Angel P, Hess J: S100A8 and S100A9 in inflammation and cancer. Biochem Pharmacol. 2006, 72: 1622-1631. 10.1016/j.bcp.2006.05.017.

    Article  CAS  PubMed  Google Scholar 

  31. Moon A, Yong HY, Song JI, Cukovic D, Salagrama S, Kaplan D, Putt D, Kim H, Dombkowski A, Kim HR: Global gene expression profiling unveils S100A8/A9 as candidate markers in H-ras-mediated human breast epithelial cell invasion. Mol Cancer Res. 2008, 6: 1544-1553. 10.1158/1541-7786.MCR-08-0189.

    Article  CAS  PubMed  Google Scholar 

  32. Eckert RL, Broome AM, Ruse M, Robinson N, Ryan D, Lee K: S100 proteins in the epidermis. J Invest Dermatol. 2004, 123: 23-33. 10.1111/j.0022-202X.2004.22719.x.

    Article  CAS  PubMed  Google Scholar 

  33. Dell'oste V, Azzimonti B, Mondini M, De Andrea M, Borgogna C, Mesturini R, Accardi R, Tommasino M, Landolfo S, Dianzani U, Gariglio M: Altered expression of UVB-induced cytokines in human papillomavirus-immortalized epithelial cells. J Gen Virol. 2008, 89: 2461-2466. 10.1099/vir.0.83586-0.

    Article  PubMed  Google Scholar 

  34. Madsen P, Rasmussen HH, Leffers H, Honore B, Celis JE: Molecular cloning and expression of a novel keratinocyte protein (psoriasis-associated fatty acid-binding protein [PA-FABP]) that is highly up-regulated in psoriatic skin and that shares similarity to fatty acid-binding proteins. J Invest Dermatol. 1992, 99: 299-305. 10.1111/1523-1747.ep12616641.

    Article  CAS  PubMed  Google Scholar 

  35. Fang LY, Wong TY, Chiang WF, Chen YL: Fatty-acid-binding protein 5 promotes cell proliferation and invasion in oral squamous cell carcinoma. J Oral Pathol Med. 2010, 39: 342-348.

    Article  CAS  PubMed  Google Scholar 

  36. Melle C, Ernst G, Winkler R, Schimmel B, Klussmann JP, Wittekindt C, Guntinas-Lichius O, von Eggeling F: Proteomic analysis of human papillomavirus-related oral squamous cell carcinoma: identification of thioredoxin and epidermal-fatty acid binding protein as upregulated protein markers in microdissected tumor tissue. Proteomics. 2009, 9: 2193-2201. 10.1002/pmic.200800882.

    Article  CAS  PubMed  Google Scholar 

  37. Kusakari Y, Ogawa E, Owada Y, Kitanaka N, Watanabe H, Kimura M, Tagami H, Kondo H, Aiba S, Okuyama R: Decreased keratinocyte motility in skin wound on mice lacking the epidermal fatty acid binding protein gene. Mol Cell Biochem. 2006, 284: 183-188. 10.1007/s11010-005-9048-8.

    Article  CAS  PubMed  Google Scholar 

  38. Mhawech P: 14–3–3 proteins–an update. Cell Res. 2005, 15: 228-236. 10.1038/

    Article  CAS  PubMed  Google Scholar 

  39. Gasco M, Sullivan A, Repellin C, Brooks L, Farrell PJ, Tidy JA, Dunne B, Gusterson B, Evans DJ, Crook T: Coincident inactivation of 14-3-3sigma and p16INK4a is an early event in vulval squamous neoplasia. Oncogene. 2002, 21: 1876-1881. 10.1038/sj.onc.1205256.

    Article  CAS  PubMed  Google Scholar 

  40. Hermeking H, Lengauer C, Polyak K, He TC, Zhang L, Thiagalingam S, Kinzler KW, Vogelstein B: 14-3-3 sigma is a p53-regulated inhibitor of G2/M progression. Mol Cell. 1997, 1: 3-11. 10.1016/S1097-2765(00)80002-7.

    Article  CAS  PubMed  Google Scholar 

  41. Lee MH, Lozano G: Regulation of the p53-MDM2 pathway by 14-3-3 sigma and other proteins. Semin Cancer Biol. 2006, 16: 225-234. 10.1016/j.semcancer.2006.03.009.

    Article  CAS  PubMed  Google Scholar 

  42. Wang Z, Trope CG, Suo Z, Troen G, Yang G, Nesland JM, Holm R: The clinicopathological and prognostic impact of 14-3-3 sigma expression on vulvar squamous cell carcinomas. BMC Cancer. 2008, 8: 308-10.1186/1471-2407-8-308.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Dellambra E, Golisano O, Bondanza S, Siviero E, Lacal P, Molinari M, D'Atri S, De Luca M: Downregulation of 14-3-3sigma prevents clonal evolution and leads to immortalization of primary human keratinocytes. J Cell Biol. 2000, 149: 1117-1130. 10.1083/jcb.149.5.1117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Nakajima T, Shimooka H, Weixa P, Segawa A, Motegi A, Jian Z, Masuda N, Ide M, Sano T, Oyama T, et al: Immunohistochemical demonstration of 14-3-3 sigma protein in normal human tissues and lung cancers, and the preponderance of its strong expression in epithelial cells of squamous cell lineage. Pathol Int. 2003, 53: 353-360. 10.1046/j.1440-1827.2003.01481.x.

    Article  CAS  PubMed  Google Scholar 

  45. Bhatia K, Siraj AK, Hussain A, Bu R, Gutierrez MI: The tumor suppressor gene 14-3-3 sigma is commonly methylated in normal and malignant lymphoid cells. Cancer Epidemiol Biomarkers Prev. 2003, 12: 165-169.

    CAS  PubMed  Google Scholar 

  46. Pastorekova S, Zatovicova M, Pastorek J: Cancer-associated carbonic anhydrases and their inhibition. Curr Pharm Des. 2008, 14: 685-698. 10.2174/138161208783877893.

    Article  CAS  PubMed  Google Scholar 

  47. Brzozowski Z, Slawinski J, Innocenti A, Supuran CT: Carbonic anhydrase inhibitors. Regioselective synthesis of novel 1-substituted 1,4-dihydro-4-oxo-3-pyridinesulfonamides and their inhibition of the human cytosolic isozymes I and II and transmembrane cancer-associated isozymes IX and XII. Eur J Med Chem. 2010, 45: 3656-3661. 10.1016/j.ejmech.2010.05.011.

    Article  CAS  PubMed  Google Scholar 

  48. Delva E, Tucker DK, Kowalczyk AP: The desmosome. Cold Spring Harb Perspect Biol. 2009, 1: a002543-10.1101/cshperspect.a002543.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Chidgey M, Dawson C: Desmosomes: a role in cancer?. Br J Cancer. 2007, 96: 1783-1787. 10.1038/sj.bjc.6603808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Chang KP, Wu CC, Chen HC, Chen SJ, Peng PH, Tsang NM, Lee LY, Liu SC, Liang Y, Lee YS, et al: Identification of candidate nasopharyngeal carcinoma serum biomarkers by cancer cell secretome and tissue transcriptome analysis: potential usage of cystatin A for predicting nodal stage and poor prognosis. Proteomics. 2010, 10: 2644-2660. 10.1002/pmic.200900620.

    Article  CAS  PubMed  Google Scholar 

  51. Cormio G, Loizzi V, Carriero C, Cazzolla A, Putignano G, Selvaggi L: Groin recurrence in carcinoma of the vulva: management and outcome. Eur J Cancer Care (Engl). 2010, 19: 302-307. 10.1111/j.1365-2354.2008.01011.x.

    Article  CAS  Google Scholar 

  52. Oonk MH, van Hemel BM, Hollema H, de Hullu JA, Ansink AC, Vergote I, Verheijen RH, Maggioni A, Gaarenstroom KN, Baldwin PJ, et al: Size of sentinel-node metastasis and chances of non-sentinel-node involvement and survival in early stage vulvar cancer: results from GROINSS-V, a multicentre observational study. Lancet Oncol. 2010, 11: 646-652. 10.1016/S1470-2045(10)70104-2.

    Article  PubMed  Google Scholar 

  53. Origoni M, Sideri M, Garsia S, Carinelli SG, Ferrari AG: Prognostic value of pathological patterns of lymph node positivity in squamous cell carcinoma of the vulva stage III and IVA FIGO. Gynecol Oncol. 1992, 45: 313-316. 10.1016/0090-8258(92)90311-6.

    Article  CAS  PubMed  Google Scholar 

  54. Van der Zee AG, Oonk MH, De Hullu JA, Ansink AC, Vergote I, Verheijen RH, Maggioni A, Gaarenstroom KN, Baldwin PJ, Van Dorst EB, et al: Sentinel node dissection is safe in the treatment of early-stage vulvar cancer. J Clin Oncol. 2008, 26: 884-889. 10.1200/JCO.2007.14.0566.

    Article  PubMed  Google Scholar 

  55. Oonk MH, de Hullu JA, van der Zee AG: Current controversies in the management of patients with early-stage vulvar cancer. Curr Opin Oncol. 2010, 22: 481-486. 10.1097/CCO.0b013e32833c06da.

    Article  PubMed  Google Scholar 

  56. Tantipalakorn C, Robertson G, Marsden DE, Gebski V, Hacker NF: Outcome and patterns of recurrence for International Federation of Gynecology and Obstetrics (FIGO) stages I and II squamous cell vulvar cancer. Obstet Gynecol. 2009, 113: 895-901.

    Article  PubMed  Google Scholar 

  57. Julian TB, Blumencranz P, Deck K, Whitworth P, Berry DA, Berry SM, Rosenberg A, Chagpar AB, Reintgen D, Beitsch P, et al: Novel intraoperative molecular test for sentinel lymph node metastases in patients with early-stage breast cancer. J Clin Oncol. 2008, 26: 3338-3345. 10.1200/JCO.2007.14.0665.

    Article  PubMed  Google Scholar 

  58. Funasako Y, Uenosono Y, Hirata M, Arigami T, Yanagita S, Arima H, Ehi K, Kijima Y, Yoshinaka H, Natsugoe S: Utility of the GeneSearch breast lymph node assay for the rapid evaluation of sentinel lymph nodes in breast cancer. Cancer. 2010, 116: 4450-4455. 10.1002/cncr.25479.

    Article  PubMed  Google Scholar 

  59. Mansel RE, Goyal A, Douglas-Jones A, Woods V, Goyal S, Monypenny I, Sweetland H, Newcombe RG, Jasani B: Detection of breast cancer metastasis in sentinel lymph nodes using intra-operative real time GeneSearch BLN Assay in the operating room: results of the Cardiff study. Breast Cancer Res Treat. 2009, 115: 595-600. 10.1007/s10549-008-0155-6.

    Article  PubMed  Google Scholar 

Pre-publication history

Download references


The authors are grateful to Drs. R. Nowak and M. Chechlinska for the critical reading of the manuscript. This work was supported by the research grants No. P05E 11624 and NN407 125 937 from the Ministry of Science and Higher Education. K.G. is supported by the Ministry of Science and Higher Education grant No. NN518331736.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Magdalena Kowalewska.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MK and JR designed and coordinated the study. MBidzinski, JAS and MJ participated in its coordination. MK wrote the manuscript. KG and MBujko performed the statistical analysis. MO-W and MK carried out the microarray and qRT-PCR experiments, respectively. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S2. Genes differentially expressed in LN(+) and LN(−) sample obtained from VC patient No 15. Abbreviations: P - p value; cP - Benjamini-Hochberg corrected p value; FC - fold change of gene expression ratio in LN(+) and LN(−) sample; exp. - expression level. (DOC 628 KB)


Additional file 2: Table S3. Genes differentially expressed in the four LN(+) and LN(−) sample pairs obtained from VC patients No 15, 21, 46 and 61. Abbreviations: P - p value; NA - not applicable. (DOC 5 MB)


Additional file 3: Table S1. Mean levels of expression of PERP, S100A8, FABP5, SFN, CA12, JUP and CSTA in LN(+) and LN(−) samples, measured by qRT-PCR. Abbreviations: ME – mean expression level in all involved [LN(+)] and uninvolved [LN(−)] lymph node samples included in the qRT-PCR analysis. (DOC 30 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, 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

Kowalewska, M., Radziszewski, J., Goryca, K. et al. Estimation of groin recurrence risk in patients with squamous cell vulvar carcinoma by the assessment of marker gene expression in the lymph nodes. BMC Cancer 12, 223 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: