Skip to main content
Download PDF

HPV-positive murine oral squamous cell carcinoma: development and characterization of a new mouse tumor model for immunological studies

Journal of Translational Medicine volume 21, Article number: 376 (2023) Cite this article

Abstract

Background

Infection with high-risk human papillomavirus (HPV) strains is one of the risk factors for the development of oral squamous cell carcinoma (OSCC). Some patients with HPV-positive OSCC have a better prognosis and respond better to various treatment modalities, including radiotherapy or immunotherapy. However, since HPV can only infect human cells, there are only a few immunocompetent mouse models available that enable immunological studies. Therefore, the aim of our study was to develop a transplantable immunocompetent mouse model of HPV-positive OSCC and characterize it in vitro and in vivo.

Methods

Two monoclonal HPV-positive OSCC mouse cell lines were established by inducing the expression of HPV-16 oncogenes E6 and E7 in the MOC1 OSCC cell line using retroviral transduction. After confirming stable expression of HPV-16 E6 and E7 with quantitative real-time PCR and immunofluorescence staining, the cell lines were further characterized in vitro using proliferation assay, wound healing assay, clonogenic assay and RNA sequencing. In addition, tumor models were characterized in vivo in C57Bl/6NCrl mice in terms of their histological properties, tumor growth kinetics, and radiosensitivity. Furthermore, immunofluorescence staining of blood vessels, hypoxic areas, proliferating cells and immune cells was performed to characterize the tumor microenvironment of all three tumor models.

Results

Characterization of the resulting MOC1-HPV cell lines and tumor models confirmed stable expression of HPV-16 oncogenes and differences in cell morphology, in vitro migration capacity, and tumor microenvironment characteristics. Although the cell lines did not differ in their intrinsic radiosensitivity, one of the HPV-positive tumor models, MOC1-HPV K1, showed a significantly longer growth delay after irradiation with a single dose of 15 Gy compared to parental MOC1 tumors. Consistent with this, MOC1-HPV K1 tumors had a lower percentage of hypoxic tumor area and a higher percentage of proliferating cells. Characteristics of the newly developed HPV-positive OSCC tumor models correlate with the transcriptomic profile of MOC1-HPV cell lines.

Conclusions

In conclusion, we developed and characterized a novel immunocompetent mouse model of HPV-positive OSCC that exhibits increased radiosensitivity and enables studies of immune-based treatment approaches in HPV-positive OSCC.

Background

Oral squamous cell carcinoma (OSCC) includes malignant neoplasms of the lip, oral cavity and oropharynx. In 2020, approximately 480,000 cases of OSCC were diagnosed worldwide, representing close to 2.5% of all diagnosed malignancies [1]. Standard approaches for the treatment of locally confined OSCC include resection, radiotherapy, and chemoradiotherapy. For recurrent or metastatic disease, systemic approaches such as chemotherapy and immunotherapy with cetuximab or immune checkpoint inhibitors are used [2]. The main risk factors for the development of OSCC are alcohol and tobacco consumption [3, 4]. Additionally, several epidemiological studies have shown that persistent infection with high-risk human papillomavirus (HPV) strains correlates with an increased risk for the development of oropharyngeal OSCC. High-risk HPV DNA was also detected with increased frequency in malignant lesions of the oral epithelium; however, the role of HPV in the development of non-oropharyngeal OSCC is not yet clear [5,6,7,8,9,10]. Moreover, clinical data indicate that patients with HPV-positive OSCC can have a better prognosis [9, 11, 12] with a better response to treatment with chemo-radiotherapy [13], radiotherapy alone [14, 15] or immunotherapy with immune checkpoint inhibitors [16].

HPVs are separated into low-risk and high-risk strains based on their oncogenic potential. Out of more than 200 identified HPV strains, 12 have been recognized as high risk [17]. Among those, the most important for OSCC development are HPV-16 and HPV-18, which contribute to 85% of all HPV-positive head and neck cancers [18]. The main factors associated with HPV carcinogenesis are viral oncogenes E6 and E7, which induce ubiquitination and proteasomal degradation of two key tumor-suppressing proteins, p53 and retinoblastoma protein (pRb), respectively [19, 20], leading to uncontrolled DNA replication, limited DNA damage response, and consequently genome instability [21,22,23]. Importantly, high-risk HPV DNA can also integrate into the host cell genome through integration breakpoints that commonly occur in the viral genes E1 and E2, which are responsible for the regulation of viral oncogene expression, resulting in elevated expression levels and stabilization of E6 and E7 mRNA [24,25,26].

Although a significant proportion of OSCCs are associated with HPV infection, there is a lack of preclinical models suitable for studies of these mechanisms and immunological studies. In recent years, great progress has been made in the development of innovative in vitro approaches for 3D modelling of HPV-positive and HPV-negative OSCC (reviewed in [27, 28]) which have several advantages over conventional adherent 2D cell cultures, such as better replication of the tumor tissue architecture and physical properties of tumor microenvironment. Culturing of human OSCC cell lines or primary cells on collagen-based 3D scaffolds enables the crosstalk between cancer cells and extracellular matrix, which mimics the tumor microenvironment [29]. Tumor tissue has also been cultured in the form of histocultures, in which primary tumor tissue is excised and cultured ex vivo. Although this method enables co-culturing of cancer cells in their original microenvironment with immune, stromal and vascular cells, the short lifespan of such cultures is a limiting factor. Engelmann et al. [30] used supportive scaffolds of healthy human-derived fibroblasts and viscose fibers to prolong the proliferation of the tumor cells in histocultures to up to 21 days. In addition, microfluidic chips have been used to co-culture OSCC cells with immune cells to study the effect of immunotherapies [31, 32]. Nevertheless, immunocompetent in vivo mouse models are still commonly used and needed in research of novel immune-based treatment approaches. Since HPV can only infect human cells, there is a lack of immunocompetent mouse models of HPV-positive OSCC. The majority of the in vivo models used in preclinical research are xenograft mouse models using cells originating from patient-derived HPV-positive tumors transplanted into immunodeficient athymic nude mice, which are not suitable for immunologic studies [33,34,35,36,37].

Therefore, the goal of our study was to develop and characterize a novel transplantable immunocompetent mouse model of HPV-positive OSCC that would resemble the clinical features of HPV-positive OSCC. Such a model would enable studies of the immunological components of currently employed treatments and evaluation of novel treatment approaches that are based on mechanisms of adaptive immunity in both HPV-positive and HPV-negative models of OSCC. The proposed model was developed through stable integration of the E6 and E7 oncogenes into the genome of a mouse OSCC cell line using retroviral transduction and further characterized in vitro and in vivo.

Methods

Establishment of HPV-positive OSCC cell line

An HPV-positive OSCC cell line was established by transduction of the murine OSCC cell line MOC1 with LXSN16E6E7 retroviral particles encoding the HPV-16 E6 and E7 genes along with a gene for neomycin resistance. After antibiotic selection (600 µg/mL G418 disulfate solution (AppliChem, Darmstadt, Germany)) for transduced cells, monoclonal cell lines were established. Monoclonal cell lines used in this study were named MOC1-HPV K1 and MOC1-HPV K3 and were cultured under standard conditions (37 °C, 5% CO2 humidified atmosphere) in MOC1 medium (Iscove's Modified Dulbecco's Medium (IMDM)/Ham’s F12 Nutrient Mixture (Gibco, Thermo Fisher Scientific, Waltham, MA, USA), in a ratio of 2:1, supplemented with 5% (v/v) fetal bovine serum (FBS, Gibco), 100 U/mL penicillin, 100 µg/mL streptomycin (100 × penicillin‒streptomycin, Sigma Aldrich, Darmstadt, Germany), 5 µg/mL insulin (Sigma Aldrich), 40 ng/mL hydrocortisone (Sigma Aldrich) and 5 ng/mL epidermal growth factor (EGF, Gibco)) with added 200 µg/mL G418 disulfate. To confirm the stable expression of HPV-16 E6 and E7, RNA was extracted using a peqGOLD Total RNA Kit (VWR, West Chester, PA, USA) following the manufacturer’s instructions, and qRT‒PCR for E6, E7, and housekeeping genes (GAPDH and β-actin) was performed. Immunofluorescence staining was used to confirm the expression of E6 and E7 at the protein level.

In vitro characterization of the newly established cell lines

Cell proliferation assay

To determine the proliferation kinetics of the cell lines, PrestoBlue™ Cell Viability Reagent (Invitrogen, Thermo Fisher Scientific) was performed on four consecutive days according to the manufacturers’ instructions.

Wound healing assay

Cell migration capacity was assessed by wound-healing assay using 24-well plates with two-well silicone inserts (Ibidi, Gräfelfing, Germany).

Radiosensitivity evaluation in vitro

To determine the radiosensitivity of the cell lines, cells were irradiated, and a clonogenic assay was performed.

Transcriptome analysis

For transcriptome analysis RNA-sequencing was performed on four biological replicates of each cell line. In addition, to compare the transcriptomic profile of murine MOC1-HPV cell lines with patient data, transcriptome analysis of publicly available data from four HPV-positive (UM-SCC-104, UM-SCC-47, and UPC-SCC-090) and one HPV-negative cell line (FaDu), as well as HPV-positive and HPV-negative patient tumor samples (NCBI Gene Expression Omnibus (GEO) series accession number GSE211322) [38] was carried out. NCBI-generated counts data of GSE211322 samples was downloaded and analyzed using free online tool iDEP.951 [39].

