Investigation of the molecular profile of basal cell carcinoma using whole genome microarrays

Background Skin cancer accounts for 1/3 of all newly diagnosed cancer. Although seldom fatal, basal cell carcinoma (BCC) is associated with severe disfigurement and morbidity. BCC has a unique interest for researchers, as although it is often locally invasive, it rarely metastasises. This paper, reporting the first whole genome expression microarray analysis of skin cancer, aimed to investigate the molecular profile of BCC in comparison to non-cancerous skin biopsies. RNA from BCC and normal skin specimens was analysed using Affymetrix whole genome microarrays. A Welch t-test was applied to data normalised using dCHIP to identify significant differentially-expressed genes between BCC and normal specimens. Principal component analysis and support vector machine analysis were performed on resulting genelists, Genmapp was used to identify pathways affected, and GOstat aided identification of areas of gene ontology more highly represented on these lists than would be expected by chance. Results Following normalisation, specimens clustered into groups of BCC specimens and of normal skin specimens. Of the 54,675 gene transcripts/variants analysed, 3,921 were differentially expressed between BCC and normal skin specimens. Of these, 2,108 were significantly up-regulated and 1,813 were statistically significantly down-regulated in BCCs. Conclusion Functional gene sets differentially expressed include those involved in transcription, proliferation, cell motility, apoptosis and metabolism. As expected, members of the Wnt and hedgehog pathways were found to be significantly different between BCC and normal specimens, as were many previously undescribed changes in gene expression between normal and BCC specimens, including basonuclin2 and mrp9. Quantitative-PCR analysis confirmed our microarray results, identifying novel potential biomarkers for BCC.


Background
It is estimated that the incidence of cutaneous basal cell carcinoma is increasing worldwide by up to 10% per year [1] and it currently accounts for approximately 80% of all non-melanoma skin cancer -with highest rates in elderly men and increasing incidence in young women [2]. Several sub-types of BCC have been identified. These include nodular-ulcerated BCC (the most frequently occurring type; often with ulceration ("rodent ulcer")); superficial BCC (often multiple); sclerosing BCC (cancer cells surrounded by dense fibrosis and so resemble scars; highest recurrence rate of BCCs after treatment); cystic BCC (uncommon; tumour undergoes central degradation to form a cystic lesion); linear BCC (recently recognised clinical entity with increased risk for aggressive histopathology); and micronodular BCC (small tumour nests; often with subclinical growth). Although BCC only occasionally (0.003-0.55% of cases [3][4][5]), results in metastasis and is seldom fatal, BCC is often locally invasive with destructive growth and may be associated with severe disfigurement and morbidity as a result of local tissue destruction or due to necessary surgery. Furthermore, people with BCC are at higher risks of developing further BCCs and other malignancies, including squamous cell carcinomas, malignant melanomas, and possibly also non-cutaneous malignancies [1].
Current selection of best treatment for BCC is based on accurate diagnosis and sub-classification of these cancers, mainly on histomorphology/pathology of H&E stained sections [6]. Features associated with recurrence and metastasis are also considered -including tumour diameter >2 cm, location in the central part of the face or ear, present for long duration, incomplete excision, aggressive growth pattern (based on histology) and perinuclear or perivascular involvement [2]. Genes reported to be associated with susceptibility to BCC include CYP2D6, GST-T1, vitamin D receptor, and TNF; with UVB irradiaton known to cause mutations in the p53 tumour suppressor gene, leading to the development of this cancer [1].
With the exception of a single study of BCCs using a small cDNA microarray -representing 1,718 genes [7] -investigations aimed at identifying BCC biomarkers and understanding the molecular events involved in this disease have, in general, been limited to one-at-a-time studies, built on chance analyses of proteins or mRNAs. Examples of such protein analysis have identified CD10 [8,9], p63 [10], low expression levels of CD44 [11] to be associated with the presence of BCC, generally with absence of ICAM-1 and LFA-3 adhesion molecule expression [12] and with Ki67 expression levels differing between BCCs that recur, compared to those that do not recur [13]. RT-PCR analysis has indicated 1, 25-dihydroxyvitamin D(3) receptor mRNA levels to be increased in BCCs compared to normal skin [14], while qPCR quantification of gli1 transcripts has been found to discriminate BCC (and trichoepithelioma) from other skin cancers [15].
While such studies have indicated the importance and relevance of gene expression analyses in BCC, the number of gene products simultaneously analysed have been very limited. In order to increase our understanding of the molecular events involved in the development/expression of BCC, here we report our findings from whole genome microarray analysis of BCC and normal skin specimens.

