Purpose: Imiquimod represents a synthetic local immune response modifier that has demonstrated efficacy in clearing basal cell carcinoma. Via interaction with Toll-like receptor 7 on immune cells, imiquimod induces local production of cytokines, such as interferon (IFN)-α.

Experimental Design: To more closely define and elucidate mechanisms leading to basal cell carcinoma clearance in vivo, we examined gene expression profiles of skin basal cell carcinoma before and after treatment with 5% imiquimod cream (Aldara) by using high-density oligonucleotide arrays.

Results: We show that imiquimod predominantly induces genes involved in different aspects of immune response. In addition to effects on immunity, imiquimod treatment modulates the expression of genes involved in the control of apoptosis and oncogenesis. Array data indicated that imiquimod treatment induces expression of opioid growth factor receptor, a molecule recently reported to be a target for antitumor antibody responses. Immunohistochemistry revealed in vivo up-regulation of opioid growth factor receptor protein on tumor and on infiltrating cells after treatment. By using basal cell carcinoma cell lines treated with IFN-α or imiquimod, we show that opioid growth factor receptor up-regulation is IFN-α-mediated, rather then directly imiquimod-mediated. By using tissue microarray containing 52 basal cell carcinomas, we demonstrate opioid growth factor receptor expression in almost half of the cases. Expression of opioid growth factor receptor correlated with a longer recurrence-free period in basal cell carcinoma that recurred after radiotherapy (Kaplan-Meier analysis, P = 0.041).

Conclusions: In addition to its immunomodulatory and antiproliferative activity, opioid growth factor receptor seems to have a prognostic significance in basal cell carcinoma patients. Our data add to the growing list of basal cell carcinoma-associated tumor antigens.

Basal cell carcinoma is the most common malignancy in the white human population worldwide, and its incidence has been increasing in the last decade (1). Basal cell carcinoma tends to grow slowly and rarely metastasizes. Basal cell carcinoma can be successfully treated by intralesional injections of interferons (IFNs) or by topical application of imiquimod (1). Imiquimod has been tested and found to be efficacious in clearing basal cell carcinoma [basal cell carcinoma (2, 3, 4)]. Imiquimod [1-(2-methylpropyl)-1H-imidazo[4,5-c]quinolin-4-amine] belongs to the group of imidazoquinolones, synthetic local immune response modifiers that have demonstrated potent antiviral and antitumor activity (5). Imiquimod activates macrophages and other monocyte-derived immune cells via interaction with Toll-like receptor-7 (6) inducing the local production of cytokines such as IFN-α, tumor necrosis factor α, and interleukin (IL)-12, resulting in an enhanced innate immune response (5, 7). In addition, imiquimod induces migration and activation of skin Langerhans cells (8), all of which promote the biasing of naïve T cells to Th1 cells, which are potent producers of IFN-γ, thus leading to the enhancement of adaptive immunity. The stimulation of innate and adaptive immunity may finally lead to regression of viral-induced lesions and neoplasms (5, 7).

A functional genomic approach offers a unique possibility to define and analyze gene expression profiles in response to a specific drug treatment. The implementation of high-density oligonucleotide arrays containing more than 10,000 genes on a single array has enabled such a genome-wide assessment of changes in the overall expression profile (9). Using this approach, we have analyzed in vivo gene expression profiles of basal cell carcinomas treated with 5% imiquimod cream. In addition, we established basal cell carcinoma cell lines from skin lesions, which we treated with imiquimod or IFN-α, and their gene expression profiles were then compared with the ones from skin basal cell carcinoma after imiquimod therapy in vivo.

Patient Samples

This substudy used material from a 3M Pharmaceuticals-supported study (10) investigating infiltrate alterations in skin basal cell carcinoma after treatment with 5% imiquimod cream (Aldara) approved by Institutional Ethnics Committee. Before entering the study, all patients provided written informed consent for both studies. At the screening visit, a biopsy specimen was taken for histological confirmation of skin basal cell carcinoma. Patients applied the cream once daily 5 times/week for a maximum of 6 weeks. At the point when the tumor began to show signs of erosion, the tumor was surgically excised. The surplus material was available from five patients, one of whom had a second basal cell carcinoma lesion, which was treated on signed informed consent after completion of the study. An additional patient with M. Bowen was treated according to the same scheme, on signed informed consent.

Cell Lines

Establishment of Primary Basal Cell Carcinoma Cell Cultures.

Four primary basal cell carcinoma cell cultures were established according to the previously described methodology (11). These basal cell carcinoma cell lines exhibited a low culture-doubling rate, ranging from 14 to 30 days. Basal cell carcinoma cell lines expressed keratin markers detected by two antibodies: Lu-5 (BMA Biomedicals AG, Augst, Switzerland) and NCL-L-AE1/AE3 (Medite AG, Nunningen, Switzerland). They were negative for expression of vimentin, neuron-specific enolase, and S-100. To remove contaminating keratinocytes, the primary basal cell carcinoma cultures were exposed to “calcium shift” (11). The increase in calcium concentration in growth medium led to their cornification in the next 96 h, whereas basal cell carcinoma cells flattened, elongated, and stratified but neither cornified nor died (11).

Immunomodulatory Treatment.

The HCl salt form of imiquimod was provided by 3M Pharmaceuticals (St. Paul, MN). Primary basal cell carcinoma cell lines and human keratinocytes were stimulated with 15 μg/ml imiquimod or with 1000 units/ml IFN-α2a (Roferon-A; Roche) for 24 h. Primary human keratinocytes were obtained from foreskins as described elsewhere (12) and grown in defined keratinocyte serum-free medium (Gibco; Invitrogen AG, Basel, Switzerland) containing bovine pituitary extract and recombinant epidermal growth factor.

To evaluate the viability of the cells after the aforementioned treatments, the cells were stained with propidium iodine (Fluka Chemie GmbH, Buchs, Switzerland) and counted under an inverted epifluorescence microscope. Neither imiquimod nor IFN-α treatment induced a significant decrease in cell numbers after 24 h (37.25 ± 8.18 cells/cm2 before treatment versus 40.50 ± 11.09 cells/cm2 after imiquimod treatment or 38.75 ± 14.5 cells/cm2 after IFN-α treatment).

Gene Expression Profiling

Total RNA was isolated using TRIzol reagent (Invitrogen AG). Double-stranded cDNA was generated using a Superscript cDNA synthesis kit (Invitrogen AG) using an oligo(dT)24 primer containing a T7 RNA polymerase promoter at the 3′ end (Microsynth, Balgach, Switzerland). Labeled cRNA was prepared from double-stranded cDNA by in vitro transcription using a T7 polymerase [MEGAscript T7 kit; Ambion (Europe) Ltd.] in the presence of biotin-11-CTP and biotin-16-UTP (Enzo Diagnostics, Farmingdale, NY) and purified using RNeasy columns (Qiagen AG, Basel, Switzerland). Fifteen μg of biotinylated cRNA were fragmented and hybridized to HG-U95Av2 GeneChip arrays (Affymetrix, Santa Clara, CA) that contain probe sets representing ∼12,000 genes. Chip hybridization, washing, and staining were performed according to the Affymetrix-recommended protocols.