In vivo characterization of MOC1-HPV tumor models

Animals and tumor induction

All experiments were approved by the Ministry of Agriculture, Forestry and Food of the Republic of Slovenia (permission no. U34401-35/2020/8). Up to six mice were housed per cage in specific-pathogen-free conditions in a carousel mouse IVC rack system (Animal Care Systems Inc., Revere Parkway, USA) with a 12 h light–dark cycle, food and water ad libitum, and cage enrichment (VWR). To induce tumor growth 1 × 106 MOC1, MOC1-HPV K1, or MOC1-HPV K3 cells were injected subcutaneously into the flank of 8- to 12-week-old female C57Bl/6NCrl mice (Charles River, Lecco, Italy).

Determination of HPV-16 E6 and E7 in vivo

Expression of HPV-16 E6 and E7 was confirmed by quantitative real-time PCR (qRT‒PCR) in tumors at volumes of 50–60 mm3.

Tumor histology

Tumors were collected at 50–60 mm3 or 100 mm3, fixed in 10% neutral buffered formalin and embedded in paraffin. Paraffin sections (2-µm-thick) were stained with hematoxylin and eosin (HE), and evaluated by an experienced pathologist.

Tumor microenvironment assessment

Characteristics of the tumor microenvironment were assessed on tumors at 50–60 mm3. Three hours before the animals were sacrificed, EdU (Abcam, Cambridge, United Kingdom) and EF5 (EF5 Hypoxia Detection Kit, Cyanine 3; EMD Millipore, CA, USA) were injected intraperitoneally to mark the proliferating cells and hypoxic areas, respectively. Collected tumors were fixed in 4% PFA and dehydrated in 30% sucrose before embedding in Tissue-Tek O.C.T. compound (Sakura Finetek, VWR). Fourteen µm thick frozen tissue sections were then immunofluorescently labeled for EdU, EF5, CD31 (R&D Systems, MN, USA), CD4 (Abcam), CD8 (Abcam) and F4/80 (Thermo Fisher Scientific). Nuclei were counterstained with Hoechst 33343 pentahydrate (Thermo Fisher Scientific).

Radiosensitivity evaluation in vivo

To evaluate the radiosensitivity of tumors in vivo, mice were randomly assigned to the control or irradiated group when the tumors reached 50–60 mm3 using a Graphpad QuickCalcs random number generator. The irradiated group received a single irradiation dose of 15 Gy. Tumor growth was evaluated by determining the tumor volume three times weekly using the formula for ellipsoid (V = a × b × c × π/6, where a, b, and c are perpendicular tumor diameters). Animal weight was monitored three times a week as a sign of wellbeing. Growth delay was calculated by determining the time irradiated tumors needed to grow from 50–60 mm3 to 100 mm3 (time to 100 mm3) and subtracting the time to 100 mm3 of control tumors.

Statistical analysis

Statistical analysis and graph plotting were performed using GraphPad Prism 9 (La Jolla, CA, USA). Datasets were tested for normal distribution using the D’Agostino-Pearson normality test. Statistical significance was evaluated using ordinary one-way ANOVA with Tukey’s multiple comparisons test for normally distributed data and nonparametric one-way ANOVA (Kruskal‒Wallis test) for data without a normal distribution. A p < 0.05 was considered statistically significant. Throughout the manuscript, the following symbols indicate statistical significance: *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001. The sample size (n) for each experiment is presented in the figure legends and represents the number of biological replicates unless otherwise stated. All experiments were repeated at least two times.

A detailed description of the methods is included in the Additional file 2: Supplementary methods.

Results

Establishment of HPV-16 E6- and E7-expressing cell lines MOC1-HPV

To establish an HPV-positive OSCC cell line, we inserted the HPV-16 E6 and E7 oncogenes into the genome of MOC1 cells by retroviral transduction, as expression of the E6 and E7 oncogenes was shown to be sufficient for cell transformation [40, 41]. Using qRT‒PCR, we confirmed that insertion of HPV-16 E6 and E7 oncogenes into the cell genome was stable in both monoclonal cell lines (MOC1-HPV K1 and MOC1-HPV K3), resulting in stable expression of E6 and E7 mRNA across at least 20 cell passages (Fig. 1A, B). Furthermore, stability of the overall gene expression in MOC1-HPV cell lines is also evident from RNA sequencing results (Additional file 1: Fig. S1). Expression of E6 and E7 proteins was also confirmed by immunofluorescence staining (Fig. 1C, D). E6 (Fig. 1E, F) and E7 (Fig. 1G, H) were detected in both the nuclei and cytoplasm of the cells, with significantly higher mean fluorescence intensity in both MOC1-HPV cell lines than in the parental MOC1 cell line.

Fig. 1
figure 1

MOC1-HPV monoclonal cell lines express HPV-16 E6 and E7 at the RNA and protein levels. A, B Expression of HPV-16 E6 (A) and E7 (B) at the mRNA level remains stable with cell passaging as determined by qRT‒PCR. Housekeeping genes: GAPDH and β-actin. The data are presented as the mean ± SD. n = 2. C, D Representative images of cells immunofluorescently stained for HPV-16 E6 (C) or HPV-16 E7 (D). Scale bar: 10 µm. E–H Mean fluorescence intensity of E6 (E, F) and E7 (G, H) in the nucleus and cytoplasm. Dot plots present all data with the black line representing the mean. Up to 100 cells of each cell line were analyzed. *p < 0.05, ***p < 0.001, ****p < 0.0001, one-way ANOVA

Full size image

Morphological examination, cell growth and migration capacity in vitro

For in vitro characterization of the cell lines, we first focused on cell morphology, cell growth kinetics, and cell migration ability. Cell morphology was examined by confocal and bright-field microscopy. MOC1 and MOC1-HPV K3 cells grew in close contact, whereas MOC1-HPV K1 cells showed less pronounced cell-to-cell contacts (Fig. 2A), which was also reflected in looser MOC1-HPV K1-cell colonies in the clonogenic assay (Fig. 2B). In addition, MOC1-HPV K1 cells appeared larger than MOC1 and MOC1-HPV K3 cells (Fig. 2A) with significantly larger cell diameters when attached (Fig. 2C) but not when in suspension after trypsinization (Fig. 2D), also suggesting differences in cell adhesion properties.

Fig. 2
figure 2

In vitro characterization of MOC1 and MOC1-HPV cell lines. A Representative immunofluorescent images of MOC1 and MOC1-HPV cells after staining of the plasma membrane with WGA (red) and nuclei with Hoechst 33342 (blue). Scale bar: 40 µm. B Two representative cell colonies of MOC1, MOC1-HPV K1, and MOC1-HPV K3 cells six days after seeding. Scale bar: 200 µm. C Average cell diameter of adherent cells. Dot plots present all data with the black line representing the mean. Between 50 and 65 cells of each cell line were analyzed. ****p < 0.0001, ns: not significant, nonparametric one-way ANOVA (Kruskal‒Wallis test). D Cell diameter of trypsinized cells in suspension. In total, up to 11,000 cells of each cell line from 10 individual experiments were analyzed. The bar plot presents the mean of each experiment ± SD. **p < 0.01, ****p < 0.0001, ns not significant, ordinary one-way ANOVA. E Proliferation of MOC1 and MOC1-HPV cell lines. Presented are the mean + SD. n = 3. F Cell front velocity of MOC1 and MOC1-HPV cell lines in the wound-healing assay. Presented are the mean ± SD. ****p < 0.0001, ordinary one-way ANOVA. n = 6. G, H Representative bright-field images of cell migration in the wound-healing assay. Scale bar: 500 µm (G) or 250 µm (H)

Full size image

Although MOC1-HPV K3 cells proliferated slower in vitro, the differences in cell growth between the cell lines were not statistically significant (Fig. 2E). To assess the migration potential of the cell lines in vitro, a wound-healing assay was performed (Fig. 2F–H). The cell front velocity was significantly lower in the MOC1-HPV K3 cells (Fig. 2F). While the MOC1 and MOC1-HPV K1 cell lines closed the 500 µm wide cell-free gap approximately 12 h after removal of the inserts, the MOC1-HPV K3 cells required more than 18 h (Fig. 2G). Cells of all three cell lines showed signs of collective migration, with most cells maintaining intercellular connections during migration and wound closure, although looser cell-to-cell contacts were still observed in MOC1-HPV K1 (Fig. 2H). We have also confirmed that the observed phenotype in MOC1-HPV cells with regards to the differences in the migration kinetics is independent of cell passaging, by repeating the wound-healing assay after prolonged culturing of MOC1-HPV cells to reach passages between 17 and 20 (Additional file 1: Fig. S2).

Transcriptome analysis

