Estimation of groin recurrence risk in patients with squamous cell vulvar carcinoma by the assessment of marker gene expression in the lymph nodes
© Kowalewska et al.; licensee BioMed Central Ltd. 2012
Received: 9 October 2011
Accepted: 6 June 2012
Published: 6 June 2012
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.
KeywordsVulvar carcinoma Lymph node Microarray Expression marker Real-time RT-PCR
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 . 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) . 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 . 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 . 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% . 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  need to be identified. In this study, to identify specific mRNA markers of VC dissemination into LNs, expression microarrays were employed.
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) .
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(−)].
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).
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  (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  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 . Differentially expressed probe sets were annotated with Gene Ontology (GO) terms (GO.db version 2.2.5) using the Bioconductor  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; http://www.r-project.org). 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 .
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).
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 . 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(−).
Significant Gene Ontology categories overrepresented among differentially expressed genes
multicellular organismal development
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.  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.
Correlations between gene expression levels and time to groin recurrence in VC patients
Impact of gene expression levels in LN samples on TTR and its significance
The results of the comparison of Kaplan-Meier estimates for TTR for VC patients according to their LN gene expression levels
TTR for patients with high expression level
TTR for patients with low expression level
Median TTR (months)
Median TTR (months)
95% CI, 0.06-0.86
95% CI, 0.04-0.61
95% CI, 0.05-0.74
95% CI, 0.04-0.66
Lymph node status is the most important prognostic factor in patients with VC . 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 . As stated by Regauer , 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  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. , 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  and constitutes a proapoptotic transcriptional target of TP53 . PERP-deficiency promotes cancer by enhancing cell survival, desmosome loss, and inflammation . S100 calcium binding protein A8, S100A8, belongs to S100 proteins, a family of EF-hand signalling proteins . 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 . However, Dell'oste et al.  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 . In cancer cells, FABP5 increases cell proliferation and invasiveness, as recently demonstrated in oral squamous cell carcinoma  where its expression may be HPV-related . In non-cancer keratinocytes FABP5 elevation may be necessary for the activation of cell motility during epidermal wound healing . Stratifin, SFN (14-3-3 sigma), plays a role in various cellular processes being the most cancer-associated 14-3-3 isoform . 14-3-3 sigma expression is lost in numerous tumours , including VC , while its increased expression could be associated with a loss of TP53 function [40, 41]. However, Wang et al.  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 . 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 . The majority of studies suggest that plakoglobin has a tumour suppressor role . 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 .
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” . The size of LN metastases correlates with survival in patients with early  and advanced  stage VC. On the other hand, “inguinofemoral lymphadenectomy can be avoided when the sentinel node is negative for disease” . However, although Oonk and colleagues  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% . 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 [57–59]. As sentinel node procedure was confirmed to be safe in the early-stage VC patients , 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.
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.
- Cancer incidence in five continents vol. VIII. Edited by: Parkin DM, Ferlay J, Teppo L, Thomas DB. 2002, IARC Scientific Publications No. 155 IARC,
- National Cancer Registry. available at: http://epid.coi.waw.pl/krn/
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ihaka R, Gentleman R: R: a language for data analysis and graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314.Google Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- 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-View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Mhawech P: 14–3–3 proteins–an update. Cell Res. 2005, 15: 228-236. 10.1038/sj.cr.7290291.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Pastorekova S, Zatovicova M, Pastorek J: Cancer-associated carbonic anhydrases and their inhibition. Curr Pharm Des. 2008, 14: 685-698. 10.2174/138161208783877893.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Delva E, Tucker DK, Kowalczyk AP: The desmosome. Cold Spring Harb Perspect Biol. 2009, 1: a002543-10.1101/cshperspect.a002543.View ArticlePubMedPubMed CentralGoogle Scholar
- Chidgey M, Dawson C: Desmosomes: a role in cancer?. Br J Cancer. 2007, 96: 1783-1787. 10.1038/sj.bjc.6603808.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/223/prepub