After scanning of the probe arrays, digitalized image data were processed using Affymetrix Microarray Suite 5.0 software. In addition to determining the expression level for each gene (the so-called “signal”), this software assigns a detection call [absent (not expressed), present (expressed), or marginal (marginally expressed)] for each gene measured. The overall fluorescence intensity was scaled to a global intensity of 500 to enable the comparison between chips.

Data Analysis

The expression data for all genes were exported to Microsoft Excel. The following criteria were used to select differentially expressed genes on imiquimod (or IFN-α) treatment: (a) paired Student’s t test significance level of P < 0.05; (b) fold change cutoff of 1.5 in both directions; and (c) the present call(s) had to coincide with the regulation in question, e.g., if a candidate gene were to be called “up-regulated” after immunomodulatory treatment, at least one sample in the treated group had to be called present; vice versa, if a gene were to be called “down-regulated” after immunomodulatory treatment, at least one sample in the untreated group had to have a present call. The selected differentially expressed genes were normalized to a mean of 0 and a SD of 1, log2-transformed, and subjected to average linkage hierarchical clustering using the uncentered Pearson similarity matrix (13). Clustering analysis was performed using the GENE CLUSTER program, and the figures were generated with the TREE VIEW program (13).4 The selected genes were annotated using the NetAffx analysis system offered by Affymetrix (14).

Real-Time Polymerase Chain Reaction

Approximately 1 μg of total RNA was reverse transcribed using oligo-p(dT)15 priming and avian myeloblastosis virus reverse transcriptase [1st Strand cDNA Synthesis Kit for Reverse Transcription-PCR (avian myeloblastosis virus); Roche Molecular Biochemicals, GmbH, Mannheim, Germany] at 42°C for 1 h. PCR amplifications were carried out using HotStart system (LightCycler-Faststart DNA Master SYBR Green I; Roche Molecular Biochemicals) in the LightCycler thermocycler, as described previously (15). The primer sets used and cycling conditions are shown in Supplementary Table 1. To generate external standards for subsequent quantification, PCR products were purified and cloned into the pCRII-TOPO vector using TOPO TA Cloning Kit (Invitrogen AG), according to the manufacturer’s recommendations. The incorporation of the insert into the vector was confirmed by sequencing and BLAST comparison.

Protein Isolation and Immunoblotting Analysis

After the RNA extraction with TRIzol, proteins were isolated from the intermediate fraction according to the manufacturer’s recommendation. For immunoblotting analysis, protein extracts were electrophoresed using Novex Tris-Glycine gels (Invitrogen AG) and the buffer system of Laemmli and subsequently transferred to nitrocellulose in 12.5 mm Tris base, 100 mm glycine in 20% methanol transfer buffer. The membranes were blocked with 10 mm Tris, 150 mm NaCl, 5% nonfat powdered milk, and 0.25% Tween 20 and incubated with primary antibody for 60 min. Primary antibodies were used at the following concentrations: opioid growth factor receptor, 0.25 μg/ml; and actin, 0.2 μg/ml (clone I-19-R; rabbit polyclonal IgG; Santa Cruz Biotechnology Inc., Heidelberg, Germany). After incubation with the primary antibody, the membrane was washed extensively with 10 mm Tris-HCl, 150 mm NaCl, and 0.25% Tween 20 and incubated for 45 min with the secondary antibody coupled to horseradish peroxidase. Finally, the blots were developed with chemiluminescence (ECL; Amersham Biosciences, Freiburg, Germany). To reprobe the blots with other antibody, the membrane was stripped in 62.5 mm Tris-HCl, 2% SDS, and 100 mm β-mercaptoethanol for 30 min at 50°C. To remove any remaining SDS, the membrane was washed three times with PBS and blocked again.

Tissue Microarray and Immunohistochemistry

A tissue microarray containing 52 basal cell carcinoma tissue samples was constructed as described previously (16). Paraffin-embedded tissue and tissue microarray sections, stained with H&E, were evaluated by an experienced dermatopathologist (R. D.) and subsequently stained with antihuman opioid growth factor receptor antibody raised against the tandem repeat region of opioid growth factor receptor, a kind gift from Dr. Joseph A. Mollick (Department of Medicine, Division of Oncology, Stanford University Hospital, Stanford, CA). Immunohistochemistry was performed with a 1:30 working antibody dilution using the alkaline phosphatase-anti-alkaline phosphatase technique, as described previously (10).

We generated gene expression profiles for six basal cell carcinomas and one M. Bowen sample before and after treatment with 5% imiquimod cream and for four basal cell carcinoma cell lines and two keratinocyte controls before and after treatment with IFN-α or imiquimod. Based on our filtering criteria, we found the expression of 1305 genes (∼10% of the genes on the HG-U95Av2 chip) to be differentially regulated on treatment with imiquimod of the skin basal cell carcinoma. Of these genes, 372 (28.5%) genes were up-regulated, whereas 933 (71.5%) genes were down-regulated. Two-way clustering analysis, a mathematical approach that essentially organizes the data by grouping genes with a similar expression pattern, yielded two major array clusters (i.e., cell lines and skin samples) with this set of 1305 differentially expressed genes and, within each group, a subcluster relating to a particular type of the treatment (Fig. 1, A and B).

Topical Imiquimod Treatment Preferentially Modulates Genes Involved in Immune Response.

The vast majority of up-regulated genes consisted of genes involved in different aspects of immune response. Given the previous reports showing that imiquimod is a potent IFN-α inducer (17, 18, 19, 20, 21), we first focused on this group of genes. Using hierarchical clustering, we compared the expression profiles of the skin samples before and after imiquimod treatment with basal cell carcinoma cell lines treated with IFN-α or imiquimod. The genes that group together according to this method are likely to be coregulated and functionally related (22). Fig. 1,C shows the results of the clustering analysis, revealing that genes up-regulated in basal cell carcinoma cell lines and keratinocytes by IFN-α treatment solely represent the same group of genes up-regulated by imiquimod treatment of skin basal cell carcinoma. In vivo imiquimod induced the expression of different 2′,5′-oligoadenylate synthetase isoforms, IFN-induced proteins including MxA and MxB protein, as well as STAT1 and STAT2 transcription factors. By comparing the genes that were down-regulated in basal cell carcinoma cell lines treated with IFN-α (and cell line controls) and basal cell carcinoma skin lesions treated with imiquimod, we could define 79 genes expressed in both evaluated groups. Clustering analysis demonstrated that the majority of these genes grouped together, confirming that they belong to the same functional category (Fig. 1 E).

Apart from IFN-α and IFN-γ, imiquimod induced a whole spectrum of cytokines (IL-1β), chemokines (RANTES, ecalectin, MIP-1α, and CXCL11), and cognate receptors [MIP-1α/RANTES receptor (CCR1), CCR7 receptor, IL-3 receptor α, IL-15 receptor α], some of which are presented in Fig. 2 A.