Next, we used RNA sequencing to analyze the changes in gene expression in vitro after integration and expression of HPV-16 E6 and E7. The three cell lines showed different gene expression patterns with between 2000 and 4000 differentially expressed genes (Fig. 3A). With Gene Ontology (GO) analysis of the differentially expressed genes, we found that in the MOC1-HPV K1 cell line compared to the MOC1 or MOC1-HPV K3 cell line, the majority of the 20 most enriched GO terms within the “Biological Processes” compartment were related to angiogenesis, taxis, cell migration and motility, cell adhesion, and extracellular matrix organization (Fig. 3B, C), which correlates well with the observed differences in cell morphology and migration. When comparing the MOC1-HPV K3 cell line to the parental MOC1 cell line, the 20 most enriched GO biological process terms were mainly related to developmental processes and morphogenesis but also to cell taxis and migration (Fig. 3D). In terms of GO cellular component terms, the 20 most enriched in the MOC1-HPV K1 cell line both when compared to MOC1 (Fig. 4A) or MOC1-HPV K3 cell line (Fig. 4B) included extracellular matrix and cell junctions, which is in line with morphologic differences observed in vitro. When comparing MOC1-HPV K3 cells to the parental MOC1 cell line, only six significantly enriched GO cellular component terms were identified, which were related to the extracellular matrix, recycling endosome membrane, and synapse (Fig. 4C).

Fig. 3
figure 3

Transcriptome analysis of MOC1 and MOC1-HPV cell lines. A Number of differentially expressed genes in MOC1 and MOC1-HPV cell lines. Genes with an adjusted p value < 0.05 and fold change > 1 were considered differentially expressed. B–D The twenty most significantly enriched GO biological process terms with an adjusted p value < 0.05 comparing MOC1-HPV K1 to MOC1 (B), MOC1-HPV K1 to MOC1-HPV K3 (C), and MOC1-HPV K3 to MOC1 (D) are shown. The position of the dots relates to the GeneRatio, which is the fraction of all differentially expressed genes that are found in the respective GO term. The size and color of the dots represent the number of differentially expressed genes in the pathway and adjusted p value, respectively. Color-coding of the squares next to the biological terms refers to the observed cell phenotype—red: cell migration, blue: angiogenesis, green: cell‒cell contacts

Full size image
Fig. 4
figure 4

Enriched GO cellular compartment terms in MOC1 and MOC1-HPV cell lines. A–C The twenty most significantly enriched GO cellular compartment terms comparing MOC1-HPV K1 to MOC1 (A), MOC1-HPV K1 to MOC1-HPV K3 (B), and MOC1-HPV K1 to MOC1 (C) are shown. The position of the dots relates to the GeneRatio, which is the fraction of differentially expressed genes that are found in the respective GO term. The size and color of the dots represent the number of differentially expressed genes in the pathway and adjusted p value, respectively. Color-coding of the squares next to the biological terms refers to the observed cell phenotype—red: cell migration, blue: angiogenesis, green: cell‒cell contacts

Full size image

To evaluate the transcriptional resemblance of MOC1-HPV cell lines to transcriptional characteristics of human OSCC tumors, we analyzed publicly available transcriptome data of human patient-derived cell lines of HPV-positive (UM-SCC-104, UM-SCC-47, UPCI-SCC-090) and HPV-negative (FaDu) OSCC and of patient tumor samples (Additional file 1: Tables S1–S4). The most enriched GO pathways when comparing human HPV-positive and HPV-negative cell lines included pathways related to extracellular matrix, angiogenesis, and cell migration, which largely overlapped with the results of the comparison of MOC1-HPV and MOC1 cell lines. On the other hand, the enriched pathways in HPV-positive patient tumors did not overlap with those enriched in MOC1-HPV cell lines, and also to a lesser extent with enriched pathways in human HPV-positive OSCC cell lines.

MOC1 and MOC1-HPV cell lines induce tumors in C57Bl/6 mice

Next, we aimed to characterize the HPV-positive OSCC models in vivo. Subcutaneous injection of 1 × 106 MOC1 or MOC1-HPV cells into the flank of C57Bl/6 mice induced subcutaneous tumors. Palpable tumors were obtained approximately one week after tumor cell injection in approximately 88% of the mice. MOC1 tumors reached a volume of 50–60 mm3 in 26.3 ± 4.0 days (min 15, max 46 days), MOC1-HPV K1 tumors in 21.4 ± 7.2 days (min 8, max 48 days), and MOC1-HPV K3 tumors in 15.7 ± 1.8 days (min 10, max 34 days) with comparable tumor growth kinetics between the experiments. Using qRT‒PCR, we confirmed the expression of HPV-16 E6 and E7 mRNA in MOC1-HPV K1 and MOC1-HPV K3 tumors at 50–60 mm3 (Fig. 5A). Histological analysis of paraffin sections of MOC1 tumors stained with hematoxylin and eosin showed typical features of differentiated squamous cell carcinoma, such as keratinization and desmoplasia. Parts of the tumor core showed necrotic areas (Fig. 5B). MOC1-HPV K1 tumors exhibited a high degree of immune cell infiltration both at the tumor margin and in the tumor core (Fig. 5C), with few necrotic areas, few apoptotic cells, and some mitotic cells (Fig. 5D). On the other hand, MOC1-HPV K3 tumors showed larger necrotic areas and less immune cell infiltration (Fig. 5E). In addition, several cells with large nuclei and multinuclear cells were found in MOC1-HPV K3 tumors (Fig. 5F). In MOC1-HPV K1 and MOC1-HPV K3 tumors, several cells had a perinuclear halo (Fig. 5F).

Fig. 5
figure 5

MOC1-HPV tumors express HPV oncogenes E6 and E7 and show histological features of squamous cell carcinomas. A Expression of the HPV-16 E6 and E7 oncogenes in tumor tissue was determined at a tumor volume of 50 mm3 using qRT‒PCR. Housekeeping genes: GAPDH and β-actin. Presented is the mean ± SD. n = 4–5. B–F Representative images of hematoxylin and eosin-stained paraffin sections of MOC1 (B), MOC1-HPV K1 (C, D), and MOC1-HPV K3 (E, F) tumor cores. Tumors were isolated when they reached volumes of 50–100 mm3. Dotted lines denote necrotic areas, asterisks denote desmoplasia, black denote arrows keratinization, white arrows denote apoptotic cells, red arrows denote cells with large nuclei, and yellow arrows denote cells with a perinuclear halo. Images B–E were acquired using a 10 × objective, and image F was acquired with a 20 × objective. Scale bar: 100 µm. n = 4

Full size image

Tumor microenvironment of MOC1 and MOC1-HPV tumors

Next, we were interested in the tumor microenvironment of the three tumor models, particularly in the extent of hypoxia and proliferating cells. To this end, we immunofluorescently labeled proliferating cells with EdU, hypoxic areas with EF5, and blood vessels with CD31 (Fig. 6A). By quantifying the fluorescence signal, we showed that MOC1-HPV K1 tumors had significantly more proliferating cells than MOC1 and MOC1-HPV K3 tumors (Fig. 6B). Moreover, MOC1-HPV K1 tumors were significantly less hypoxic than MOC1 and MOC1-HPV K3 tumors (Fig. 6C), even though the percentage of blood vessel area was comparable in all three tumor models (Fig. 6D). However, when examining the area of each vessel, the MOC1-HPV K1 tumors had significantly smaller vessels than the MOC1 and MOC1-HPV K3 tumors (Fig. 6E).

Fig. 6
figure 6

Tumor microenvironment of MOC1 and MOC1-HPV tumors. A Representative images of immunofluorescently stained frozen sections of the tumor core. Proliferating cells are marked with EdU (Green), hypoxic regions with EF5 (Orange), blood vessels with CD31 (Magenta), and cell nuclei with Hoechst 33342 (Blue). Scale bar: 200 µm. B–E The fluorescent signal was quantified and expressed as the percentage of proliferating cells (B), percentage of hypoxic area (C), percentage of tumor blood vessel area (D), and average vessel size (E). Bar plots present the mean ± SD. The black line in the dot plot presents the mean. n = 4. ns not significant, **p < 0.01, ***p < 0.001, ****p < 0.0001, ordinary one-way ANOVA (B–D) and nonparametric one-way ANOVA (Kruskal‒Wallis test) (E)

Full size image

Since there are reports that HPV-positive OSCC tumors have higher infiltration of immune cells than HPV-negative tumors [42, 43], we were also interested if E6 and E7 expression induced changes in immune cell infiltration in our tumor models. Therefore, we immunofluorescently stained for CD8 T lymphocytes, CD4 T lymphocytes, macrophages expressing the surface marker F4/80 and blood vessels expressing CD31 (Fig. 7A, B). A comparable number of CD4 T lymphocytes per mm2 of tumor area was detected in the tumor core (Fig. 7C) as well as in the tumor edge (Fig. 7D) in all three tumor models. A similar trend was observed for CD8 T lymphocytes (Fig. 7E, F). Although the average CD4/CD8 ratio was slightly higher in MOC1-HPV K1 and MOC1-HPV K3 tumors, the differences between the three tumor models were not statistically significant (Fig. 7G, H). Similarly, macrophages were present at comparable frequencies in all three tumor models in both the tumor core (Fig. 7I) and tumor edge (Fig. 7J).

Fig. 7
figure 7

