- Research article
- Open Access
Establishment and characterization of patient-derived xenografts as paraclinical models for head and neck cancer
BMC Cancer volume 20, Article number: 316 (2020)
We investigated whether head and neck squamous cell carcinoma (HNSCC) patient-derived xenografts (PDXs) reaffirm patient responses to anti-cancer therapeutics.
Tumors from HNSCC patients were transplanted into immunodeficient mice and propagated via subsequent implantation. We evaluated established PDXs by histology, genomic profiling, and in vivo anti-cancer efficacy testing to confirm them as the authentic in vivo platform.
From 62 HNSCCs, 15 (24%) PDXs were established. The primary cancer types were tongue (8), oropharynx (3), hypopharynx (1), ethmoid sinus cancer (1), supraglottic cancer (1), and parotid gland (1); six PDXs (40%) were established from biopsy specimens from advanced HNSCC. PDXs mostly retained donor characteristics and remained stable across passages. PIK3CA (H1047R), HRAS (G12D), and TP53 mutations (H193R, I195T, R248W, R273H, E298X) and EGFR, CCND1, MYC, and PIK3CA amplifications were identified. Using the acquisition method, biopsy showed a significantly higher engraftment rate when compared with that of surgical resection (100% [6/6] vs. 16.1% [9/56], P < 0.001). Specimens obtained from metastatic sites showed a significantly higher engraftment rate than did those from primary sites (100% [9/9] vs. 11.3% [6/53], P < 0.001). Three PDX models from HPV-positive tumors were established, as compared to 12 from HPV-negative (15.8% [3/19] and 27.9% [12/43] respectively, P = 0.311), suggesting that HPV positivity tends to show a low engraftment rate. Drug responses in PDX recapitulated the clinical responses of the matching patients with pan-HER inhibitors and pan-PI3K inhibitor.
Genetically and clinically annotated HNSCC PDXs could be useful preclinical tools for evaluating biomarkers, therapeutic targets, and new drug discovery.
Head and neck squamous cell carcinoma (HNSCC) is the sixth common malignancy worldwide . Most patients experience local/regional recurrences, and systemic metastases occurred in about 20% of these patients after front-line treatment. Therapeutic options for recurrent or metastatic HNSCC patients include platinum-based combined or single cytotoxic chemotherapy, targeted therapy, and immunotherapy. Despite aggressive treatment regimens, survival outcomes in patients with metastatic disease remains dismal, showing median overall survival of less than one year . Recently, personalized therapy targeting specific genetic alterations has been introduced in HNSCC, and clinical trials are underway .
Although preclinical studies presented favorable responses to numerous novel anticancer agents, actual clinical trials in patients have not been as successful, with only 5% of novel anti-cancer agents showing efficacy. These disappointing results emphasize the necessity for novel therapeutic strategies to improve prediction of efficacy on the preclinical stage [3, 4]. Increasing attention has been engrossed on the development and characterization of patient-derived tumor xenograft (PDX) . PDX models are established by direct engraftment of the patient’s tumor into immunodeficient mice. These models recapitulate the principal histologic and genetic properties of corresponding patients’ tumor, and retain more tumor heterogeneity than other preclinical models [5, 6]. Therefore, we can anticipate patient’s clinical outcome by referring to PDXs, and thus are useful for clinical translational research, drug screening, and biomarker discovery and validation. PDXs are in vivo models to test clinically guided hypotheses, capable of evaluating drug activity and novel drug combined strategies, and also to elucidate their predictive biomarkers [7, 8].
p63 is a p53-related DNA-binding protein that helps regulate differentiation and proliferation in epithelial progenitor cells . The expression of p63 is increased in squamous cell cancer of various cancer type including the lung, esophagus, and head and neck but not expressed in adenocarcinoma. In this study, we conducted immunohistochemical staining panel including p16, p63, Ki67 as surrogate marker for diagnosis of squamous cell carcinoma.
Several groups have established PDX models of HNSCC, characterized their genomic fidelity and therapeutic response profiles [10,11,12,13,14,15]. Herein, we established PDXs using surgically resected or biopsied tumor tissues from HNSCC patients and characterized both histological and genomic fidelities of the established PDXs. We aim to evaluate the clinical factors affecting the establishment of HNSCC PDX models and whether established PDX models can reflect actual patients. Established PDX models in this study faithfully replicated the histologic, genomic, and responses with novel agents (pan-HER inhibitor and pan PI3K inhibitor) observed in the matched HNSCC patients. Our HNSCC PDX repository, which is clinically and genomically annotated in a comprehensive manner, may serve as a genuine platform for testing novel anti-cancer drugs.
The current research work was approved by the Severance Institutional Review Board (#4–2013-0526), and patients’ tumor tissues were from the Yonsei University Severance Hospital. All patients agreed with written informed consent. We collected tumors and its paired peripheral blood species to establish PDX and further conduct genetic analysis.
Establishment of PDX models
We established the HNSCC PDX model using a similar method as that used to generate lung cancer PDX . PDX models were generated using 6- to 8-week-old female nude (nu/nu) and severe combined immunodeficient (NOG, NOD/Shi-scid/IL-2Rγnull) mice (OrientBio, Seoul, Korea). We carefully removed the necrotic lesion from biopsied specimens with a surgical blade and subcutaneously implanted small pieces (3 × 3 × 3-mm3) obtained from each patient into 1–2 mice. When a tumor became a diameter of approximately 1.5 cm, it was dissected and chopped into small specimens (3 × 3× 3-mm3) and engrafted into other mice set. We defined the patient-originated specimen as the F0 generation, whereas subsequent generations were numbered consecutively by the number of re-implantations (e.g., F1, F2, or F3) . We engrafted patient-derived tumor (F0) into the NOG mice and after F1 generation in nude mice. We usually expanded F3 for drug responses of in vivo PDX models. We assigned Yonsei Human In Mouse (YHIM) identifiers in tumors and related PDXs corresponding to the original patient-derived tumors. After completion of experiments, we sacrificed mice by inhalation of anesthetics with CO2. All mice models were maintained in the specific pathogen-free facility of the Avison BioMedical Research Center (ABMRC) Animal Research Center at Yonsei University College of Medicine. We performed all procedures according the Animal Research Committee’s Guidelines at Yonsei University College of Medicine.
Tumor fragments from original primary tumors and PDX models were sliced into 3 × 3× 3-mm3 fragments using a surgical blade in a laminar flow cabinet. For storage, three fragments were placed in RPMI 1640 medium (supplemented with streptomycin, penicillin, gentamicin, and dimethyl sulfoxide), flash frozen, and kept at − 80 °C. The stored fragments were used for genetic profiling and/or conserved in RNAlater RNA stabilization reagent (Qiagen). Moreover, we fixed fragments in 10% neutral-buffered formalin and made paraffin embedded tissue samples for further additional pathological analysis .
Tumor growth measurement
Tumor size was measured once or twice every week in two dimensions with a caliper. Tumor volume was calculated using following formula: (width2 × length)/2. When the tumor was measured more than 1 cm3 or the mice reached an endpoint written in the Dutch Code of Practice for animal experiments in cancer research , the tumor was harvested and stored in media for either preservation or propagation into a following generation. The time of latency was calculated as the duration of time from engraftment until when tumor became measurable (about 70 mm3) .
We fixed harvested tissues from all PDX models in 10% buffered formalin in 30 min of resection. After 24 h of tissue fixation, routine procedures were followed for further tissue processing. Hematoxylin and eosin- and immunohistochemistry (IHC) staining with p63 antibody (Abcam: ab735, clone A4A), Ki67 antibody (Cell Signaling Technology, clone D2H10) stained slides were prepared, and the histopathological characteristics of each tumor including differentiation and tumor necrosis were reviewed pathologist (S.O.Y) in our institution. Moreover, p16 IHC (CINtec Histology, clone E6H4) was performed using a Ventana Bench Mark XT Autostainer to confirm the matched PDX F2 tissue originated from HPV positive patients (YHIM-3001, 3007, 3011). Regarding HPV positivity, we used conventionally accepted criteria and positivity was defined as the presence of strong and diffuse nuclear and cytoplasmic staining in more than 70% of tumor cells .
Targeted deep sequencing
We performed targeted deep sequencing on nine PDX and primary tumor paired specimens.. The SureSelect custom capture probe for the targeted deep sequencing was designed to detect 244 genes which are based on reporting of cancer related journals and literature. The target coverage of the captured region was 1000x, and to accomplish this, we used Illumina HiSeq 2500 platform as paired-end reads. We analyzed the targeted deep sequencing data based on the previous our report . In the general case, PDTX sequences of raw FASTQ were mixed human and mouse genome. For that reason, it is necessary to remove mouse reads from FASTQ. To perform this process, the human(hg19) and mouse(mm10) reference genome were integrated into one reference genome. Then Burrows-Wheeler Aligner  was run to identify the best similar read against the human genome refer to the primary alignments. We mapped reads to the most similar regions in both human and mouse genomes and obtained best-matched reads of the human genome . The Genome Analysis Tool Kit ver. 3.5 was used to conduct mark duplication, local realignment, and base quality score recalibration for the obtained reads. Germline mutations were identified using a GATK Haplotype caller. We filtered out variants with poor quality through GATK VariantFiltration and discarded according to the followings: QD < 2.0, DP < 20, FS > 60.0, MQ < 40.0, MQRankSum<− 12.5, ReadPosRankSum<− 8.0, and mutant allele frequency (MAF) < 0.2. Somatic mutations were called by MuTect v.1.1.7  using a default criteria, and databases including CIVIC  and DoCM were used to annotate these mutations.
Oncomine Cancer panel analysis
We used the Oncomine Cancer Panel to identify additional gene fusions and to cross-validate the genetic variations with the PDTX targeted deep sequencing data as described in our previous reports [8, 24]. Fourteen HNSCC and PDX samples (YHIM-3001, YHIM-3014) were sequenced but YHIM-3015 was not because of sample insufficiency. For amplified libraries, the Ion PGM Template kit on a OneTouch v2 was performed using the multiplexed DNA template according to the manufacturer’s protocol. The Torrent Suite Server v4.4 was utilized for the alignment process using the human genome reference(hg19). The variant calling and annotation were performed by self-owned variant analysis system (Ion Reporter Workflow v5.0, Oncomine Panel Annotations set v1.1, Oncomine Panel Baseline v2.0, Oncomine Variant annotator plugin v2.0) which is contained in Torrent Suite Server v4.4. To figure out the pattern of the variants over all sequenced data, we used Oncoprint provided by ComplexHeatmap package .
In vivo drug treatment
As described in previous reports , mice were set apart and randomized into each group (six to seven mice per group), when tumor volumes became 200–250 mm3. We administered vehicle (5 mM citrate buffer) to the control group. We administered afatinib (15 mg/kg) and BKM120 (35 mg/kg) methotrexate (30 mg/kg) via oral gavage or intraperitoneal injection mixed with 5 mM citrate buffer in deionized water into each model. We measured dimensions of tumor twice weekly using a digital caliper and calculated tumor volume using the following formula: tumor volume = [length×width2]/2. We defined percentage tumor growth inhibition as [%TGI = 1 − (change of tumor volume in treatment group/change of tumor volume in control group) × 100] to evaluate antitumor efficacy as primary experimental outcomes.
In vivo drug test, the sample size was six to seven in each group. The data are presented as the mean ± SD. Statistical differences between vehicle and afatinib treatment group were determined using ANOVA test. A p-value< 0.05 was considered statistically significant.
Establishment of the HNSCC PDX models
We collected 62 tumor specimens from patients and engrafted these specimens into mice between May 2013 and March 2017. Among these, 15 HNSCC PDXs were successfully generated, which include, 8 tongue, 3 oropharyngeal, 1 hypopharyngeal cancers, 1 ethmoid sinus cancer, 1 supraglottic cancer, and 1 parotid gland. Baseline characteristics for the established 15 PDX models are summarized in Supplementary Table 1. Patients were median 55 years old (range: 36–72), and tumor stages ranged from II to IV. Six PDXs (40%) were successfully generated using biopsied tumor tissues from advanced-stage HNSCC patients. Tissue samples were obtained from surgical resection (60%, 9/15) and core biopsy (40%, 6/15). The donor tissue for the 15 established PDX models were derived from primary (6/15, 40%) and metastatic sites (9/15; 60%).
All PDXs were generated through 3–5 sequential transfers with varying latency times from 10 to 41 days. We established a HNSCC PDX repository, consisting of the original patient tissues and the PDX tissues from each transfer, including formalin-fixed paraffin-embedded (FFPE) blocks, viable cryopreserved tissue specimens, and RNAlater-treated and/or fresh-frozen tissues .
PDX models preserve histological stability compared to primary cancers
In order to determine the biological stability of the PDX tumors, the histological properties from original patient tumors (F0) and matching second-generation PDX tumors (F2) were compared. F0 and F2 generation tumors had close histopathological characteristics including differentiation and tumor necrosis, as well as similar p63 expression (Fig. 1 and Supplementary Figure 1). Detailed pathologic information between F0 and F2 were summarized in Supplementary Table 2. Additionally, we performed the Ki67 IHC staining as surrogate marker of proliferation (Supplementary Figure 2). Taken together, the F2 tumors preserved comparable histological characteristics to those of F0 tumor specimens, implicating morphological fidelity with their matched patient tumors. Moreover, we validate HPV positivity in corresponding PDX tumor tissues by IHC staining of p16 protein in tumor cells as shown in Supplementary Figure 3.
Engraftment and tumor growth rates of HNSCC PDX models
We further analyzed successful engraftment according to tumor acquisition method, tumor acquisition site, stage at tumor acquisition, and human papillomavirus (HPV) infection (Table 1). We defined engraftment success as completion of the transfer to the F3, third generation. We achieved a successful overall engraftment rate of 24% (15/62).
When analyzing engraftment rates using the acquisition method, biopsy showed a significantly higher engraftment rate than that of surgical resection (100% [6/6] vs. 16.1% [9/56], P = 0.010). Establishment using specimens from metastatic sites demonstrated a significantly higher engraftment rate compared to those from primary sites (100% [9/9] vs. 11.3% [6/53], P = 0.002). A higher stage (stage IV) at sample acquisition tended to show a higher engraftment rate than that of stage I–III (32.4% [11/34] and 14.3% [4/28], respectively, P = 0.106). Of the 62 tumors classified, 19 were HPV-positive, based on p16 protein expression in the major part of cancer cells. We successfully generated three PDX models from HPV-positive tumors, as compared to 12 from HPV-negative tumors (15.8% 3/19 and 27.9% 1243, respectively, P = 0.311), suggesting that HPV positivity tends to show a low engraftment rate. Engraftment rates were comparable across primary cancer types.
If we performed the multivariate analysis with several variables including “Method of tumor acquisition”, “Site of sampling”, and “Stage at tumor acquisition”, “Site of sampling” was observed the significant factor related with the success of PDX establishment (OR = 0.016, 95%CI = < 0.001–0.355, P = 0.009) (Supplementary Table 3).
Figure 2 illustrates the tumor growth curves in the established HNSCC PDXs for three generations (F1–F3). F3 exhibited relatively faster growth than that of F1 and F2. Moreover, the growth rate stabilized after the third generation. The median latency time was 27 days, but it varied between tumor acquisition method, tumor acquisition site, stage at tumor acquisition, and HPV infection; moreover, samples acquired at stage IV showed a significantly longer latent time than did other sources of variation (Supplementary Table 4).
Genomic fidelity of the PDX tumors
We conducted targeted next-generation sequencing to figure out single-nucleotide variants and small insertions/deletions using nine pairs of the original patient (F0) and the PDX F2 tumor. The goal of current study was to discover identical mutations in F0 and F2, which would indicate that F0 and F2 are identical if most of the germline mutations coincided. The Jaccard similarity score was used to measure the similarity; the score of most samples was > 82%. The highest scoring models were YHIM-3002 and -3009, with respective Jaccard similarity scores of 93.6 and 95.2% (Fig. 3a). The MAF values for common mutations between F0 and F2 demonstrated overall concordance (r2 = 0.93 and 0.95 for YHIM-3002 and -3009, respectively), implicating that most germline mutations in the F2 sample corresponded with those in the F0 sample (Fig. 3b). HRAS G12D and PIK3CA H1047R mutations were observed simultaneously in F0 and F2 (YHIM-3003, − 3013). F0 and F2 showed similar copy numbers for CCND1 amplification.
Genetic alteration based on Oncomine Cancer panel
We next performed a genomic overview of the most recurring somatic mutations of interest in 14 PDX tumor specimens using the Oncomine Cancer Panel (Fig. 3c, Supplementary Table 5). After filtering the predefined Oncomine variants, we found an average of 0.64 relevant somatic point mutations and 1.1 high-level CNAs per specimen. An integrative heat map showed that the prioritized alterations across the YHIM cohort and the copy number profiles for all samples (Fig. 3c). All Oncomine-derived relevant genetic alterations in the RNA and DNA components of the 14 PDX specimens are shown in the heat map. PIK3CA, HRAS, and TP53 mutations and EGFR, CCND1, MYC, and PIK3CA amplifications were identified in established PDXs. Gene amplifications included CCND1 (36% [5/14]), EGFR (29% [4/14]), MYC (29% [4/14]), and PIK3CA (14% [2/14]). Regarding somatic mutations, TP53 mutation (36% [5/14], H193R, I195T, R248W, R273H, E298X), PIK3CA mutation (7% [1/14], 1047R), and HRAS mutation (7% [1/14], G12D) were observed.
PDXs faithfully recapitulated the anti-cancer efficacy of their matched patients
We then assessed whether the generated PDXs can recapitulate the anti-cancer efficacy of matched patients and thus work for a genuine platform for novel anti-cancer drug efficacy testing. We had two PDX models (YHIM-3006, and YHIM-3011) which were involved in co-clinical trial with afatinib, a pan-HER inhibitor. YHIM-3006 was established using surgical resection from tongue cancer patients. In YHIM-3006, afatinib induce growth delay (TGI = 44.9%; P < 0.001) (Fig. 4a), which was similar with the modest clinical outcome in the actual corresponding patient. This patient showed stable disease as best response for 5 months, but the tumor was abruptly increased after 5 months of afatinib treatment. Next, YHIM-3011 was generated using small biopsied tissue from a patient with recurrent oropharyngeal cancer after surgery and radiation therapy and platinum-based chemotherapy. Corresponding patient of YHIM-3011 represented marked tumor shrinkage showing partial response for 8 months, which was recapitulated with the PDX model which showed remarkable tumor shrinkage with afatinib treatment (TGI = 105.7%; P < 0.001) (Fig. 4b). Moreover, these two patients were previously treated with MTX as anti-cancer therapy and showed strong resistance to MTX (Supplementary Figure 4). These findings were correlated with those of MTX response in the matched PDX models. Moreover, YHIM-3004, PDX model were treated with BKM120 as pan-PI3K inhibitor as co-clinical trial. This model showed strong resistance to BKM120 monotherapy, which is recapitulated with response of corresponding patient. This patient showed rapid progression after 1 month to BKM120 monotherapy showing abrupt increase of multiple metastatic LN (Fig. 4 c). The weights of mice in each group were not significantly reduced (data not shown).
We need to interpret cautiously due to small sample size, PDXs seem to recapitulate the drug responses of the corresponding patients and may serve as a clinically relevant platform for novel drug testing.
The present study successfully established 15 HNSCC PDXs by directly engrafting tumor tissue samples from patients into immunodeficient mice. These PDXs maintained features of the donor patients’ cancers as determined by histology as well as genetic aberrations. These established PDXs may be useful as preclinical platforms to test newly developed anti-cancer drug efficacy. Of the 15 established PDXs, most were derived from tongue cancer; six were made using a small amount of biopsied tissue originated from patients with advanced HNSCC. The established PDXs harbored some clinically implicated druggable targets, including copy number amplifications (CCND1, PIK3CA, EGFR) and somatic mutations (PIK3CA H1047R, HRAS G12D).
Our data suggest that HNSCC PDX establishment has a 24% (15/62) success rate. We comprehensively evaluated which factors, including tumor acquisition method, tumor acquisition site, stage, and HPV infection, affected the engraftment rate of PDX models. PDXs from a biopsy sample, from a metastatic site, or of a higher stage showed a significantly higher engraftment rate. Additionally, PDXs from HPV-negative tumors tended to show a higher engraftment rate compared with those from HPV-positive tumors. These results of our study are consistent with prior reports .
We conducted targeted next-generation sequencing to figure out single-nucleotide variants and small insertions/deletions using nine pairs of the F0 and the PDX F2 tumor. We observed that genetic mutations or amplifications of PDX matched to those of original tumors. However, it could be possible that tumors with stable mutation/amplification were selected during the course of PDX establishment. To find out the identification of PDX models and their matched patients, we exam histology and IHC staining using p63 and p16 in this study. Other reports also evaluate the histology and IHC staining including p16 and CD31, to identify whether established PDX models retain the histopathologic characteristics of original patients [26, 27].
The essential strength of the PDX platform is that it is generated from unselected and uncultured patient samples and has the homogeneous genetic background as its matched patient tumor . In our study, 83–95% of genetic aberrations and the morphology from the corresponding donor patient cancers were maintained in the established PDXs. This suggests that PDXs maintain relatively stable genomes without significant accumulation of DNA structure rearrangements, but do have some of enrichment for PDX-specific single nucleotide variants . Therefore, we conclude that anti-cancer efficacy research works performed using PDXs will be similar to what actually happens in patient tumors. Previous studies have also shown great concordance of response between donor patient tumors and PDXs [30, 31]. The concordance between patient and xenograft responses despite the absence of the immune system is encouraging and consistent with previous literature on targeted therapy [10, 11]. In this study, drug responses in PDX recapitulated the clinical responses of the corresponding patients with afatinib. Since clinical studies are examining the combination of Afatinib with RT, future studies examining the combination in their PDX panel would be relevant. In this regard, previous studies have shown concordance of radiation response of HPV+ and HPV- PDX models of HNSCC . In this study, the responses of PDXs to pan-HER inhibitor mimicked the clinical responses of the matched patients, indicating that established PDXs represent a clinically relevant platform for predicting drug response.
HPV infection has become as a major risk factor, specifically in the oropharynx, accounting for more than 20% of all HNSCCs . While the overall HNSCC incidence continues to decline, there has been an observable and highly increase in the incidence of HPV-associated HNSCC, predominantly among young patients [33,34,35]. Of the 62 tumors classified, 19 were HPV-positive. We successfully established three PDX models (15.8% [3/19] engraftment rate) from HPV-positive tumors, suggesting that HPV positivity decreases the engraftment rate (HPV-negative tumors: 27.9%, [12/43]). Our HPV-positive PDXs could help identify novel and more effective therapeutic strategies and transition patients with these tumors to personalized therapies. Interestingly, YHIM-3001 and YHIM-3011 models were originated from HPV positive and never smoker tongue and oropharyngeal cancer patients. Additionally, these two models have I195T and H193R mutation in TP53 gene, which is unusual. These models (YHIM-3001 and 3007) were originated from patients with recurrence who were previously treated with chemotherapy. Therefore, the presence of TP53 mutation could be explained to be chemotherapeutic agents induced mutations.
The Cancer Genome Atlas (TCGA) with comprehensive genetic profiling of 279 HNSCC specimens demonstrated that EGFR, CCND1, PIK3CA, and FGFR1 amplifications and TP53, CDKN2A, and HRAS mutations are common genetic alterations in HNSCC . These genetic alterations are potential therapeutic and druggable targets in HNSCC patients. EGFR as a cell surface receptor member of the ErbB family has been extensively studied in HNSCC. Genetic alterations in EGFR pathway are potential therapeutic and druggable targets in HNSCC. Currently, EGFR targeting strategies include EGFR antibody (i.e., cetuximab) and EGFR tyrosine kinase inhibitors including erlotinib, gefitinib, or afatinib . Although several prospective clinical trials were conducted, only a minor portion of R/M HNSCC responded to EGFR targeting agents; predictive molecular biomarkers have not yet been identified . Furthermore, genetic alterations in the PI3K signaling pathway are commonly observed in HNSCC, and pan-PI3K inhibitor such as buparlisib with paclitaxel improved efficacy in R/M HNSCC patients compared with paclitaxel alone, suggesting that PI3K inhibition plays an important anti-cancer role. However, predictive markers to determine the clinical efficacy for PI3K inhibitors are unknown . To maximize the treatment outcome of targeted therapy in HNSCC patients, there is a need to identify molecular determinants for drug efficacy. In this study, clinically relevant therapeutic targets such as CCND1, PIK3CA, and EGFR amplifications and PIK3CA H1047R, HRAS G12D mutations were observed in several PDX models. These PDX models that retain the potential therapeutic target can help to guide efficient targeted therapy and overcome resistance mutations. Recently, several umbrella clinical trials, in which patients were allocated to targeted therapy according to genetic alterations based on next-generation sequencing, are underway for advanced HNSCC. The MOSCATO-01 clinical trial demonstrated that genetic analyses of advanced solid cancers patients improved treatment outcomes through matched targeted therapeutics . The ongoing NCI-MATCH trial and TRIUMPH trial (NCT03292250) are evaluating whether molecular biomarkers can predict response to target therapy in advanced cancer patients [24, 40]. This process might render greater success in identifying more potent treatment strategies and transitioning patients to personalized therapies.
Notably, the PDX models have several advantages in precision medicine in translational cancer research. Conventional cell lines are usually cultured in vitro for a long time, which leads to the development of extra genetic alterations that differ from the original cancer. Thus, it is difficult to predict clinical anti-cancer efficacy only by relying on the genomic information of cell lines (i.e., the Genomics of Drug Sensitivity in Cancer  and the Cancer Cell Line Encyclopedia  . Moreover, PDXs may aid in identifying drug–tumor relationships on the basis of inter- or intra-patient tumoral heterogeneity. Prior report showed that intertumoral angiogenic heterogeneity was observed in HNSCC PDX models and the therapeutic response anti-VEGF therapy was dependent on their angiogenic heterogeneity .
In addition, it is feasible to predict and overcome acquired resistance to novel anti-cancer drugs by investigating resistance mechanisms. Thus, investigators can synchronize preclinical experiments and clinical trials and systematize study design and data analysis protocols for PDX . Recent development of high-throughput technologies provided important molecular insights in genetically characterized PDX models.
On the other hand, there are some disadvantages of the PDX models . First, it routinely takes 4–5 months to generate and stabilize the first generation of PDX (F1), but only 2–5 weeks for subsequent passages. This long time lapse means that PDX platform can only be used as a research tool, especially for cancer patients with a fast disease progression and a short survival time [46,47,48]. Second, PDX platforms lack an intact immune system and may be of limited use for immunotherapeutic drug testing. These limitations could be solved by establishing humanized in vivo models, that were used to evaluate anti-cancer drug efficacy and elucidate predictive biomarkers and resistant mechanisms to immunotherapies.
In summary, we have developed and characterized an HNSCC PDX models that is useful for screening novel targeted drugs or combination drugs and identifying resistance mechanisms. Our data illustrate that PDXs may be a robust preclinical platform for analyzing biomarkers, therapeutic targets, and novel anti-cancer drug development.
Availability of data and materials
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Association of Assessment and Accreditation of Laboratory Animal Care
Head and neck squamous cell carcinoma
Patient-derived tumor xenograft
Schoenfeld JD. Immunity in head and neck cancer. Cancer Immunol Res. 2015;3(1):12–7.
Parsel SM, Grandis JR, Thomas SM. Nucleic acid targeting: towards personalized therapy for head and neck cancer. Oncogene. 2016;35(25):3217–26.
Abaan OD, Polley EC, Davis SR, Zhu YJ, Bilke S, Walker RL, Pineda M, Gindin Y, Jiang Y, Reinhold WC, et al. The exomes of the NCI-60 panel: a genomic resource for cancer biology and systems pharmacology. Cancer Res. 2013;73(14):4372–82.
Shoemaker RH. The NCI60 human tumour cell line anticancer drug screen. Nat Rev Cancer. 2006;6(10):813–23.
Fritsch M, Schmidt N, Groticke I, Frisk AL, Keator CS, Koch M, Slayden OD. Application of a patient derived Xenograft model for predicative study of uterine fibroid disease. PLoS One. 2015;10(11):e0142429.
Griffith M, Rovira J, Torres R, Calzada J, Victoria C, Caceres L. knowledge, attitudes and practices regarding malaria in the indigenous Guna population of the Madungandi region, Panama, 2012. Biomedica. 2015;35(4):480–95.
Williams JA. Using PDX for Preclinical Cancer Drug Discovery: The Evolving Field. J Clin Med. 2018;7(3):41.
Kang HN, Choi JW, Shim HS, Kim J, Kim DJ, Lee CY, Hong MH, Park SY, Park AY, Shin EJ, et al. Establishment of a platform of non-small-cell lung cancer patient-derived xenografts with clinical and genomic annotation. Lung Cancer (Amsterdam, Netherlands). 2018;124:168–78.
Westfall MD, Pietenpol JA. p63: molecular complexity in development and cancer. Carcinogenesis. 2004;25(6):857–64.
Peng S, Creighton CJ, Zhang Y, Sen B, Mazumdar T, Myers JN, Lai SY, Woolfson A, Lorenzi MV, Bell D, et al. Tumor grafts derived from patients with head and neck squamous carcinoma authentically maintain the molecular and histologic characteristics of human cancers. J Transl Med. 2013;11:198.
Klinghammer K, Raguse JD, Plath T, Albers AE, Joehrens K, Zakarneh A, Brzezicha B, Wulf-Goldenberg A, Keilholz U, Hoffmann J, et al. A comprehensively characterized large panel of head and neck cancer patient-derived xenografts identifies the mTOR inhibitor everolimus as potential new treatment option. Int J Cancer. 2015;136(12):2940–8.
Karamboulas C, Bruce JP, Hope AJ, Meens J, Huang SH, Erdmann N, Hyatt E, Pereira K, Goldstein DP, Weinreb I, et al. Patient-Derived Xenografts for Prognostication and Personalized Treatment for Head and Neck Squamous Cell Carcinoma. Cell Rep. 2018;25(5):1318–1331.e1314.
Facompre ND, Sahu V, Montone KT, Harmeyer KM, Nakagawa H, Rustgi AK, Weinstein GS, Gimotty PA, Basu D. Barriers to generating PDX models of HPV-related head and neck cancer. Laryngoscope. 2017;127(12):2777–83.
Kimple RJ, Harari PM, Torres AD, Yang RZ, Soriano BJ, Yu M, Armstrong EA, Blitzer GC, Smith MA, Lorenz LD, et al. Development and characterization of HPV-positive and HPV-negative head and neck squamous cell carcinoma tumorgrafts. Clin Cancer Res. 2013;19(4):855–64.
Campbell KM, Lin T, Zolkind P, Barnell EK, Skidmore ZL, Winkler AE, Law JH, Mardis ER, Wartman LD, Adkins DR, et al. Oral cavity squamous cell carcinoma Xenografts retain complex genotypes and Intertumor molecular heterogeneity. Cell Rep. 2018;24(8):2167–78.
Sun S, Liu S, Duan SZ, Zhang L, Zhou H, Hu Y, Zhou X, Shi C, Zhou R, Zhang Z. Targeting the c-met/FZD8 signaling axis eliminates patient-derived cancer stem-like cells in head and neck squamous carcinomas. Cancer Res. 2014;74(24):7546–59.
Alkema NG, Tomar T, Duiker EW, Jan Meersma G, Klip H, van der Zee AG, Wisman GB, de Jong S. Biobanking of patient and patient-derived xenograft ovarian tumour tissue: efficient preservation with low and high fetal calf serum based methods. Sci Rep. 2015;5:14495.
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.
Li H, Durbin R. Fast and accurate long-read alignment with burrows-Wheeler transform. Bioinformatics. 2010;26(5):589–95.
Tso KY, Lee SD, Lo KW, Yip KY. Are special read alignment strategies necessary and cost-effective when handling sequencing reads from patient-derived tumor xenografts? BMC Genomics. 2014;15:1172.
Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, Gabriel S, Meyerson M, Lander ES, Getz G. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31(3):213–9.
Griffith M, Spies NC, Krysiak K, McMichael JF, Coffman AC, Danos AM, Ainscough BJ, Ramirez CA, Rieke DT, Kujan L, et al. CIViC is a community knowledgebase for expert crowdsourcing the clinical interpretation of variants in cancer. Nat Genet. 2017;49(2):170–4.
Ainscough BJ, Griffith M, Coffman AC, Wagner AH, Kunisaki J, Choudhary MN, McMichael JF, Fulton RS, Wilson RK, Griffith OL, et al. DoCM: a database of curated mutations in cancer. Nat Methods. 2016;13(10):806–7.
Lim SM, Cho SH, Hwang IG, Choi JW, Chang H, Ahn MJ, Park KU, Kim JW, Ko YH, Ahn HK, et al. Investigating the feasibility of targeted next-generation sequencing to guide the treatment of head and neck squamous cell carcinoma. Cancer Res Treat. 2019;51(1):300–12.
Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–9.
Wennerberg J, Willen R, Biorklund A, Trope C. Histopathological characteristics predictive for growth of squamous cell carcinomas of the head and neck heterotransplanted to nude mice. Anticancer Res. 1986;6(5):1165–70.
Folaron M, Merzianu M, Duvvuri U, Ferris RL, Seshadri M. Profiling the Stromal and Vascular Heterogeneity in Patient-derived Xenograft Models of Head and Neck Cancer: Impact on Therapeutic Response. Cancers. 2019;11(7):951.
Pan CX, Zhang H, Tepper CG, Lin TY, Davis RR, Keck J, Ghosh PM, Gill P, Airhart S, Bult C, et al. Development and characterization of bladder Cancer patient-derived Xenografts for molecularly guided targeted therapy. PLoS One. 2015;10(8):e0134346.
Li S, Shen D, Shao J, Crowder R, Liu W, Prat A, He X, Liu S, Hoog J, Lu C, et al. Endocrine-therapy-resistant ESR1 variants revealed by genomic characterization of breast-cancer-derived xenografts. Cell Rep. 2013;4(6):1116–30.
Zhang X, Claerhout S, Prat A, Dobrolecki LE, Petrovic I, Lai Q, Landis MD, Wiechmann L, Schiff R, Giuliano M, et al. A renewable tissue resource of phenotypically stable, biologically and ethnically diverse, patient-derived human breast cancer xenograft models. Cancer Res. 2013;73(15):4885–97.
Marangoni E, Vincent-Salomon A, Auger N, Degeorges A, Assayag F, de Cremoux P, de Plater L, Guyader C, De Pinieux G, Judde JG, et al. A new model of patient tumor-derived breast cancer xenografts for preclinical assays. Clin Cancer Res. 2007;13(13):3989–98.
D'Souza G, Kreimer AR, Viscidi R, Pawlita M, Fakhry C, Koch WM, Westra WH, Gillison ML. Case-control study of human papillomavirus and oropharyngeal cancer. N Engl J Med. 2007;356(19):1944–56.
Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. 2008;26(4):612–9.
Ryerson AB, Peters ES, Coughlin SS, Chen VW, Gillison ML, Reichman ME, Wu XC, Chaturvedi AK, Kawaoka K. Burden of potentially human papillomavirus-associated cancers of the oropharynx and Oral cavity in the US, 1998-2003. Cancer. 2008;113(10):2901–9.
Shiboski CH, Schmidt BL, Jordan RCK. Tongue and tonsil carcinoma - increasing trends in the US population ages 20-44 years. Cancer. 2005;103(9):1843–9.
The Cancer Genome Atlas Network. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature. 2015;517(7536):576–82.
Kim HS, Kwon HJ, Jung I, Yun MR, Ahn MJ, Kang BW, Sun JM, Kim SB, Yoon DH, Park KU, et al. Phase II clinical and exploratory biomarker study of dacomitinib in patients with recurrent and/or metastatic squamous cell carcinoma of head and neck. Clin Cancer Res. 2015;21(3):544–52.
Soulieres D, Faivre S, Mesia R, Remenar E, Li SH, Karpenko A, Dechaphunkul A, Ochsenreither S, Kiss LA, Lin JC, et al. Buparlisib and paclitaxel in patients with platinum-pretreated recurrent or metastatic squamous cell carcinoma of the head and neck (BERIL-1): a randomised, double-blind, placebo-controlled phase 2 trial. Lancet Oncol. 2017;18(3):323–35.
Massard C, Michiels S, Ferte C, Le Deley MC, Lacroix L, Hollebecque A, Verlingue L, Ileana E, Rosellini S, Ammari S, et al. High-throughput genomics and clinical outcome in hard-to-treat advanced cancers: results of the MOSCATO 01 trial. Cancer Discov. 2017;7(6):586–95.
Mullard A. NCI-MATCH trial pushes cancer umbrella trial paradigm. Nat Rev Drug Discov. 2015;14(8):513–5.
McDermott U, Sharma SV, Dowell L, Greninger P, Montagut C, Lamb J, Archibald H, Raudales R, Tam A, Lee D, et al. Identification of genotype-correlated sensitivity to selective kinase inhibitors by using high-throughput tumor cell line profiling. Proc Natl Acad Sci U S A. 2007;104(50):19936–41.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehar J, Kryukov GV, Sonkin D, et al. The Cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483(7391):603–7.
Johnson JI, Decker S, Zaharevitz D, Rubinstein LV, Venditti J, Schepartz S, Kalyandrug S, Christian M, Arbuck S, Hollingshead M, et al. Relationships between drug activity in NCI preclinical in vitro and in vivo models and early clinical trials. Br J Cancer. 2001;84(10):1424–31.
Hasina R, Whipple ME, Martin LE, Kuo WP, Ohno-Machado L, Lingen MW. Angiogenic heterogeneity in head and neck squamous cell carcinoma: biological and therapeutic implications. Lab Invest. 2008;88(4):342–53.
Hidalgo M, Amant F, Biankin AV, Budinska E, Byrne AT, Caldas C, Clarke RB, de Jong S, Jonkers J, Maelandsmo GM, et al. Patient-derived Xenograft models: an emerging platform for translational Cancer research. Cancer Discovery. 2014;4(9):998–1013.
Tentler JJ, Tan AC, Weekes CD, Jimeno A, Leong S, Pitts TM, Arcaroli JJ, Messersmith WA, Eckhardt SG. Patient-derived tumour xenografts as models for oncology drug development. Nat Rev Clin Oncol. 2012;9(6):338–50.
Siolas D, Hannon GJ. Patient-derived tumor Xenografts: transforming clinical samples into mouse models. Cancer Res. 2013;73(17):5315–9.
Houghton JA, Taylor DM. Growth characteristics of human colorectal tumours during serial passage in immune-deprived mice. Br J Cancer. 1978;37(2):213–23.
We thank to the Yonsei Genomic Institute for supporting Oncomine and sequencing.
This work was supported by a grant from the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (HI15C1592 to BCC), by Basic Science Research Program through the NRF funded by the Ministry of Science and ICT (NRF-2016R1A2B3016282 to BCC), and by a Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (HI13C2162 to SP). The funding bodies had not involved in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.
All applicable international, national, and/or institutional guidelines for the care and use of animals were followed. All methods were approved by the AAALAC (Association of Assessment and Accreditation of Laboratory Animal Care). All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee. The study was approved by the appropriate Institutional Review Board (4–2013-0526), and patient tumor specimens were obtained from the Yonsei University Severance Hospital.
All patients to participate in this study provided written informed consent based on Institutional Review Board (4–2013-0526).
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.
: Supplementary Table 1. Baseline characteristics of patients (N = 15). Supplementary Table 2. Detailed pathologic information of F0 and F2 generation. Supplementary Table 3. Multivariate analysis to determine the association between covariates and patients- derived xenograft establishment. Supplementary Table 4. Latency time of 15 established PDX models. Supplementary Table 5. Detailed information of established head and neck cancer patient-derived xenografts
: Supplementary Figure 1. Comparison of histopathologic and immunohistochemistry findings between patients and their matched patient-derived xenografts (PDXs, F2 generation) of head and neck cancer squamous cell carcinoma (HNSCC). Left columns show hematoxylin and eosin-staining and right columns show immunohistochemical staining of p63, a marker of squamous cell carcinoma, from patients and the PDX of each model. Representative stained sections are shown (magnification: 200× in patient samples; scale bars = 100 μm).Supplementary Figure 2. P16 immunohistochemistry staining of tumor tissue between patient and PDX model. Supplementary Figure 3. Ki67 immunohistochemistry staining on tumor tissue of patient derived xenograft model. Supplementary Figure 4. Tumor growth rate to anti-cancer therapy in YHIM-3006 and 3011 models.
About this article
Cite this article
Kang, H.N., Kim, J., Park, A. et al. Establishment and characterization of patient-derived xenografts as paraclinical models for head and neck cancer. BMC Cancer 20, 316 (2020). https://doi.org/10.1186/s12885-020-06786-5
- Patient-derived xenograft
- Head and neck cancer
- Squamous cell cancer