In addition to the IFN-α cluster, we could define a IFN-γ cluster, namely, genes that are up-regulated by IFN-γ as a result of Th1 shift and the activation of T cells. Fig. 1 D shows parts of the “infiltrate” cluster containing some of the IFN-γ-induced genes involved in antigen processing (TAP-1 and MECL-1) and the formation of immunoproteasome (PA28), function of activated T cells (granzyme B, LAG-3, and T-cell immune regulator 1), and autoregulatory apoptosis in T cells (tumor necrosis factor receptor 2).

This cytokine boom is likely responsible for the influx of immune cells observed in histology (Fig. 4, E−G), which is reflected in the up-regulation of genes representing receptors (e.g., killer cell inhibiting receptors and killer cell activating receptors, Toll-like receptors, and HLA molecules) pertinent to these cells (summarized in Fig. 2,B). Topical imiquimod treatment resulted in an enhanced expression of molecules related to dendritic cell development, such as CD1d, flt-3 ligand, and Toll-like receptor-3 (Fig. 2, B and C). Moreover, imiquimod seems to be a very potent inducer of costimulatory molecules such as CD86, CD40, ICAM-2 (CD102), or SLAM (CDw150; Fig. 2,C). Infiltrating T cells seem to be activated and fully armed to kill, which is reflected in the up-regulation of cytotoxic (cytolytic) molecules such as granzyme B and granulysin (Fig. 2,D). There was also an up-regulation of genes involved in respiratory burst and production of reactive oxygen species in phagocytes (neutrophils and macrophages), i.e., neutrophil cytosolic factors 1 and 4 (Fig. 2 D).

Imiquimod Treatment Up-Regulates Opioid Growth Factor Receptor Expression: An Effect of Interferon-α?

One molecule within the infiltrate cluster with no clear functional relation to this cluster drew our attention. We observed posttreatment induction of opioid growth factor receptor gene expression in skin basal cell carcinomas (2.5-fold increase; P = 0.005). Quantitative real-time PCR analysis confirmed microarray results, demonstrating up-regulation of opioid growth factor receptor expression in six of seven treated skin samples (one remained unchanged), as shown in Fig. 3,A. Normalization of opioid growth factor receptor to PTCH gene expression (defined as basal cell carcinoma-specific) again revealed an increase in opioid growth factor receptor expression in five of six tested skin samples (Fig. 3,B). Immunohistochemistry using the antibody raised against the immunogenic opioid growth factor receptor tandem repeat region revealed an impressive opioid growth factor receptor up-regulation in all skin basal cell carcinoma after imiquimod treatment (Table 1 and Fig. 4), including one M. Bowen sample. In addition to tumor cells, infiltrating cells also up-regulated opioid growth factor receptor expression (Fig. 4, D–G). Before treatment, basal cell carcinoma tumor cells demonstrated either weak or no opioid growth factor receptor immunoreactivity (Fig. 4, B and D). Normal epidermis demonstrated a nuclear pattern of opioid growth factor receptor expression in the cells from the basal to the granular cell layer (Fig. 4,A). It is of note that the irregular nests of atypical, neoplastic keratinocytes in M. Bowen exhibited opioid growth factor receptor immunoreactivity before treatment (Fig. 4 C).

To further differentiate whether imiquimod induces opioid growth factor receptor directly or via IFN-α, we analyzed opioid growth factor receptor expression in basal cell carcinoma cell lines (keratinocytes) treated with imiquimod and IFN-α for 24 h by real-time PCR. IFN-α up-regulated opioid growth factor receptor expression in two basal cell carcinoma cell lines (L1 and L2) and keratinocyte controls (t test, P < 0.05), whereas two cell lines demonstrated no significant change in opioid growth factor receptor expression (Fig. 3,C). After treatment with imiquimod, there was no significant change in opioid growth factor receptor expression in any of the cell lines investigated. Immunoblotting analysis of basal cell carcinoma cell culture extracts demonstrated an increase in opioid growth factor receptor protein after IFN-α treatment (Fig. 3 D). We observed no changes in opioid growth factor receptor protein expression after imiquimod treatment of basal cell carcinoma cell lines.

Opioid Growth Factor Receptor Expression in Basal Cell Carcinomas Treated with Radiotherapy Is Associated with Longer Recurrence-Free Time.

Using a tissue microarray containing 27 basal cell carcinomas of aggressive, sclerosing type and 25 basal cell carcinomas of nonaggressive type (14 nodular and 11 superficial basal cell carcinomas), we demonstrate opioid growth factor receptor expression by tumor cells in 23 of 52 (44.2%) evaluated cases (Fig. 4,I). Opioid growth factor receptor tissue microarray immunohistochemistry results are presented in Table 2A. Opioid growth factor receptor was predominantly expressed in aggressive, sclerosing basal cell carcinoma (16 of 23 opioid growth factor receptor-positive cases; Spearman’s ρ = −0.294; P = 0.035). Nonaggressive basal cell carcinomas displayed opioid growth factor receptor expression in seven cases (four nodular and three superficial basal cell carcinomas, respectively). In a group of 36 primary basal cell carcinomas, 13 of which have recurred after radiotherapy, opioid growth factor receptor expression was detected in 12 of 23 basal cell carcinomas that showed no recurrence (Table 2B). The absence of opioid growth factor receptor expression in the primary tumor correlated with the history of recurrence after radiotherapy (Spearman’s ρ = −0.457; P = 0.005). Moreover, primary basal cell carcinomas that did not recur during the follow-up period exhibited stronger opioid growth factor receptor expression than the ones that have recurred, irrespective of the histology type (Mann-Whitney U test, P = 0.020). In patients whose basal cell carcinomas recurred, opioid growth factor receptor expression correlated with a longer recurrence-free period (Kaplan-Meier analysis; log-rank test, P = 0.041).

Imiquimod Treatment Affects the Expression of Various Apoptosis Genes.

In the skin, imiquimod treatment induced a spectrum of genes implicated in apoptotic signaling and the execution of apoptosis (Fig. 2 F). Clustering analysis revealed that the up-regulation of these genes was an infiltrate-specific phenomenon encountered exclusively on imiquimod treatment of the skin basal cell carcinoma. Namely, none of the cell lines investigated (treated with either IFN-α or imiquimod) exhibited a similar pattern of expression.

Apart from induction, imiquimod-treated skin basal cell carcinomas showed down-regulation of additional apoptosis genes, among others, the decrease in bcl-2 homolog expression (Fig. 2,F). Some of these genes organized together according to the two-way hierarchical clustering algorithm (Fig. 1 G). Down-regulation of these genes was not observed in any of the treated cell lines.

The Clearance of Skin Basal Cell Carcinoma Is Associated with the Decrease in Expression of Well-Defined and Potentially Novel Genes Implicated in Basal Cell Carcinoma Oncogenesis.