Immune cell infiltration of MOC1 and MOC1-HPV tumors. A, B Representative images of immunofluorescently stained frozen sections of the tumor core (A) or tumor edge (B). T lymphocytes were marked with CD4 or CD8 (orange), macrophages with F4/80 (green), blood vessels with CD31 (magenta), and cell nuclei with Hoechst 33,342 (blue). Scale bar: 100 µm. C, D Numbers of CD4 T lymphocytes per mm2 of tumor core (C) and tumor edge (D). E, F Numbers of CD8 T lymphocytes per mm2 of tumor core (E) and edge (F). G, H CD4/CD8 ratio in the tumor core (G) and tumor edge (H). I, J Percentage of F4/80-positive area in the tumor core (I) and tumor edge (J). Bar plots present the mean ± SD. n = 4. ns not significant, ordinary one-way ANOVA

Full size image

Response of MOC1 and MOC1-HPV cell lines and tumors to radiotherapy

Clinical reports suggest that patients with HPV-positive OSCC respond better to radiotherapy [14, 15]; therefore, we were interested in the effect of E6 and E7 expression on the radiosensitivity of MOC1-HPV cells. First, we determined the radiosensitivity of the three cell lines in vitro using a clonogenic assay. The reproductive integrity after irradiation with a single irradiation dose between 0 and 12 Gy did not differ significantly among all tested cell lines, resulting in comparable IC50 values calculated with the linear quadratic model—6.0 ± 1.3 Gy for MOC1, 5.8 ± 2.2 Gy for MOC1-HPV K1, and 4.9 ± 0.5 Gy for MOC1-HPV K3, as well as IC90 values of 9.9 ± 1.5 Gy for MOC1, 9.5 ± 2.5 Gy for MOC1-HPV K1, and 9.1 ± 0.7 Gy for MOC1-HPV K3 (Fig. 8A). Next, we examined the sensitivity of the subcutaneous tumor models to irradiation with a single dose of 15 Gy. By measuring tumor volumes after irradiation, we showed that irradiation with a single dose of 15 Gy slowed the growth of all three tumor models compared to the unirradiated controls (Fig. 8B), but the growth delay was significantly longer in MOC1-HPV K1 tumors (15.5 days compared to 8.5 days in MOC1 and 9 days in MOC1-HPV K3) (Fig. 8C).

Fig. 8
figure 8

Radiosensitivity of MOC1, MOC1-HPV K1, and MOC1-HPV K3. A Radiosensitivity of MOC1 and MOC1-HPV cell lines as assessed with a clonogenic assay. Presented are the mean + SD. N = 3. B Tumor growth after irradiation with a single dose of 15 Gy. Presented are the mean + SEM. n = 9–14. C Tumor growth delay after irradiation with a single dose of 15 Gy. Presented are the mean ± SEM. n = 9–14. ns not significant, *p < 0.05, ordinary one-way ANOVA. Ctrl control, IR irradiated

Full size image

Discussion

In this study, we developed a syngeneic transplantable mouse tumor model of HPV-positive OSCC suitable for immunological studies in immunocompetent C57Bl/6 mice. The MOC1-HPV K1 tumor model exhibits increased radiosensitivity in vivo, reduced hypoxia, and increased proliferation compared to its HPV-negative MOC1 counterpart. Tumors present with CD8 and CD4 T lymphocytes in the tumor core and edge, which is crucial for the efficacy of immunotherapies. RNA sequencing of the MOC1 and MOC1-HPV cell lines demonstrated transcriptomic differences in angiogenesis, cell migration and cell adhesion-related pathways, which are reflected in functional in vitro assays and in in vivo tumor model characteristics.

Mouse models are crucial tools for preclinical research. Early studies of HPV-induced carcinogenesis used immunocompetent FVB/N or C57Bl/6 mouse strains genetically engineered to express HPV-16 E6 and E7 oncogenes under the control of the human keratin 14 (K14) promoter to direct transgene expression into stratified squamous epithelium [44,45,46,47,48]. Additionally, mutant oncogenes such as the mutated c-H-Ras gene or carcinogen treatment were included to increase the primary tumor formation rate [48]. To direct the expression of E6 and E7 to the head and neck region, inducible transgenic mouse models were developed using the Cre recombinase system, in which the expression of E6 and E7 can be induced by local intralingual administration of tamoxifen or doxycycline [49,50,51], resulting in OSCC that recapitulates the histological and molecular characteristics of human HPV-positive OSCC [49]. However, genetically engineered mouse tumor models are time-consuming and have variable tumor formation rates and locations. Therefore, transplanted mouse tumor models are easier to handle in the context of 3Rs due to the reduced variability of the experiments. Most available transplanted mouse models of HPV-positive OSCC are patient-derived xenografts (PDXs) that are transplanted subcutaneously into immunocompromised mice. There are several reports on establishing PDXs of HPV-positive and HPV-negative head and neck SCC, which retain the tumor characteristics of the donor and in some cases recapitulate patient clinical responses [33,34,35,36, 52]. The main obstacles in developing HPV-positive PDX models are low engraftment rates and often the formation of Epstein‒Barr virus-positive lymphomas derived from cotransfer of stroma from the tonsil or base of the tongue [34, 52], both of which contribute to higher costs. Moreover, because immunocompromised mouse strains are used for PDXs, their use in preclinical immunological studies is limited.

Because of the high tumor induction success rate and time effectiveness of transplantable immunocompetent mouse models, we aimed to develop an HPV-positive murine cell line of OSCC that would resemble the clinical features of HPV-positive OSCC and would enable both in vitro studies and in vivo studies in immunocompetent mice. Since HPV can only infect human cells, there are very few murine cell lines of HPV-positive OSCC. Hoover et al. and Williams et al. reported the establishment of an HPV-positive cell line, mEER, by transducing mouse pharyngeal epithelial cells isolated from C57Bl/6 mice with a retroviral vector encoding the H-ras and HPV-16 E6 and E7 oncogenes, capable of forming transplantable tumors [53, 54]. Similarly, MOC2-E6E7 cells were established by transduction of highly metastatic murine OSCC MOC2 cells with a retrovirus encoding HPV-16 E6 and E7 [55]. Furthermore, Paolini et al. established an HPV-positive cell line, AT-84-E7, by stably transfecting the AT-84 cell line with a plasmid containing the HPV-16 E7 oncogene. AT-84-E7 cells can be injected into the floor of the mouth of immunocompetent C3H mice to obtain an orthotopic mouse model of HPV-positive OSCC [56]. The characteristics of all available HPV-positive murine cell lines are described in Table 1.

Table 1 Comparison of the newly established MOC1-HPV tumor models with the existing models of HPV-positive OSCC
Full size table

Expression of the HPV oncogenes E6 and E7 can lead to several phenotypic changes in HPV-positive cells. Studies show that HPV oncogenes increase migration capacity in some cell types, including cervical cancer cells and trophoblastic cells [57,58,59], while the results of studies on HPV-positive head and neck SCC cells are inconclusive. Nagel et al. tested the migratory capacity of several OSCC cell lines using a wound-healing assay and showed that the migratory capacity of HPV-positive cell lines varies and is not significantly different from that of HPV-negative SCC cell lines [60]. Similarly, Kahue et al. showed that the migratory capacity of HPV-positive OSCC cell lines does not significantly differ from the migratory capacity of HPV-negative OSCC cell lines; however, when the expression of E6 and E7 is abrogated in HPV-positive cell lines, their migratory capacity decreases, indicating a potential role of E6 and E7 in migration properties [61]. In our study, differences in cell migration were detected both in in vitro functional assays and at the transcriptome level. However, the expression of HPV-16 E6 and E7 did not change the migration velocity in MOC1-HPV K1 cells and even decreased the migration velocity of MOC1-HPV K3 cells in comparison to the parental MOC1 cell line. Therefore, it seems that the effect on migration properties is multifactorial and cannot be attributed solely to HPV oncogenes. Furthermore, E6 and E7 impair cell adhesion in some cell types proposedly through E6-mediated downmodulation of focal adhesion [59, 62], which may contribute to the observed differences in migration and adhesion properties of the MOC1-HPV K1-cell line compared to the parental MOC1 cell line.

Studies on patient tumor samples observed differences in the tumor microenvironment of HPV-positive and HPV-negative tumors, namely, greater T lymphocyte infiltration with a lower CD4/CD8 ratio in HPV-positive OSCC tumors, which correlated with better overall survival [42, 43, 63,64,65]. Greater infiltration of CD8 T lymphocytes was also observed in flow cytometric evaluation of the orthotopic MOC2-E6E7 mouse tumor model [55] but not in the immunofluorescence staining of our subcutaneous MOC1-HPV tumor models, even though MOC1-HPV K1 tumors presented with higher lymphocyte and neutrophilic granulocyte infiltration on HE-stained paraffin-embedded tumor tissue sections. Therefore, the focus of further studies will be the precise characterization of the immune profile of MOC1-HPV tumors with flow cytometry, where multiple immune cell populations can be detected simultaneously in the whole tumor. Moreover, additional parameters could be assessed to better evaluate the immune phenotype of MOC1-HPV tumors. A study by Mytilineos et al. [66] investigated and compared the levels of circulating cytokines in serum of patients with HPV-negative and HPV-positive OSCC, which would be interested to investigate and compare in MOC1-HPV mouse model as well. Nevertheless, we have shown that MOC1-HPV tumors have an immune-inflamed phenotype [67, 68] with T lymphocyte infiltration both in the tumor edge and in the tumor core, making them suitable mouse models for immunotherapy studies.