Quality Control of Microarray Data
As indicated in Table 1, quality control (Q.C.) analysis of all 25 microarray data sets (from 20 BCC and 5 normal skin specimens) indicated an average percentage present call of 42.68% (+/-5.79 SD). This would be of the order expected for high quality RNA from cell lines from many origins (Affymetrix Inc. "Genechip ® Expression Analysis Data Analysis Fundamentals" [16], indicating that these results are acceptable for further analysis. Again, based on cell line Q.C. parameters, the accepted background levels are <100, while 52.31 +/-4.08 was found in this study; acceptable noise levels are <3, here we report 1.67 +/-0.4; and the acceptable scaling factor is <3 fold between data sets being compared. This was generally, but not always, achieved. The acceptable 3'/M ratio of <3 was achieved in 18/25 cases.

Data Analysis
Approximately 7% (3,921/54,675) of the probe sets representing transcripts on the microarray were significantly differentially expressed between BCC and normal skin specimens (Tables 2 &3; [see Additional Files 1 &2] for further information). The scatter plot of data differentially expressed ≥1.2 fold (Fig. 1) indicates an even, normal distribution of data. As shown in Table 4, of the 2,108 upregulated by ≥1.2 fold, genes involved in many crucial aspects of cellular biology, including metabolism, transcription, cell cycle regulation, cell adhesion, cell migration, cell proliferation and cell motility were amongst the largest groups of genes affected, while oxidative phosphorylation, lipid metabolism, translation, and apoptosis were among the main categories down-regulated ≥1.2 fold in BCCs compared to normal skin (Table 5). Of the 748 probesets representing transcripts up-regulated by ≥2 fold, approximately 11 were described as cloned cDNAs, 122 were ESTs, 6 were described as hypothetical genes, and 46 encoded hypothetical proteins. Of the 484 transcripts showing ≥2 fold down-regulation, twenty-five represented hypothetical proteins, 11 cloned cDNAs/RIKENS and 49 ESTs.
In order to identify specimen similarity/diversity in our group of 25 skin specimens, condition tree clustering (using Pearson's correlation coefficient as similarity metric) and principal component analysis (PCA) were performed using GeneSpring software. As indicated in Fig.  2(A), while Pearson's correlation coefficient does not suggest a significant difference (i.e. no values <0.9) between the 25 specimens analysed, the 5 normal specimens form a discrete cluster in relation to the BCC specimens. (A similar clustering pattern was obtained when Spearman's correlation coefficient was investigated). Two of the BCC specimens, T19 and T22, are apparently different (but not significantly different) to the other 18 BCCs. While the 3'/ M ratios were greater than expected (at least in comparison to cell line data) for these specimens, this was also the case for BCC26, T16, T24, T25 and T28, which did not group with T19 and T22. The higher scaling factor resulting from analysis of T19 and T22, compared to all other specimens, may be responsible for/contribute to their apparently somewhat different behaviour as represented on the condition tree. For the purpose of investigating the effects of T19 and T22 on the overall dataset, a re-analysis was performed excluding data relating to those specimens. This resulted in a reduction from 3,921 (i.e. 7.17%) to 3,865 (7.06%) of significantly differentially expressed transcripts between BCC and normal specimens, indicat-ing that the vast majority of these transcripts are unbiased by the slightly different behaviour of T19 and T22 compared to all other BCCs (as described above).
To further identify the relatedness of the BCC samples to each other, we performed principal components analysis [17] on the entire data set. PCA was carried out on logtransformed data, using mean centering and scaling. As can be seen in Fig. 2(B), while our PCA analysis has divided the specimens into two groups, the results indicate that the BCC sub-group is much more varied than the normal skin group and that the BCC and normal skin specimens do not differ greatly. Support Vector Machine analysis [18] is a machine learning classification approach which is suitable for application to the dimensionality of microarray data. It operates by examining the expression information of a set of data points whose classification is known (referred to as the "training set"), from which a defined number of classification predictor genes are identified. This predictor genelist can then be applied to a separate set of genes which are known not to be members of the functional class (referred to as the "test set"). The predictor genelist is user-defined, so it is more beneficial to the user for validation if a smaller number of genes comprise the predictor list.    Biological processes that are significantly enriched in our set of 2,108 transcripts found to be significantly up-regulated in BCC compared to normal skin specimens. Shown here is a sub-set of 37 representative significant GO annotations.
For the purposes of this study, we carried out SVM analysis repeatedly on the tumour vs. normal dataset, in order to identify the lowest number of probesets that discriminated 100% of the time between the two classifications. The gene selection method used was Fisher's Exact Test, the kernel function was Polynomial Dot Product (Order 1) with zero Diagonal Scaling Factor. A minimum number of six identifiers, shown in Table 6, were required for 100% classification of every specimen as either tumour or normal in origin. Attempts to further prioritise these genes resulted in a decrease in the efficiency of the classification, as further analysis aimed at classifying the specimens using less than six transcripts resulted in the mis-classification of normal specimen N5 as a BCC specimen.

qPCR Validation of Microarray Data
Quantitative-PCR (qPCR) analysis for 10 potential endogenous controls was performed on a random selection of four BCC and four normal skin specimens, to select a suitable endogenous control(s) amplifiable in all specimens and where levels of expression did not differ greatly between specimens being analysed. All 10 transcripts studied were found to be expressed in all specimens analysed and, as indicated in Fig. 3 [Additional File 3], expression levels and standard deviation results suggested that Biological processes that are significantly enriched in our set of 1,813 transcripts found to be significantly down-regulated in BCC compared to normal skin specimens. Shown here is a sub-set of 29 representative significant GO annotations.
the 10 transcripts are likely to be of similar suitability as endogenous controls. Two controls, β-actin and GAPDH, were subsequently selected for amplification in all specimens, in parallel with a selection of 5 target transcripts of interest. The relative expression levels of PTCH1, gli2, Frizzled D2, basonuclin 2, and chromagranin A were analysed by qPCR. As indicated in Table 7, while the fold differences detected by microarray and qPCR methods differed to some extent, the trend (i.e. absent/low expression in normal skin and induced/increased expression in BCCs) was always found.

Discussion
Basal cell carcinoma (BCC), the most common skin cancer in humans, is locally aggressive and relentlessly invasive, but generally does not metastasise [5,19]. Despite this, studies aimed at investigating the molecular mechanisms associated with -possibly responsible for -BCCs, have been very limited.  [7] were not reported in their manuscript and information is not publicly available on transcripts that were present on their microarray, but were not significantly changed, further comparisons between these studies cannot be performed.
The development of BCC is known to be associated with dys-regulation of the hedgehog and Wnt pathways [2,20]. Lack of expression and/or suppressed activation of patched homologue 1 (PTCH1), a tumour suppressor gene that forms part of the hedgehog signaling network [21], has been reported to be fundamental to the development of BCC [22]. Disruption of this tumour suppressor gene results in up-regulated cell proliferation [6]. The accepted mechanism of PTCH1's action is via its binding to another transmembrane molecule, smoothed (SMO), thus suppressing intracellular signaling. Following binding of sonic hedgehog (shh) to PTCH1, this suppressor activity is, however, quenched, resulting in uninterrupted signal transduction by SMO, via GLI transcription factors, and subsequent constitutive activation of target genes, Scatter Plot of the 3,921 gene transcripts identified as signifi-cantly differentially expressed (by ≥1.2 fold; ≥100 difference in expression intensity; P < 0.05) between BCC and normal skin specimens Figure 1 Scatter Plot of the 3,921 gene transcripts identified as significantly differentially expressed (by ≥1.2 fold; ≥100 difference in expression intensity; P < 0.05) between BCC and normal skin specimens. Transcripts significantly up-regulated are shown in red; those down-regulated are shown as green.
(A) Condition tree distribution of the 20 BCC (red) and 5 normal skin (yellow) specimens (following dCHIP normalisation, sim-ilarity measure, Pearson's correlation, clustering algorithm, average linkage) including members of the Wnt pathway [23] and PTCH1, itself [21]. From our analysis of BCC compared to normal skin tissue, while we found no significant changes in shh expression levels, we report an approximately 11 Table 7). It is important to note that the likely reason for the difference in fold expression detected by microarray and qPCR methods is due to their differing baseline sensitivity. Low expression levels detected by microarrays -e.g. gli2 in normal skinresulting in large fold differences are considered as induction from absent in normal skin compared to present in BCC specimens, when analysed by qPCR. However, although PTCH1 mRNA levels have previously been reported as enhanced in nodular BCC but undetectable in superficial BCC [24], here we report PTCH1 to be detectable in both of these histological types of BCC, with no significant difference in their respective expression values (ttest: p = 0.637).
While the increased expression of gli2 detected in BCCs compared to normal skin may be expected and associated with the development/presence of BCC, the lack of tumour suppressor activity by PTCH1 -despite its increased mRNA levels -may be due to lack of expression of its corresponding protein and/or lack of binding to SMO (levels of which were not significantly different between BCC and normal skin). As PTCH1 may shuttle between the cell membrane and endocytotic vesicles in response to active hedgehog ligand, it is obvious that not only its mRNA expression, but also its protein expression (at the relevant location, binding of SMO, etc.) is necessary to exert its tumour suppressor activity [21]. Furthermore, as loss-of-function mutations of PTCH1 have been identified in 30-40% of sporadic cases of BCC, it may be that the mRNA over-expressed in the BCCs is not coding for a functional protein.
Wnt signaling may be able to regulate a number of the aspects of the biology of tumour cells and thus contribute in several ways to the tumour phenotypes. The strongest link is to the control of proliferation. Knockouts of Wnt signal transduction components, including Wnt5A, can result in proliferative failure [25] while up-regulation of Wnt5A mRNA expression been associated with a range of cancer types, including breast, lung, prostate and malignant melanomas [26,27]. In our study, the involvement of Wnt signaling pathway in BCC is suggested by the significantly increased expression of a number of Wnt family members. These include Wnt5A (  Note: 1 = Expression levels, Affymetrix arbitrary units; 2 = fold increased expression in BCC compared to normal; 3 = following normalisation of data on mean of β-actin + GAPDH expression and calibrating to a pooled skin specimen control (as described in Materials and Methods); P = induced in BCC, while undetected in normal skin In the "canonical" Wnt signaling pathway, secreted ligands bind to Frizzled receptors and regulate the stability of β-catenin. (Given the large number of mammalian Frizzleds and Wnts, considerable ligand-receptor specificity might be expected; however, redundancy of function seems to be the rule [25]). The subsequent accumulation of β-catenin -the central player in the Wnt pathway [28] -in the nucleus, results in its participation in transcriptionally active complexes with members of the LEF/TCF family of transcription factors [29]. While we did not find levels of β-catenin to be significantly different between BCC and normal skin, decreased (-2.  [33], were also found to be increased in BCCs compared to normal skin. Differential expression of other transcription factors associated with cancer has also been found in this study. These include CHES1 (checkpoint suppressor 1) which is apparently involved in repressing expression of genes important for tumorigenesis [34]. CHES1 mRNA has been reported as down-regulated in oral squamous cell carcinoma [35] and in hepatocellular carcinoma [36].
Here we found CHES1 mRNA levels to be significantly (-2.03 fold; P = 0.045) [expression values 920.32 vs. 452.95], down-regulated in BCC compared to normal skin. Not unexpectedly, mRNAs encoding proteins involved in inducing apoptosis were also found to be down-regulated. These include CIDE [37] and CARD15 [38]  Using support vector machine analysis we have identified 6 transcripts that, as a group, enable the accurate classification of all 25 specimens as BCC or normal skin. These include matrix metalloproteinase 11/mmp11 (previously associated with the presence of other cancer types, including oral [39,40], lung [41] and breast [42]); hydroxysteroid (17-beta) dehydrogenase 4/hsd17b4 (changes in expression of which have been associated with breast cancer ( [43]); B-cell receptor-associated protein 29/bap29 (a member of the B cell receptor-associated family of proteins [44][45][46]); 1acylglycerol-3-phosphate O-acyltransferase 3/agpat3 (which catalyses the acylation of lysophosphatidic acid to form phosphatidic acid, the precursor of all glycerolipids [47]); as well as zinc finger protein 548 and full-length cDNA clone CS0DK008YI09 of HeLa cells Cot 25-normalized of homo sapiens, on which no literature has previously been published. Future studies involving the co-analysis of this group of 6 transcripts in larger cohort of BCC and normal skin specimens should enable validation of their diagnostic relevance.
Molecular events responsible for the quite unique invasive, but non-metastatic, nature of BCCs are not known. However, it is interesting to note that expression of certain genes believed to be involved in malignant invasion and metastasis of another form of skin cancer, i.e. malignant melanoma, apparently differ in their expression patterns in BCCs. AP-2 transcription factor is not expressed in malignant melanoma cells [48], but it is significantly upregulated (by 10 Differential expression, between BCC and normal skin, of many transcripts not previous associated with the presence of BCC was also found during the course of our study. Basonuclin (now termed Basonuclin 1; was first discovered in cultured human epidermal cells [51]) and the more recently discovered basonuclin 2 [52] are zinc finger proteins. Basonuclin 1 is expressed at high levels in proliferating keratinocytes of stratified squamous epithelium. During terminal differentiation of squamous epithelium, basonuclin mRNA and protein disappear from the suprabasal cells [53,54]. In normal cultured human keratinocytes, basonuclin 1 is the predominant transcript -with expression levels approximately 10 fold that of basonuclin 2 [52]. In this study we found that basonuclin 1 levels did not differ significantly between BCC and normal skin, although a previous study of 3 BCC and 2 normal skin specimens indicated increased expression of basonuclin 1 in BCC, dependent on Gli protein expression [55]. In contrast, we report basonuclin 2 levels to be significantly increased in BCCs. This finding was observed with all 4 basonuclin 2 probe sets present on the microarray, indicating a 6.7-9.6 fold increase level (P < 0.00005) of expression [greatest change in expression values being 85.63 vs. 828.01], in BCC. This increased expression of basonuclin 2 in BCC compared to normal skin was confirmed by qPCR analysis. This, we believe, is the first study indicating an association between expression of basonuclin 2 and BCC.
ABCC12/mrp9, identified in 2001 [56] is one of a superfamily of 9 ATP-binding cassette (ABC) multiple drug resistant proteins [57]. Mrp9 encodes an approximately 100 kDa protein detectable in breast cancer, normal breast tissue and testis, while an alternative mrp9 transcript -encoding an approximately 25 kDa protein -has been detected in normal brain, skeletal muscle and ovary tissues. Due to the differential levels of expression of mrp9 transcripts in breast tumour and normal tissue, MRP9 has been proposed as an immunotherapy target for breast cancer [58]. The functional relevance of our observation of approximately 8.7 fold greater levels of mrp9 [expression values 21.12 vs. 184.71] in BCC compared to normal skin has yet to be determined; its presence is unlikely to be involvement in drug resistance, as all of the BCCs included in this study were chemotherapy-naive (indeed, except in a limited number of advanced cases of BCC, chemotherapy is not used as a therapy for this disease. In these exceptional cases, excellent response rates have been reported with cisplatin in combination with either 5fluorouracil or doxorubicin [59]). However, as for breast cancer, mrp9 mRNA may be useful as a member of a panel of BCC biomarkers or as an immunotherapy.
Chromagranin A (ChgA; parathyroid secretory protein 1) is an established tissue marker associated with neuroendocrine differentiation -and indicative of outcome -in non-small cell lung carcinoma [60]. Increased levels of ChgA in serum have been associated with poor prognosis/ shortened survival for prostate cancer patients [61]. ChgA protein levels have been proposed to assist in the diagnosis of Merkel cell carcinoma patients who may benefit from oncological therapy [62][63][64]. Although described as relatively uncommon -using analysis of markers including ChgA -neuroendocrine differentiation in BCC has been reported [65]. In this study we present the first whole genome expression microarray analysis of skin cancer, aimed at investigating the molecular profile of BCC in comparison to non-cancerous skin biopsies. This investigation has not only confirmed previous findings from analyses of limited numbers of transcripts, but it has also identified changes in expression of mRNAs that had never previously been associated with this disease. The results from this work are interesting and exciting, but it is necessary to recognise their preliminary nature. Further analyses, building on our findings, should include independent replication studies so that the true relevance of these findings may be realised.

Conclusion
The success of this study indicates the feasibility and relevance of using whole genome microarrays to study BCC specimens. In addition to confirming previous findings, this work has increased our understanding of molecular differences between BCC and normal skin and has identified a number of novel potential biomarkers for BCC. Future studies including BCC tissue and normal skin tissue from the same individual, thus lowering inter-individual variability and ruling out genetic influences; analyses of age-and gender-matched cases; studies focusing on molecular profiling and comparisons of sub-types of BCC, with due consideration given to disease duration (as early tumours may have a different gene expression profile to prolonged tumours); as well as analysing BCCs that metastasise compared to those that do not, should further increase our understanding of this disease and assist in management of the individual BCC patients. Further-more, as gene expression may be independent of protein levels, future confirmatory analysis at the protein level would complement these findings.

Patient Characteristics
This study involved analysis of 20 BCC biopsies from both male and female patients aged between 47 years and 83 years (mean age = 65 +/-11 years; median = 67 years) at the time of diagnosis. In order to gain an understanding of the most common types of BCC, we elected to include a range of BCC sub-types in this study, rather than to focus on any particular sub-type. For this reason, BCC sub-classifications included were nodular/micronodular, superficial and sclerosing. Tissue specimens from these twenty cases of BCC were procured at Blackrock Clinic and the Bons Secours Hospital, Dublin, examined macroscopically, immediately snap-frozen in liquid nitrogen, and were subsequently stored at -80°C until required for analysis. Five normal skin specimens (from consenting male and female volunteers of a similar age range who do not/ never had skin cancer) were also included in these studies.

RNA Extraction
For RNA analyses, dissected tumours that had been snapfrozen in liquid nitrogen and then stored at -80°C until required were homogenised, on ice, in 1 ml TriReagent (Sigma; Poole, England) and total RNA was subsequently isolated according to the manufacturer's instructions. RNA quantity and purity were assessed at 260 nm and 280 nm using a Nanodrop (ND-1000; Labtech. International); an Agilent bioanalyser (Agilent 2100; Agilent Technologies) was used to assess RNA qualitatively after isolation and, subsequently, after biotin-labelling and after fragmentation. and stained with streptavidin-phycoerythrin (Stain Buffer, 2 mg/ml acetylated BSA and 10 µg/ml streptavidin R-phycoerythrin; Molecular Probes, Inc., Eugene, OR), and were scanned on an Affymetrix GCS GeneChip GeneArray scanner. Resulting data was analysed using GCOS, dCHIP, and GeneSpring (Agilent Technologies).

Normalisation and Filtering
Cel files obtained from the GCOS server were processed and normalized by dCHIP [67] algorithm. In this normalisation procedure, an array with median overall intensity is chosen as the baseline array against which other arrays are normalized at probe intensity level. Subsequently, a subset of PM probes, with small within-subset rank difference in the two arrays, serves as the basis for fitting a normalization curve. A filter was designed to include a fold change of at least 1.2 fold between normal and BCC specimens, a difference of at least 100 Affymetrix arbitrary units between normal and BCC average values, and a t-test between normal and BCC (with a p-value cut-off <0.05).

Gene Ontology and Pathway Analysis Analysis
In order to establish which gene ontologies (GO) are overrepresented in our lists of 2,108 significantly up-regulated and 1,813 significantly down-regulated (in BCCs compared to normal skin) transcripts, we compared these to the list of all human genes from the EBI [68], using Gostat [69]. In brief, for all of the gene transcripts analysed, GOstat determines the associated annotated GO terms and all branches/splits on their connection path. The program then counts the number of appearances of each GO term for the gene transcripts in the list being analysed, as well as in the reference list. For each GO term, a Chi-squared pvalue is calculated representing the probability that the observed numbers of counts could have resulted from randomly distributing this GO term between the tested and the reference lists. If the expected value for any analysis is <5, the Chi-squared approximation is considered to be inaccurate. Genmapp [70] was used to identify pathways affected by the differentially-expressed genelist.

Real-time PCR (qPCR)
Following priming with oligo (dT) at 65°C for 5 minutes, followed by 1 minute incubation on ice, cDNA was synthesised from 100 ng total RNA, using Superscript III RNase H-(with increased thermal stability; Invitrogen), RNase OUT Ribonuclease (active against RNase A, B and C; Invitrogen) and a cocktail of dNTPs, by incubating at 50°C for 1 hour, followed by 70°C for 15 minutes, in a 40 µl reaction volume. The cDNA (diluted 1:10), was amplified in 25 µl reactions, by qPCR, using an ABI 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA). Following evaluation of 11 potential endogenous controls in a random selection* of 4 BCC and 4 normal specimens, this study involved evaluation, in all 20 BCC and 5 normal skin specimens, of 5 target transcripts normalised to 2 suitable endogenous controls -β-actin and GAPDH -and calibrated against a pooled cDNA of BCC and normal skin specimens, the relative quantity of which was set to 1. The temperature profile of all reactions was 50°C for 2 minutes 95°C for 10 minutes, 40 cycles of 95°C and 60°C for 1 minute. Individual specimens were analysed in triplicate. [* Note: these 4 BCC and normal specimens were analysed by microarrays with all of the other specimens included in this study].