In the group of down-regulated genes, 63 genes were assigned to be involved in oncogenesis according to the known or predicted function (Fig. 2,E). In the skin, imiquimod treatment resulted in a decrease in expression of Patched (PTCH), GLI1, and GLI3 genes, which are implicated in basal cell carcinoma pathogenesis (23). Real-time PCR analysis showed down-regulation of PTCH in five of six tested samples, GLI1 down-regulation of four of six tested samples, and down-regulation of GLI3 in five of six tested samples (Fig. 5). In addition to these genes, we observed a down-regulation of genes reported to be expressed in other human epithelial cancers, such as frizzled (FZD)-7 (24), and c-erb-B2 (HER2). Because the expression of C-erb-B2 (HER2) has been reported previously in basal cell carcinoma (25), we assessed the expression of FZD7 by real-time PCR. All basal cell carcinoma cases investigated demonstrated FZD7 expression. After imiquimod treatment, there was a down-regulation of FZD7 gene expression in five of six treated skin lesions (Fig. 5).

On the other hand, a group of genes (tyrosinase, Pmel-17/gp100, Melan-A/MART-1, and p15/MAAT-1) involved in melanin biosynthesis and known as melanoma-associated antigens was singled out (Fig. 2,E). The expression of these genes, although not detectable in any of the basal cell carcinoma cell lines or keratinocytes, was down-regulated in skin lesions on imiquimod treatment (Fig. 1 H). Immunohistochemical staining showed that the melanocytes were the cells accounting for tyrosinase, gp100, and MART-1 expression in basal cell carcinoma tissue sections. After imiquimod treatment, the melanocyte density decreased, which is probably the reason for the decrease in expression of these markers (data not shown).

Immune response modifiers are substances that directly or indirectly (through the modification of one or more components of the immunoregulatory network) affect a specific immune function (26). Imidazoquinolone compounds, of which imiquimod is the best characterized to date, are such agents. Imiquimod is found to be efficacious in clearing superficial basal cell carcinoma (2, 3, 4). However, its exact mechanism of action has yet to be fully understood. To more closely define and elucidate the mechanism(s) leading to tumor clearance in vivo, we examined gene expression profiles of basal cell carcinoma of the skin before and after treatment with 5% imiquimod cream by using Affymetrix oligonucleotide arrays. Of a total of 12,000 probe sets, we identified 1,305 that were differentially regulated in the skin basal cell carcinoma on imiquimod treatment. Of these genes, one third were up-regulated after imiquimod application, comprising mostly genes that were functionally involved in the immune response. We first focused on genes modulated by IFN-α, knowing that the regression of viral warts or high-grade intraepithelial neoplasia was associated with increased production of IFN-α, IFN-γ, and 2′,5′-oligoadenylate synthetase(18, 19, 21). In addition to these, imiquimod treatment of the skin basal cell carcinoma induced many other IFN-α-regulated genes [various IFN-induced proteins, STAT1, STAT2, MxA, and MxB protein(27, 28)]. Topical imiquimod treatment resulted in an IFN-γ increase as well and in the induction of molecules preferentially regulated by this IFN type.

Apart from IFNs, imiquimod treatment induced the expression of various cytokines, of which IL-1β and MIP-1α were recently shown to be up-regulated in macrophages treated with S-28463, an imiquimod analog (29). In patients with genital warts, the erythema at the treatment site correlated with IFN-γ and RANTES mRNA expression (18), which were up-regulated in our samples as well.

Early infiltrate developing within 3–5 days after treatment initiation appears to be of lymphoid/myeloid origin, which is reflected in the up-regulation of genes encoding for different receptors and cell surface markers expressed in these cells. According to the expression of genes (e.g., IFN-γ, tumor necrosis factor receptor 2, CD97, ICAM-2, and granzyme B) previously described to be expressed in differentiated Th1 cells (30), imiquimod treatment effectively induced biasing of naïve T cells into Th1 cells. Imiquimod activates macrophages via interaction with Toll-like receptor-7 resulting in the secretion of IL-12 (6, 31), favoring the development of IFN-γ-producing T cells. Exposure to IL-12 results not only in the increased expression of Th1-specific genes but also in the induction of additional genes such as CCR1 (MIP-1α/RANTES receptor), LAG-3, and SLAM(30), whose up-regulation we noted as well on imiquimod treatment. A particularly interesting finding was the up-regulation of CD1d, whose expression was recently shown on dermal dendritic cells and monocyte-derived dendritic cells (32).

Apart from monocytes/macrophages, plasmacytoid dendritic cells represent another important target for Toll-like receptor-7-mediated imiquimod activation (6). Activation of the Toll-like receptors on plasmacytoid dendritic cells results in the expression of costimulatory molecules (e.g., CD40 and CD86) necessary for the appropriate activation of antigen-specific T cells (5). Plasmacytoid cells are characterized by a particular set of cell surface molecules, some of which are quite specific for these cells (IL-3 receptor α) and have been up-regulated after imiquimod treatment in our samples (CD32, CD38, CD40, and CD43). IL-3 promotes survival but does not drive the plasmacytoid dendritic cells to proliferate (33). Flt3 ligand and granulocyte colony-stimulating factor, however, have the capability to increase the number of plasmacytoid dendritic cells in vivo(33, 34). Our data show the up-regulation of flt-3 ligand and granulocyte colony-stimulating factor receptor after imiquimod treatment. In view of a recent report demonstrating that imidazoquinolones can directly induce plasmacytoid dendritic cell maturation (35), it is highly suggestive that in addition to activating these cells, imiquimod treatment induces a microenvironment promoting survival and maturation of plasmacytoid dendritic cells. The arrays that we used (HG-U95Av2) contain only Toll-like receptor-1, -2, -3, -5, and 6 genes. Nevertheless, we observed the increase in Toll-like receptor-2 (myelomonocytic cells, neutrophils, lymphocytes) and -3 [predominantly dendritic cells (36)] transcripts that correlated with the infiltrate characteristics observed by histology.