Clinical data show increased radiosensitivity of HPV-positive OSCC, leading to improved patient survival [14, 15]. Subcutaneous mouse tumors derived from E6E7/Ras pharyngeal cells showed a better response to irradiation with a single dose between 8 and 32 Gy than tumors derived from HPV-negative shPTPN13/Ras pharyngeal cells [69]. However, the results of in vitro studies are inconclusive. Some report on higher radiosensitivity of HPV-positive cell lines in clonogenic or cell viability assays [70,71,72], while others observed no significant differences [60] or even lower radiosensitivity compared with HPV-negative cell lines [69]. Similarly, in this study, all three cell lines showed comparable radiosensitivity in vitro, but one of the tumor models, MOC1-HPV K1, presented with significantly higher sensitivity to irradiation with a single dose of 15 Gy in vivo, indicating that the intrinsic radiosensitivity of the cells is not the main factor determining the response of these tumor models to radiotherapy. The irradiation dose of 15 Gy was selected as the dose that would not cure the tumors, which would enable us to perform the growth delay assay and compare the radiosensitivity of the tumor models without generating dose response curves for each tumor model. The latter would require a significantly higher number of animals [73] and would not provide significantly more information on the radiosensitivity of the models.

One of the factors compromising the response to irradiation is hypoxia, which reduces the amount of radiation-induced free radicals and DNA damage [74]. In line with this, we observed significantly less tumor hypoxia in the more radiosensitive MOC1-HPV K1 tumor model. Based on literature data, more than half of patients with head and neck cancer present with tumor hypoxia [75, 76], with no statistically significant differences in the level or distribution of hypoxia among HPV-negative and HPV-positive tumors [76,77,78]. Interestingly, differences in hypoxia correlated with in vitro transcriptional differences in the pathways related to blood vessel morphogenesis and angiogenesis. Moreover, despite all three tumor models showing a comparable percentage of tumor blood vessel area, vessels in MOC1-HPV K1 tumors were significantly smaller, which could lead to better perfusion of these tumors [79]. Additional experiments to determine whether the functionality of the vessels in the three tumor models differ would therefore be of interest.

Differences in the characteristics of MOC1-HPV K1 and MOC1-HPV K3 tumor models could also be due to the different E6 and E7 integration sites in the established monoclonal cell lines. During persistent HPV infection, parts of the viral genome, including the E6 and E7 oncogenes, frequently integrate into the host genome. Based on currently available studies, HPV integration into the genome is random; however, more integration sites seem to be located in or in close proximity to transcriptionally active and genic regions [26, 80, 81]. In addition, a recent study by Mima et al. [38] suggests that epigenetic regulation of genes in the proximity of HPV integration sites may play an important role in the biology of HPV-related OSCC. Parfenov et al. reported on integration sites resulting in disruption of the DNA repair gene RAD51 homolog 2 (RAD51B) or tumor suppressor gene ETS2, elevated expression of the tumor oncogene NR4A2 or elevated expression of a PD-L1 splicing form, which could all contribute to HPV-driven carcinogenesis [26]. Therefore, determination of E6 and E7 integration sites in the monoclonal cell lines MOC1-HPV K1 and MOC1-HPV K3 through whole-genome sequencing would offer insight into the genomic changes induced by E6 and E7 retroviral integration and whether they are similar to those after HPV infection in patients.

This study and the developed HPV-positive mouse OSCC models have some limitations. First, the anatomical origin of the parental MOC1 cell line is not the oropharynx, which is clinically associated with better prognosis and increased radiosensitivity in patients [11, 13,14,15]. However, mice do not have the tonsillar crypt epithelium [82, 83] where most human HPV-positive oropharyngeal tumors arise [84]. Moreover, tumors in the study were not grown orthotopically, but rather a flank model was used due to the ease of tumor induction, measurement, handling, and therapy application. Furthermore, the MOC1 cell line originates from carcinogen-induced tumors and therefore presents with different driver mutations than HPV-positive OSCC tumors in patients. For example, MOC1 cells present with mutated p53, which is an important target of HPV oncogenes in HPV-induced carcinogenesis [85, 86]. Nevertheless, the developed MOC1-HPV K1 tumor model still presents with increased radiosensitivity similar to that observed in patients. To our knowledge, all available murine OSCC cell lines are derived from carcinogen-induced tumors, therefore we have chosen MOC1 for parental cell line since it has been extensively used in various studies of OSCC biology (reviewed in [87]). Another limitation of the study is that the transcriptomic analysis was performed only on the monoclonal cell lines in vitro. Transcriptomic differences at the in vitro level correspond well with the characteristics of the tumors induced in C57Bl/6 mice as well as with transcriptomic characteristics of human patient-derived OSCC cell lines. On the other hand, the most enriched GO Biological process pathways in MOC1-HPV cell lines as well as in patient-derived OSCC cell lines mainly do not overlap with those enriched in patient tumors. Therefore, transcriptomic profile characterization of MOC1-HPV tumors would be of interest for further studies in order to better determine the transcriptomic resemblance of MOC1-HPV tumors to HPV-positive OSCC tumors in patients. Moreover, to further increase the resemblance of MOC1-HPV models to patient tumors, different co-culture setups with naïve murine immune cells prior to tumor implantation could be investigated, similarly to the Winn Assay, where tumor and immune cells are mixed prior to tumor implantation [88,89,90]. This way, immune hot or immune cold MOC1-HPV tumors could potentially be established. Another promising approach for developing the MOC1-HPV model further would also be the co-culture of MOC1-HPV tumor cells with stromal or immune cells in 2D or 3D in vitro cell model. This could increase the complexity of the in vitro model and conform with the 3R animal welfare guidelines.

Conclusions

Altogether, in this study, we established and systematically characterized two murine OSCC cell lines that stably express the HPV-16 oncogenes E6 and E7 and can be used for the induction of subcutaneous tumors in immunocompetent C57Bl/6 mice. One of the resulting mouse models of HPV-positive OSCC, MOC1-HPV K1, shows increased radiosensitivity similar to patient HPV-positive OSCC. To our knowledge, our developed syngeneic transplantable mouse tumor models are the best characterized tumor models described to date that allow us to study responses to immune-mediated treatment approaches in HPV-negative and HPV-positive OSCC, which is currently limited due to very few syngeneic HPV-positive murine OSCC cell lines.

Availability of data and materials

The data that support the findings of this study are available from the corresponding authors upon reasonable request. RNA-sequencing data obtained and discussed in this publication have been deposited to the NCBI Gene Expression Omnibus (GEO) and are accessible through GEO Series accession number GSE224448 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE224448).

Abbreviations

OSCC:

Oral squamous cell carcinoma

HPV:

Human papillomavirus

pRb:

Retinoblastoma protein

mRNA:

Messenger ribonucleic acid

qRT-PCR:

Quantitative real-time polymerase chain reaction

HE:

Hematoxylin and eosin

GO:

Gene ontology

PDX:

Patient-derived xenograft

SCC:

Squamous cell carcinoma

IC50:

Half maximal inhibitory concentration

IMDM:

Iscove’s Modified Dulbecco’s Medium

FBS:

Fetal bovine serum

EGF:

Epidermal growth factor

padj:

Adjusted p value

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.

    Article  PubMed  Google Scholar 

  2. Johnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR. Head and neck squamous cell carcinoma. Nat Rev Dis Primers. 2020;6:92.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Franceschi S, Talamini R, Barra S, Barón AE, Negri E, Bidoli E, et al. Smoking and drinking in relation to cancers of the oral cavity, pharynx, larynx, and esophagus in northern Italy. Cancer Res. 1990;50:6502–7.

    CAS  PubMed  Google Scholar 

  4. Morse DE, Psoter WJ, Cleveland D, Cohen D, Mohit-Tabatabai M, Kosis DL, et al. Smoking and drinking in relation to oral cancer and oral epithelial dysplasia. Cancer Causes Control. 2007;18:919–29.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Gillison ML, Koch WM, Capone RB, Spafford M, Westra WH, Wu L, et al. Evidence for a causal association between human papillomavirus and a subset of head and neck cancers. J Natl Cancer Inst. 2000;92:709–20.

    Article  CAS  PubMed  Google Scholar 

  6. Miller CS, Johnstone BM. Human papillomavirus as a risk factor for oral squamous cell carcinoma: a meta-analysis, 1982–1997. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2001;91:622–35.

    Article  CAS  PubMed  Google Scholar 

  7. Herrero R, Castellsagué X, Pawlita M, Lissowska J, Kee F, Balaram P, et al. Human papillomavirus and oral cancer: the International Agency for Research on Cancer multicenter study. J Natl Cancer Inst. 2003;95:1772–83.

    Article  PubMed  Google Scholar 

  8. Dalla Torre D, Burtscher D, Soelder E, Offermanns V, Rasse M, Puelacher W. Human papillomavirus prevalence in a Mid-European oral squamous cell cancer population: a cohort study. Oral Dis. 2018;24:948–56.

    Article  PubMed  Google Scholar 

  9. Kaminagakura E, Villa LL, Andreoli MA, Sobrinho JS, Vartanian JG, Soares FA, et al. High-risk human papillomavirus in oral squamous cell carcinoma of young patients. Int J Cancer. 2012;130:1726–32.

    Article  CAS  PubMed  Google Scholar 

  10. Smith EM, Rubenstein LM, Haugen TH, Pawlita M, Turek LP. Complex etiology underlies risk and survival in head and neck cancer human papillomavirus, tobacco, and alcohol: a case for multifactor disease. J Oncol. 2012. https://doi.org/10.1155/2012/571862.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Ritchie JM, Smith EM, Summersgill KF, Hoffman HT, Wang D, Klussmann JP, et al. Human papillomavirus infection as a prognostic factor in carcinomas of the oral cavity and oropharynx. Int J Cancer. 2003;104:336–44.

    Article  CAS  PubMed  Google Scholar 

  12. Sugiyama M, Bhawal UK, Kawamura M, Ishioka Y, Shigeishi H, Higashikawa K, et al. Human papillomavirus-16 in oral squamous cell carcinoma: clinical correlates and 5-year survival. Br J Oral Maxillofac Surg. 2007;45:116–22.

    Article  PubMed  Google Scholar 

  13. Rischin D, Young RJ, Fisher R, Fox SB, Le QT, Peters LJ, et al. Prognostic significance of p16INK4A and human papillomavirus in patients with oropharyngeal cancer treated on TROG 0202 phase III trial. J Clin Oncol. 2010;28:4142–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Fakhry C, Westra WH, Li S, Cmelak A, Ridge JA, Pinto H, et al. Improved survival of patients with human papillomavirus-positive head and neck squamous cell carcinoma in a prospective clinical trial. J Natl Cancer Inst. 2008;100:261–9.

    Article  CAS  PubMed  Google Scholar 

  15. Lindel K, Beer KT, Laissue J, Greiner RH, Aebersold DM. Human papillomavirus positive squamous cell carcinoma of the oropharynx: a radiosensitive subgroup of head and neck carcinoma. Cancer. 2001;92:805–13.

    Article  CAS  PubMed  Google Scholar 

  16. Galvis MM, Borges GA, de Oliveira TB, de Toledo IP, Castilho RM, Guerra ENS, et al. Immunotherapy improves efficacy and safety of patients with HPV positive and negative head and neck cancer: a systematic review and meta-analysis. Crit Rev Oncol Hematol. 2020;150:102966.

    Article  PubMed  Google Scholar 

  17. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, el Ghissassi F, et al. A review of human carcinogens—Part B: biological agents. Lancet Oncol. 2009;10:321–2.

    Article  PubMed  Google Scholar 

  18. Castellsagué X, Alemany L, Quer M, Halec G, Quirós B, Tous S, et al. HPV involvement in head and neck cancers: comprehensive assessment of biomarkers in 3680 patients. J Natl Cancer Inst. 2016;108:6.

    Article  Google Scholar 

  19. Scheffner M, Werness BA, Huibregtse JM, Levine AJ, Howley PM. The E6 oncoprotein encoded by human papillomavirus types 16 and 18 promotes the degradation of p53. Cell. 1990;63:1129–36.

    Article  CAS  PubMed  Google Scholar 

  20. Boyer SN, Wazer DE, Band V. Protein of human papilloma virus-16 induces degradation of retinoblastoma protein through the ubiquitin-proteasome pathway. Cancer. 1996;56:4620–4.

    CAS  Google Scholar 

  21. Cheng S, Schmidt-Grimminger DC, Murant T, Broker TR, Chow LT. Differentiation-dependent up-regulation of the human papillomavirus E7 gene reactivates cellular DNA replication in suprabasal differentiated keratinocytes. Genes Dev. 1995;9:2335–49.

    Article  CAS  PubMed  Google Scholar 

  22. Zerfass K, Schulze A, Spitkovsky D, Friedman V, Henglein B, Jansen-Dürr P. Sequential activation of cyclin E and cyclin A gene expression by human papillomavirus type 16 E7 through sequences necessary for transformation. J Virol. 1995;69:6389–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Moody CA, Laimins LA. Human papillomavirus oncoproteins: pathways to transformation. Nat Rev Cancer. 2010;10:550–60.

    Article  CAS  PubMed  Google Scholar 

  24. Jeon S, Lambert PF. Integration of human papillomavirus type 16 DNA into the human genome leads to increased stability of E6 and E7 mRNAs: implications for cervical carcinogenesis. Proc Natl Acad Sci USA. 1995;92:1654–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Rautava J, Syrjänen S. Biology of human papillomavirus infections in head and neck carcinogenesis. Head Neck Pathol. 2012;6(S3):15.

    Google Scholar 

  26. Prfenov M, Pedamallu CS, Gehlenborg N, Freeman SS, Danilova L, Bristow CA, et al. Characterization of HPV and host genome interactions in primary head and neck cancers. Proc Natl Acad Sci USA. 2014;111:15544–9.

    Article  Google Scholar 

  27. Chaves P, Garrido M, Oliver J, Pérez-Ruiz E, Barragan I, Rueda-Domínguez A. Preclinical models in head and neck squamous cell carcinoma. Br J Cancer. 2023;128:1819–27.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Demers I, Donkers J, Kremer B, Speel EJ. Ex vivo culture models to indicate therapy response in head and neck squamous cell carcinoma. Cells. 2020;9:2527.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Miserocchi G, Cocchi C, De Vita A, Liverani C, Spadazzi C, Calpona S, et al. Three-dimensional collagen-based scaffold model to study the microenvironment and drug-resistance mechanisms of oropharyngeal squamous cell carcinomas. Cancer Biol Med. 2021;18:502–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Engelmann L, Thierauf J, Laureano NK, Stark HJ, Prigge ES, Horn D, et al. Organotypic co-cultures as a novel 3D model for head and neck squamous cell carcinoma. Cancers. 2020;12:2330.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Al-Samadi A, Poor B, Tuomainen K, Liu V, Hyytiäinen A, Suleymanova I, et al. In vitro humanized 3D microfluidic chip for testing personalized immunotherapeutics for head and neck cancer patients. Exp Cell Res. 2019;383:111508.

    Article  CAS  PubMed  Google Scholar 

  32. Sieviläinen M, Saavalainen J, Adnan-Awad S, Salo T, Al-Samadi A. IDO1 Inhibition reduces immune cell exclusion through inducing cell migration while PD-1 blockage increases IL-6 and -8 secretion from T cells in head and neck cancer. Front Immunol. 2022. https://doi.org/10.3389/fimmu.2022.812822.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Kimple RJ, Harari PM, Torres AD, Yang RZ, Soriano BJ, Yu M, et al. Development and characterization of HPV-positive and HPV-negative head and neck squamous cell carcinoma tumorgrafts. Clin Cancer Res. 2013;19:855–64.

    Article  CAS  PubMed  Google Scholar 

  34. Kang HN, Kim JH, Park AY, Choi JW, Lim SM, Kim J, et al. Establishment and characterization of patient-derived xenografts as paraclinical models for head and neck cancer. BMC Cancer. 2020;20:316.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Facompre ND, Rajagopalan P, Sahu V, Pearson AT, Montone KT, James CD, et al. Identifying predictors of HPV-related head and neck squamous cell carcinoma progression and survival through patient-derived models. Int J Cancer. 2020;147:3236–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Keysar SB, Astling DP, Anderson RT, Vogler BW, Bowles DW, Morton JJ, et al. A patient tumor transplant model of squamous cell cancer identifies PI3K inhibitors as candidate therapeutics in defined molecular bins. Mol Oncol. 2013;7:776–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lei ZG, Ren XH, Wang SS, Liang XH, Tang YL. Immunocompromised and immunocompetent mouse models for head and neck squamous cell carcinoma. Onco Targets Ther. 2016;9:545–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Mima M, Okabe A, Hoshii T, Nakagawa T, Kurokawa T, Kondo S, et al. Tumorigenic activation around HPV integrated sites in head and neck squamous cell carcinoma. Int J Cancer. 2023;152:1847–62.

    Article  CAS  PubMed  Google Scholar 

  39. Ge SX, Son EW, Yao R. iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data. BMC Bioinform. 2018;19:534.

    Article  CAS  Google Scholar 

  40. Barbosa MS, Schlegel R. The E6 and E7 genes of HPV-18 are sufficient for inducing two-stage in vitro transformation of human keratinocytes. Oncogene. 1989;4:1529–32.

    CAS  PubMed  Google Scholar 

  41. Hawley-Nelson P, Vousden KH, Hubbert NL, Lowy DR, Schiller JT. HPV16 E6 and E7 proteins cooperate to immortalize human foreskin keratinocytes. EMBO J. 1989;8:3905–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Nordfors C, Grün N, Tertipis N, Ährlund-Richter A, Haeggblom L, Sivars L, et al. CD8+ and CD4+ tumour infiltrating lymphocytes in relation to human papillomavirus status and clinical outcome in tonsillar and base of tongue squamous cell carcinoma. Eur J Cancer. 2013;49:2522–30.

    Article  CAS  PubMed  Google Scholar 

  43. Partlová S, Bouček J, Kloudová K, Lukešová E, Zábrodský M, Grega M, et al. Distinct patterns of intratumoral immune cell infiltrates in patients with HPV-associated compared to non-virally induced head and neck squamous cell carcinoma. Oncoimmunology. 2015;4:1–12.

    Article  Google Scholar 

  44. Song S, Pitot HC, Lambert PF. The human papillomavirus type 16 E6 gene alone is sufficient to induce carcinomas in transgenic animals. J Virol. 1999;73:5887–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Herber R, Liem A, Pitot H, Lambert PF. Squamous epithelial hyperplasia and carcinoma in mice transgenic for the human papillomavirus type 16 E7 oncogene. J Virol. 1996;70:1873–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Melero I, Singhal MC, Mcgowan P, Haugen HS, Blake J, Hellstro¨m KE, et al. Immunological ignorance of an E7-encoded cytolytic T-lymphocyte epitope in transgenic mice expressing the E7 and E6 oncogenes of human papillomavirus type 16. J Virol. 1997;71:3998–4004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Jabbar S, Strati K, Shin MK, Pitot HC, Lambert PF. Human papillomavirus type 16 E6 and E7 oncoproteins act synergistically to cause head and neck cancer in mice. Virology. 2010;407:60–7.

    Article  CAS  PubMed  Google Scholar 

  48. Schreiber K, Cannon RE, Karrison T, Beck-Engeser G, Huo D, Tennant RW, et al. Strong synergy between mutant ras and HPV16 E6/E7 in the development of primary tumors. Oncogene. 2004;23:3972–9.

    Article  CAS  PubMed  Google Scholar 

  49. Carper MB, Troutman S, Wagner BL, Byrd KM, Selitsky SR, Parag-Sharma K, et al. An immunocompetent mouse model of HPV16(+) head and neck squamous cell carcinoma. Cell Rep. 2019;29:1660–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Tan MT, Wu JG, Callejas-Valera JL, Schwarz RA, Gillenwater AM, Richards-Kortum RR, et al. A PIK3CA transgenic mouse model with chemical carcinogen exposure mimics human oral tongue tumorigenesis. Int J Exp Pathol. 2020;101:45–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhong R, Pytynia M, Pelizzari C, Spiotto M. Bioluminescent imaging of HPV-positive oral tumor growth and its response to image-guided radiotherapy. Cancer Res. 2014;74:2073–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Facompre ND, Sahu V, Montone KT, Harmeyer KM, Nakagawa H, Rustgi AK, et al. Barriers to generating PDX models of HPV-related head and neck cancer. Laryngoscope. 2017;127:2777–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Williams R, Wook Lee D, Elzey BD, Anderson ME, Hostager BS, Lee JH. Preclinical models of HPV+ and HPV− HNSCC in mice: an immune clearance of HPV+ HNSCC. Head Neck. 2009;31:911–8.

    Article  PubMed  Google Scholar 

  54. Hoover AC, Spanos WC, Harris GF, Anderson ME, Klingelhutz AJ, Lee JH. The role of human papillomavirus 16 E6 in anchorage-independent and invasive growth of mouse tonsil epithelium. Arch Otolaryngol Head Neck Surg. 2007;133:495–502.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Dharmaraj N, Piotrowski SL, Huang C, Newton JM, Golfman LS, Hanoteau A, et al. Anti-tumor immunity induced by ectopic expression of viral antigens is transient and limited by immune escape. Oncoimmunology. 2019;8:1–13.

    Article  Google Scholar 

  56. Paolini F, Massa S, Manni I, Franconi R, Venuti A. Immunotherapy in new pre-clinical models of HPV-associated oral cancers. Hum Vaccin Immunother. 2013;9:534–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Wang Q, Song R, Zhao C, Liu H, Yang Y, Gu S, et al. HPV16 E6 promotes cervical cancer cell migration and invasion by downregulation of NHERF1. Int J Cancer. 2019;144:1619–32.

    Article  CAS  PubMed  Google Scholar 

  58. Hu D, Zhou J, Wang F, Shi H, Li Y, Li B. HPV-16 E6/E7 promotes cell migration and invasion in cervical cancer via regulating cadherin switch in vitro and in vivo. Arch Gynecol Obstet. 2015;292:1345–54.

    Article  CAS  PubMed  Google Scholar 

  59. Boulenouar S, Weyn C, van Noppen M, Ali MM, Favre M, Delvenne PO, et al. Effects of HPV-16 E5, E6 and E7 proteins on survival, adhesion, migration and invasion of trophoblastic cells. Carcinogenesis. 2010;31:473–80.

    Article  CAS  PubMed  Google Scholar 

  60. Nagel R, Martens-De Kemp SR, Buijze M, Jacobs G, Braakhuis BJM, Brakenhoff RH. Treatment response of HPV-positive and HPV-negative head and neck squamous cell carcinoma cell lines. Oral Oncol. 2013;49:560–6.

    Article  PubMed  Google Scholar 

  61. Kahue CN, Jerrell RJ, Parekh A. Expression of human papillomavirus oncoproteins E6 and E7 inhibits invadopodia activity but promotes cell migration in HPV-positive head and neck squamous cell carcinoma cells. Cancer Rep. 2018;1:3.

    Google Scholar 

  62. Khalifa Y, Teissier S, Tan MKM, Phan QT, Daynac M, Wong WQ, et al. The human papillomavirus E6 oncogene represses a cell adhesion pathway and disrupts focal adhesion through degradation of TAp63β upon transformation. PLoS Pathog. 2011;7:9.

    Google Scholar 

  63. Krupar R, Robold K, Gaag D, Spanier G, Kreutz M, Renner K, et al. Immunologic and metabolic characteristics of HPV-negative and HPV-positive head and neck squamous cell carcinomas are strikingly different. Virchows Arch. 2014;465:299–312.

    Article  CAS  PubMed  Google Scholar 

  64. Oguejiofor K, Hall J, Slater C, Betts G, Hall G, Slevin N, et al. Stromal infiltration of CD8 T cells is associated with improved clinical outcome in HPV-positive oropharyngeal squamous carcinoma. Br J Cancer. 2015;113:886–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Wansom D, Light E, Worden F, Prince M, Urba S, Chepeha DB, et al. Correlation of cellular immunity with human papillomavirus 16 status and outcome in patients with advanced oropharyngeal cancer. Arch Otolaryngol Head Neck Surg. 2010;136:1267–73.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Mytilineos D, Ezić J, von Witzleben A, Mytilineos J, Lotfi R, Fürst D, et al. Peripheral cytokine levels differ by HPV status and change treatment-dependently in patients with head and neck squamous cell carcinoma. Int J Mol Sci. 2020;21:1–16.

    Article  Google Scholar 

  67. Chen L, Flies DB. Molecular mechanisms of T cell co-stimulation and co-inhibition. Nat Rev Immunol. 2013;13:227–42.

    Article  PubMed  PubMed Central  Google Scholar 

  68. de Ochoa Olza M, Navarro Rodrigo B, Zimmermann S, Coukos G. Turning up the heat on non-immunoreactive tumours: opportunities for clinical development. Lancet Oncol. 2020;21:9.

    Google Scholar 

  69. Spanos WC, Nowicki P, Lee DW, Hoover A, Hostager B, Gupta A, et al. Immune response during therapy with cisplatin or radiation for human papillomavirus-related head and neck cancer. Arch Otolaryngol Head Neck Surg. 2009;135:1137–46.

    Article  PubMed  Google Scholar 

  70. Arenz A, Ziemann F, Mayer C, Wittig A, Dreffke K, Preising S, et al. Increased radiosensitivity of HPV-positive head and neck cancer cell lines due to cell cycle dysregulation and induction of apoptosis. Strahlenther Onkol. 2014;190:839–46.

    Article  PubMed  Google Scholar 

  71. Reid P, Staudacher AH, Marcu LG, Olver I, Moghaddasi L, Brown MP, et al. Intrinsic radiosensitivity is not the determining factor in treatment response differences between HPV negative and HPV positive head and neck cancers. Cells. 2020;9:8.

    Article  Google Scholar 

  72. Todorovic V, Groselj B, Cemazar M, Prevc A, Zakelj MN, Strojan P, et al. Expression of DNA-damage response and repair genes after exposure to DNA-damaging agents in isogenic head and neck cells with altered radiosensitivity. Radiol Oncol. 2022;56:173–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Ciecior W, Ebert N, Borgeaud N, Thames HD, Baumann M, Krause M, et al. Sample-size calculation for preclinical dose-response experiments using heterogeneous tumour models. Radiother Oncol. 2021;158:262–7.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Toustrup K, Sørensen BS, Lassen P, Wiuf C, Alsner J, Overgaard J. Gene expression classifier predicts for hypoxic modification of radiotherapy with nimorazole in squamous cell carcinomas of the head and neck. Radiother Oncol. 2012;102:122–9.

    Article  CAS  PubMed  Google Scholar 

  75. Rischin D, Hicks RJ, Fisher R, Binns D, Corry J, Porceddu S, et al. Prognostic significance of [18F]-misonidazole positron emission tomography-detected tumor hypoxia in patients with advanced head and neck cancer randomly assigned to chemoradiation with or without tirapazamine: a substudy of Trans-Tasman Radiation Oncology Group Study 98.02. J Clin Oncol. 2006;24:2098–104.

    Article  PubMed  Google Scholar 

  76. Mortensen LS, Johansen J, Kallehauge J, Primdahl H, Busk M, Lassen P, et al. FAZA PET/CT hypoxia imaging in patients with squamous cell carcinoma of the head and neck treated with radiotherapy: results from the DAHANCA 24 trial. Radiother Oncol. 2012;105:14–20.

    Article  PubMed  Google Scholar 

  77. Trinkaus ME, Hicks RJ, Young RJ, Peters LJ, Solomon B, Bressel M, et al. Correlation of p16 status, hypoxic imaging using [18F]-misonidazole positron emission tomography and outcome in patients with loco-regionally advanced head and neck cancer. J Med Imaging Radiat Oncol. 2014;58:89–97.

    Article  PubMed  Google Scholar 

  78. Kong CS, Narasimhan B, Cao H, Kwok S, Erickson JP, Koong A, et al. The relationship between human papillomavirus status and other molecular prognostic markers in head and neck squamous cell carcinomas. Int J Radiat Oncol Biol Phys. 2009;74:553–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Bernabeu MO, Köry J, Grogan JA, Markelc B, Beardo A, d’Avezac M, et al. Abnormal morphology biases hematocrit distribution in tumor vasculature and contributes to heterogeneity in tissue oxygenation. Proc Natl Acad Sci USA. 2020;117:27811–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Bodelon C, Untereiner ME, Machiela MJ, Vinokurova S, Wentzensen N. Genomic characterization of viral integration sites in HPV-related cancers. Int J Cancer. 2016;139:2001–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Lawrence MS, Sougnez C, Lichtenstein L, Cibulskis K, Lander E, Gabriel SB, et al. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature. 2015;517:576.

    Article  CAS  Google Scholar 

  82. Roberts S, Evans D, Mehanna H, Parish JL. Modelling human papillomavirus biology in oropharyngeal keratinocytes. Philos Trans R Soc Lond B Biol Sci. 2019;374:1773.

    Article  Google Scholar 

  83. Wei T, Buehler D, Ward-Shaw E, Lambert PF. An infection-based murine model for papillomavirus-associated head and neck cancer. MBio. 2020;11:3.

    Article  Google Scholar 

  84. Haeggblom L, Ramqvist T, Tommasino M, Dalianis T, Näsman A. Time to change perspectives on HPV in oropharyngeal cancer. A systematic review of HPV prevalence per oropharyngeal sub-site the last 3 years. Papillomavirus Res. 2017;4:1–11.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Onken MD, Winkler AE, Kanchi KL, Chalivendra V, Law JH, Rickert CG, et al. A surprising cross-species conservation in the genomic landscape of mouse and human oral cancer identifies a transcriptional signature predicting metastatic disease. Clin Cancer Res. 2014;20:2873–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Judd NP, Winkler AE, Murillo-Sauca O, Brotman JJ, Law JH, Lewis JS, et al. ERK1/2 regulation of CD44 modulates oral cancer aggressiveness. Cancer Res. 2012;72:365–74.

    Article  CAS  PubMed  Google Scholar 

  87. Kono M, Saito S, Egloff AM, Allen CT, Uppaluri R. The mouse oral carcinoma (MOC) model: a 10-year retrospective on model development and head and neck cancer investigations. Oral Oncol. 2022;132:106012.

    Article  CAS  PubMed  Google Scholar 

  88. Predina J, Eruslanov E, Judy B, Kapoor V, Cheng G, Wang LC, et al. Changes in the local tumor microenvironment in recurrent cancers may explain the failure of vaccines after surgery. Proc Natl Acad Sci USA. 2013;110:E415–24.

    Article  CAS  PubMed  Google Scholar 

  89. Predina JD, Kapoor V, Judy BF, Cheng G, Fridlender ZG, Albelda SM, et al. Cytoreduction surgery reduces systemic myeloid suppressor cell populations and restores intratumoral immunotherapy effectiveness. J Hematol Oncol. 2012;5:34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Pagura L, Cáceres JM, Cardinale A, Scharovsky OG, Di Masso RJ, Zacarías-Fluck MF, et al. A mammary adenocarcinoma murine model suitable for the study of cancer immunoediting. J Biomed Sci. 2014;21:52.

    Article  PubMed  PubMed Central  Google Scholar 

  91. Miyauchi S, Sanders PD, Guram K, Kim SS, Paolini F, Venuti A, et al. HPV16 E5 mediates resistance to PD-L1 blockade and can be targeted with rimantadine in head and neck cancer. Cancer Res. 2020;80:732–46.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Teja Valant (Institute of Oncology Ljubljana, Ljubljana, Slovenia) and Mira Lavric (Institute of Oncology Ljubljana, Ljubljana, Slovenia) for their technical help. We would also like to thank Maja Ota (Institute of Oncology Ljubljana, Ljubljana, Slovenia) for her help with preparation and staining of histological samples. We would like to thank Simona Miceska (Institute of Oncology Ljubljana, Ljubljana, Slovenia) for her help with preparation of cytospin slides, and Nina Mele (Biotechnical faculty, Ljubljana, Slovenia) for her help with immunofluorescence staining of cytospin slides and cells in vitro.

