A novel human in vitro papillomavirus type 16 positive tonsil cancer cell line with high sensitivity to radiation and cisplatin
BMC Cancer volume 19, Article number: 265 (2019)
Human papillomavirus (HPV) is an established risk factor for oropharyngeal squamous cell carcinoma (OSCC). The aim was to establish cell lines from HPV-positive tonsil carcinomas to be used for treatment development.
Fresh samples from 23 HPV-positive tonsil carcinomas were cultivated in vitro. The established cell line was analyzed for viral characteristics, cell karyotype, TP53 status, and growth capabilities in nude mice. In vitro studies of sensitivities to radiation, cisplatin and cetuximab were performed.
After 19 months (eight passages), one cell line, LU-HNSCC-26, was established in vitro and also grew as xenografts. The tumor was from a 48 year old non-smoking man with non-keratinizing, p16 positive tonsil OSCC, stage T2N0M0 with HPV16. It contained 19.5 (CV% 3.7) HPV16 copies/cell (passage 8). The complete HPV16 genome sequence was obtained. Episomal HPV16 was present with an E2/E7 ratio of 1.1 (CV% 2.6). In addition, HPV16 mRNA specific for the intact E2 gene was detected. The viral expression manifested 1.0 (CV% 0.1) E7 mRNA copies per HPV16 genome. The karyotype was determined and the cell line demonstrated wild type TP53. The ID50 for radiation was 0.90 Gy and the IC50 for cisplatin was 0.99 μmol/L. The cell line was inhibited to a maximum of 18% by cetuximab.
We established an in vitro tonsil carcinoma cell line containing episomal HPV16. This is an important step towards efficient treatment development.
In recent years, there has been an increase in the incidence of oropharyngeal squamous cell carcinoma (OSCC) including cancers in the tonsils and the base of the tongue [1, 2]. In Sweden, the increase has been close to 5% per year between 1970 and 2016  which represents a doubling in incidence every 14 years. The patients in this group are generally younger and often lack the traditional risk factors for head and neck squamous carcinoma (HNSCC) such as a history of smoking and alcohol abuse .
There is strong evidence that this rise in OSCC incidence is caused by human papillomavirus (HPV) infection  and that this also explains an observed shift to higher survival rates for this cancer type  as patients with HPV-positive OSCC have a better overall survival and better treatment response compared to patients with HPV-negative OSCC .
In many respects, including incidence, survival and treatment response, HPV-positive and HPV-negative OSCC behave differently, and it has become clear that the treatment strategies should differ significantly between the groups and that new regimes should be developed specifically for HPV-positive OSCC.
An important tool for testing new drugs or combinations of drugs as well as for the elucidation of the biology of HPV-induced transformation is representative, in vitro and in vivo, HPV-positive OSCC cell lines. So far, a small number of HPV16-positive HNSCC cell lines have been established in vitro but none originating from tonsil tumors and all from patients with a history of smoking and alcohol use [7, 8]. As the outcome is different for tumors linked to HPV compared to those linked to the classical risk factors, and for HPV-positive tumors in the oropharynx (especially in the tonsils) compared to those originating in other sites [9,10,11], it is very important to establish cell lines from tonsil tumors in order to perform valid pre-clinical therapy development.
For uterine cervical cancer, the integration process has been shown to be of great importance for the carcinogenic process, though not being an exclusive prerequisite as episomal HPV16 was present in 47% of cervical HPV16-positive cancers, either as pure episomal (27%) or as a mixture of integrated and episomal HPV16 (20%) . Measurements of the physical state of the HPV genome in OSCC have shown very different frequencies of integration, possibly as a consequence of a combination of different methods used to assess integration status and varying sample sizes. But it has been reported that, among HPV16 positive cases, 76% of OSCCs and 93% of tonsil carcinomas manifest episomal HPV16 either as pure episomal or as the mixed status [13,14,15], thus implicating that the frequency of the episomal state of HPV16 is more common in OSCC than in cervical cancer. To the best of our knowledge, no cell line with episomal HPV16 derived from tonsil carcinoma has been established in vitro.
The aim of the presented work was to establish cell lines in vitro from HPV-positive OSCC tumors that could be used for therapy development and basic research purposes.
Materials and methods
The study was approved by Lund University regional ethical review board (LU 376–01), and the samples were collected, after informed consent was obtained, at Skåne University Hospital, Lund, Sweden.
Establishment of cell line in vitro
Tumor samples were rinsed consecutively in phosphate buffered saline (PBS) and R10 medium (RPMI 1640 with stable glutamine supplemented with 1 mmol/L sodium pyruvate, 1 × MEM non-essential amino acids, 20 μg/mL gentamicin and 10% fetal bovine serum (FBS), all from GE Healthcare (Piscataway, NJ, USA). The tumor samples were cut in pieces and transferred to cell culture flasks with R10 medium and left to attach and grow at 37 °C under a humidified atmosphere with 5% CO2. The cells used in the described experiments were passaged less than 18 times after establishment. Mycoplasma tests were performed on a regular basis by GATC Biotech AG (Konstanz, Germany) and were negative.
For cell doubling determination, cells were seeded in eight 96-well plates. The plates were analyzed by the sulphorhodamine B (SRB) method as previously described  consecutively at the indicated time points. Error bars indicate standard error of the mean. The data points were fitted to a logarithmic equation using the GraphPad Prism (5.04) software package (GraphPad Software, La Jolla, CA, USA).
Single tandem repeat analysis
Single tandem repeat (STR) analysis was performed by Uppsala Genome Center using the AmpFlSTR Identifiler polymerase chain reaction (PCR) amplification kit (Applied Biosystems, Carlsbad, CA, USA). Separation of DNA fragments was performed with capillary electrophoresis on the ABI3730XL DNA analyzer (Applied Biosystems, Foster City, CA, USA) and allele scoring was performed with the GeneMapper 4.0 software package.
Establishment of xenograft
Experiments on mice were approved by the local ethics committee (488/2017). We used in-house bred, athymic 5- to 8-week-old BALB/c nude (nu/nu) mice. The mice were given food and water ad libitum and treated according to regulations issued by the Swedish Board of Agriculture (SJVFS 2017:40). Cells were injected subcutaneously into the flanks of non-anesthetized animals. To determine the doubling time, tumors were retransplanted to 10 nude mice under ketamine anesthesia, two tumors per mouse. The length and width of the tumors were measured with calipers three times per week and the relative tumor size calculated in relation to day 0 (set as the day when all tumors had shown growth over 3 consecutive time points). The data points were fitted to a logarithmic equation in GraphPad Prism. After the measuring period, the mice were euthanized using cervical dislocation and the tumors removed for histological sectioning.
Xenografts and cell pellets were sectioned in 3 μm slices and stained with hematoxylin and eosin, and the following antibodies: the CINtec p16-kit (E6H4, cat. no. 725–4713), anti-pan keratin, (AE1/AE3/PCK26, cat. no. 760–2595), and anti-p63 (4A4, cat. no. 790–4509) all from Ventana Medical Systems (Tucson, AZ, USA), anti-cytokeratin 5 (XM26, cat. no. NCL-L-CK5) from Leica Biosystems (Nussloch, Germany), and p40 (BC28, cat. no. ACI3066C) from Biocare Medical (Pacheco, CA, USA), All stainings were performed using the BenchMark Ultra platform (Ventana Medical Systems).
Three HNSCC cell lines, previously established in our laboratory, were used for comparison in several experiments: LU-HNSCC-4, LU-HNSCC-5 and LU-HNSCC-6 . They were propagated in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% FBS, 100 units/mL penicillin, and 100 units/mL streptomycin sulfate, and used within 15 passages after thawing. Single tandem repeat analysis was performed before the experiments as detailed above showing no cross-contamination between the cell lines or with other common contaminants. Tests for mycoplasma infection prior to the experiments were negative.
A 2–3 mm biopsy was immersed in 1 mL saline solution and immediately transferred to the laboratory. The saline was removed and the biopsy was digested in lysis buffer and the DNA extracted with the Total NA-kit using MagNA Pure LC  (Roche, Pleasanton, CA, USA). Sample adequacy was assessed by PCR for the human beta-globin gene . Identification of 39 genital HPV types was carried out by modified general primer polymerase chain reaction (MGP-PCR) and subsequent Luminex analysis (Luminex, Austin, TX, USA) [18, 20, 21].
Samples from subsequent passages in cell cultures were immersed in 1 mL GITS (4 mol/L guanidinium thiocyanate, 22 mmol/L sodium citrate, and 5% N-lauroylsarcosine sodium salt) and stored at − 20 °C. After thawing, 200 μL was used for DNA extraction with the Total NA-kit using MagNA Pure LC (100 μL output).
Briefly, the number of viral genomes per cell was quantified by carrying out two separate real-time PCR tests, including appropriate standard curves, to amplify the HPV16 E7 gene and the human beta-globin gene as previously described . The samples were analyzed in duplicate.
Physical status of HPV16
As a marker for the presence of mixed or episomal forms of HPV16, samples were analyzed for the ratio of E2/E7 gene copy numbers. HPV16 was classified as mixed status when complete E2 was evident from complete sequencing of the HPV16 genome and the E2/E7 ratio was < 0.9, and pure episomal status in the presence of complete E2 open reading frame (ORF) and E2/E7 ratio of > 0.9. The E7 copy numbers were retrieved from the viral load measurement as described above. For quantification of the E2 copy numbers we used a modified version of the PCR described by Peitsaro et al. , where E2/E6 ratios greater than 1 indicate predominance of the episomal form . Quantification was extrapolated from a linear regression standard curve that was included in each batch assay, as previously described . The samples were analyzed in duplicate.
Long range HPV16 PCR
The complete HPV16 genome was amplified in a 25 μL solution with 2.5 μL (15.5 ng) of extracted nucleic acid from passage 2, 0.2 μmol/L of each primer (HPV16 7465F, 5′-ATGCTTTTTGGCACAAAATGTG, HPV16 7464R 5′-GCAACCGAATTCGGTTGAAG), 200 μmol/L of each dNTP, 1 × PrimeStar GXL Buffer, and 0.625 U PrimeStar GXL DNA Polymerase (TaKaRa Bio, Shiga, Japan). The PCR was performed at 98 °C for 10 s, 60 °C for 10 s, and 68 °C for 8 min for 45 cycles. Amplified DNA was separated by electrophoresis in a 1.0% agarose gel (agarose Type I-A, A0169, Sigma-Aldrich, Stockholm, Sweden) with 1 × Gel RedTM Nucleic Acid Gel Stain (VWR International, Lund, Sweden) in 0.5 × TBE buffer  and visualized by UV light.
DNA sequence of the complete HPV16 genome
The complete HPV16 genome was obtained by next-generation sequencing (NGS) in the MiSeq system (Illumina, San Diego, CA, USA) from passage 2.
Briefly, extracted DNA (5.6 ng/μL in a total volume of 20 μL) was whole genome amplified (WGA) with multiple displacement amplification using the Illustra™ Ready-To-Go™ GenomiPhi™ DNA Amplification Kit (GE Healthcare) based on the manufacturer’s protocol, but with some modifications. Five microliter of extracted DNA was diluted with 20 μL PCR-grade water and 25 μL denaturation buffer, incubated at 95 °C for 3 min, cooled on ice, added to the lyophilized cake, incubated at 30 °C for 7 h and inactivated at 65 °C for 10 min. The WGA product was diluted 1:2 in PCR-grade water and quantified with QuantiFluor-ST (Promega, Madison, WI, USA). A DNA library with 2 unique indices was prepared from 50 ng of the WGA product using the Nextera DNA Sample Preparation kit, user guide revision B (Illumina) and validated by fragment size using 2100 Bioanalyzer High sensitivity DNA assay (Agilent Technologies, Santa Clara, CA, USA) and quantified with QuantiFluor-ST (Promega). The library was normalized to 4 nmol/L in EB-buffer (Qiagen, Hilden, Germany), denatured, diluted to 1.8 pmol/L, spiked with 1% PhiX control and paired-end sequenced by 151 plus 151 cycles using NextSeq 500 High Output reagent kit (Illumina) as described in the user guides Denaturing and Diluting Libraries for the NextSeq 500, revision A and NextSeq 500 kit Reference Guide, revision F.
Papillomavirus-related contigs having > 90% identity with each other were clustered. The total number of sequencing reads from each cluster was calculated and the longest contig was selected as a representative to remove redundancy. A de novo HPV-sequence was obtained by the use of clustering analysis, conducted separately for contigs longer and shorter than 700 bp length. For the HPV type, the total number of reads was calculated. All analyses were performed using in-house R (http://www.r-project.org/), python (http://www.python.org) and bash (http://www.gnu.org/software/bash/) scripts that ran on a high performance (40 core, 2 TB RAM) Linux server.
The obtained sequence was then analyzed by the use of BioEdit v7.0  for presence of the expected ORFs and in cases of ambiguities, re-sequencing from the passage 2 was performed to resolve the ambiguity.
The primary de novo sequence had two segments of N-ambiguities of 72 bases and of 129 bases at positions 4633–4704 and 7150–7278, respectively. In order fill these gaps, two PCRs were performed, each with 50 μL PCR mixture that contained 5 μL of template of extracted DNA from passage 2, 0.3 μmol/L of each forward and reverse primer, 200 μmol/L of each dNTP (Roche), 1.5 U AmpliTaq Gold, 1 × reaction buffer and 3.0 mmol/L MgCl2 (all from Applied Biosystems). The following primer pairs were used to generate amplicons of 887 bp with HPV16 forward primer at position 3992 (5′- TATTGTGGATAACAGCAGCCTCTG) and reverse primer at position 4878 (5′- GGTGTGCTACTAGTTACTGTGTTAGGG), and of 316 bp with HPV16 forward primer at position 7048 (5′-AAGCAGGATTGAAGGCCAAAC) and reverse primer at position 7363R (5′- AATAACCACAACACAATTAGTAGGTGTTG). We also separately sequenced a fragment of 250 bp in order to verify the guanosine at position 1376, by a HPV16 forward primer at position 1250 (5′- GCGAAGACAGCGGGTATGG) and a reverse primer at position 1489 (5′- TGCTAACATTGCTGCCTTTGC).
PCR was carried out in an automated thermocycler (Eppendorf) programmed for tube setting: with the following parameters: 94 °C for 10 min; 45 × (95 °C for 15 s, 55 °C for 30 s, and 72 °C for 1 s). Five microlitres of the amplicons were separated by electrophoresis in a 2.0% agarose gel (agarose Type I-A, A0169, Sigma-Aldrich) with 1 × Gel RedTM Nucleic Acid Gel Stain (VWR International, Lund Sweden) in 1 × TBE buffer  and visualized by UV light. The remaining amplified material (45 μL) was purified by Illustra MicroSpin S-300 HR Columns (GE Healthcare) and sequenced with the HPV16 primers described above (Eurofins, Ebersberg, Germany). The sequences were then used to fil the gaps of the de novo sequence and the final sequence was compiled.
mRNA extraction and quantitative mRNA RT-PCR
RNA extraction and RT-PCR
Total RNA was extracted using TRI Reagent (Sigma-Aldrich) and Direct-zol RNA MiniPrep (ZYMO Research, Irvine, CA, USA) according to the manufacturer’s protocol. One μg of total RNA was reverse transcribed in a 20 μL reaction at 37 °C by using M-MLV Reverse Transcriptase and random primers (both from Invitrogen, Carlsbad, CA, USA) according to the protocol of the manufacturer. One microliter of cDNA was subjected to PCR amplification with primers indicated in Fig. 3 and listed in Table 1.
Protein extraction and Western blotting
Cells were lysed in radioimmunoprecipitation assay buffer (50 mmol/L Tris, pH 7.8, 150 mmol/L NaCl, 1% sodium deoxycholate, 0.1% sodium dodecyl sulphate, 1% Triton X-100, 1 mmol/L dithiothreitol and protease inhibitors), subjected to SDS-PAGE, transferred to nitrocellulose membrane, and stained with antibodies against HPV16 L1 (ab30908, Abcam, Cambridge, UK), L2 (sc-65,708, Santa Cruz Biotechnology, Dallas, TX, USA), and epidermal growth factor receptor (EGFR) (no. 2232, Cell Signaling Technology, Danvers, MA, USA). Proteins were detected using chemiluminescence. For quantification of EGFR, loading control was performed using Coomassie R-350 as previously described .
Screening of mutations in TP53
Primers for all exons of TP53 were designed using the Primer3 software (http://primer3.ut.ee/); primer sequences are available on request. PCR was done on genomic DNA according to standard methods and amplified fragments were sequenced using the BigDye v1.1 Cycle Sequencing Kit on an ABI-3130 Genetic Analyzer (Applied Biosystems). ChromasLite 2.1.1 free online software (http://technelysium.com.au/) was used to analyze the data.
For cytogenetic analysis, cells were cultured in DMEM/F12 medium with 10% fetal bovine serum and antibiotics until the mitotic index was high, and then exposed to colchicine (0.03 μg/mL) for 4–5 h. The cells were detached by trypsin-EDTA, treated with 0.06 mol/L KCl for 35 min, and then fixed three times with methanol/acetic acid 3:1. Chromosome preparations were incubated at 60 °C overnight and then treated with 0.3 mol/L sodium chloride, 30 mmol/L sodium citrate, pH 7.0, at 60 °C for 4 h. G-banding was done using Wright’s stain and karyotypes were described according to ISCN 2009 .
Determination of treatment sensitivity
For cisplatin, cells were seeded and, after 48 h, treated with increasing concentrations of cisplatin (Sandoz AS, Copenhagen, Denmark) for one hour. After 5 days at 37 °C under a humidified atmosphere with 5% CO2, the cell amounts were measured by the SRB assay.
For determination of cetuximab sensitivity, cells were treated with increasing concentrations of cetuximab (Merck AB, Solna, Sweden) for 5 days, after which the cell amounts were determined by the SRB assay.
Radiosensitivity was analyzed by irradiating the cells in a Gamma Cell-3000 Elan from MDS Nordion (Ottawa, ON, Canada) equipped with a Cs-137 source. After irradiation, the cells were seeded and left to grow until the untreated controls were approximately 90% confluent, and then analyzed by the SRB assay.
For all experiments the data were fitted to a sigmoidal dose-response equation with variable slope using non-linear regression, and the half maximal inhibitory concentration (IC50), inhibitory dose (ID50) or maximum inhibition (Emax) calculated in GraphPad Prism.
Establishment of cell line
A total of 27 tumor samples (23 HPV-positive), were collected and treated as detailed in the Materials and Methods section. In a single case, we were able to establish a cell line. The original tumor tissue was obtained from a 48 year old non-smoking man with moderate daily consumption of alcohol, and, besides exercise induced asthma, no other intercurrent diseases. The initial symptom was left-sided pain in the throat. He had no weight loss. Three days after his first symptom he was examined by an ENT-specialist who noted a prominent left sided tonsil without mucosal ulceration. No nodes could be felt. A fine-needle aspirate was performed, with negative result. Two and a half weeks later he was examined under anesthesia and biopsied through intact mucosa. This biopsy confirmed a non-keratinizing, p16 (CDKN2A) positive squamous cell tonsil carcinoma with HPV16 detected in the original biopsy. After further investigation with CT and CT-PET scan, the tumor was classified as stage T2N0M0. The patient was treated with full dose radiotherapy (Intensity-modulated radiation therapy) without concomitant chemotherapy leading to full remission, and was relapse-free after four and a half years. The primary tumor received 64 Gy (2.0 Gy/fraction) and the neck 54.4 Gy (1.6 Gy/fraction).
Approximately 19 months (8 passages) after the initial seeding, the cell line was visibly free from fibroblasts and displayed a stable growth rate (Fig. 1a) and was denominated LU-HNSCC-26. The doubling time was determined to be 84 h (Fig. 2a).
For future identification and authentication of the cell line, STR profiling was performed using 16 loci (Table 2). No exact matches were found when comparing with other analyzed cell lines in the DSMZ repository .
Cell pellets were sectioned and stained for p16 and epithelial markers. The cells were p16-positive (cytoplasmic and nuclear staining) (Fig. 3a). They were also positive for CKAE1/AE3 and CK5 (cytoplasmic staining) (Fig. 3b and c), and p63 and p40 (strong nuclear staining) (Fig. 3d and e), establishing the squamous epithelial origin of the cells.
The LU-HNSCC-26 cells were injected subcutaneously in nude mice resulting in xenografts that could be successfully retransplanted to new mice. The doubling time of tumor growth was determined to be 177 h (Fig. 2b). The histology of the resulting tumors was compared with that of the original tumor. Both tumors showed non-keratinizing morphology, growing in solid sheets and trabeculae. Scattered mitoses and apoptotic cells was seen in both, whereas prominent tumor necrosis was only is observed in the xenograft (Fig. 1b and c). Both xenograft and original tumor were strongly positive for p16 (Fig. 1d and e).
Detection and characterization of HPV
At passage two and eight the viral DNA load was 14.8, coefficient of variation (CV%) 6.6, and 19.5 (CV% 3.7) copies per cell, respectively. Long range PCR for detection of the complete HPV16 genome from passage two and eight demonstrated amplicons of approximately 8 kb (Fig. 4). Episomal HPV16 was present in the cell line with E2/E7 ratios of 0.8 (CV% 1.7) and 1.1 (CV% 2.6) from passage two and eight, respectively. The viral expression manifested 1.0 (CV% 0.1) E7 mRNA copies per HPV16 genome of passage two.
To investigate if correctly spliced HPV16 mRNAs were produced, we extracted total RNA (Fig. 5b) and performed RT-PCR with various HPV16 specific primers (Fig. 5a). We found production of mRNAs encoding full-length E6, E6*I, E6*II and E7 proteins. The E7 mRNA spliced from SD226 to SA409 was the most abundant of the early mRNAs (Fig. 5c).). Both E2 mRNAs and mRNA spliced from SD880 to SA3358 were detected (Fig. 5c), confirming that episomal copies of the HPV16 genome were present and transcriptionally active. We detected mRNAs spliced from SD880 to the exclusively late 3′-splice site SA5639, suggesting that L1i mRNAs were produced (Fig. 5c), and L2 mRNAs were also detected. These results suggested that HPV16 late proteins L1 and L2 could be produced in the LU-HNSCC-26 cell. We therefore performed a Western blot which confirmed the presence of L1 and L2 proteins in the LU-HNSCC-26 cells (Fig. 5d), albeit at very low levels.
From passage two, the complete DNA sequence of HPV16 was obtained by NGS with a minimum nucleotide coverage of 134 (Fig. 6). The HPV16 genome comprised 7904 bp, GenBank accession KY994539. The genomic organization showed the characteristic features shared by HPV types. The HPV16 genome showed closest identity (99%, 7900/7906) to the HPV16 strain CU4 (GenBank accession FJ610149), originally isolated from a uterine cervical cytology sample with CINI. In comparisons to this isolate, our HPV16 demonstrated a deletion of two bp in the second non-coding region between the E5 and L2 ORFs. Furthermore, within the E1 ORF at nucleotide position 1376, a guanosine was present instead of adenosine, translating to cysteine instead of tyrosine in the putative E1 protein.
Genetic characterization of the cell line
The karyotype was determined to be 43–47,XY,del(2)(q33), del(3)(p21),der(4)t(4;5)(q35;q31), add(5) (q35),-7,-14,add (15) (p11),?i(16)(q10),-21,+mar[cp6]/44–47,idem,i(8)(q10)[cp3] (see karyogram in Additional file 1).
The gene TP53 was wild-type but homozygous for a single nucleotide polymorphism chr17:7579472G > C, giving a proline to arginine substitution at position 72 (see partial sequencing results in Additional file 2).
The sensitivities to commonly used treatments were assessed in vitro and compared to three HPV-negative HNSCC cell lines. LU-HNSCC-26 was considerably more sensitive to cisplatin than the other HNSCC cell lines tested (Fig. 7a and Table 3). The IC50 for LU-HNSCC-26 (0.99 μmol/mL) was 10–20 times lower than those of the other cell lines.
For cetuximab, we compared the maximum inhibition (Emax) of the drug instead of the IC50. LU-HNSCC-26 was only inhibited to a maximum of 18% (Fig. 7b and Table 3), which was less than the other cell lines . This difference could not be correlated with the EGFR expression as assessed by Western blotting showing a similar expression in all four cell lines (Fig. 7d).
As for cisplatin, LU-HNSCC-26 was considerably more sensitive to radiation compared to the other cell lines (Fig. 7c and Table 3). The ID50 for LU-HNSCC-26 (0.90 Gy) was 18–26 times lower than those of the other cell lines.
In the current work we demonstrated the in vitro establishment of an HPV16-positive cell line, LU-HNSCC-26, derived from a tonsil carcinoma. The cell line manifested episomal HPV16 with active expression of viral mRNA. The TP53 gene was wild type and the cell line displayed a relatively simple karyotype including loss of 3p, 4q, and 7p as well as isochromosomes, which are common findings in HPV-positive HNSCC-cases, that often have less complex karyotypes with fewer signs of chromosomal instability compared to HPV-negative ones . The doubling time in vitro was 84 h which is approximately 3 times slower than most other HNSCC cell lines established in our laboratory. The cell line was also grown as a xenograft in nude mice and the tumors displayed a histology that closely matched that of the original tumor. The cell line showed high sensitivity to cisplatin and to radiation but was less sensitive to treatment with cetuximab, compared with that of three HPV-negative HNSCC cell lines previously described by us .
It is widely accepted that HPV infection plays an etiological role in many oropharyngeal cancers and there is ample evidence that HPV-positivity is coupled to increased survival and better prognosis. As a result of this, there are currently several clinical trials addressing the question whether it is possible to de-intensify the treatment for this patient group to reduce the toxicities while preserving the efficacy . However, as the carcinogenic process most certainly is different for HPV-driven tumors, there is also an opportunity to increase the efficacy of treatment for these patients by targeting mechanisms specific to HPV-positive cancers.
Well characterized cellular and animal models are valuable tools for efficient treatment development. Such models can provide fast and efficient testing of novel drugs and combinations of drugs or other treatment modalities. So far, no in vitro cell line has been established from HPV-positive tonsil cancer and only few from other head and neck locations [7, 8]. Importantly, there are clear indications that, in contrast to OSCC, the outcome for HPV-positive cases in sites other than oropharynx is not favorable compared to HPV-negative ones [9,10,11]. Cell lines derived from other sites, though probably mirroring the properties of the originating tumors , might therefore not correctly represent the properties of HPV-positive OSCC. Thus, the presented cell line provides an important step forward in the quest for improved therapies against HPV-positive OSCC and seems to be representative for this cancer type in the measured variables.
The presence of the episomal state of the HPV16 in the LU-HNSCC-26 cell line, as demonstrated by the complete circular genome of HPV 16 and as well as the intact HPV16 E2 mRNA, is similar to that of our previously established HPV16-positive tonsil carcinoma xenograft, which manifested episomal HPV16 and integrated HPV16 simultaneously . Overall, the episomal form of HPV16, either as pure episomal or as a mixed status of both integrated and episomal HPV16, is in predominance in tonsil carcinomas [13, 14, 32]. The relatively high proportion of intact episomal HPV16 in these lesions may lead to repression of the viral oncogene transcription via expression of intact HPV16 E2-protein, which binds to conserved sites of the upper regulatory region of HPV16 genome . Interestingly, the latency period between HPV infection and oropharyngeal cancer development has been estimated to be up to 30 years whereas the corresponding period for cancer of the cervix uteri is around 20 years . Speculatively, the longer latency period for HPV-positive tonsil carcinomas compared to that of cervical cancer might partly be explained by a predominance of the episomal HPV16 among tonsil carcinomas.
The HPV16 positive LU-HNSCC-26 cell line manifested wild-type TP53, homozygous for arginine at codon 72 of the TP53. Interestingly, this arginine form of the gene product p53 has been shown to be more susceptible to HPV16 E6 mediated degradation than the proline form . In addition, the arginine form of p53 was found to be associated with HNSCC, but not with worse prognosis of the disease . Furthermore, Habbous et al. reported in a meta-analysis that the arginine variant of p53 was associated with progression of squamous intraepithelial lesions to cervical cancer only in the presence of HPV positivity . Further studies are required to investigate the link between the arginine variant of p53 and susceptibility for development of tonsil carcinoma among HPV-infected individuals.
The p53 protein is inactivated by the HPV16 E6 protein . Here we demonstrated HPV16 mRNAs encoding full-length E6 protein as well as E6*I, E6*II and E7, strongly supporting the idea that p53 is inactivated by HPV16 E6 protein in LU-HNSCC-26 cells. The splicing pattern of the various E6 and E7 mRNAs is reminiscent of the splicing pattern observed in HPV16-positive cervical cancer cell lines with relatively low levels of the mRNA encoding the full-length E6 protein, while the most abundant early mRNA is coding for E6*I and E7. The relative levels of the various HPV16 E6 and E7 mRNAs are determined by cis-acting RNA elements on the HPV16 mRNAs, and by cellular proteins that bind to these elements to control splicing efficiency [39, 40]. One may speculate that the levels of these cellular splicing regulatory proteins are similar in cervical cancer and tonsil cancer cells. Interestingly we detected HPV16 L1 and L2 proteins from the LU-HNSCC-26 cell line. However, initial attempts to identify virus particles from the cell line have been unsuccessful (data not shown, personal communication H. Faust).
In this work we attempted to set up cell lines from 27 tonsil tumors of which 23 were HPV positive but only one of these actually resulted in an established cell line. This success rate is considerably lower than that for HNSCC in general, previously reported to be around 33% in our lab . In our hands, the establishment of HPV-positive tonsil cancer cell lines has thus been much more difficult than that of other HNSCC tumor types. The low number of HPV-negative tumors included in this study precludes any conclusion about the take rate of these cells.
There are around 14 tonsil cell lines described in the literature  of which, as far as we know, all are negative for HPV, though one in vitro cell line from tonsil fossa with only integrated HPV16 was recently reported . Overall, this implies that HPV-positive tonsil cancers are indeed less prone to grow in vitro than the HPV-negative ones, which is surprising considering the ease with which HPV-positive cell lines have been established from cervical cancers.
Conflicting results regarding the sensitivity of HPV-positive HNSCC cell lines to cisplatin have previously been reported indicating either similar sensitivities between HPV-positive and negative cell lines  or higher sensitivities in HPV-positive cell lines . This is remarkable as several of the cell lines used were identical between the two studies, but indicates that the different results depend on methodological issues. In our work, LU-HNSCC-26 was clearly more sensitive to cisplatin than the three HPV-negative HNSCC cell lines used for comparison, with the IC50 being 10 to 20 times lower (Table 3). In combination with the low sensitivity of the cell line to cetuximab (maximum 18% inhibition), this is interesting considering the current discussion on deescalating therapies by omitting cisplatin or exchanging with targeting drugs such as cetuximab; an alternative not suggested by the present results.
To the best of our knowledge, we have established the first in vitro tonsil carcinoma cell line containing episomal HPV16. The cell line demonstrated high sensitivity to radiation and cisplatin treatment. The cell line provides an important experimental system of episomal HPV16-positive tonsil carcinomas in an era of increasing incidence of tonsil cancer in several countries.
coefficient of variation
Dulbecco’s modified Eagle’s medium
epidermal growth factor receptor
- Emax :
fetal bovine serum
head and neck squamous cell carcinoma
- IC50 :
half maximal inhibitory concentration
- ID50 :
half maximal inhibitory dose
modified general primer polymerase chain reaction
open reading frame
oropharyngeal squamous cell carcinoma
phosphate buffered saline
polymerase chain reaction
single tandem repeat
- Td :
whole genome amplified
Chaturvedi AK, Engels EA, Pfeiffer RM, Hernandez BY, Xiao W, Kim E, Jiang B, Goodman MT, Sibug-Saber M, Cozen W, et al. Human papillomavirus and rising oropharyngeal cancer incidence in the United States. J Clin Oncol. 2011;29(32):4294–301.
Hammarstedt L, Lindquist D, Dahlstrand H, Romanitan M, Dahlgren LO, Joneberg J, Creson N, Lindholm J, Ye W, Dalianis T, et al. Human papillomavirus as a risk factor for the increase in incidence of tonsillar cancer. Int J Cancer. 2006;119(11):2620–3.
Socialstyrelsens statistikdatabas. https://www.socialstyrelsen.se/statistik/statistikdatabas.
Gillison ML, Koch WM, Capone RB, Spafford M, Westra WH, Wu L, Zahurak ML, Daniel RW, Viglione M, Symer DE, et al. Evidence for a causal association between human papillomavirus and a subset of head and neck cancers. J Natl Cancer Inst. 2000;92(9):709–20.
Ang KK, Harris J, Wheeler R, Weber R, Rosenthal DI, Nguyen-Tan PF, Westra WH, Chung CH, Jordan RC, Lu C, et al. Human papillomavirus and survival of patients with oropharyngeal cancer. N Engl J Med. 2010;363(1):24–35.
Dayyani F, Etzel CJ, Liu M, Ho CH, Lippman SM, Tsao AS. Meta-analysis of the impact of human papillomavirus (HPV) on cancer risk and overall survival in head and neck squamous cell carcinomas (HNSCC). Head Neck Oncol. 2010;2:15.
Tang AL, Hauff SJ, Owen JH, Graham MP, Czerwinski MJ, Park JJ, Walline H, Papagerakis S, Stoerker J, McHugh JB, et al. UM-SCC-104: a new human papillomavirus-16-positive cancer stem cell-containing head and neck squamous cell carcinoma cell line. Head Neck. 2012;34(10):1480–91.
Olthof NC, Huebbers CU, Kolligs J, Henfling M, Ramaekers FC, Cornet I, van Lent-Albrechts JA, Stegmann AP, Silling S, Wieland U, et al. Viral load, gene expression and mapping of viral integration sites in HPV16-associated HNSCC cell lines. Int J Cancer. 2015;136(5):E207–18.
Chakravarthy A, Henderson S, Thirdborough SM, Ottensmeier CH, Su X, Lechner M, Feber A, Thomas GJ, Fenton TR. Human papillomavirus drives tumor development throughout the head and neck: improved prognosis is associated with an immune response largely restricted to the oropharynx. J Clin Oncol. 2016;34(34):4132–41.
Chung CH, Zhang Q, Kong CS, Harris J, Fertig EJ, Harari PM, Wang D, Redmond KP, Shenouda G, Trotti A, et al. p16 protein expression and human papillomavirus status as prognostic biomarkers of nonoropharyngeal head and neck squamous cell carcinoma. J Clin Oncol. 2014;32(35):3930–8.
Lohaus F, Linge A, Tinhofer I, Budach V, Gkika E, Stuschke M, Balermpas P, Rodel C, Avlar M, Grosu AL, et al. HPV16 DNA status is a strong prognosticator of loco-regional control after postoperative radiochemotherapy of locally advanced oropharyngeal carcinoma: results from a multicentre explorative study of the German Cancer consortium radiation oncology group (DKTK-ROG). Radiother Oncol. 2014;113(3):317–23.
Cullen AP, Reid R, Campion M, Lorincz AT. Analysis of the physical state of different human papillomavirus DNAs in intraepithelial and invasive cervical neoplasm. J Virol. 1991;65(2):606–12.
Deng Z, Hasegawa M, Kiyuna A, Matayoshi S, Uehara T, Agena S, Yamashita Y, Ogawa K, Maeda H, Suzuki M. Viral load, physical status, and E6/E7 mRNA expression of human papillomavirus in head and neck squamous cell carcinoma. Head Neck. 2013;35(6):800–8.
Faust H, Eldenhed Alwan E, Roslin A, Wennerberg J, Forslund O. Prevalence of human papillomavirus types, viral load and physical status of HPV16 in head and neck squamous cell carcinoma from the south Swedish health care region. J Gen Virol. 2016;97(11):2949–56.
Lim MY, Dahlstrom KR, Sturgis EM, Li G. Human papillomavirus integration pattern and demographic, clinical, and survival characteristics of patients with oropharyngeal squamous cell carcinoma. Head Neck. 2016;38(8):1139–44.
Yamatodani T, Ekblad L, Kjellen E, Johnsson A, Mineta H, Wennerberg J. Epidermal growth factor receptor status and persistent activation of Akt and p44/42 MAPK pathways correlate with the effect of cetuximab in head and neck and colon cancer cell lines. J Cancer Res Clin Oncol. 2009;135(3):395–402.
Henriksson E, Kjellén E, Baldetorp B, Bendahl PO, Borg Å, Brun E, Mertens F, Ohlsson T, Rennstam K, Wennerberg J, et al. Comparison of cisplatin sensitivity and the 18F fluoro-2-deoxy 2 glucose uptake with proliferation parameters and gene expression in squamous cell carcinoma cell lines of the head and neck. J Exp Clin Cancer Res. 2009;28:17.
Letsolo BT, Faust H, Ekblad L, Wennerberg J, Forslund O. Establishment and characterization of a human papillomavirus type 16-positive tonsillar carcinoma xenograft in BALB/c nude mice. Head Neck. 2016;38(3):417–25.
Sturegård E, Johansson H, Ekström J, Hansson BG, Johnsson A, Gustafsson E, Dillner J, Forslund O. Human papillomavirus typing in reporting of condyloma. Sex Transm Dis. 2013;40(2):123–9.
Schmitt M, Bravo IG, Snijders PJ, Gissmann L, Pawlita M, Waterboer T. Bead-based multiplex genotyping of human papillomaviruses. J Clin Microbiol. 2006;44(2):504–12.
Soderlund-Strand A, Carlson J, Dillner J. Modified general primer PCR system for sensitive detection of multiple types of oncogenic human papillomavirus. J Clin Microbiol. 2009;47(3):541–6.
Peitsaro P, Hietanen S, Johansson B, Lakkala T, Syrjanen S. Single copy heterozygote integration of HPV 33 in chromosomal band 5p14 is found in an epithelial cell clone with selective growth advantage. Carcinogenesis. 2002;23(6):1057–64.
Sambrook JE, Fritsch F, Maniatis T: Molecular cloning: a laboratory manual., 2nd edn. Cold Spring Harbor, N.Y.: Cold Spring Harbor Laboratory Press; 1989.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98. Nucl Acids Symp Ser. 1999;41:95–8.
Welinder C, Ekblad L. Coomassie staining as loading control in Western blot analysis. J Proteome Res. 2011;10(3):1416–9.
Schaffer LG, Slovak ML, L.J. C (eds.): ISCN 2009: An international system for human cytogenetic nomenclature. Basel: S. Karger; 2009.
DSMZ Online STR Analysis [https://www.dsmz.de/services/services-human-and-animal-cell-lines/online-str-analysis.html].
Oshima G, Wennerberg J, Yamatodani T, Kjellen E, Mineta H, Johnsson A, Ekblad L. Autocrine epidermal growth factor receptor ligand production and cetuximab response in head and neck squamous cell carcinoma cell lines. J Cancer Res Clin Oncol. 2012;138(3):491–9.
Riaz N, Morris LG, Lee W, Chan TA. Unraveling the molecular genetics of head and neck cancer through genome-wide approaches. Genes Dis. 2014;1(1):75–86.
Hay A, Ganly I. Targeted therapy in oropharyngeal squamous cell carcinoma: the implications of HPV for therapy. Rare Cancers Ther. 2015;3:89–117.
Kalu NN, Mazumdar T, Peng S, Shen L, Sambandam V, Rao X, Xi Y, Li L, Qi Y, Gleber-Netto FO, et al. Genomic characterization of human papillomavirus-positive and -negative human squamous cell cancer cell lines. Oncotarget. 2017;8(49):86369–83.
Syrjanen S. HPV infections and tonsillar carcinoma. J Clin Pathol. 2004;57(5):449–55.
Doeberitz M, Vinokurova S. Host factors in HPV-related carcinogenesis: cellular mechanisms controlling HPV infections. Arch Med Res. 2009;40(6):435–42.
Gillison ML, Chaturvedi AK, Anderson WF, Fakhry C. Epidemiology of human papillomavirus-positive head and neck squamous cell carcinoma. J Clin Oncol. 2015;33(29):3235–42.
Storey A, Thomas M, Kalita A, Harwood C, Gardiol D, Mantovani F, Breuer J, Leigh IM, Matlashewski G, Banks L. Role of a p53 polymorphism in the development of human papillomavirus-associated cancer. Nature. 1998;393(6682):229–34.
Pinheiro UB, Fraga CA, Mendes DC, Farias LC, Cardoso CM, Silveira CM, D'Angelo MF, Jones KM, Santos SH, de Paula AM, et al. Fuzzy clustering demonstrates that codon 72 SNP rs1042522 of TP53 gene associated with HNSCC but not with prognoses. Tumour Biol. 2015;36(12):9259–65.
Habbous S, Pang V, Eng L, Xu W, Kurtz G, Liu FF, Mackay H, Amir E, Liu G. p53 Arg72Pro polymorphism, HPV status and initiation, progression, and development of cervical cancer: a systematic review and meta-analysis. Clin Cancer Res. 2012;18(23):6407–15.
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(6):1129–36.
Johansson C, Schwartz S. Regulation of human papillomavirus gene expression by splicing and polyadenylation. Nat Rev Microbiol. 2013;11(4):239–51.
Schwartz S. Papillomavirus transcripts and posttranscriptional regulation. Virology. 2013;445(1–2):187–96.
Lansford CD, Grenman R, Bier H, Somers KD, Kim SY, Whiteside TL, Clayman GL, Welkoborsky H-J, Carey TE: Head and neck cancers. In: Human Cell Culture Volume 2, edn. Edited by Masters JRW, Palsson B: Kluwer Academic Publishers; 1999: 185–255.
Akagi K, Li J, Broutian TR, Padilla-Nash H, Xiao W, Jiang B, Rocco JW, Teknos TN, Kumar B, Wangsa D, et al. Genome-wide analysis of HPV integration in human cancers reveals recurrent, focal genomic instability. Genome Res. 2014;24(2):185–99.
Busch CJ, Becker B, Kriegs M, Gatzemeier F, Kruger K, Mockelmann N, Fritz G, Petersen C, Knecht R, Rothkamm K, et al. Similar cisplatin sensitivity of HPV-positive and -negative HNSCC cell lines. Oncotarget. 2016;7(24):35832–42.
Ziemann F, Arenz A, Preising S, Wittekindt C, Klussmann JP, Engenhart-Cabillic R, Wittig A. Increased sensitivity of HPV-positive head and neck cancer cell lines to x-irradiation +/− cisplatin due to decreased expression of E6 and E7 oncoproteins and enhanced apoptosis. Am J Cancer Res. 2015;5(3):1017–31.
We thank Catarina Blennow, Camilla Ekenstierna, and Katarzyna Kawka for excellent technical assistance.
This work was supported by grants from the Swedish Research Council-Medicine, the Swedish Cancer Society, the Laryngology Fund, the King Gustaf V Jubilee Fund, Governmental funding of clinical research within the National Health System, Region of Scania R&D funding, the Foundations of the University Hospital of Lund, and the Berta Kamprad Foundation for Investigation and Control of Cancer Diseases. Funding sources had no influence on the study design or the interpretation of results.
Availability of data and materials
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study. Any raw data not included in the article is available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by Lund University regional ethical review board (LU 376–01), and the samples were collected, after a written informed consent was obtained, at Skåne University Hospital, Lund, Sweden.
Experiments on mice were approved by the local ethics committee (488/2017).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Karyogram from LU-HNSCC-26. (PDF 73 kb)
Partial sequencing results for TP53. The LU-HNSSC-26 is homozygous for the single nucleotide polymorphism chr17:7579472G > C (rs1042522). Partial sequencing results for TP53 mapped against GRCh37 (top) and Sanger sequencing electropherogram (bottom) with the rs1042522 indicated by arrows. (PDF 172 kb)
About this article
Cite this article
Forslund, O., Sugiyama, N., Wu, C. et al. A novel human in vitro papillomavirus type 16 positive tonsil cancer cell line with high sensitivity to radiation and cisplatin. BMC Cancer 19, 265 (2019). https://doi.org/10.1186/s12885-019-5469-8