The posttreatment up-regulation of opioid growth factor receptor in the skin (on both the mRNA and protein level) represents a novel and intriguing finding. Opioid growth factor belongs to the group of endogenous opioid peptides and acts as negative regulator of cell proliferation, in addition to its role in pain modulation/transmission (37). Opioid growth factor action is targeted to DNA synthesis, with no effect on the apoptotic pathway (38). Opioid growth factor exerts its action through opioid growth factor receptor, which differs markedly from classical opioid receptors in terms of its function (growth), distribution (neural and nonneural), subcellular location (nucleus), and competitive inhibitory profile (37). Zagon et al.(39) demonstrated constitutive expression of opioid growth factor receptor in normal skin, suggesting its regulatory role in DNA synthesis in the basal cell layer. In addition, McLaughlin et al.(40) showed decreased opioid growth factor receptor protein expression in squamous cell carcinoma of the head and neck as compared with normal epithelium, unleashing tumor cells from growth control. Opioid growth factor receptor transcriptional levels, however, remained unaffected, which led the authors to conclude that the changes responsible for decreased opioid growth factor receptor expression are most likely due to translational/posttranslational defects (40). Using real-time PCR, we could demonstrate opioid growth factor receptor mRNA expression in skin basal cell carcinomas. Our findings confirmed the opioid growth factor receptor expression pattern described previously by Zagon et al.(39) in normal skin and, furthermore, revealed weak or absent immunoreactivity in basal cell carcinoma tumor cells. After immunomodulatory treatment with imiquimod, we describe up-regulation of opioid growth factor receptor protein in tumor by using the antibody raised against the opioid growth factor receptor tandem repeat region. In addition, we show in vitro that opioid growth factor receptor up-regulation in keratinocytes and basal cell carcinoma cell lines is IFN-α-mediated. These findings correspond to in situ up-regulation of opioid growth factor receptor shown by immunohistochemistry. Two recent reports by Wu et al.(41) and Mollick et al.(42) imply that opioid growth factor receptor protein itself could be immunogenic by demonstrating anti-opioid growth factor receptor responses in cancer patients after immunotherapy. Opioid growth factor receptor protein structure is characterized by the presence of tandem repeats that share homology with the extracellular domain of MUC1 (42). It is conceivable that through up-regulation of opioid growth factor receptor, imiquimod might not only be suppressing tumor growth but also increasing tumor immunogenicity. Induction of opioid growth factor receptor in infiltrating cells might be a part of the local immune response control as well as an additional signal to produce specific cytokines (Th1/Th2 switch) or antibodies analogous to other opioid receptors (43). An additional positive side affect of opioid growth factor receptor up-regulation might be analgesia despite substantial tissue damage that develops locally. In a large series of more than 600 patients with anogenital warts treated with imiquimod, patients preferred imiquimod 5% cream over other therapies due to the lack of associated pain (44). Using a tissue microarray containing a panel of different basal cell carcinoma subtypes, we demonstrate opioid growth factor receptor expression in almost half of the evaluated tumors, primarily those with an aggressive phenotype. On the other hand, basal cell carcinomas that did not recur after radiotherapy show strong opioid growth factor receptor expression irrespective of their histology type. Even in the case of recurrence, opioid growth factor receptor expression had an impact on the length of the recurrence-free time. In addition to its putative immunomodulatory as well as antiproliferative activity, opioid growth factor receptor seems to have prognostic significance in basal cell carcinoma patients.

According to our array data, imiquimod treatment had an impact on the expression of various genes involved in apoptotic pathways. The expression pattern of some of these genes appeared to be part of the control of locally induced immune response. Noteworthy was the posttreatment down-regulation of a bcl-2 homolog in skin samples. Aberrant expression of bcl-2 in basal cell carcinoma was implicated in extended tumor cell survival without influencing cell proliferation (10, 45, 46). A recent study by Schön et al.(47) showed that imiquimod induces apoptosis through bcl-2-dependent translocation of cytochrome c, which is independent of membrane-bound death receptors (Fas, tumor necrosis factor-related apoptosis-inducing ligand, and tumor necrosis factor). Accordingly, we have shown that imiquimod induces down-regulation of bcl-2 protein in tumor cells in vivo, thereby not affecting tumor Fas expression (10).

Imiquimod treatment resulted in the reduction of the tumor volume, thus decreasing the expression of various genes implicated in oncogenesis. We detected and confirmed the down-regulation of genes involved in the hedgehog signaling pathway (PTCH, GLI1, and GLI3) implicated in the pathogenesis of sporadic and hereditary basal cell carcinoma (23). Moreover, we demonstrated expression and down-regulation of another gene involved in Wnt signaling, FZD7, which is oncogenic in other epithelial cancer types (24, 48). Apart from alterations in the sonic hedgehog pathway and some members of the Wnt cascade (23), there is little available information on other genes involved in basal cell carcinoma tumor formation. Therefore, high-throughput microarray technology might be a good tool to define and investigate novel basal cell carcinoma tumor antigens.

In conclusion, our study using high-density oligonucleotide arrays provides for the first time a detailed insight into a complex set of in vivo changes after basal cell carcinoma treatment with imiquimod cream. We have summarized and confirmed the existing literature evidence on imiquimod mechanism of action in one go. Our approach revealed for the first time the up-regulation of opioid growth factor receptor by imiquimod, implying that this effect is likely IFN-α-mediated. Future studies will provide more evidence on the role of opioid growth factor receptor as a tumor antigen in basal cell carcinoma and other cancer types.