Funding

This work was supported by the Slovenian Research Agency (ARRS), Grant number P3-0003, and with bilateral research project funding, Grant number BT-AT/23-24-027.

Author information

Authors and Affiliations

  1. Department of Experimental Oncology, Institute of Oncology Ljubljana, Zaloska cesta 2, Ljubljana, Slovenia

    Ziva Modic, Maja Cemazar, Bostjan Markelc, Gregor Sersa, Simona Kranjc Brezar & Tanja Jesenko

  2. Faculty of Medicine, University of Ljubljana, Vrazov trg 2, Ljubljana, Slovenia

    Ziva Modic, Simona Kranjc Brezar & Tanja Jesenko

  3. Faculty of Health Sciences, University of Primorska, Polje 42, Izola, Slovenia

    Maja Cemazar

  4. Faculty of Health Sciences, University of Ljubljana, Zdravstvena pot 5, Ljubljana, Slovenia

    Bostjan Markelc & Gregor Sersa

  5. Department of Research, Valdoltra Orthopedic Hospital, Jadranska cesta 31, Ankaran, Slovenia

    Andrej Cör

  6. Faculty of Education, University of Primorska, Cankarjeva pot 5, Koper, Slovenia

    Andrej Cör

Authors
  1. Ziva Modic
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Maja Cemazar
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Bostjan Markelc
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Andrej Cör
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Gregor Sersa
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Simona Kranjc Brezar
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Tanja Jesenko
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

