<title>Background</title> <p>Human papillomaviruses (HPVs) are double-stranded DNA tumor viruses that infect keratinocytes of the anogenital tract epithelium [##REF##1651607##1##]. Persistent infection by high risk HPV types (e.g., HPV-16, -18, -31, and -45) is the main risk factor for the development of cervical intraepithelial neoplasia and cervical cancer [##REF##14585766##2##,##REF##11754676##3##]. High-risk HPV DNA is detected in more than 90% of cervical carcinomas worldwide [##REF##10451482##4##] and it has been shown that HPV types 16 and 18 can immortalize normal cells in culture, a function that is attributed to E6 and E7 oncogenes [##REF##15479788##5##]. These are the only HPV genes consistently retained and expressed in cervical carcinomas. Besides, their continued expression is required to maintain the malignant phenotype [##REF##16609012##6##, ####REF##15047823##7##, ##REF##11070078##8####11070078##8##]. The proteins encoded by these genes disturb cell proliferation and differentiation by physical and functional interaction with several cellular factors involved in cell cycle regulation [##REF##10793105##9##]. E6 is best known for its ability to bind to p53 and induce its ubiquitin-dependent degradation [##REF##2157286##10##,##REF##8221889##11##], whereas E7 was initially recognized by its ability to interact with members of the retinoblastoma protein family, namely pRb, p107 and p130 [##REF##2537532##12##] and its capacity of enhancing their degradation [##REF##8840974##13##].</p> <p>Persistence of HPV infections and development of neoplasia is influenced by local cell-mediated immune response [##UREF##0##14##]. Tumor necrosis factor-alpha (TNF) is one of the main mediators of skin and mucosa inflammation and has a potent antiproliferative effect on normal primary human keratinocytes (PHKs). This cytokine is a key regulator of diverse inflammatory and immune processes in human epithelia and its expression by keratinocytes is enhanced in response to tissue injury, inflammation, viral infection, and UV radiation [##REF##1712219##15##, ####UREF##1##16##, ##REF##12088639##17####12088639##17##]. Furthermore, TNF has been identified as a key mediator for the regression of HPV-induced lesions [##REF##3925977##18##, ####REF##2829077##19##, ##UREF##2##20##, ##REF##3543677##21####3543677##21##]. Previous studies from our group had shown that TNF exerts a potent cytostatic effect on normal and HPV16 immortalized keratinocytes. On the other hand, keratinocytes immortalized by HPV18 or SV40, as well as HPV16 or HPV18-positive cervical tumor-derived cell lines continue to proliferate normally in the presence of this cytokine [##REF##8625327##22##,##REF##1323969##23##]. In addition, it has been observed that continuous HPV18-gene expression in malignant HeLa-fibroblasts hybrids, as well as increased tumorigenicity of HPV16-transformed human keratinocytes is associated with TNF resistance [##REF##8112896##24##,##REF##8138998##25##]. These observations underscore the importance of TNF-resistance acquisition in HPV mediated pathogenesis and suggest that this event could be an important factor in HPV-associated neoplasia outcome. However, the molecular basis of HPV-mediated TNF resistance has not been elucidated.</p> <p>The aim of the present study was to characterize and compare the global transcription profile of normal and HPV-immortalized keratinocytes. Furthermore, we sought to analyze their response to TNF in order to identify differences that contribute to explain their divergent response to this cytokine. For this purpose, we used microarray analysis to determine transcriptional changes upon 3 and 60 hours after TNF exposure. The 3 hours treatment would favor the identification of immediate early TNF regulated genes. On the other hand, the 60 hours treatment was used because the cytostatic effect exerted by this cytokine on normal and HPV16-immortalized keratinocytes reaches its maximum at this time-point [##REF##8625327##22##,##REF##1323969##23##]. Our experimental setting allowed us to: 1) identify genes that are differentially expressed between TNF-sensitive and TNF-resistant cells; 2) identify genes that are differentially modulated by TNF at two-time points (3 and 60 hours); 3) analyze the effect of HPV-induced immortalization on TNF-regulated genes and, 4) find genes that are differentially expressed between cells immortalized by two different high-risk HPV types. Using this approach, we identified differentially expressed genes that are involved in different cell processes such as immune and inflammatory responses, cell differentiation, cell death, proliferation, extracellular matrix remodeling and DNA repair. The implications of these results are discussed.</p>
<title>Methods</title> <title>Cell Culture and TNF treatment</title> <p>Cultures of primary human keratinocytes (PHK), recovered from newborn foreskin (Cambrex, Walkersville, MD, USA), were maintained in keratinocyte serum-free medium – KSFM (Life Technology, Inc., Gaithersburg, MD, USA) for 3 to seven passages [##REF##2460337##26##]. HF698 and HF18Nco are cell lines obtained from human keratinocytes immortalized by HPV16 and HPV18 whole genome, respectively. These cell lines (from now on referred as HPV16 and HPV18, respectively) were kindly provided by R. Schlegel, Georgetown University Medical Center, Washington, DC [##REF##2556677##27##], were grown in 3+1 medium, consisting of a mixture of 3 parts KSFM and 1 part DMEM, supplemented with 10% fetal calf serum. Cells were grown in 100-mm tissue culture dishes to 30% confluence and treated with 2 nM of human recombinant TNF (Boehringer Mannheim, Germany), for 3 or 60 hours. Cells were then trypsinized, washed 3 times with PBS and frozen until RNA extraction. For all time points, RNA was obtained from two independent experiments, including the control plates.</p> <title>RNA extraction, amplification, labeling, and hybridizations</title> <p>For each sample, total RNA was extracted using TRIzol Reagent (Life Technologies, Inc., Grand Island, NY, USA) following the procedure recommended by the manufacturer. Three micrograms of target and reference (a pool of RNA from all control conditions) total RNA were linearly amplified using T7-based protocol, converted to cyanine-modified cDNA, and labeled as described previously [##REF##14511690##28##].</p> <p>Hybridizations were performed in duplicate, using dye-swap, on a cDNA platform of ORESTES representing 4,600 unique genes with known full-length sequence selected from the clone collection derived from the Human Cancer Genome Project [##REF##15843092##29##]. cDNA amplification, purification, identity verification and printing were performed as previously described [##REF##14511690##28##]. A detailed description of the cDNA microarray platform used and the raw data of this study are available at the GEO website under the accession numbers GPL1930 and GSE4524, respectively [##UREF##3##30##]. Slides were scanned on a confocal laser scanner (Arrayexpress; Packard Bioscience, USA) and, for each spot, signal and background intensities were measured using histogram method of Quantarray software (version 3.0, Packard BioScience, BioChip Technologies LLC, USA).</p> <title>Statistical Analysis</title> <p>Data analysis was performed with R project for statistical computing [##UREF##4##31##] and tools of the associated project, Bioconductor [##UREF##5##32##]. Prior to analysis, signal intensity was corrected by background subtraction, and data normalized by loess method, using span = 0.4 and degree = 2. For the identification of differentially expressed genes, we used ANOVA model when just one variable was considered. For the identification of differentially expressed genes in a pair-wise manner, we used <italic>t</italic>-test and determined the nominal p-value for each individual gene. Those nominal p-values can be conservatively adjusted for multiple testing with the Bonferroni correction by multiplying them by the number of genes in our chip. For clustering samples on the basis of their expression profile, we applied hierarchical clustering based on correlation distance and complete linkage.</p> <title>Northern Blotting</title> <p>For Northern blot analysis, 15 μg of total RNA was fractionated through a 1% denaturing agarose gel and transferred by capillarity onto Hybond N filters (GE Healthcare BioSciences, NJ, USA). Prehybridization, hybridization, and washes were performed as described by Church and Gilbert [##REF##6326095##33##]. The KLK7 and SOD2 cDNA probes were the same used for immobilization in the array. The human GAPDH cDNA probe was used as control for ensuring equal RNA loading. Probes were labeled by random priming, using Ready-To-Go Labeling Beads (GE Healthcare Bio-Sciences, NJ, USA) and [α-<sup>32</sup>P]dCTP (3000 Ci/mmol). Nylon filters were exposed to Kodak Hyperfilm (GE Healthcare BioSciences, NJ, USA) with intensifying screen.</p> <title>Electrophoretic mobility shift assay (EMSA)</title> <p>Nuclear extracts were obtained from monolayer cultures of PHK, and from cell lines HPV16 and HPV18 treated with 2 nM of human recombinant TNF for 1 h. Briefly, cell plates were washed with ice-cold PBS and cells were scraped in 5 ml of PBS. Cells were transferred to a 15 ml Falcon tube and centrifuged at 3000 rpm for 3 min. Cell pellets were ressuspended in 4 ml of lysis buffer (10 mM HEPES pH 7,9, 10 mM KCl, 0,2 mM EDTA, 1 mM DTT), incubated on ice for 5 min, centrifuged and ressuspended in 4 ml of lysis buffer. Nuclei obtained were centrifuged at 2000 rpm for 2 min, ressuspended in 100 μl of extraction buffer (20 mM HEPES pH7,9, 0,42 M NaCl, 2 mM EDTA, 1 mM DTT, 1 mM PMSF, 2 μM pepstatin, 0,6 μM leupeptin, 25 mU/ml aprotinin) and incubated on ice for 30 min. Finally, the samples were centrifuged at 12000 rpm for 15 min at 4°C. The supernatants were stored at -80°C. The protein concentration was determined by the Bradford method (Bio-Rad, CA, USA).</p> <p>For gel retardation the following double-stranded oligonucleotide, corresponding to the NF-κB binding sequence, was used: forward-5'-GCCTGGGAAAGTCCCCTCAACT-3' (Invitrogen, CA, USA) was used. The annealed oligonucleotide was labeled with [γ-<sup>32</sup>]ATP (Amersham, Buckinghamshire, UK; 3,000 Ci/mmol) using TK polynucleotide kinase according to the manufacturer instructions (Biolabs, MA, USA) and purified using Sephadex G50 columns followed by phenol:chloroform extraction and precipitation using 10 μg of salmon sperm DNA as a carrier (Invitrogen, CA, USA). DNA pellets were ressuspended in binding buffer (20 mM HEPES pH 7,9, 20% glycerol, 0,1 M KCl, 2 mM EDTA, 1 mM PMSF, 2 μM pepstatin, 0,6 μM leupeptin, 25 mU/ml aprotinin) to a final concentration of 2,5 fmol/μl. The incorporated radioactivity was quantitated using a LS6500 scintillation counter (Beckman Coulter, CA, USA).</p> <p>The binding of NF-κB was performed in a reaction containing 5 μg of protein extract, 5 μg of BSA, 5 μg of salmon sperm DNA and binding buffer to a final volume of 32 μl on ice. After 10 min, 8 μl of the [γ-32]ATP 5'-end-labeled double-stranded oligonucleotide probe was added, and the incubation was continued for an additional 15 min at 30°C. The DNA-protein complexes were resolved on 4% nondenaturing polyacrylamide gels (29:1 cross-linking ratio), dried, and exposed overnight to X-ray films (Amersham, Buckinghamshire, UK).</p>
<title>Results</title> <p>Glass arrays containing 4.800 cDNA sequences were used in order to determine the effects of HPV infection in human keratinocytes as well as the impact of TNF treatment on global gene expression, in HPV negative or positive cells.</p> <p>In order to identify differentially expressed genes as a function of a unique variable (cell type or TNF-treatment) our dataset was first analyzed by one way ANOVA. The comparisons performed allowed us to determine 1) genes that are differentially expressed between TNF-sensitive and TNF-resistant cells; 2) identify genes that are differentially modulated by TNF at two-time points (3 and 60 hours); 3) analyze the effect of HPV-induced immortalization on TNF-regulated genes and, 4) find genes that are differentially expressed between cells immortalized by two different high-risk HPV types (Figure ##FIG##0##1##).</p> <title>Differentially expressed genes as a function of cell type or TNF treatment</title> <p>The identification of differentially expressed genes, with statistical significance, as a function of cell type was performed by ANOVA. Samples and differentially expressed genes (cutoff p-value &lt;10<sup>-10</sup>) were grouped hierarchically, using correlation distance and complete linkage (Figure ##FIG##1##2##). As it can be observed, normal (PHK) and HPV16-immortalized keratinocytes (HPV16), which are sensitive to TNF cytostatic effect, grouped together while TNF-resistant HPV18-immortalized cell line (HPV18) formed an independent branch. This indicates that TNF-sensitive cell lines share a group of genes which are regulated in a way that clearly differentiate them from the TNF-resistant one. Samples were further clusterized by the time in culture after the last medium change (3 or 60 hs) and finally separated as a function of TNF treatment. This clusterization pattern may reflect differences in cell density and other cultures variables such as nutrients availability or medium conditioning. Initially, all treatments were performed using 30% cell density cultures. As expected, due to TNF cytostatic effect on normal and HPV16-immortalized keratinocytes, cell density at the end of the 60 hours period was different between treated (40–50%) and control cells (70–80%) for these cell lines. On the other hand, both cytokine treated and control HPV18-immortalized cells reached 80–90% cell density by the end of the 60 hours period. Flow-cytometry analysis revealed that the TNF effect on sensitive cells was characterized by the accumulation of cells in the G1-phase of the cell cycle. Conversely, TNF-induced G1-arrest was not observed in HPV18-immortalized keratinocytes [[##REF##8625327##22##,##REF##1323969##23##] and data not shown]. Finally, no differences in cell density were observed for cultures corresponding to 3 hours-treatment group.</p> <p>Among the differentially regulated genes we found some related with inflammatory response (<italic>SOD2, TGFB1, CD44, INHBA, OAS1, SIMP</italic>), epidermal development, differentiation and proliferation (<italic>ADAMST1, RARRES1, CREG, HBP17, MCM2, PRSS1, S100P, CREG1</italic>), proteolysis regulation (<italic>KLK7, PI3, LXN</italic>), and cell adhesion (<italic>CD44, PARVA, PROS1</italic>). The name and function of the genes described are listed in Table ##TAB##0##1##. We next determined the global changes in gene expression as a function of TNF treatment. The name and annotated function of the identified genes that best distinguish samples based on TNF treatment (cutoff p-value &lt;10<sup>-2,9</sup>) are listed in Table ##TAB##1##2##. As expected, many of these genes are involved in the inflammatory response and/or are direct targets of TNF e.g. CCL20, CD44, HLA-F, IL1F9, NFKBIA, INHBA, SOD2, MARCKS, RFX5. Samples and genes were hierarchically clusterized on the basis of their correlation distance using complete linkage (Figure ##FIG##2##3##). Samples from TNF-treated normal keratinocytes (PHK) grouped together and apart from the others. HPV-positive samples exhibited a complex clusterization pattern suggesting that the presence of either HPV16 or 18 has an impact on TNF-regulated gene expression. Furthermore, the grouping of treated PHKs apart from the other samples could reflect the fact that PHKs are the only normal cells used in this study and, as such, the only cell type expected to have an unaltered TNF-signaling network. This could contribute to explain the differences in gene expression upon TNF treatment observed between normal and HPV-immortalized keratinocytes.</p> <p>Since our experimental setting included the comparative analysis of global gene expression at two time points (3 or 60 hours), we searched for genes that best differentiate our samples as a function of time. We found 48 genes that clearly differentiate samples from the analyzed time points. The name and annotated function of the identified genes (cutoff p-value &lt;10<sup>-9</sup>) are listed in additional file ##SUPPL##0##1##. Hierarchical clusterization divided samples in two main branches (additional file ##SUPPL##1##2##). Each branch was exclusively composed of samples from the same time point, namely, 3 or 60 hours. Samples from the 3 hours-time point formed a secondary branch that divided normal from HPV-immortalized keratinocytes. On the other hand, samples from the 60 hours-time point formed a secondary branch that divided normal and HPV16-immortalized keratinocytes (TNF-sensitive samples) from HPV18-immortalized keratinocytes (TNF-resistant samples).</p> <title>Differentially expressed genes between TNF-sensitive and TNF-resistant cells</title> <p>In order to identify differentially expressed genes between specific samples we performed a series of pair-wise comparisons. For each pair-wise comparison, we generated a list of differentially expressed genes with p-value lower than 0,01. The complete list of all pair-wise comparisons performed is presented in additional file ##SUPPL##2##3##. We next aimed to characterize genes that were differentially expressed between TNF-sensitive (PHK and HPV16) and TNF-resistant cells (HPV18). To achieve this goal we selected the thirty genes with the lowest p-values that best distinguish both PHK and HPV16 from HPV18 and the thirty genes with lowest p-values that best distinguish both PHK+TNF and HPV16+TNF from HPV18+TNF (considering treatment with TNF for 3 h). Twelve genes were common to both lists giving a total of 48 different genes identified (Table ##TAB##2##3##). Using the expression profile of these 48 genes, samples were grouped hierarchically, based on their correlation distance and complete linkage (Figure ##FIG##3##4##).</p> <p>Using this approach we observed that genes involved with cell cycle control (CCNA2, CDCA2, CDK2AP1), epidermis development, differentiation and proliferation (KLK7, ALDH3A2, PI3, APG12L, BCLAF1, DEK, MAPRE1, S100P, RRAGA, SFRP1), protein ubiquitination (APPBP1) and cell adhesion (BOC, PROS1, SDCBP, THBS1, JPH3), among others, were differentially expressed between TNF-sensitive and TNF-resistant cells (Table ##TAB##2##3##). These analyses were also performed considering TNF treatment for 60h (available as additional files ##SUPPL##3##4## and ##SUPPL##4##5##).</p> <title>Validation on KLK7 and SOD2 as differentially expressed genes</title> <p>We identified a group of genes whose differential expression could be associated with the differential response to TNF of the cell lines studied, namely: KLK7 (kallikrein 7), SOD2 (<italic>superoxide dismutase 2</italic>), S100P (<italic>S100 calcium binding protein P</italic>), PI3 (<italic>protease inhibitor 3, skin-derived</italic>), CSTA (<italic>cystatin A</italic>), RARRES1 (<italic>retinoic acid receptor responder 1</italic>), and LXN (<italic>latexin</italic>). Based on the reported function as well as the expression profile observed, KLK7 and SOD2 genes were selected for further analysis. The expression pattern of these genes observed by microarray was confirmed by Northern Blot in control and TNF-treated (60 hours) samples from all cell lines used (Figures ##FIG##4##5A## and ##FIG##4##5B##). As it can be observed, KLK7 is equally expressed in TNF-treated or untreated HPV18-immortalized cells but is not detected in PHK or HPV16-immortalized cells, even after cytokine treatment. On the other hand, we observed that SOD2 expression is up-regulated by TNF in both PHK and HPV16-immortalized cells but not in HPV18-immortalized cells, confirming the data obtained by microarray (Figures ##FIG##4##5A## and ##FIG##4##5B##).</p> <title>NF-κB is differentially activated in HPV-16- and HPV-18-infected cells</title> <p>It has been reported that NF-κB activation plays an important role in SOD2 induction by TNF. So we hypothesized that the differential expression of SOD2 could be due to the presence of different levels of activated NF-κB after TNF treatment between TNF-sensitive and TNF resistant cells. In order to address this hypothesis NF-κB activation was determined by electrophoretic mobility shift assay (EMSA) using specific oligonucleotides and nuclear protein extracts. Interestingly, we observed that normal as well as HPV16-immortalized keratinocytes exhibited a clear activation of NF-κB as shown by the increase of this factor levels in nuclear protein extracts after TNF treatment (Figure ##FIG##5##6##). On the other hand, NF-κB activation in TNF-resistant HPV18-immortalized cells was below the level of detection (Figure ##FIG##5##6##, lanes 9 and 10). This prompted us to analyze if NF-κB activation was also altered in other HPV-positive cell lines previously reported to be resistant to TNF cytostatic effect [[##REF##8625327##22##,##REF##1323969##23##], and data not shown]. To address this issue we performed EMSA using nuclear protein extracts obtained from HPV16-positive (SiHa) or HPV18-positive (HeLa and SW756) cervical cancer derived cell lines cultures. We observed that TNF-resistant cells exhibited reduced NF-κB activation when compared to normal PHK (additional file ##SUPPL##5##6##). Altogether, these observations suggest that alteration of TNF-signaling pathway leading to NF-κB activation is a common event in HPV-positive cell lines resistant to this cytokine.</p>
<title>Discussion</title> <p>Production and secretion of inflammatory cytokines are among the main events that take place upon viral infection. These molecules coordinate host cell-mediated immune response by recruiting cellular elements from the immune system and by regulating gene expression on target cells [##REF##17892846##34##,##REF##17291188##35##]. The pleiotropic cytokine TNF is a key regulator of inflammation of the epithelia with a well-documented capacity to induce growth arrest in normal or HPV16-immortalized keratinocytes, mainly in the G0/G1 phase of the cell cycle [##REF##11306566##36##]. Conversely, we have previously reported that HPV18-immortalized and both HPV16 or HPV18-transformed cell lines are resistant to TNF-induced growth arrest [##REF##8625327##22##,##REF##1323969##23##].</p> <p>In order to address the yet unknown molecular bases of this difference we applied cDNA microarray technology to compare the global gene expression profiles of TNF-sensitive normal and HPV16-immortalized keratinocytes with that of TNF-resistant HPV18-immortalized ones. Some limitations of this study are the use of a reduced number of samples and the existence of differences in cell culture conditions which are inherent to our experimental setting, i.e. the differences in cell density at the different time-points described above. However, using this approach we identified a group of genes that clearly distinguish both cells groups (Figure ##FIG##1##2## and Table ##TAB##0##1##). This indicates that TNF-sensitive cell lines share a group of genes which are regulated in a way that clearly differentiate them from the TNF-resistant one.</p> <p>On the other hand, when we analyzed changes in global gene expression as a function of TNF treatment we observed that HPV16 and HPV18 samples could not be distinguished from each other while normal keratinocytes could be readily discriminated (Figure ##FIG##2##3##). This observation suggests that the presence of either HPV16 or 18 has an impact on TNF-regulated gene expression. In line with these observations, several studies have shown that HPV positive cells exhibit impaired TNF pathways [##REF##11934887##37##,##REF##11739693##38##]. Moreover, it has been reported that the effects of TNF on HPV-harboring cells depends on variables as cell type studied, the virus type present and culture conditions (i.e., growth factors availability). This cytokine is capable of inducing the proliferation of HPV16-immortalized human cervical epithelial cells cultures in the absence of growth factors through an autocrine, EGF receptor-dependent, pathway [##REF##7708734##39##]. Besides, TNF can upregulate E6/E7 RNA expression and cyclin-dependent kinase activity in these cells [##REF##10657902##40##]. Conversely, it has been reported that TNF exerts a potent cytostatic effect on HPV16-immortalized keratinocytes while HPV18-immortalized as well as cervical carcinoma-derived HPV-positive cell lines remain unaffected [##REF##8625327##22##,##REF##1323969##23##]. Furthermore, it has been observed that increased tumorigenicity of human keratinocytes transformed by HPV16 is associated with resistance to TNF cytostatic effect [##REF##8112896##24##]. Finally, it was demonstrated that TNF downregulates HPV18 transcription in non-malignant HeLa-fibroblasts hybrids, while viral expression in tumorigenic hybrids segregants as well as in parental HeLa cells remained undisturbed [##REF##8138998##25##]. On the other hand, it has been consistently observed that TNF negatively regulates normal keratinocytes proliferation in monolayer [##REF##8625327##22##,##REF##1323969##23##,##REF##11306566##36##] as well as in organotypic cell cultures [##REF##15464843##41##,##REF##7887441##42##]. Altogether, these data support the notion that acquisition of resistance to TNF by HPV-infected cells may represent an important step towards malignancy.</p> <p>Despite the existence of similarities between the two high-risk HPV types used to generate the cell lines studied, the fact that HPV16 and HPV18 are different viruses that exhibit clear differences in their biological activities must be highlighted. For instance, epidemiological studies have shown that HPV18 is more associated to cervical adenocarcinomas while HPV16 is more prevalent in squamous cell carcinomas [##REF##12556961##43##, ####REF##14648656##44##, ##UREF##6##45####6##45##]. Furthermore, compared to other HPV types HPV18 has been associated with increased transforming potential in cell culture systems and with poorer cancer prognosis at the clinical level [##REF##2460337##26##,##REF##2556677##27##,##REF##1850048##46##,##REF##1847260##47##]. On the other hand, HPV16 exhibits a greater potential to establish persistent infections that can progress to high-grade lesions [##REF##17230525##48##,##REF##12838308##49##]. Although we cannot explain the molecular bases of the differences in gene expression between these cell lines, we believe that this may reflect the divergences that exist between these HPV types.</p> <p>We next searched for genes that best distinguish between TNF-sensitive and TNF-resistant cells by pair-wise comparison both before and after cytokine treatment for 3 or 60 hours. By this means we identified 48 and 52 different genes, respectively, that set apart TNF-sensitive from TNF resistant cells (Figure ##FIG##3##4##, Table ##TAB##2##3##, additional files ##SUPPL##3##4## and ##SUPPL##4##5##). The functional characterization of these genes shows that they are involved in critical cellular processes such as regulation of proliferation, differentiation and cell adhesion. Altogether, the differential expression of these genes may contribute to the differential response to the cytostatic effect of TNF observed in these cells.</p> <p>Two genes, namely KLK7 and SOD2, were selected for further analysis based on their reported function and expression profile (Figure ##FIG##4##5A##). KLK7 expression pattern was validated by Northern blot and showed that it is equally expressed in TNF-treated or untreated HPV18-immortalized cells but is not detected in PHK or HPV16-immortalized cells (Figure ##FIG##4##5B##). Kallikreins are a sub-group of serine proteases with different physiological functions. In humans, kallikreins are encoded by 15 structurally similar, steroid hormone-regulated genes that co-localize to chromosome 19q13.4, representing the largest cluster of contiguous protease genes in the entire genome [##REF##15192120##50##, ####REF##12142373##51##, ##REF##11294823##52####11294823##52##]. These proteins mediate the proteolytic degradation of cohesive intracellular structures associated to epithelial differentiation. Recent data also suggest that kallikreins may be causally involved in carcinogenesis, particularly in tumor metastasis and invasion, and, thus, may represent attractive drug targets to consider for therapeutic intervention [##REF##15192120##50##]. Consistent with our findings, it has been observed that KLK7 expression is up-regulated in cervical tumors as well as in cells lines derived from them. On the other hand, normal keratinocytes express low levels of this protein [##REF##15492784##53##,##REF##15297163##54##]. Furthermore, KLK7 expression has been found up-regulated in breast [##REF##14691584##55##] and ovary tumors [##REF##12738725##56##] and is being considered a new tumor progression marker.</p> <p>The superoxide dismutase 2 (SOD2) expression pattern was also validated by Northern blot (Figure ##FIG##4##5B##). This gene is up-regulated in TNF-sensitive but not in TNF-resistant cells. The superoxide dismutase 2 (SOD2) belongs to a family of enzymes involved in the conversion of superoxide radicals in molecular oxygen. Reactive oxygen metabolites have multifactorial effects on the regulation of cell growth and malignant invasion. Furthermore, numerous in vivo studies have shown that the superoxide dismutases can be highly expressed in aggressive human solid tumors [##REF##14990352##57##, ####REF##15581626##58##, ##REF##12032830##59####12032830##59##].</p> <p>Previous reports have shown that activation of the transcription factor NF-κB is essential for the induction of SOD2 by TNF and IL-1β [##REF##11115408##60##,##REF##11912921##61##]. Here we show that TNF-sensitive cells exhibit higher levels of activated NF-κB than TNF-resistant ones after cytokine treatment (Figure ##FIG##5##6## and additional file ##SUPPL##5##6##). Several studies have shown that NF-κB is a negative regulator of keratinocytes proliferation in the epidermis, and that it plays an important role in cell differentiation and tissue homeostasis [##REF##9482881##62##, ####REF##10887174##63##, ##REF##11287960##64####11287960##64##]. In stratified epithelia NF-κB is found in the cytoplasm of proliferating cells from the basal layer while it is detected in the nuclei of non-proliferating cells from the upper layers. Furthermore, it has been observed that NF-κB superexpression is associated with epidermal hypoplasia while its down-regulation promotes hyperplasia [##REF##9482881##62##]. Overall, these data suggest that alterations in TNF-mediated NF-κB activation pathways can play a role in the development and progression of HPV-associated epithelial and mucosal lesions.</p>
<title>Conclusion</title> <p>Progression of HPV-associated lesions depends on the many alterations caused by this virus in the infected cells. We have identified multiple genes differentially regulated by TNF in HPV16 and HPV18 immortalized keratinocytes. Among them we found KLK7 (<italic>kallikrein 7</italic>), SOD2 (<italic>superoxide dismutase 2</italic>), S100P (<italic>S100 calcium binding protein P</italic>), PI3 (<italic>protease inhibitor 3, skin-derived</italic>), CSTA (<italic>cystatin A</italic>), RARRES1 (<italic>retinoic acid receptor responder 1</italic>), and LXN (<italic>latexin</italic>). The differential expression of the <italic>KLK7 </italic>and <italic>SOD2 </italic>transcripts was further confirmed at the RNA level. Moreover, we present evidence that differential <italic>SOD2 </italic>expression correlates with the levels of NF-κB activation exhibited by TNF-sensitive and TNF-resistant cells.</p> <p>This is the first time that the effect of TNF on global gene expression of normal and HPV-immortalized keratinocytes is addressed at two time points. The thorough analysis of the expression pattern of the identified genes may contribute to the understanding of critical differences between transient and chronic events. Furthermore, it may provide insights of the molecular mechanisms of HPV-induced TNF resistance, contribute to the identification of key functions and pathways associated to specific HPV types and, finally, lead to the identification of new cervical tumor progression markers.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Persistent infection by high risk HPV types (e.g. HPV-16, -18, -31, and -45) is the main risk factor for development of cervical intraepithelial neoplasia and cervical cancer. Tumor necrosis factor (TNF) is a key mediator of epithelial cell inflammatory response and exerts a potent cytostatic effect on normal or HPV16, but not on HPV18 immortalized keratinocytes. Moreover, several cervical carcinoma-derived cell lines are resistant to TNF anti-proliferative effect suggesting that the acquisition of TNF-resistance may constitute an important step in HPV-mediated carcinogenesis. In the present study, we compared the gene expression profiles of normal and HPV16 or 18 immortalized human keratinocytes before and after treatment with TNF for 3 or 60 hours.</p> <title>Methods</title> <p>In this study, we determined the transcriptional changes 3 and 60 hours after TNF treatment of normal, HPV16 and HPV18 immortalized keratinocytes by microarray analysis. The expression pattern of two genes observed by microarray was confirmed by Northern Blot. NF-κB activation was also determined by electrophoretic mobility shift assay (EMSA) using specific oligonucleotides and nuclear protein extracts.</p> <title>Results</title> <p>We observed the differential expression of a common set of genes in two TNF-sensitive cell lines that differs from those modulated in TNF-resistant ones. This information was used to define genes whose differential expression could be associated with the differential response to TNF, such as: <italic>KLK7 </italic>(<italic>kallikrein 7</italic>), <italic>SOD2 </italic>(<italic>superoxide dismutase 2</italic>), <italic>100P </italic>(<italic>S100 calcium binding protein P</italic>), <italic>PI3 </italic>(<italic>protease inhibitor 3, skin-derived</italic>), <italic>CSTA </italic>(<italic>cystatin A</italic>), <italic>RARRES1 </italic>(<italic>retinoic acid receptor responder 1</italic>), and <italic>LXN </italic>(<italic>latexin</italic>). The differential expression of the <italic>KLK7 </italic>and <italic>SOD2 </italic>transcripts was confirmed by Northern blot. Moreover, we observed that <italic>SOD2 </italic>expression correlates with the differential NF-κB activation exhibited by TNF-sensitive and TNF-resistant cells.</p> <title>Conclusion</title> <p>This is the first in depth analysis of the differential effect of TNF on normal and HPV16 or HPV18 immortalized keratinocytes. Our findings may be useful for the identification of genes involved in TNF resistance acquisition and candidate genes which deregulated expression may be associated with cervical disease establishment and/or progression.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>LT conceived the study, participated in its design, was involved in cell culture, microarray studies, gene expression validation and manuscript preparation. EB was involved in co-ordination, EMSA assays and manuscript preparation. GHE, RHJ and EJN were involved in microarray data analyses and statistical modeling. WKM and AELC were involved in RNA amplification and microarray studies. LLV participated in the study design and co-ordination. LFLR conceptualized the study format. All authors have read and approved the final version of the manuscript.</p> <title>Pre-publication history</title> <p>The pre-publication history for this paper can be accessed here:</p> <p><ext-link ext-link-type="uri" xlink:href="http://www.biomedcentral.com/1755-8794/1/29/prepub"/></p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported in part by Fundação de Amparo à Pesquisa do Estado de São Paulo grant 98/14335-2. L. Termini had a PhD fellowship from Fundação de Amparo à Pesquisa do Estado de São Paulo (grant 01/01006-5).</p> <p>We would like to thank Dr. Alex Fiorini, Chamberlein Neto, Mariana Santos and Aline Pacífico for technical assistance, Ana Carolina Quirino Simões and Lucas Fahham for the review of the statistical analysis, and Dr. Ana Paula Lepique and members of the Laboratory of Functional Genomics and the Laboratory of Virology of Ludwig Institute for Cancer research for helpful discussions.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Experimental setup for the analysis of HPV and TNF effects on keratinocytes gene expression</bold>. In order to characterize and compare the global transcription profile of normal and HPV-immortalized keratinocytes and to analyze their response to TNF we used an experimental setting that allowed us to: 1) identify differential expressed genes between normal PHK, HPV16 and HPV18-immortalized keratinocytes (comparisons represented by dashed arrows); 2) identify genes modulated by TNF upon treatment for three and sixty hours (comparisons represented by solid arrows) and; 3) compare the effect of TNF between normal PHK and cells immortalized by two different high-risk HPV types (comparisons represented by round dot arrows).</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Hierarchical grouping based on differentially expressed genes as a function of cell type</bold>. These genes where identified by the ANOVA method and the samples where grouped considering the correlation distance and complete linkage. After sample grouping the genes (p values &lt;10<sup>-10</sup>) were hierarchically grouped by their correlation distances. High gene expression is shown in red, low gene expression is shown in green and black indicates non-differential gene expression. Samples: Primary human keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (PHK_3H, PHK_60H, PHK_3H.TNF, PHK_60H.TNF); HPV16-immortalized keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (HPV16_3H, HPV16_60H, HPV16_3H.TNF, HPV16_60H.TNF); HPV18-immortalized keratinocytes: controls and treated with 3 or 60 hours for TNF, respectively (HPV18_3H, HPV18_60H, HPV18_3H.TNF, HPV18_60H.TNF).</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Hierarchical grouping based on differentially expressed genes as a function of TNF treatment</bold>. These genes where identified by the ANOVA method and the samples where grouped considering the correlation distance and complete linkage. After sample grouping the genes (p values &lt;10<sup>-2,9</sup>) were hierarchically grouped by their correlation distances. High gene expression is shown in red, low gene expression is shown in green and black indicates non-differential gene expression. Samples: Primary human keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (PHK_3H, PHK_60H, PHK_3H.TNF, PHK_60H.TNF); HPV16-immortalized keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (HPV16_3H, HPV16_60H, HPV16_3H.TNF, HPV16_60H.TNF); HPV18-immortalized keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (HPV18_3H, HPV18_60H, HPV18_3H.TNF, HPV18_60H.TNF).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Supervised hierarchical grouping based on differentially expressed genes between normal/HPV16-immortalized keratinocytes and HPV 18-immortalized ones after treatment with TNF for 3 hours</bold>. High gene expression is shown in red, low gene expression is shown in green and black indicates non-differential gene expression. Samples: Primary human keratinocytes: controls and treated for 3 hours with TNF, respectively (PHK_3H, PHK_3H.TNF); HPV16-immortalized keratinocytes: controls and treated for 3 hours with TNF, respectively (HPV16_3H, HPV16_3H.TNF); HPV18-immortalized keratinocytes: controls and treated for 3 hours with TNF, respectively (HPV18_3H, HPV18_3H.TNF).</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Differential expression of KLK7 and SOD2 transcripts</bold>. <bold>A</bold>. Detail of the supervised hierarchical grouping based on differentially expressed genes between normal/HPV16-immortalized keratinocytes and HPV 18-immortalized ones, after treatment with TNF for 60 hours. <bold>B</bold>. Northern blot analysis of KLK7 and SOD2 transcription levels. Arrows indicate the two alternative splicing products of KLK7 in HPV18-immortalized keratinocytes (<italic>GenBank </italic># <ext-link ext-link-type="gen" xlink:href="NM_005046">NM_005046</ext-link>); the SOD2 transcript is induced by TNF in both PHK and HPV16-immortalized cells but not in HPV18-immortalized cells (<italic>GenBank </italic># <ext-link ext-link-type="gen" xlink:href="NM_00636">NM_00636</ext-link>). A probe against GAPDH was used to monitor comparable loading between samples.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Analysis of TNF-induced NF-κB activation in normal and HPV16 or 18-immortalized keratinocytes</bold>. Subconfluent cultures of normal and HPV16 or 18-immortalized keratinocytes treated with 2 nM of TNF for 1 h were used to obtain nuclear protein extracts. For each EMSA reaction, 5 μg of nuclear protein were incubated with 50 fmol of [γ-32P]ATP-labeled double-stranded oligonucleotide and a 50X excess of unlabeled oligonucleotide (lanes 3 and 7). Specificity of binding was further demonstrated by incubation of 1 μg of nuclear protein with the described amount of labeled consensus oligonucleotide and a 50X excess of a labeled oligonucleotide carrying a single-base mutation at the NF-κB binding site (lane 4), and incubation of nuclear extract in the absence of any labeled probe (lane 8). NF-κB DNA binding reactions were carried out as described under "Material and Methods". DNA binding complexes are indicated.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Name and function of the differentially expressed genes that best distinguish samples by cell type variable</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>GENE</bold></td><td align="left"><bold>UniGene ID</bold></td><td align="left"><bold>GENE NAME</bold></td><td align="left"><bold>FUNCTION</bold></td></tr></thead><tbody><tr><td align="left"><bold>ADAMTS1</bold></td><td align="left">Hs.534115</td><td align="left"><italic>disintegrin-like and metalloprotease (reprolysin type)</italic></td><td align="left"><italic>negative regulation of cell proliferation</italic></td></tr><tr><td align="left"><bold>ARHGEF10</bold></td><td align="left">Hs.443460</td><td align="left"><italic>Rho guanine nucleotide exchange factor (GEF) 10</italic></td><td align="left"><italic>GTPase activator activity</italic></td></tr><tr><td align="left"><bold>CD44</bold></td><td align="left">Hs.502328</td><td align="left"><italic>CD44 antigen</italic></td><td align="left"><italic>cell adhesion</italic></td></tr><tr><td align="left"><bold>CITED1</bold></td><td align="left">Hs.40403</td><td align="left"><italic>Cbp/p300-interacting transactivator</italic></td><td align="left"><italic>transcription regulator activity</italic></td></tr><tr><td align="left"><bold>CREG</bold></td><td align="left">Hs.5710</td><td align="left"><italic>cellular repressor of E1A-stimulated genes 1</italic></td><td align="left"><italic>cell proliferation</italic></td></tr><tr><td align="left"><bold>CSTA</bold></td><td align="left">Hs.518198</td><td align="left"><italic>cystatin A (stefin A)</italic></td><td align="left"><italic>cysteine protease inhibitor activity</italic></td></tr><tr><td align="left"><bold>D4S234E</bold></td><td align="left">Hs.518595</td><td align="left"><italic>DNA segment on chromosome 4 (unique) 234 expressed sequence</italic></td><td align="left"><italic>dopamine receptor signaling pathway</italic></td></tr><tr><td align="left"><bold>DHRS3</bold></td><td align="left">Hs.289347</td><td align="left"><italic>dehydrogenase/reductase (SDR family) member 3</italic></td><td align="left"><italic>fatty acid metabolism</italic></td></tr><tr><td align="left"><bold>EPB41L1</bold></td><td align="left">Hs.437422</td><td align="left"><italic>erythrocyte membrane protein band 4.1-like 1</italic></td><td align="left"><italic>structural molecule activity</italic></td></tr><tr><td align="left"><bold>FLJ20105</bold></td><td align="left">Hs.47558</td><td align="left"><italic>FLJ20105</italic></td><td align="left"><italic>regulation of transcription</italic></td></tr><tr><td align="left"><bold>FLJ21511</bold></td><td align="left">Hs.479703</td><td align="left"><italic>FLJ21511</italic></td><td align="left"><italic>unknown function</italic></td></tr><tr><td align="left"><bold>FLJ21616</bold></td><td align="left">Hs.591836</td><td align="left"><italic>FLJ21616</italic></td><td align="left"><italic>regulation of transcription</italic></td></tr><tr><td align="left"><bold>FLJ30525</bold></td><td align="left">Hs.7962</td><td align="left"><italic>FLJ30525</italic></td><td align="left"><italic>unknown function</italic></td></tr><tr><td align="left"><bold>GALNT11</bold></td><td align="left">Hs.647109</td><td align="left"><italic>UDP-N-acetyl-alpha-D-galactosamine</italic></td><td align="left"><italic>transferase activity, transferring glycosyl groups</italic></td></tr><tr><td align="left"><bold>GC20</bold></td><td align="left">Hs.315230</td><td align="left"><italic>translation factor sui1 homolog</italic></td><td align="left"><italic>regulation of translational initiation</italic></td></tr><tr><td align="left"><bold>GLDC</bold></td><td align="left">Hs.584238</td><td align="left"><italic>glycine dehydrogenase</italic></td><td align="left"><italic>glycine metabolism</italic></td></tr><tr><td align="left"><bold>GSR</bold></td><td align="left">Hs.271510</td><td align="left"><italic>glutathione reductase</italic></td><td align="left"><italic>glutathion metabolism</italic></td></tr><tr><td align="left"><bold>HBP17</bold></td><td align="left">Hs.1690</td><td align="left"><italic>fibroblast growth factor binding protein 1</italic></td><td align="left"><italic>regulation of cell proliferation</italic></td></tr><tr><td align="left"><bold>INHBA</bold></td><td align="left">Hs.28792</td><td align="left"><italic>inhibin, beta A</italic></td><td align="left"><italic>cell cycle arrest, negative regulation of immune cell differentiation</italic></td></tr><tr><td align="left"><bold>JPH3</bold></td><td align="left">Hs.592068</td><td align="left"><italic>junctophilin 3</italic></td><td align="left"><italic>unknown function</italic></td></tr><tr><td align="left"><bold>KIAA0368</bold></td><td align="left">Hs.368255</td><td align="left"><italic>KIAA0368</italic></td><td align="left"><italic>ER-associated protein catabolism</italic></td></tr><tr><td align="left"><bold>KLK7</bold></td><td align="left">Hs.151254</td><td align="left"><italic>kallikrein 7 (chymotryptic, stratum corneum)</italic></td><td align="left"><italic>epidermis development, proteolysis and peptidolysis, chymotrypsin activity</italic></td></tr><tr><td align="left"><bold>LCN2</bold></td><td align="left">Hs.204238</td><td align="left"><italic>lipocalin 2 (oncogene 24p3)</italic></td><td align="left"><italic>transporter activity</italic></td></tr><tr><td align="left"><bold>LXN</bold></td><td align="left">Hs.478067</td><td align="left"><italic>latexin</italic></td><td align="left"><italic>enzyme inhibitor activity</italic></td></tr><tr><td align="left"><bold>MAL2</bold></td><td align="left">Hs.201083</td><td align="left"><italic>T-cell differentiation protein 2</italic></td><td align="left"><italic>Unknown function</italic></td></tr><tr><td align="left"><bold>MCM2</bold></td><td align="left">Hs.477481</td><td align="left"><italic>minichromosome maintenance deficient 2</italic></td><td align="left"><italic>cell cycle</italic></td></tr><tr><td align="left"><bold>MGC45400</bold></td><td align="left">Hs.389734</td><td align="left"><italic>transcription elongation factor A (SII)-like 8</italic></td><td align="left"><italic>translation elongation factor activity</italic></td></tr><tr><td align="left"><bold>MGEA5</bold></td><td align="left">Hs.500842</td><td align="left"><italic>meningioma expressed antigen 5 (hyaluronidase)</italic></td><td align="left"><italic>glycoprotein catabolism</italic></td></tr><tr><td align="left"><bold>MYO5B</bold></td><td align="left">Hs.200136</td><td align="left"><italic>acetyl-Coenzyme A acyltransferase 2</italic></td><td align="left"><italic>fatty acid metabolism</italic></td></tr><tr><td align="left"><bold>NMES1</bold></td><td align="left">Hs.112242</td><td align="left"><italic>normal mucosa of esophagus specific 1</italic></td><td align="left"><italic>unknown function</italic></td></tr><tr><td align="left"><bold>NMU</bold></td><td align="left">Hs.418367</td><td align="left"><italic>neuromedin U</italic></td><td align="left"><italic>neuropeptide signaling pathway, digestion</italic></td></tr><tr><td align="left"><bold>NT5E</bold></td><td align="left">Hs.153952</td><td align="left"><italic>5'-nucleotidase, ecto (CD73)</italic></td><td align="left"><italic>DNA metabolism</italic></td></tr><tr><td align="left"><bold>OAS1</bold></td><td align="left">Hs.524760</td><td align="left"><italic>2',5'-oligoadenylate synthetase 1</italic></td><td align="left"><italic>immune response to viral infections</italic></td></tr><tr><td align="left"><bold>ODC1</bold></td><td align="left">Hs.467701</td><td align="left"><italic>ornithine decarboxylase 1</italic></td><td align="left"><italic>polyamine biosynthesis</italic></td></tr><tr><td align="left"><bold>PARVA</bold></td><td align="left">Hs.607144</td><td align="left"><italic>parvin, alpha</italic></td><td align="left"><italic>cell adhesion, actin binding</italic></td></tr><tr><td align="left"><bold>PEX3</bold></td><td align="left">Hs.7277</td><td align="left"><italic>peroxisomal biogenesis factor 3</italic></td><td align="left"><italic>peroxisome organization</italic></td></tr><tr><td align="left"><bold>PI3</bold></td><td align="left">Hs.112341</td><td align="left"><italic>protease inhibitor 3, skin-derived (SKALP)</italic></td><td align="left"><italic>elastase-specific inhibitor</italic></td></tr><tr><td align="left"><bold>PLAU</bold></td><td align="left">Hs.77274</td><td align="left"><italic>plasminogen activator</italic></td><td align="left"><italic>chemotaxis</italic></td></tr><tr><td align="left"><bold>PPGB</bold></td><td align="left">Hs.517076</td><td align="left"><italic>protective protein for beta-galactosidase</italic></td><td align="left"><italic>intracellular protein transport</italic></td></tr><tr><td align="left"><bold>PROS1</bold></td><td align="left">Hs.64016</td><td align="left"><italic>protein S (alpha)</italic></td><td align="left"><italic>cell adhesion, endopeptidase inhibitor activity</italic></td></tr><tr><td align="left"><bold>PRSS11</bold></td><td align="left">Hs.501280</td><td align="left"><italic>protease, serine, 11 (IGF binding)</italic></td><td align="left"><italic>insulin-like growth facto binding, regulation of cell growth</italic></td></tr><tr><td align="left"><bold>RARRES1</bold></td><td align="left">Hs.131269</td><td align="left"><italic>retinoic acid receptor responder</italic></td><td align="left"><italic>negative regulation of cell proliferation</italic></td></tr><tr><td align="left"><bold>RDH-E2</bold></td><td align="left">Hs.170673</td><td align="left"><italic>epidermal retinal dehydrogenase 2</italic></td><td align="left"><italic>oxidoreductase activity</italic></td></tr><tr><td align="left"><bold>RPL15</bold></td><td align="left">Hs.381219</td><td align="left"><italic>ribosomal protein L15</italic></td><td align="left"><italic>protein biosynthesis</italic></td></tr><tr><td align="left"><bold>RUTBC3</bold></td><td align="left">Hs.474914</td><td align="left"><italic>RUN and TBC1 domain containing 3</italic></td><td align="left"><italic>unknown function</italic></td></tr><tr><td align="left"><bold>S100P</bold></td><td align="left">Hs.440880</td><td align="left"><italic>S100 calcium binding protein P</italic></td><td align="left"><italic>cell cycle progression and differentiation</italic></td></tr><tr><td align="left"><bold>SEPT10</bold></td><td align="left">Hs.469615</td><td align="left"><italic>septin 10</italic></td><td align="left"><italic>cell cycle</italic></td></tr><tr><td align="left"><bold>SF3B4</bold></td><td align="left">Hs.516160</td><td align="left"><italic>myotubularin related protein 11</italic></td><td align="left"><italic>inositol or phosphatidylinositol phosphatase activity</italic></td></tr><tr><td align="left"><bold>SIMP</bold></td><td align="left">Hs.475812</td><td align="left"><italic>immunodominant MHC-associated peptides</italic></td><td align="left"><italic>protein amino acid glycosylation</italic></td></tr><tr><td align="left"><bold>SMOC2</bold></td><td align="left">Hs.487200</td><td align="left"><italic>SPARC related modular calcium binding</italic></td><td align="left"><italic>calcium ion binding</italic></td></tr><tr><td align="left"><bold>SOD2</bold></td><td align="left">Hs.487046</td><td align="left"><italic>superoxide dismutase 2</italic></td><td align="left"><italic>age-dependent response to reactive oxygen species, cellular defense response</italic></td></tr><tr><td align="left"><bold>STAF65 (gamma)</bold></td><td align="left">Hs.6232</td><td align="left"><italic>SPTF-associated factor 65 gamma</italic></td><td align="left"><italic>regulation of transcription, DNA- dependent</italic></td></tr><tr><td align="left"><bold>SYTL3</bold></td><td align="left">Hs.436977</td><td align="left"><italic>synaptotagmin-like 3</italic></td><td align="left"><italic>intracellular protein transport</italic></td></tr><tr><td align="left"><bold>TFRC</bold></td><td align="left">Hs.529618</td><td align="left"><italic>transferrin receptor (p90, CD71)</italic></td><td align="left"><italic>endocytosis</italic></td></tr><tr><td align="left"><bold>TGFB1</bold></td><td align="left">Hs.645227</td><td align="left"><italic>transforming growth factor, beta 1</italic></td><td align="left"><italic>cell proloferation</italic></td></tr><tr><td align="left"><bold>TPD52</bold></td><td align="left">Hs.368433</td><td align="left"><italic>tumor protein D52</italic></td><td align="left"><italic>morphogenesis</italic></td></tr><tr><td align="left"><bold>TPX2</bold></td><td align="left">Hs.244580</td><td align="left"><italic>microtubule-associated protein homolog</italic></td><td align="left"><italic>cell proliferation</italic></td></tr><tr><td align="left"><bold>TRIM31</bold></td><td align="left">Hs.493275</td><td align="left"><italic>tripartite motif-containing 31</italic></td><td align="left"><italic>protein ubiquitination, ubiquitin ligase activity</italic></td></tr><tr><td align="left"><bold>YME1L1</bold></td><td align="left">Hs.499145</td><td align="left"><italic>YME1-like 1 (S. cerevisiae)</italic></td><td align="left"><italic>protein catabolism</italic></td></tr><tr><td align="left"><bold>ZNF198</bold></td><td align="left">Hs.644041</td><td align="left"><italic>zinc finger protein 198</italic></td><td align="left"><italic>regulation of transcription, DNA- dependent</italic></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Name and function of the differentially expressed genes that best distinguish samples by TNF treatment variable</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>GENE</bold></td><td align="left"><bold>UniGene ID</bold></td><td align="left"><bold>GENE NAME</bold></td><td align="left"><bold>FUNCTION</bold></td></tr></thead><tbody><tr><td align="left"><bold>ADORA2b</bold></td><td align="left">Hs.167046</td><td align="left"><italic>adenosine A2b receptor</italic></td><td align="left"><italic>activation of MAPK</italic></td></tr><tr><td align="left"><bold>AKAP1</bold></td><td align="left">Hs.463506</td><td align="left"><italic>A kinase (PRKA) anchor protein 1</italic></td><td align="left"><italic>RNA binding</italic></td></tr><tr><td align="left"><bold>BTG2</bold></td><td align="left">Hs.519162</td><td align="left"><italic>BTG family, member 2</italic></td><td align="left"><italic>negative regulation of cell proliferation</italic></td></tr><tr><td align="left"><bold>C3</bold></td><td align="left">Hs.529053</td><td align="left"><italic>complement component 3</italic></td><td align="left"><italic>inflammatory response</italic></td></tr><tr><td align="left"><bold>CCL20</bold></td><td align="left">Hs.75498</td><td align="left"><italic>chemokine (C-C motif) ligand 20</italic></td><td align="left"><italic>inflammatory response</italic></td></tr><tr><td align="left"><bold>CD44</bold></td><td align="left">Hs.502328</td><td align="left"><italic>CD44 antigen</italic></td><td align="left"><italic>cell adhesion</italic></td></tr><tr><td align="left"><bold>cig5</bold></td><td align="left">Hs.17518</td><td align="left"><italic>radical S-adenosyl methionine domain containing 2</italic></td><td align="left"><italic>Catalytic activity</italic></td></tr><tr><td align="left"><bold>CLCA4</bold></td><td align="left">Hs.546343</td><td align="left"><italic>chloride channel, calcium activated, family member 4</italic></td><td align="left"><italic>chloride transport</italic></td></tr><tr><td align="left"><bold>DC-UbP</bold></td><td align="left">Hs.179852</td><td align="left"><italic>dendritic cell-derived ubiquitin-like protein</italic></td><td align="left"><italic>Protein modification</italic></td></tr><tr><td align="left"><bold>FAD104</bold></td><td align="left">Hs.159430</td><td align="left"><italic>fibronectin type III domain containing 3B</italic></td><td align="left"><italic>cell differentiation</italic></td></tr><tr><td align="left"><bold>FLJ21511</bold></td><td align="left">Hs.479703</td><td align="left"><italic>FLJ21511</italic></td><td align="left"><italic>Unknown function</italic></td></tr><tr><td align="left"><bold>FMNL3</bold></td><td align="left">Hs.179838</td><td align="left"><italic>formin-like 3</italic></td><td align="left"><italic>cell organization and biogenesis</italic></td></tr><tr><td align="left"><bold>GFPT2</bold></td><td align="left">Hs.30332</td><td align="left"><italic>glutamine-fructose-6-phosphate transaminase 2</italic></td><td align="left"><italic>carbohydrate biosynthesis</italic></td></tr><tr><td align="left"><bold>HLA-F</bold></td><td align="left">Hs.519972</td><td align="left"><italic>major histocompatibility complex, class I, F</italic></td><td align="left"><italic>antigen presentation, endogenous antigen</italic></td></tr><tr><td align="left"><bold>IL1F9</bold></td><td align="left">Hs.211238</td><td align="left"><italic>interleukin 1 family, member 9</italic></td><td align="left"><italic>inflammatory response</italic></td></tr><tr><td align="left"><bold>INHBA</bold></td><td align="left">Hs.28792</td><td align="left"><italic>inhibin, beta A</italic></td><td align="left"><italic>cell cycle arrest, negative regulation of immune cell differentiation</italic></td></tr><tr><td align="left"><bold>KIAA0303</bold></td><td align="left">Hs.133539</td><td align="left"><italic>microtubule associated serine/threonine kinase family member 4</italic></td><td align="left"><italic>protein kinase activity</italic></td></tr><tr><td align="left"><bold>KIAA1279</bold></td><td align="left">Hs.279580</td><td align="left"><italic>KIAA1279</italic></td><td align="left"><italic>Unknown function</italic></td></tr><tr><td align="left"><bold>LAP3</bold></td><td align="left">Hs.479264</td><td align="left"><italic>leucine aminopeptidase 3</italic></td><td align="left"><italic>Protein metabolism</italic></td></tr><tr><td align="left"><bold>MARCKS</bold></td><td align="left">Hs.75061</td><td align="left"><italic>MARCKS-like 1</italic></td><td align="left"><italic>calmodulin binding, macrophage activation</italic></td></tr><tr><td align="left"><bold>MGAT4B</bold></td><td align="left">Hs.437277</td><td align="left"><italic>mannosyl (alpha-1,3-)-glycoprotein beta-1,4-N-acetylglucosaminyltransferase, isoenzyme B</italic></td><td align="left"><italic>cytokine activity</italic></td></tr><tr><td align="left"><bold>MGC45400</bold></td><td align="left">Hs.389734</td><td align="left"><italic>transcription elongation factor A (SII)-like 8</italic></td><td align="left"><italic>translation elongation factor activity</italic></td></tr><tr><td align="left"><bold>MMP9</bold></td><td align="left">Hs.297413</td><td align="left"><italic>matrix metalloproteinase 9</italic></td><td align="left"><italic>proteolysis and peptidolysis</italic></td></tr><tr><td align="left"><bold>NFKBIA</bold></td><td align="left">Hs.81328</td><td align="left"><italic>nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, alpha</italic></td><td align="left"><italic>cytoplasmic sequestering of NF-kappaB</italic></td></tr><tr><td align="left"><bold>NMES1</bold></td><td align="left">Hs.112242</td><td align="left"><italic>normal mucosa of esophagus specific 1</italic></td><td align="left"><italic>Unknown function</italic></td></tr><tr><td align="left"><bold>OAS1</bold></td><td align="left">Hs.524760</td><td align="left"><italic>2',5'-oligoadenylate synthetase 1</italic></td><td align="left"><italic>immune response to viral infections</italic></td></tr><tr><td align="left"><bold>PLAU</bold></td><td align="left">Hs.77274</td><td align="left"><italic>plasminogen activator</italic></td><td align="left"><italic>chemotaxis</italic></td></tr><tr><td align="left"><bold>RDH-E2</bold></td><td align="left">Hs.170673</td><td align="left"><italic>epidermal retinal dehydrogenase 2</italic></td><td align="left"><italic>oxidoreductase activity</italic></td></tr><tr><td align="left"><bold>RFX5</bold></td><td align="left">Hs.166891</td><td align="left"><italic>regulatory factor X, 5</italic></td><td align="left"><italic>inflammatory response, HLA class II expression</italic></td></tr><tr><td align="left"><bold>RIG-1</bold></td><td align="left">Hs.17466</td><td align="left"><italic>retinoic acid receptor responder (tazarotene induced) 3</italic></td><td align="left"><italic>negative regulation of cell proliferation</italic></td></tr><tr><td align="left"><bold>RIPK2</bold></td><td align="left">Hs.103755</td><td align="left"><italic>receptor-interacting serine-threonine kinase 2</italic></td><td align="left"><italic>inflammatory response</italic></td></tr><tr><td align="left"><bold>SASH1</bold></td><td align="left">Hs.193133</td><td align="left"><italic>SAM and SH3 domain containing 1</italic></td><td align="left"><italic>Negative regulation of cell cycle</italic></td></tr><tr><td align="left"><bold>SDCBP</bold></td><td align="left">Hs.200804</td><td align="left"><italic>syndecan binding protein (syntenin)</italic></td><td align="left"><italic>intracellular signaling cascade, interleukin-5 receptor binding</italic></td></tr><tr><td align="left"><bold>SEC24A</bold></td><td align="left">Hs.211612</td><td align="left"><italic>SEC24 related gene family, member A</italic></td><td align="left"><italic>intracellular protein transport</italic></td></tr><tr><td align="left"><bold>SERPINB2</bold></td><td align="left">Hs.514913</td><td align="left"><italic>encoding serine (or cysteine) proteinase inhibitor, clade B (ovalbumin), member 2</italic></td><td align="left"><italic>anti-apoptosis</italic></td></tr><tr><td align="left"><bold>SF3B4</bold></td><td align="left">Hs.412818</td><td align="left"><italic>myotubularin related protein 11</italic></td><td align="left"><italic>RNA splicing</italic></td></tr><tr><td align="left"><bold>SOD2</bold></td><td align="left">Hs.487046</td><td align="left"><italic>superoxide dismutase 2</italic></td><td align="left"><italic>age-dependent response to reactive oxygen species, cellular defense response</italic></td></tr><tr><td align="left"><bold>TMSB4</bold></td><td align="left">Hs.522584</td><td align="left"><italic>thymosin, beta 4, X-linked</italic></td><td align="left"><italic>cytoskeleton organization and biogenesis</italic></td></tr><tr><td align="left"><bold>VMP1</bold></td><td align="left">Hs.444569</td><td align="left"><italic>transmembrane protein 49</italic></td><td align="left"><italic>Unknown function</italic></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>List of differentially expressed genes that best distinguish TNF-resistant cells (HPV18) from TNF-sensitive cells (PHK and HPV16), in normal culture conditions or upon treatment with TNF for 3 hours</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="left" colspan="2"><bold>(PHK and HPV16) <italic>vs </italic>HPV18</bold></td><td align="left" colspan="2"><bold>(PHK_TNF and HPV16_TNF) <italic>vs </italic>HPV18_TNF</bold></td></tr><tr><td/><td/><td colspan="2"><hr/></td><td colspan="2"><hr/></td></tr><tr><td align="left"><bold>GENE</bold></td><td align="left"><bold>GENE NAME</bold></td><td align="left"><bold>FOLD</bold></td><td align="left"><bold>p VALUE</bold></td><td align="left"><bold>FOLD</bold></td><td align="left"><bold>p VALUE</bold></td></tr></thead><tbody><tr><td align="left"><bold>ABCE1</bold></td><td align="left"><italic>ATP-binding cassette, sub-family E (OABP), member 1</italic></td><td align="left">0.592</td><td align="left">0.001384</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>ACBD5</bold></td><td align="left"><italic>acyl-Coenzyme A binding domain containing 5</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.378</td><td align="left">0.00086</td></tr><tr><td align="left"><bold>ALDH3A2</bold></td><td align="left"><italic>encoding aldehyde dehydrogenase 3 family, member A2</italic></td><td align="left">1.623</td><td align="left">0.00122</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>APG12L</bold></td><td align="left"><italic>APG12 autophagy 12-like (S. cerevisiae)</italic></td><td align="left">1.739</td><td align="left">0.000924</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>APPBP1</bold></td><td align="left"><italic>amyloid beta precursor protein binding protein 1</italic></td><td align="left">0.540</td><td align="left">0.000283</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>ARF4L</bold></td><td align="left"><italic>ADP-ribosylation factor 4-like</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.578</td><td align="left">5.70E-05</td></tr><tr><td align="left"><bold>BCLAF1</bold></td><td align="left"><italic>BCL2-associated transcription factor 1</italic></td><td align="left">0.611</td><td align="left">0.000876</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>BOC</bold></td><td align="left"><italic>brother of CDO</italic></td><td align="left">----</td><td align="left">----</td><td align="left">2.337</td><td align="left">6.00E-06</td></tr><tr><td align="left"><bold>CCNA2</bold></td><td align="left"><italic>cyclin A2</italic></td><td align="left">0.577</td><td align="left">0.000604</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>CDCA2</bold></td><td align="left"><italic>cell division cycle associated 2</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.539</td><td align="left">7.20E-05</td></tr><tr><td align="left"><bold>CDK2AP1</bold></td><td align="left"><italic>CDK2-associated protein 1</italic></td><td align="left">----</td><td align="left">----</td><td align="left">1.523</td><td align="left">0.000148</td></tr><tr><td align="left"><bold>CPSF3</bold></td><td align="left"><italic>cleavage and polyadenylation specific factor 3</italic></td><td align="left">0.586</td><td align="left">0.000466</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>CYP1B1</bold></td><td align="left"><italic>cytochrome P450, family 1, subfamily B, polypeptide 1</italic></td><td align="left">0.499</td><td align="left">0.001386</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>DEK</bold></td><td align="left"><italic>DEK oncogene</italic></td><td align="left">0.480</td><td align="left">0.000278</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>FAM31C</bold></td><td align="left"><italic>family with sequence similarity 31, member C</italic></td><td align="left">----</td><td align="left">----</td><td align="left">2.548</td><td align="left">0.000663</td></tr><tr><td align="left"><bold><underline>FLJ20105</underline></bold></td><td align="left"><italic>hypothetical protein LOC54821</italic></td><td align="left">0.026</td><td align="left">5.00E-06</td><td align="left">0.029</td><td align="left">2.00E-06</td></tr><tr><td align="left"><bold><underline>GALNAC4S-6ST</underline></bold></td><td align="left"><italic>B cell RAG associated protein</italic></td><td align="left">1.827</td><td align="left">0.000145</td><td align="left">1.673</td><td align="left">0.000215</td></tr><tr><td align="left"><bold>H105E3</bold></td><td align="left"><italic>encoding NAD(P) dependent steroid dehydrogenase-like</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.569</td><td align="left">6.80E-05</td></tr><tr><td align="left"><bold>HLCS</bold></td><td align="left"><italic>holocarboxylase synthetase</italic></td><td align="left">----</td><td align="left">----</td><td align="left">1.551</td><td align="left">0.00011</td></tr><tr><td align="left"><bold>JPH3</bold></td><td align="left"><italic>junctophilin 3</italic></td><td align="left">0.154</td><td align="left">0.000192</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>KIAA0795</bold></td><td align="left"><italic>kelch-like 18 (Drosophila)</italic></td><td align="left">0.857</td><td align="left">0.001325</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>KIAA1023</bold></td><td align="left"><italic>IQ motif containing E</italic></td><td align="left">1.575</td><td align="left">0.000723</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>KIF1B</bold></td><td align="left"><italic>kinesin family member 1B</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.570</td><td align="left">0.000632</td></tr><tr><td align="left"><bold><underline>KLK7</underline></bold></td><td align="left"><italic>encoding kallikrein 7 (chymotryptic, stratum corneum)</italic></td><td align="left">0.421</td><td align="left">0.000416</td><td align="left">0.374</td><td align="left">2.30E-05</td></tr><tr><td align="left"><bold>LCN2</bold></td><td align="left"><italic>lipocalin 2 (oncogene 24p3)</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.216</td><td align="left">0.000686</td></tr><tr><td align="left"><bold>LOC151242</bold></td><td align="left"><italic>protein phosphatase 1, regulatory (inhibitor)</italic></td><td align="left">1.928</td><td align="left">0.000268</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>Lrp2bp</bold></td><td align="left"><italic>low density lipoprotein receptor-related protein binding protein</italic></td><td align="left">0.629</td><td align="left">0.000574</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold><underline>MAPRE1</underline></bold></td><td align="left"><italic>encoding microtubule-associated protein, RP/EB family, member 1</italic></td><td align="left">0.410</td><td align="left">0.001068</td><td align="left">0.503</td><td align="left">4.20E-05</td></tr><tr><td align="left"><bold>MBD2</bold></td><td align="left"><italic>methyl-CpG binding domain protein 2</italic></td><td align="left">0.680</td><td align="left">0.001247</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>MGC35048</bold></td><td align="left"><italic>hypothetical protein MGC35048</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.499</td><td align="left">0.000211</td></tr><tr><td align="left"><bold>MRPS6</bold></td><td align="left"><italic>mitochondrial ribosomal protein S6</italic></td><td align="left">----</td><td align="left">----</td><td align="left">1.460</td><td align="left">0.000161</td></tr><tr><td align="left"><bold><underline>MYO5B</underline></bold></td><td align="left"><italic>acetyl-Coenzyme A acyltransferase 2</italic></td><td align="left">0.353</td><td align="left">0.000104</td><td align="left">0.292</td><td align="left">4.10E-05</td></tr><tr><td align="left"><bold><underline>NMES1</underline></bold></td><td align="left"><italic>normal mucosa of esophagus specific 1</italic></td><td align="left">0.324</td><td align="left">0.000294</td><td align="left">0.244</td><td align="left">7.00E-06</td></tr><tr><td align="left"><bold>NPR2</bold></td><td align="left"><italic>encoding natriuretic peptide receptor B/guanylate cyclase B</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.435</td><td align="left">0.000234</td></tr><tr><td align="left"><bold>ODC1</bold></td><td align="left"><italic>ornithine decarboxylase 1</italic></td><td align="left">1.660</td><td align="left">0.000886</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold><underline>PI3</underline></bold></td><td align="left"><italic>protease inhibitor 3, skin-derived (SKALP)</italic></td><td align="left">0.213</td><td align="left">0.000274</td><td align="left">0.207</td><td align="left">1.80E-05</td></tr><tr><td align="left"><bold>PROS1</bold></td><td align="left"><italic>protein S (alpha)</italic></td><td align="left">1.807</td><td align="left">0.000633</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold><underline>PTP4A1</underline></bold></td><td align="left"><italic>protein tyrosine phosphatase type IVA, member 1</italic></td><td align="left">0.478</td><td align="left">0.000385</td><td align="left">0.497</td><td align="left">1.80E-05</td></tr><tr><td align="left"><bold>RRAGA</bold></td><td align="left"><italic>Ras-related GTP binding A</italic></td><td align="left">----</td><td align="left">----</td><td align="left">1.510</td><td align="left">0.000301</td></tr><tr><td align="left"><bold><underline>RUTBC3</underline></bold></td><td align="left"><italic>RUN and TBC1 domain 3</italic></td><td align="left">0.518</td><td align="left">0.000704</td><td align="left">0.450</td><td align="left">0.000131</td></tr><tr><td align="left"><bold><underline>S100P</underline></bold></td><td align="left"><italic>S100 calcium binding protein P</italic></td><td align="left">0.101</td><td align="left">0</td><td align="left">0.102</td><td align="left">0</td></tr><tr><td align="left"><bold>SDCBP</bold></td><td align="left"><italic>syndecan binding protein (syntenin)</italic></td><td align="left">0.456</td><td align="left">0.000395</td><td align="left">----</td><td align="left">----</td></tr><tr><td align="left"><bold>SFRP1</bold></td><td align="left"><italic>secreted frizzled-related protein 1</italic></td><td align="left">----</td><td align="left">----</td><td align="left">0.502</td><td align="left">5.10E-05</td></tr><tr><td align="left"><bold>SLC35B3</bold></td><td align="left"><italic>solute carrier family 35, member B3</italic></td><td align="left">----</td><td align="left">----</td><td align="left">1.990</td><td align="left">0.000304</td></tr><tr><td align="left"><bold><underline>STAF65 (gamma)</underline></bold></td><td align="left"><italic>SPTF-associated factor 65 gamma</italic></td><td align="left">0.021</td><td align="left">1.00E-06</td><td align="left">0.028</td><td align="left">0</td></tr><tr><td align="left"><bold>THBS1</bold></td><td align="left"><italic>thrombospondin 1</italic></td><td align="left">----</td><td align="left">----</td><td align="left">2.712</td><td align="left">0.000225</td></tr><tr><td align="left"><bold>VMP1</bold></td><td align="left"><italic>likely ortholog of rat vacuole membrane protein 1</italic></td><td align="left">----</td><td align="left">----</td><td align="left">1.547</td><td align="left">0.000899</td></tr><tr><td align="left"><bold><underline>YME1L1</underline></bold></td><td align="left"><italic>YME1-like 1 (S. cerevisiae)</italic></td><td align="left">0.373</td><td align="left">2.60E-05</td><td align="left">0.333</td><td align="left">5.00E-06</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Table with name and function of the differentially expressed genes that best distinguish samples by time variable. The cutoff p-value was set as &lt;10<sup>-9</sup>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Hierarchical grouping based on differentially expressed genes as a function of time. These genes where identified by the ANOVA method and the samples where grouped considering the correlation distance and complete linkage. After sample grouping the genes were hierarchically grouped by their correlation distances. High gene expression is shown in red, low gene expression is shown in green and black indicates non-differential gene expression, p values &lt;10<sup>-9</sup>. Samples: Primary human keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (PHK_3H, PHK_60H, PHK_3H.TNF, PHK_60H.TNF); HPV16-immortalized keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (HPV16_3H, HPV16_60H, HPV16_3H.TNF, HPV16_60H.TNF); HPV18-immortalized keratinocytes: controls and treated for 3 or 60 hours with TNF, respectively (HPV18_3H, HPV18_60H, HPV18_3H.TNF, HPV18_60H.TNF).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Table of differentially expressed genes obtained by pair-wise comparison of all the samples and their respective expression ratio ordered by p-value. PHK_3H (primary human keratinocytes – control/3 hours), PHK_60H (primary human keratinocytes – control/60 hours), HPV16_3H (HPV16-immortalized keratinocytes – control/3 hours), HPV16_60H (HPV16-immortalized keratinocytes – control/60 hours), HPV18_3H (HPV18-immortalized keratinocytes – control/3 hours), HPV18_60H (HPV18-immortalized keratinocytes – control/60 hours), PHK_TNF3H (primary human keratinocytes – TNF/3 hours), PHK_TNF60H (primary human keratinocytes – TNF/60 hours), HPV16_TNF3H (HPV16-immortalized keratinocytes – TNF/3 hours), HPV16_TNF60H (HPV16-immortalized keratinocytes – TNF/60 hours), HPV18_TNF3H (HPV18-immortalized keratinocytes – TNF/3 hours), HPV18_TNF60H (HPV18-immortalized keratinocytes – TNF/60 hours).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Supervised hierarchical grouping based on the 30 genes that best differentiate normal and HPV 16-immortalized keratinocytes from the HPV 18-immortalized ones after treatment with TNF for 60 hours. High gene expression is shown in red, low gene expression is shown in red and black indicates non-differential gene expression. PHK_3H (primary human keratinocytes – control/3 hours), PHK_3H TNF (primary human keratinocytes – TNF/3 hours), HPV16_3H (HPV16-immortalized keratinocytes – control/3 hours), HPV16_3H TNF (HPV16-immortalized keratinocytes – TNF/3 hours), HPV18_3H (HPV18-immortalized keratinocytes – control/3 hours), HPV18_3H TNF (HPV18-immortalized keratinocytes – TNF/3 hours).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Table of differentially expressed genes that best distinguish TNF-resistant (HPV18) from TNF-sensitive cells (PHK and HPV16), in normal culture conditions or upon treatment with TNF for 60 hours.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Analysis of TNF-induced NF-κB activation in normal and HPV16 or 18-transformed keratinocytes. Nuclear protein extracts were obtained from sub confluent cultures of normal keratinocytes (PHK), HPV16-positive (SiHa) and HPV18-positive (SW756 and HeLa) cervical cancer-derived cell lines treated with 2 nM of TNF for 1 h. For each EMSA reaction, 5 mg of nuclear protein were incubated with 50 fmol of [γ-32P]ATP-labeled double-stranded oligonucleotide and a 50X excess of consensus unlabeled oligonucleotide (lanes 4, 8 and 13). Specificity of binding was further demonstrated by incubation of nuclear extracts with the described amount of labeled consensus oligonucleotide and a 50X excess of a labeled oligonucleotide carrying a single-base mutation at the NF-κB binding site (lanes 3, 7 and 14) and incubation of nuclear extract in the absence of any labeled probe (lane 15). NF-κB-DNA binding reactions were carried out as described under "Material and Methods". HeLa nuclear extracts obtained from two different experiments were included as controls of NF-κB activation levels in HPV-positive cell lines. DNA binding complexes are indicated.</p></caption></supplementary-material>
<table-wrap-foot><p><bold>*</bold>Genes are listed in alphabetical order. The cutoff p-value was set as &lt;10<sup>-10</sup>.</p></table-wrap-foot> <table-wrap-foot><p><bold>*</bold>Genes are listed in alphabetical order. The cutoff p-value was set as &lt;10<sup>-2,9</sup>.</p></table-wrap-foot> <table-wrap-foot><p><bold>*</bold>Genes are listed in alphabetical order. Underlined genes were identified as differentially expressed between TNF sensitive and TNF resistant cells in both culture conditions.</p></table-wrap-foot>
<graphic xlink:href="1755-8794-1-29-1"/> <graphic xlink:href="1755-8794-1-29-2"/> <graphic xlink:href="1755-8794-1-29-3"/> <graphic xlink:href="1755-8794-1-29-4"/> <graphic xlink:href="1755-8794-1-29-5"/> <graphic xlink:href="1755-8794-1-29-6"/>
<media xlink:href="1755-8794-1-29-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-29-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-29-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-29-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-29-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-29-S6.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Stanley"], "given-names": ["M"], "article-title": ["Immune responses to human papillomavirus"], "source": ["Vaccine"], "year": ["2006"], "volume": ["24"], "fpage": ["16"], "lpage": ["22"], "pub-id": ["10.1016/j.vaccine.2005.09.002"]}, {"surname": ["Tizard"], "given-names": ["IR"], "article-title": ["Immunity at body surfaces"], "source": ["Vet Immunol"], "year": ["2000"], "fpage": ["222"], "lpage": ["34"]}, {"surname": ["Tay", "Jenkins", "Maddox Campion", "Singer"], "given-names": ["SK", "D", "M", "A"], "article-title": ["Subpopulation of Langerhans' cells in cervical neoplasia"], "source": ["Br J Obstet Gynecol"], "year": ["1987"], "volume": ["94"], "fpage": ["10"], "lpage": ["5"]}, {"article-title": ["Gene Expression Omnibus-NCBI"], "comment": ["under acc numbers GPL1930 and GSE4524."]}, {"article-title": ["The R Project for Statistical Computing"]}, {"article-title": ["Bioconductor software for bioinformatics"]}, {"collab": ["IARC \u2013 International Agency for Research on Cancer-World Health Organization"], "article-title": ["IARC Monographs on the evaluation of carcinogenic risks to humans: human papillomaviruses"], "source": ["Lyon: Human papillomavirus"], "year": ["2007"], "volume": ["90"]}]
BMC Med Genomics. 2008 Jun 27; 1:29