Grant support: 3M Pharmaceuticals, Gottfried and Julia Bangerter-Rhyner Foundation, and Functional Genomics Center Zurich.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Note: Supplementary data for this article are available at Clinical Cancer Research online (http://clincancerres.aacrjournals.org).

Requests for reprints: Reinhard Dummer, Department of Dermatology, University Hospital Zurich, Gloriastrasse 31, 8091 Zurich, Switzerland. Phone: 41-1-255-25-07; Fax: 41-1-255-89-88; E-mail: [email protected]

4

http://rana.lbl.gov.

Fig. 1.

Gene expression patterns determined by two-way hierarchical clustering analysis of 1305 genes. Relationships between experimental samples are summarized in the dendrogram, in which the pattern and length of the branches reflect the relatedness of the samples. Red indicates a high relative level of expression, whereas green represents a low relative level of expression. A, within the “cell lines” cluster, two major clusters aredistinguishable: basal cell carcinoma cell lines and keratinocyte controls. One group of samples initially considered as a basal cell carcinoma cell line (L4), which was eventually overgrown by keratinocytes, clustered together with keratinocyte control. Within the “basal cell carcinoma skin samples” array cluster, two major clusters are distinguishable: the samples before and after (red letters) topical imiquimod treatment. B, moreover, one pretreatment M. Bowen sample (no_MB) clustered quite peripherally to the basal cell carcinoma samples, reflecting the different biology of this entity. After topical imiquimod treatment, this difference disappeared, grouping this sample (IMIQ_MB) with other basal cell carcinomas. Interestingly, two basal cell carcinoma samples from the same patient (P5-1 and P5-2) remained grouped together irrespective of the treatment. C, group of genes up-regulated by IFN-α (“IFN-α cluster”). Although IMIQ_P4 clustered together with other basal cell carcinoma samples treated with 5% imiquimod cream, this sample failed to up-regulate IFN-α-inducible genes. D, genes up-regulated by IFN-γ belonging to a bigger “infiltrate cluster.” The infiltrate cluster is composed of genes characteristic for immune cells, which are up-regulated on topical imiquimod treatment. E, genes down-regulated by IFN-α. Most of these genes are encoding for ribosomal proteins. F, a cluster of apoptosis genes down-regulated in skin samples after imiquimod treatment. G, the down-regulation of melanocytic markers after topical imiquimod treatment.

Fig. 1.

Gene expression patterns determined by two-way hierarchical clustering analysis of 1305 genes. Relationships between experimental samples are summarized in the dendrogram, in which the pattern and length of the branches reflect the relatedness of the samples. Red indicates a high relative level of expression, whereas green represents a low relative level of expression. A, within the “cell lines” cluster, two major clusters aredistinguishable: basal cell carcinoma cell lines and keratinocyte controls. One group of samples initially considered as a basal cell carcinoma cell line (L4), which was eventually overgrown by keratinocytes, clustered together with keratinocyte control. Within the “basal cell carcinoma skin samples” array cluster, two major clusters are distinguishable: the samples before and after (red letters) topical imiquimod treatment. B, moreover, one pretreatment M. Bowen sample (no_MB) clustered quite peripherally to the basal cell carcinoma samples, reflecting the different biology of this entity. After topical imiquimod treatment, this difference disappeared, grouping this sample (IMIQ_MB) with other basal cell carcinomas. Interestingly, two basal cell carcinoma samples from the same patient (P5-1 and P5-2) remained grouped together irrespective of the treatment. C, group of genes up-regulated by IFN-α (“IFN-α cluster”). Although IMIQ_P4 clustered together with other basal cell carcinoma samples treated with 5% imiquimod cream, this sample failed to up-regulate IFN-α-inducible genes. D, genes up-regulated by IFN-γ belonging to a bigger “infiltrate cluster.” The infiltrate cluster is composed of genes characteristic for immune cells, which are up-regulated on topical imiquimod treatment. E, genes down-regulated by IFN-α. Most of these genes are encoding for ribosomal proteins. F, a cluster of apoptosis genes down-regulated in skin samples after imiquimod treatment. G, the down-regulation of melanocytic markers after topical imiquimod treatment.

Close modal
Fig. 2.

Genes modulated by imiquimod treatment of the skin basal cell carcinoma. Bars represent the posttreatment fold change of the mRNA level of a particular gene when compared with the status before treatment. Black bars indicate up-regulation, whereas white bars represent gene down-regulation. Every gene is represented by gene name and GenBank accession number.

Fig. 2.

Genes modulated by imiquimod treatment of the skin basal cell carcinoma. Bars represent the posttreatment fold change of the mRNA level of a particular gene when compared with the status before treatment. Black bars indicate up-regulation, whereas white bars represent gene down-regulation. Every gene is represented by gene name and GenBank accession number.

Close modal
Fig. 3.

The expression of opioid growth factor receptor mRNA assessed by real-time PCR (LightCycler). A, opioid growth factor receptor expression in skin basal cell carcinoma before and after topical imiquimod treatment. Results are presented as glyceraldehyde-3-phosphate dehydrogenase-normalized ratios to compensate for variations in quantity and quality of starting mRNA. Px, basal cell carcinoma sample not included in the microarray analysis due to unsatisfactory RNA quality; MB, M. Bowen sample. B, opioid growth factor receptor mRNA expression in skin basal cell carcinoma normalized to PTCH gene (basal cell carcinoma-specific) to demonstrate tumor-specific opioid growth factor receptor induction. †, due to the high variability of glyceraldehyde-3-phosphate dehydrogenase copy numbers (up to 20-fold difference), these samples were normalized to β-actin showing less variability. C, changes in opioid growth factor receptor expression in basal cell carcinoma cell lines and keratinocyte controls after IFN-α and imiquimod treatment. D, expression of opioid growth factor receptor protein in basal cell carcinoma cell lines as revealed by immunoblotting analysis. The U937 cell line represents the positive control.

Fig. 3.

The expression of opioid growth factor receptor mRNA assessed by real-time PCR (LightCycler). A, opioid growth factor receptor expression in skin basal cell carcinoma before and after topical imiquimod treatment. Results are presented as glyceraldehyde-3-phosphate dehydrogenase-normalized ratios to compensate for variations in quantity and quality of starting mRNA. Px, basal cell carcinoma sample not included in the microarray analysis due to unsatisfactory RNA quality; MB, M. Bowen sample. B, opioid growth factor receptor mRNA expression in skin basal cell carcinoma normalized to PTCH gene (basal cell carcinoma-specific) to demonstrate tumor-specific opioid growth factor receptor induction. †, due to the high variability of glyceraldehyde-3-phosphate dehydrogenase copy numbers (up to 20-fold difference), these samples were normalized to β-actin showing less variability. C, changes in opioid growth factor receptor expression in basal cell carcinoma cell lines and keratinocyte controls after IFN-α and imiquimod treatment. D, expression of opioid growth factor receptor protein in basal cell carcinoma cell lines as revealed by immunoblotting analysis. The U937 cell line represents the positive control.

Close modal
Fig. 4.

Opioid growth factor receptor immunoreactivity in skin basal cell carcinoma. A, the pattern of opioid growth factor receptor expression in normal epidermis showing clear nuclear opioid growth factor receptor immunoreactivity in the basal, spinous, and part of the granular cell layer (original magnification, ×20). B, cytoplasmic opioid growth factor receptor protein expression in skin basal cell carcinoma (patient P5-2 before treatment). Tumor islets are indicated by arrows. Note the nuclear pattern of opioid growth factor receptor expression in neighboring epidermis, which is not present in tumor islets (original magnification, ×10). C, opioid growth factor receptor expression in Bowen’s disease (original magnification, ×40). D, absence of opioid growth factor receptor expression in skin basal cell carcinoma (patient 4 before treatment; original magnification, ×10). E−G, massive up-regulation of opioid growth factor receptor immunoreactivity in skin basal cell carcinoma after imiquimod treatment (patients Px, P1, and P2, respectively; original magnification, ×10). For better discrimination of the infiltrate, tumor borders are accentuated by a dotted black line. H, nuclear opioid growth factor receptor expression in basal cell carcinoma tumor cells after treatment (original magnification, ×100). I, tissue samples from the basal cell carcinoma tissue microarray stained with antibody against the immunogenic opioid growth factor receptor tandem repeat region (original magnification, ×10).

Fig. 4.

Opioid growth factor receptor immunoreactivity in skin basal cell carcinoma. A, the pattern of opioid growth factor receptor expression in normal epidermis showing clear nuclear opioid growth factor receptor immunoreactivity in the basal, spinous, and part of the granular cell layer (original magnification, ×20). B, cytoplasmic opioid growth factor receptor protein expression in skin basal cell carcinoma (patient P5-2 before treatment). Tumor islets are indicated by arrows. Note the nuclear pattern of opioid growth factor receptor expression in neighboring epidermis, which is not present in tumor islets (original magnification, ×10). C, opioid growth factor receptor expression in Bowen’s disease (original magnification, ×40). D, absence of opioid growth factor receptor expression in skin basal cell carcinoma (patient 4 before treatment; original magnification, ×10). E−G, massive up-regulation of opioid growth factor receptor immunoreactivity in skin basal cell carcinoma after imiquimod treatment (patients Px, P1, and P2, respectively; original magnification, ×10). For better discrimination of the infiltrate, tumor borders are accentuated by a dotted black line. H, nuclear opioid growth factor receptor expression in basal cell carcinoma tumor cells after treatment (original magnification, ×100). I, tissue samples from the basal cell carcinoma tissue microarray stained with antibody against the immunogenic opioid growth factor receptor tandem repeat region (original magnification, ×10).

Close modal
Fig. 5.

Expression of basal cell carcinoma-associated tumor antigens assessed by real-time PCR (LightCycler) before and after topical imiquimod treatment. Results are presented as glyceraldehyde-3-phosphate dehydrogenase-normalized ratios. Px, basal cell carcinoma sample not included in the microarray analysis due to unsatisfactory RNA quality; MB, M. Bowen sample.

Fig. 5.

Expression of basal cell carcinoma-associated tumor antigens assessed by real-time PCR (LightCycler) before and after topical imiquimod treatment. Results are presented as glyceraldehyde-3-phosphate dehydrogenase-normalized ratios. Px, basal cell carcinoma sample not included in the microarray analysis due to unsatisfactory RNA quality; MB, M. Bowen sample.

Close modal
Table 1

OGFR protein expression in basal cell carcinomas treated with topical imiquimod

PatientOGFR TUOGFR INF
P1 before 
P1 after 5 days +++ +++ 
P2 before ± 
P2 after 3 days +++ +++ 
P3 before ± − 
P3 after 3 days +++ +++ 
P4 before − ± 
P4 after 4 days 
P5-1 before − 
P5-1 after 3 days ++ +++ 
P5-2 before ± − 
P5-2 after 3 days ++ +++ 
Px before − − 
Px after 18 days ++ +++ 
MB before 
MB after 5 days +++ +++ 
PatientOGFR TUOGFR INF
P1 before 
P1 after 5 days +++ +++ 
P2 before ± 
P2 after 3 days +++ +++ 
P3 before ± − 
P3 after 3 days +++ +++ 
P4 before − ± 
P4 after 4 days 
P5-1 before − 
P5-1 after 3 days ++ +++ 
P5-2 before ± − 
P5-2 after 3 days ++ +++ 
Px before − − 
Px after 18 days ++ +++ 
MB before 
MB after 5 days +++ +++ 

NOTE. The duration of the treatment is presented in days (time point where the tumor showed signs of erosion and was excised).

Abbreviations: OGFR, opioid growth factor receptor; TU, tumor; INF, infiltrate; −, no expression; ±, single cell positivity; +, 5–35% of the cells positive; ++, 35–70% cells positive; +++, >70% cells positive; Px, basal cell carcinoma sample not included in the microarray analysis due to unsatisfactory RNA quality; MB, M. Bowen sample.

Table 2

OGFR immunohistochemistry results on a BCC tissue microarray

A. OGFR expression with respect to BCC histology in 52 samples in the TMA
OGFR expression*Histology type
Sclerosing BCCNodular BCCSuperficial BCCΣ
No 11 10 29 
Single cell positivity 
Yes 10 15 
Σ 27 11 52 
A. OGFR expression with respect to BCC histology in 52 samples in the TMA
OGFR expression*Histology type
Sclerosing BCCNodular BCCSuperficial BCCΣ
No 11 10 29 
Single cell positivity 
Yes 10 15 
Σ 27 11 52 
B. OGFR expression in primary BCCs that recurred after radiotherapy
RecurrenceOGFR expression*Histology type
Sclerosing BCCNodular BCCSuperficial BCCΣ
No No 11 
 Single cell positivity 
 Yes 
 Σ 23 
Yes No 12 
 Single cell positivity 
 Yes 
 Σ 13 
B. OGFR expression in primary BCCs that recurred after radiotherapy
RecurrenceOGFR expression*Histology type
Sclerosing BCCNodular BCCSuperficial BCCΣ
No No 11 
 Single cell positivity 
 Yes 
 Σ 23 
Yes No 12 
 Single cell positivity 
 Yes 
 Σ 13 

NOTE. The 52 samples in the TMA included 36 primary BCCs and 16 secondary BCCs.

Abbreviations: OGFR, opioid growth factor receptor; BCC, basal cell carcinoma; TMA, tissue microarray.

*

Expression refers to tumor cells.

We acknowledge the expert technical assistance of Christa Dudli and Beatrice Müller in immunohistochemistry. We are indebted to Jessica C. Hassel for her optimized Western blot protocol. We thank Philomena Selvam for excellent work in cloning of the PCR products.

1
Urosevic M, Dummer R. Immunotherapy for nonmelanoma skin cancer: does it have a future?.
Cancer (Phila)
,
94
:
477
-85,  
2002
.
2
Geisse JK, Rich P, Pandya A, et al Imiquimod 5% cream for the treatment of superficial basal cell carcinoma: a double-blind, randomized, vehicle-controlled study.
J Am Acad Dermatol
,
47
:
390
-8,  
2002
.
3
Marks R, Gebauer K, Shumack S, et al Imiquimod 5% cream in the treatment of superficial basal cell carcinoma: results of a multicenter 6-week dose-response trial.
J Am Acad Dermatol
,
44
:
807
-13,  
2001
.
4
Beutner KR, Geisse JK, Helman D, et al Therapeutic response of basal cell carcinoma to the immune response modifier imiquimod 5% cream.
J Am Acad Dermatol
,
41
:
1002
-7,  
1999
.
5
Stanley MA. Imiquimod and the imidazoquinolones: mechanism of action and therapeutic potential.
Clin Exp Dermatol
,
27
:
571
-7,  
2002
.
6
Hemmi H, Kaisho T, Takeuchi O, et al Small anti-viral compounds activate immune cells via the TLR7 MyD88-dependent signaling pathway.
Nat Immunol
,
3
:
196
-200,  
2002
.
7
Hengge UR, Benninghoff B, Ruzicka T, Goos M. Topical immunomodulators: progress towards treating inflammation, infection, and cancer.
Lancet Infect Dis
,
1
:
189
-98,  
2001
.
8
Suzuki H, Wang B, Shivji GM, et al Imiquimod, a topical immune response modifier, induces migration of Langerhans cells.
J Investig Dermatol
,
114
:
135
-41,  
2000
.
9
Lipshutz RJ, Fodor SP, Gingeras TR, Lockhart DJ. High density synthetic oligonucleotide arrays.
Nat Genet
,
21
:
20
-4,  
1999
.
10
Urosevic M, Maier T, Benninghoff B, et al Mechanisms underlying imiquimod-induced regression of basal cell carcinoma in vivo.
Arch Dermatol
,
139
:
1325
-32,  
2003
.
11
Grando SA, Schofield OM, Skubitz AP, et al Nodular basal cell carcinoma in vivo vs in vitro. Establishment of pure cell cultures, cytomorphologic characteristics, ultrastructure, immunophenotype, biosynthetic activities, and generation of antisera.
Arch Dermatol
,
132
:
1185
-93,  
1996
.
12
Kuroki T, Chida K, Hosomi J, Kondo S. Use of human epidermal cells in the study of carcinogenesis.
J Investig Dermatol
,
92
:
271S
-4S,  
1989
.
13
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns.
Proc Natl Acad Sci USA
,
95
:
14863
-8,  
1998
.
14
Liu G, Loraine AE, Shigeta R, et al NetAffx: Affymetrix probesets and annotations.
Nucleic Acids Res
,
31
:
82
-6,  
2003
.
15
Trojan A, Urosevic M, Dummer R, Nestle FO, Stahel RA. Real-time polymerase chain reaction monitoring of epithelial cell adhesion molecule-induced T-cell stimulation in patients with lung cancer and healthy individuals using LightCycler technology.
J Immunother
,
25
:
264
-8,  
2002
.
16
Zagrodnik B, Kempf W, Seifert B, et al Superficial radiotherapy for patients with basal cell carcinoma: recurrence rates, histologic subtypes, and expression of p53 and Bcl-2.
Cancer (Phila)
,
98
:
2708
-14,  
2003
.
17
Weeks CE, Gibson SJ. Induction of interferon and other cytokines by imiquimod and its hydroxylated metabolite R-842 in human blood cells in vitro.
J Interferon Res
,
14
:
81
-5,  
1994
.
18
Tyring SK, Arany I, Stanley MA, et al A randomized, controlled, molecular study of condylomata acuminata clearance during treatment with imiquimod.
J Infect Dis
,
178
:
551
-5,  
1998
.
19
Diaz-Arrastia C, Arany I, Robazetti SC, et al Clinical and molecular responses in high-grade intraepithelial neoplasia treated with topical imiquimod 5%.
Clin Cancer Res
,
7
:
3031
-3,  
2001
.
20
Arany I, Tyring SK, Brysk MM, et al Correlation between pretreatment levels of interferon response genes and clinical responses to an immune response modifier (imiquimod) in genital warts.
Antimicrob Agents Chemother
,
44
:
1869
-73,  
2000
.
21
Arany I, Tyring SK, Stanley MA, et al Enhancement of the innate and cellular immune response in patients with genital warts treated with topical imiquimod cream 5%.
Antiviral Res
,
43
:
55
-63,  
1999
.
22
Lockhart DJ, Winzeler EA. Genomics, gene expression and DNA arrays.
Nature (Lond)
,
405
:
827
-36,  
2000
.
23
Saldanha G, Fletcher A, Slater DN. Basal cell carcinoma: a dermatopathological and molecular biological update.
Br J Dermatol
,
148
:
195
-202,  
2003
.
24
Frierson HF, Jr, El-Naggar AK, Welsh JB, et al Large scale molecular analysis identifies genes with altered expression in salivary adenoid cystic carcinoma.
Am J Pathol
,
161
:
1315
-23,  
2002
.
25
Krahn G, Leiter U, Kaskel P, et al Coexpression patterns of EGFR, HER2, HER3 and HER4 in non-melanoma skin cancer.
Eur J Cancer
,
37
:
251
-9,  
2001
.
26
Straten MV, Lee P, Weitzul S, Cockerell CJ, Tyring SK. Advances in the treatment of basal cell carcinoma: the promise of pharmacologic therapy.
Adv Dermatol
,
16
:
299
-318,  
2000
.
27
Der SD, Zhou A, Williams BR, Silverman RH. Identification of genes differentially regulated by interferon alpha, beta, or gamma using oligonucleotide arrays.
Proc Natl Acad Sci USA
,
95
:
15623
-8,  
1998
.
28
de Veer MJ, Holko M, Frevel M, et al Functional classification of interferon-stimulated genes identified using microarrays.
J Leukocyte Biol
,
69
:
912
-20,  
2001
.
29
Buates S, Matlashewski G. Identification of genes induced by a macrophage activator, S-28463, using gene expression array analysis.
Antimicrob Agents Chemother
,
45
:
1137
-42,  
2001
.
30
Rogge L, Bianchi E, Biffi M, et al Transcript imaging of the development of human T helper cells using oligonucleotide arrays.
Nat Genet
,
25
:
96
-101,  
2000
.
31
Ito T, Amakawa R, Kaisho T, et al Interferon-alpha and interleukin-12 are induced differentially by Toll-like receptor 7 ligands in human blood dendritic cell subsets.
J Exp Med
,
195
:
1507
-12,  
2002
.
32
Gerlini G, Hefti HP, Kleinhans M, et al Cd1d is expressed on dermal dendritic cells and monocyte-derived dendritic cells.
J Investig Dermatol
,
117
:
576
-82,  
2001
.
33
Grouard G, Rissoan MC, Filgueira L, et al The enigmatic plasmacytoid T cells develop into dendritic cells with interleukin (IL)-3 and CD40-ligand.
J Exp Med
,
185
:
1101
-11,  
1997
.
34
Siegal FP, Kadowaki N, Shodell M, et al The nature of the principal type 1 interferon-producing cells in human blood.
Science (Wash DC)
,
284
:
1835
-7,  
1999
.
35
Gibson SJ, Lindh JM, Riter TR, et al Plasmacytoid dendritic cells produce cytokines and mature in response to the TLR7 agonists, imiquimod and resiquimod.
Cell Immunol
,
218
:
74
-86,  
2002
.
36
Medzhitov R. Toll-like receptors and innate immunity.
Nat Rev Immunol
,
1
:
135
-45,  
2001
.
37
Zagon IS, Verderame MF, McLaughlin PJ. The biology of the opioid growth factor receptor (OGFr).
Brain Res Brain Res Rev
,
38
:
351
-76,  
2002
.
38
Zagon IS, McLaughlin PJ. Opioids and the apoptotic pathway in human cancer cells.
Neuropeptides
,
37
:
79
-88,  
2003
.
39
Zagon IS, Wu Y, McLaughlin PJ. The opioid growth factor, [Met5]-enkephalin, and the zeta opioid receptor are present in human and mouse skin and tonically act to inhibit DNA synthesis in the epidermis.
J Investig Dermatol
,
106
:
490
-7,  
1996
.
40
McLaughlin PJ, Stack BC, Jr, Levin RJ, Fedok F, Zagon IS. Defects in the opioid growth factor receptor in human squamous cell carcinoma of the head and neck.
Cancer (Phila)
,
97
:
1701
-10,  
2003
.
41
Wu CJ, Yang XF, McLaughlin S, et al Detection of a potent humoral response associated with immune-induced remission of chronic myelogenous leukemia.
J Clin Investig
,
106
:
705
-14,  
2000
.
42
Mollick JA, Hodi FS, Soiffer RJ, Nadler LM, Dranoff G. MUC1-like tandem repeat proteins are broadly immunogenic in cancer patients.
Cancer Immun
,
3
:
3
2003
.
43
McCarthy L, Wetzel M, Sliker JK, Eisenstein TK, Rogers TJ. Opioids, opioid receptors, and the immune response.
Drug Alcohol Depend
,
62
:
111
-23,  
2001
.
44
O’Mahony C, Law C, Gollnick HP, Marini M. New patient-applied therapy for anogenital warts is rated favourably by patients.
Int J STD AIDS
,
12
:
565
-70,  
2001
.
45
Verhaegh ME, Sanders CJ, Arends JW, Neumann HA. Expression of the apoptosis-suppressing protein Bcl-2 in non-melanoma skin cancer.
Br J Dermatol
,
132
:
740
-4,  
1995
.
46
Delehedde M, Cho SH, Sarkiss M, et al Altered expression of bcl-2 family member proteins in nonmelanoma skin cancer.
Cancer (Phila)
,
85
:
1514
-22,  
1999
.
47
Schön M, Bong AB, Drewniok C, et al Tumor-selective induction of apoptosis and the small-molecule immune response modifier imiquimod.
J Natl Cancer Inst (Bethesda)
,
95
:
1138
-49,  
2003
.
48
http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=OMIM.