TJ and MC conceptualized, designed, and supervised the study; ZM, TJ, and BM optimized the methodology; ZM, TJ, and SKB conducted the experiments; ZM, TJ, MC, BM, and AC analyzed and interpreted the data; ZM prepared the original manuscript draft; TJ, MC, BM, SKB, AC, and GS reviewed and edited the manuscript; GS and SKB acquired the funding. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Maja Cemazar or Tanja Jesenko.

Ethics declarations

Ethics approval and consent to participate

All animal experiments were approved by the Ministry of Agriculture, Forestry and Food of the Republic of Slovenia (permission no. U34401-35/2020/8).

Consent for publication

Not applicable.

Competing interests

The authors confirm that there are no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Figure S1. MOC1-HPV cell lines show stable overall gene expression with cell culturing. Heatmap of differentially expressed genes with hierarchical clustering shows stable expression of differentially expressed genes between cells in passage 10 and passage 15. Figure S2. Cell migration kinetics of MOC1-HPV cell lines after extensive passaging. Cell front velocity of MOC1-HPV cells in the wound-healing assay after passaging for minimally 17 passages. Presented are mean ± SD. ***: p < 0.001, one-way ANOVA. Table S1. Most enriched GO Biological process pathways when comparing UM-SCC-104 (HPV-positive) and FaDu (HPV-negative) cell lines. Gene ontology analysis of differentially expressed genes. Pathways with an adjusted p value below 0.05 were considered to be enriched. Table S2. Most enriched GO Biological process pathways when comparing UM-SCC-47 (HPV-positive) and FaDu (HPV-negative) cell lines. Gene ontology analysis of differentially expressed genes. Pathways with an adjusted p value below 0.05 were considered to be enriched. Table S3. Most enriched GO Biological process pathways when comparing UM-SCC-090 (HPV-positive) and FaDu (HPV-negative) cell lines. Gene ontology analysis of differentially expressed genes. Pathways with an adjusted p value below 0.05 were considered to be enriched. Table S4. Most enriched GO Biological process pathways when comparing RNA-sequencing data from HPV-positive and HPV-negative patient tumor samples. Gene ontology analysis of differentially expressed genes. Pathways with an adjusted p value below 0.05 were considered to be enriched.

Additional file 2:

Supplementary methods.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Modic, Z., Cemazar, M., Markelc, B. et al. HPV-positive murine oral squamous cell carcinoma: development and characterization of a new mouse tumor model for immunological studies. J Transl Med 21, 376 (2023). https://doi.org/10.1186/s12967-023-04221-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-04221-4

Keywords

Journal of Translational Medicine

ISSN: 1479-5876

Contact us