Skip to main content

NSD1 mutations by HPV status in head and neck cancer: differences in survival and response to DNA-damaging agents



Compared to HPV-negative head and neck squamous cell carcinomas (HNSCCs), HPV-positive HNSCCs are associated with a favorable prognosis in part due to their improved treatment sensitivity. Inactivating mutations in NSD1 were shown to be a favorable prognostic biomarker in laryngeal cancers. Here, we characterize NSD1 mutations from the expanded The Cancer Genome Atlas (TCGA) HNSCC cohort (n = 522) and examine their prognostic implications based on HPV status of the tumor. We also begin to examine if NSD1 regulates response to platinum-based drugs and other DNA-damaging agents.


TCGA HNSCC samples were segregated by HPV and NSD1 mutations using cBioPortal and patient survival was determined. Pathogenicity of mutations was predicted using UMD-Predictor. NSD1-depleted cell lines were established by transfection with control or shRNAs against NSD1, followed by puromycin selection, and confirmed by qRT-PCR. Cell sensitivity to DNA damaging agents was assessed using short-term proliferation and long-term clonogenic survival assays.


Among 457 HPV(−) tumors, 13% contained alterations in the NSD1 gene. The majority (61.3%) of NSD1 gene alterations in HPV(−) specimens were truncating mutations within or before the enzymatic SET domain. The remaining alterations included homozygous gene deletions (6.7%), missense point mutations (30.7%) and inframe deletions (1.3%). UMD-Predictor categorized 18 of 23 missense point mutations as pathogenic. For HPV(+) HNSCC (n = 65), 6 NSD1 mutations, comprised of two truncating (33%) and 4 missense point (66%) mutations, were identified. Three of the 4 missense point mutations were predicted to be pathogenic or probably pathogenic by UMD-Predictor. Kaplan-Meier survival analysis determined significantly improved survival of HPV(−) HNSCC patients whose tumors harbored NSD1 gene alterations, as compared to patients with wild-type NSD1 tumors. Interestingly, the survival effect of NSD1 mutations observed in HPV-negative HNSCC was reversed in HPV(+) tumors. Proliferation and clonogenic survival of two HPV(−) cell lines stably expressing control or NSD1 shRNAs showed that NSD1-depleted cells were more sensitive to cisplatin and carboplatin, but not to other DNA damaging drugs.


Genetic alterations in NSD1 hold potential as novel prognostic biomarkers in HPV(−) head and neck cancers. NSD1 mutations in HPV(+) cancers may also play a prognostic role, although this effect must be examined in a larger cohort. NSD1 downregulation results in improved sensitivity to cisplatin and carboplatin, but not to other DNA-damaging agents, in epithelial cells. Increased sensitivity to platinum-based chemotherapy agents associated with NSD1 depletion may contribute to improved survival in HPV(−) HNSCCs. Further studies are needed to determine mechanisms through which NSD1 protects HPV(−) HNSCC cells from platinum-based therapy, as well as confirmation of NSD1 effect in HPV(+) HNSCC.


Head and neck squamous cell carcinomas (HNSCCs) comprise a diverse group of tumors, which are classified into various anatomical sites (e.g. oral cavity, oropharynx, hypopharynx, larynx, and nasopharynx), with unique epidemiology, anatomy, clinical behavior, and association with human papilloma virus (HPV) or Epstein-Barr Virus infection [1, 2]. HPV-negative (HPV(−)) head and neck cancers are strongly associated with tobacco and alcohol use and are most commonly found in the oral cavity and larynx [3]. Rates of HPV(−) cancers have declined in recent decades with the decreasing rates of smoking, but during this same period, rates of oropharyngeal cancers continued to increase [3,4,5,6,7]. The majority of oropharyngeal squamous cell carcinomas (OPSCCs) are associated with HPV and occur in younger patients with minimal or no tobacco exposure [8,9,10,11,12,13]. HPV-positive (HPV(+)) OPSCCs carry a favorable prognosis compared to HPV-negative cancers. Five-year survival rates for patients with HPV(+) OPSCC are 60–90%, versus survival rates of 20–25% among patients with HPV(−) OPSCCs (survival rates for HPV(−) tumors of all sites is 65%) [3]. The improved survival of patients with HPV(+) tumors can in part be attributed to their remarkable treatment sensitivity, as HPV(+) tumors have been shown to respond better to chemotherapy and radiation than HPV(−) tumors [13, 14].

For HPV(−) cancers, which are significantly more deadly than HPV(+) cancers, understanding their molecular diversity is critical to developing and selecting appropriate treatments to improve mortality. Comprehensive genomic profiling of head and neck squamous cell carcinomas by The Cancer Genome Atlas (TCGA) has identified a distinct subgroup of HPV(−) cancers with inactivating mutations in the nuclear receptor-binding SET domain protein 1 (NSD1), a methyltransferase and chromatin modifier [15, 16]. In a recent study, loss-of-function mutations in NSD1 were shown to be a favorable prognostic biomarker in laryngeal cancers, which are almost uniformly HPV(−) and strongly associated with tobacco and alcohol use [17]. However, no studies have examined the differences of NSD1 mutations between HPV-positive and HPV-negative head and neck cancers. Certainly, much remains to be elucidated regarding the role of NSD1 in oncogenesis, and knowledge of NSD1 direct targets is currently very limited. As for the distinct subgroup within HPV(−) head and neck cancers with inactivating NSD1 mutations, the clinical implications of this genetic alteration have yet to be fully explored [18].

In this study, we characterize NSD1 mutations in HPV(+) and HPV(−) cohorts from the expanded TCGA HNSCC dataset and examine prognostic implications of NSD1 mutations in the TCGA dataset. Finally, we investigate response of NSD1-depleted cells to platinum-based chemotherapy drugs and a variety of agents that have different mechanisms of DNA damage.


NSD1 mutations analysis in TCGA cohort

TCGA HNSCC database was mined for HPV(−) and HPV(+) samples (HPV(+) defined as ≥100 HPV reads per million human reads in whole-genome sequencing. Of the initial group of HPV(+) cases detected, those with p53 mutations were redesignated as HPV(−) cases, as the HPV DNA detected was likely not a driver of carcinogenesis. NSD1 mutations were identified with cBioPortal (available at: Figure 1 and Fig. 2 A were downloaded and adapted from the portal.

Fig. 1
figure 1

NSD1 mutations in the expanded TCGA HNSCC cohort. Alterations in NSD1 gene found in (a) HPV-negative (n = 457) and (b) HPV-positive (n = 65) HNSCC. Each tumor is represented as a column. The green arrow indicates tumors without alterations that were cut to fit the figure. Adapted from [41, 42]

Fig. 2
figure 2

NSD1 mutations in HPV(−) and HPV(+) HNSCC are inactivating mutations. (a) Schematic representation of missense point (green) and truncating (black) mutations in NSD1 found in HPV-negative HNSCC; functional NSD1 domains are indicated in colors. Adapted from [41, 42]. b Pathogenicity of HPV(−) NSD1 missense point mutations. (c) Pathogenicity of HPV(+) NSD1 mutations. Pathogenicity determined with UMD-predictor [21]

Kaplan-Meier survival analysis

Kaplan-Meier curves illustrating the survival of patients with HPV(−) and HPV(+) HNSCC with and without mutations in NSD1 were produced and downloaded from cBioPortal, and the log-rank (Mantel-Cox) test was used for statistical analyses.

Cell lines, constructs and chemicals

FaDu (human HPV(−) hypopharyngeal HNSCC) and SCC35 (human HPV(−) pharyngeal HNSCC) cells were cultured in DMEM/F12 medium supplemented with 0.4 μg/mL hydrocortisone. HaCat (spontaneously transformed non-tumorigenic human keratinocyte) cells were grown in DMEM with nonessential amino acids. All media was supplemented with 10% FBS (Invitrogen), 50 μg/mL penicillin, and 50 μg/mL streptomycin (Invitrogen). All cell lines tested negative for Mycoplasma and microsatellites authenticated.

Establishment of NSD1-depleted cells

FaDu [16], SCC35, and HaCat [19] cells were transfected with control (SHC016) or 2 different shRNAs against NSD1 (TRCN0000238372 and TRCN0000238371) obtained from Sigma. Stably transfected cells were selected using puromycin resistance. Clones were tested for NSD1 knockdown using quantitative RT-PCR with primer pairs listed below; GAPDH was used as a reference gene.

Primer pairs:













RNA extraction and quantitative RT-PCR

Total RNA was extracted by Qiagen RNA extraction kit and cDNA was synthesized using iScript cDNA Synthesis Kit (Bio-Rad) according to the manufacturer’s instructions. Quantitative real-time reverse transcription (qRT-PCR) was done using iQ SYBR Green Supermix (Bio-Rad) and indicated primer pairs on the iCycler iQ Real-Time PCR Detection System (Bio-Rad). Each qRT-PCR reaction was done in duplicate at least, and the ΔΔCt method was used to analyze the data; relative expression (% from GPDH) is presented.

Short-term proliferation assays

Cells (3000 cells/well) were plated in 96-well plates and treated with increasing doses of cisplatin (0–16 μM), carboplatin (0–20 μg/ml), etoposide (0–25 μM), zeocin (0–480 μg/ml), paclitaxel (0–160 nM), or mirin (0–80 μM) in triplicate. Proliferation was determined 72 h after treatment using the Cell Titer Glo reagent (Promega).

Long-term clonogenic survival assays

Cells (1000 cells/well) were plated in 6-well plates and treated with carboplatin (0–1 μg/ml for FaDu, 0–2 μg/ml for HaCat) in duplicates. After cell growth for 10 days for HaCat cells and 14 days for FaDu cells, colonies were fixed and stained with methylene blue in methanol for 1 h. Colonies consisting of at least 50 cells were counted.


Distinct subgroups of HPV(+) and HPV(−) tumors harbor inactivating NSD1 mutations

In the previously published TCGA HNSCC cohorts (n = 273 [15] or n = 421 [20]), NSD1 was identified as a significantly mutated gene. We found NSD1 gene alterations in 13% of HPV-negative (n = 457) and 9% of HPV-positive (n = 65) tumors (Fig. 1). The majority (61.3%) of NSD1 gene alterations in HPV(−) specimens were truncating mutations that occurred within or before the enzymatic SET domain of the NSD1 protein (Fig. 2A). The remaining alterations included homozygous gene deletions (6.7%), missense point mutations (30.7%) and a single in-frame deletion (1.3%). UMD-Predictor assigned 18 out of the 23 missense point mutations as pathogenic (Fig. 2B) [21]. Furthermore, mutant NSD1 was found to be underexpressed in HPV(−) tumors compared to wild-type NSD1 (Fig. 3A). Likewise, of the six HPV(+) NSD1 mutations, which were comprised of truncating (33%) and missense point (66%) mutations (Fig. 1), 3 of the 4 missense point mutations were predicted to be pathogenic or probably pathogenic by UMD-Predictor (Fig. 2C). However, in contrast to HPV(−) tumors, mutant NSD1 was not underexpressed in our limited set (n = 6) of HPV(+) head and neck cancer specimens (Fig. 3B), though this finding must be examined in a larger cohort before any definitive conclusions can be drawn.

Fig. 3
figure 3

Mutant NSD1 is underexpressed in HPV(−), but not in HPV(+) tumors. Scatterplot showing expression of wild-type and mutant NSD1 in (a) HPV(−) tumors and (b) HPV(+) tumors. Adapted from [41, 42]

Clinical characteristics of HPV(−) and HPV(+) tumors comparing those with and without NSD1 mutations

Statistical analysis of clinical characteristics of the HPV(−) cohort demonstrated that the NSD1 mutant and wild-type groups did not differ significantly in terms of average age at diagnosis or sex, race, alcohol use, or clinical stage (Table 1). Differences in smoking history were detected, with the cohort harboring NSD1 mutations having higher proportions of current and former smokers compared to those without NSD1 mutations (p < 0.001). Although only a small number of HPV(+) tumors were available for analysis, no differences between the NSD1 mutant and wild-type groups were found in average age at diagnosis, sex, race, smoking history, alcohol use, or stage (Table 2). Thus, NSD1 mutations are not associated with any clinical characteristics, except of smoking status in HPV-negative tumors.

Table 1 Clinical characteristics of TCGA HPV(−) cohort by NSD1 mutation status
Table 2 Clinical characteristics of TCGA HPV(+) cohort by NSD1 mutation status

NSD1 mutations are enriched in laryngeal HPV(−) tumors

Classification of the TCGA HPV(−) cohort by site revealed that laryngeal cancers were overrepresented amongst those with mutant NSD1: 25% of HPV(−) tumors arose in the larynx but almost 60% of HPV(−) tumors with NSD1 mutations occurred in the larynx (Fig. 4). In contrast, oral tongue tumors were underrepresented in the NSD1 mutant cohort, comprising 27% of HPV(−) tumors but only 11% of HPV(−) tumors with NSD1 mutations. In agreement with prior reports [17], NSD1 mutations are enriched in laryngeal HPV(−) tumors.

Fig. 4
figure 4

NSD1 mutations are enriched in HPV(−) laryngeal tumors. (a) Percentage by subsite of HPV(−) tumors (n = 457) in expanded TCGA cohort. (b) Percentage of NSD1 mutations (n = 60) in HPV(−) tumors, grouped by subsite

NSD1 mutations correlate with improved survival in HPV(−) HNSCC but worse survival in HPV(+) HNSCC

We next performed Kaplan-Meier survival analysis of head and neck patients in the expanded TCGA cohort based on HPV and NSD1 mutation status. Interestingly, patients with HPV-negative HNSCC whose tumors contain alterations in the NSD1 gene (n = 60), showed significantly improved overall (p = 0.006) and disease-free (p = 0.007) survival compared to patients with tumors that harbor wild-type NSD1 (n = 397) (Fig. 5A). Improved survival associated with NSD1 mutations in HPV-negative HNSCC is consistent with previously published work showing that NSD1 mutations are a favorable prognostic biomarker in head and neck tumors [17, 20]. In contrast, survival analysis on HPV(+) tumors with and without NSD1 mutations surprisingly revealed that NSD1 mutations were associated with significantly worse overall (p = 0.0019) and disease/progression-free survival (p = 0.002) (Fig. 5B). These data suggest that NSD1 mutations may have different impact on patients’ survival based on HPV status, although confirmation is needed in a larger HPV-positive HNSCC cohort.

Fig. 5
figure 5

NSD1 mutations correlate with improved survival in HPV(−) HNSCC and worse survival in HPV(+) HNSCC. Kaplan–Meier curves showing survival of patients with (a) HPV(−) tumors and (b) HPV(+) tumors, with or without alterations in NSD1 gene. Overall survival and disease/progression-free survival are shown

NSD1 mutations are associated with increased sensitivity to cisplatin

To begin exploring the factors underlying the survival differences between patients harboring tumors with wild-type and mutant NSD1, we used, a database in which > 1000 genetically characterized human cancer cell lines have been screened against a wide range of anti-cancer drugs [22]. Although the database did not include enough head and neck cancer cell lines for separate cohort analysis, the pan-cancer cohort showed significantly increased sensitivity (p = 0.0198) to cisplatin for cell lines with NSD1 mutations (n = 18) compared to those with wild-type NSD1 (n = 798) (Fig. 6A). Cisplatin, a DNA crosslinking agent, is the standard chemotherapy drug used as a radiation sensitizer for locally advanced HNSCC [23]. Increased sensitivity to cisplatin may in part explain why NSD1 mutations confer a survival advantage compared to wild-type NSD1 in HPV(−) tumors. In contrast to improved sensitivity to cisplatin, NSD1 mutant cell lines had the same sensitivity as NSD1 wild-type cells to several other types of DNA damaging agents: bleomycin (Fig. 6B), a radiomimetic drug that produces DNA single and double strand breaks [24]; etoposide (Fig. 6C), a topoisomerase II inhibitor that generates DNA double strand breaks (DSBs) [25, 26]; and paclitaxel (Fig. 6D), a DNA damage-inducing antimicrotubular drug.

Fig. 6
figure 6

Cells with NSD1 mutations show increased sensitivity to cisplatin, but not to other DNA-damaging agents. Pan-cancer cohort showing sensitivity of cells with or without NSD1 mutations to treatment with (a) cisplatin (18 mutant and 798 wild-type cell lines); (b) bleomycin (18 mutant and 817 wild-type cell lines); (c) etoposide (18 mutant and 834 wild-type cell lines); and (d) paclitaxel (8 mutant and 377 wild-type cell lines). Adapted from the Wellcome Sanger Institute website

Inactivation of ERCC5 has been shown to result in increased sensitivity to cisplatin in a variety of cancers [27,28,29]. The TCGA cohort was queried to determine if NSD1 mutations correlated with ERCC5 expression. Expression of ERCC5 was significantly lower in HPV(−) tumors with NSD1 mutations, compared to HPV(−) tumors with wild-type NSD1 (q = 0.0314) (Fig. 7A). However, for the limited number of HPV(+) tumors (n = 6 for NSD1 mutants), there was no difference in ERCC5 expression between those tumors with wild-type and mutated NSD1 (q = 0.975) (Fig. 7B).

Fig. 7
figure 7

Relative to NSD1 wild-type tumors, ERCC5 has lower expression in NSD1 mutant HPV(−) tumors, but not in NSD1 mutant HPV(+) tumors. Scatterplots showing ERCC5 expression in NSD1 wild-type and mutant HPV-negative (a) and HPV-positive (b) tumors. Adapted from [41, 42]

NSD1 knockdown increases cell sensitivity to DNA-crosslinking agents, but not to other DNA-damaging agents

With indications that NSD1 mutant and wild-type cells respond differently to cisplatin (Fig. 6A), we depleted NSD1 to determine the effect of NSD1 on survival of HNSCC cells or immortalized keratinocytes following treatment with platinum or other types of DNA-damaging agents. We first depleted NSD1 with specific shRNAs in FaDu, SCC35 (both are human head and neck cancer cell lines) cells, and HaCat (human immortalized, non-tumorigenic epithelial) cells, and selected knockdown clones using different shRNAs for each line that had > 50% depletion of NSD1 expression compared to parental cells or cells expressing control shRNA (Fig. 8).

Fig. 8
figure 8

Development of NSD1-depleted cell lines. FaDu and HaCat cells were transfected with control or NSD1 shRNAs. Puromycin-resistant clones were selected. NSD1 expression was determined by RT-PCR in wild-type FaDu (a) and HaCat (b) cells or cells expressing control or NSD1 shRNA. (csh = control shRNA; NSD1sh = NSD1 shRNA knockdown)

Cell viability assays of NSD1 wild-type (control) and NSD1-depleted (NSD1 sh) FaDu cells showed that NSD1 knockdown rendered cells more sensitive to increasing concentrations of carboplatin and cisplatin (Fig. 9A top and bottom). FaDu, HaCat, or SCC35 cells expressing another NSD1 shRNA (NSD1 sh a) were also more sensitive to cisplatin than cells expressing wild-type NSD1 (Fig. 9B, C, and D). Clonogenic survival of FaDu and HaCat cells demonstrated that increasing concentrations of carboplatin also inhibited survival of NSD1-depleted cells (Fig. 9E). Thus, NSD1 knockdown caused increased sensitivity of epithelial cells to platinum-based chemotherapy drugs in both short-term cell viability and long-term cell survival assays.

Fig. 9
figure 9

NSD1 knockdown increases cell sensitivity to DNA crosslinking agents. (a) Short-term viability assay using CellTiterGlo (Promega) of FaDu cells expressing control (black) or NSD1 (red, NSD1sh) shRNA, treated or not with carboplatin or cisplatin (as indicated) for 72 h. Short-term viability of FaDu (b), HaCat (c), or SCC35 (d) cells expressing control (black) or another NSD1 (red, NSD1sh a) shRNA, treated or not with cisplatin. (e) Long-term clonogenic survival of FaDu (top) and HaCat (bottom) cells expressing control (black) or NSD1 (red) shRNA, treated or not with carboplatin. * Indicates P < 0.05

As opposed to platinated drugs, cell viability assays showed that depletion of NSD1 did not sensitize cells to etoposide, zeocin (a chemical analog of radiomimetic drug bleomycin), paclitaxel, or mirin, an MRE11 inhibitor that prevents activation of a DNA damage sensor molecule Ataxia Telangiectasia Mutated (ATM) and DNA damage signaling (Fig. 10 and B). Conversely, NSD1 knockdown cells were relatively resistant to treatment with increasing concentrations of zeocin, paclitaxel, and mirin. Etoposide had a concentration-dependent effect on FaDu cells, with NSD1 knockdowns being more resistant to etoposide at low concentrations (5 μM and 10 μM), but more sensitive to etoposide at higher concentrations (20 μM and 25 μM).

Fig. 10
figure 10

NSD1 knockdown does not increase cell sensitivity to other types of DNA-damaging agents. Short-term viability assay using CellTiterGlo (Promega) of FaDu (a) or SCC35 (b) cells expressing control (black) or NSD1 (red) shRNA, treated or not with etoposide, zeocin, paclitaxel, and mirin. * Indicates P < 0.05

Expression of ERCC5 was decreased in HPV(−) tumors with NSD1 mutations, compared to HPV(−) tumors with wild-type NSD1 (Fig. 7A). ERCC5 mRNA levels were also lower in HaCat and FaDu, but not in SCC35 cells expressing NSD1 shRNAs, as compared to control shRNA expressing cells (Fig. 11).


A histone H3 lysine 36 (H3K36) methyltransferase and chromatin modifier, NSD1 plays a role in epigenetic DNA regulation and is essential in normal human development [30]. NSD1 is also a powerful transcriptional co-factor with the ability to repress and activate expression of many genes [31]. Germline inactivating mutations in NSD1 are responsible for Sotos syndrome, a childhood overgrowth condition in which patients are at increased risk of cancer and display an NSD1 mutation-specific DNA methylation signature [32,33,34]. Loss-of-function mutations in NSD1 are associated with global genome hypomethylation in head and neck cancer [15, 16]. Inactivating NSD1 mutations have also been described in several other malignancies, including lung cancer, clear cell renal carcinoma, and skin cancer, suggesting a role of NSD1 as a tumor suppressor [35,36,37,38]. In contrast, in acute myeloid leukemia, a recurring translocation that fuses NSD1 to nucleoporin-98 (NUP98) leads to leukemogenesis through dysregulated H3K36 methylation at certain genetic loci [39].

Alterations in NSD1 expression have already shown promise as a prognostic biomarker in various cancers. In neuroblastomas and gliomas, NSD1 promoter methylation-associated gene silencing predicts worse survival [40]. In laryngeal cancers, loss-of-function mutations in NSD1 were shown to be a favorable prognostic biomarker [17]. While NSD1 mutations were also present in oral, oropharyngeal, and hypopharyngeal cancers, NSD1 mutations at non-laryngeal sites were not associated with better overall survival or recurrence-free survival, suggesting molecular differences between anatomical subsites [17].

Unique epidemiological, molecular, biological and clinical differences have led to the increasing recognition of HPV-positive head and neck cancer as distinct from HPV-negative HNSCC. Our data examined NSD1 mutations in both HPV(+) and HPV(−) tumors in the expanded TCGA HNSCC dataset and identified a subgroup of inactivating NSD1 mutations in both the HPV(+) and HPV(−) cohorts (Figs. 1 and 2). NSD1 mutations correlated with improved survival in HPV(−) tumors (Fig. 5 A), consistent with previously published work in laryngeal tumors [17, 20]. The worse survival of HPV(−) patients with NSD1 wild-type tumors did not correlate with worse clinical characteristics; in fact these patients smoked less than the HPV(−) patients with NSD1 mutations (Table 1). While this manuscript was under revision, another study demonstrated similar results with the frequency of NSD1 mutations in HPV-negative heavy smokers (> 20 pack year) of 18.5%, 4.55% in HPV-negative light smokers (≤ 20 pack year), and 2.22% in HPV-negative never smokers [43].

Interestingly, NSD1 mutations were associated with decreased survival in patients with HPV(+) HNSCCs (Fig. 5 B). No differences in clinical characteristics were found that could explain the survival difference between HPV(+) patients with and without NSD1 mutations (Table 2). Although the TCGA is the largest molecularly characterized cohort of HPV(+) HNSCC publicly available, these statistically significant survival results are based on a small number of NSD1 mutant tumors that are HPV(+) (n = 6) and require further validation in a separate larger cohort.

NSD1 mutations were associated with differences in gene expression of ERCC5 based on HPV status (Fig. 7). We demonstrated that ERCC5 expression was lower in HPV(−) but not HPV(+) tumors that harbored NSD1 mutations (compared to HPV(−) or HPV(+) NSD1 wild-types, respectively) (Fig. 7). Moreover, ERCC5 mRNA levels were lower in two of three NSD1 shRNA expressing cell lines tested (Fig. 11). ERCC5 is a single strand-specific nuclease that participates in nucleotide excision repair. Inactivation of ERCC5 has been shown to result in increased sensitivity to cisplatin in a variety of cancers, including non-small cell lung cancer, epithelial ovarian cancer, and uroepithelial cancer [27,28,29]. Given the strong involvement of ERCC5 in response to cisplatin treatment, whether or not lower expression of ERCC5 contributes to the increased cisplatin sensitivity of NSD1 mutant HPV(−) HNSCCs deserves further investigation. Given opposite correlation of NSD1 mutations with survival and ERCC5 expression in HPV(+) versus HPV(−) HNSCC, the concept that NSD1 may play different roles in epigenetic DNA regulation based on HPV status should be additionally explored.

Fig. 11
figure 11

Relative expression of NSD1 and ERCC5 in HaCat (top), Fadu (middle), and SCC35 (bottom) cells expressing control or NSD1 shRNAs as determined by RT-PCR

Using isogenic cell lines, our data showed that even partial downregulation of NSD1 (25–50%, Figs. 8 and 11) sensitized HPV(−) HNSCC cells and immortalized keratinocytes to cisplatin and carboplatin (Fig. 9), which may, at least in part, explain the survival advantage of NSD1 mutations in HPV(−) HNSCC (Fig. 5 A). In contrast to platinated compounds, FaDu and SCC35 cells with shRNA-mediated depletion of NSD1 became somewhat more resistant to several DNA damaging drugs, including zeocin, etoposide and paclitaxel (Fig. 10), while available databases of various cancer cell lines demonstrated that NSD1 mutation status did not alter response to the same treatment (Fig. 6 B). The most likely explanation for that is a partial NSD1 depletion in isogenic cells in our experiments versus 100% NSD1 enzymatic inactivation due to mutation in cells derived from different origins of cancer.

The improved treatment response of NSD1 knockdown was limited to the DNA crosslinking agents cisplatin and carboplatin and was not seen with other types of DNA damaging agents tested: etoposide (topoisomerase II inhibitor), zeocin (radiomimetic), paclitaxel (microtubule inhibitor, also used in head and neck cancer treatment), and mirin (MRE11 inhibitor). Together, these data suggest that platinum-based chemotherapy better targets HPV(−) tumors with NSD1 mutations, while those without NSD1 mutations (which also have a worse prognosis) may require alternative treatments.


In this study, we identify in the expanded TCGA HNSCC database a subgroup of inactivating NSD1 mutations in both the HPV(+) and HPV(−) cohorts. NSD1 mutations appear to be associated with opposite survival effects in HPV(+) and HPV(−) tumors, with NSD1 mutations correlating with improved survival in HPV(−) tumors but worse survival in HPV(+) tumors. However, the prognostic role of NSD1 mutations in HPV(+) tumors must currently be taken with caution, as it has yet to be examined in a larger cohort. In HPV(−) tumors, the survival advantage of NSD1 mutations can in part be explained by increased sensitivity to the DNA-crosslinking agents cisplatin and carboplatin.

Availability of data and materials

This review article is original and has not been published elsewhere.



Ataxia-telangiectasia mutated serine/threonine kinase


Control shRNA


Double strand breaks


Head and neck squamous cell carcinoma


Human papillomavirus


Nuclear receptor-binding SET domain protein 1


NSD1 shRNA knockdown




Oropharyngeal squamous cell carcinoma


The Cancer Genome Atlas


  1. 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  Google Scholar 

  2. Götz C, Drecoll E, Straub M, Bissinger O, Wolff K-D, Kolk A. Impact of HPV infection on oral squamous cell carcinoma. Oncotarget. 2016;7:76704–12.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Rettig EM, D’Souza G. Epidemiology of head and neck Cancer. Surg Oncol Clin N Am. 2015;24:379–96.

    Article  PubMed  Google Scholar 

  4. Shiboski CH, Schmidt BL, Jordan RCK. Tongue and tonsil carcinoma: increasing trends in the U.S. population ages 20-44 years. Cancer. 2005;103:1843–9.

    Article  Google Scholar 

  5. Frisch M, Hjalgrim H, Jaeger AB, Biggar RJ. Changing patterns of tonsillar squamous cell carcinoma in the United States. Cancer Causes Control. 2000;11:489–95.

    Article  CAS  Google Scholar 

  6. Sturgis EM, Cinciripini PM. Trends in head and neck cancer incidence in relation to smoking prevalence. Cancer. 2007;110:1429–35.

    Article  PubMed  Google Scholar 

  7. Mehta V, Yu G-P, Schantz SP. Population-based analysis of oral and oropharyngeal carcinoma: changing trends of histopathologic differentiation, survival and patient demographics. Laryngoscope. 2010;120:2203–12.

    Article  PubMed  Google Scholar 

  8. Pytynia KB, Dahlstrom KR, Sturgis EM. Epidemiology of HPV-associated oropharyngeal cancer. Oral Oncol. 2014;50:380–6.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Gillison ML. Current topics in the epidemiology of oral cavity and oropharyngeal cancers. Head Neck. 2007;29:779–92.

    Article  Google Scholar 

  10. Goon PKC, Stanley MA, Ebmeyer J, Steinsträsser L, Upile T, Jerjes W, et al. HPV & head and neck cancer: a descriptive update. Head Neck Oncol. 2009;1:36.

    Article  Google Scholar 

  11. 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  Google Scholar 

  12. Dahlstrom KR, Calzada G, Hanby JD, Garden AS, Glisson BS, Li G, et al. An evolution in demographics, treatment, and outcomes of oropharyngeal cancer at a major cancer center. Cancer. 2013;119:81–9.

    Article  PubMed  Google Scholar 

  13. 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:612–9.

    Article  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  Google Scholar 

  15. Cancer Genome Atlas NetworkComprehensive genomic characterization of head and neck squamous cell carcinomas. Nature. 2015;517:576–82.

  16. Papillon-Cavanagh S, Lu C, Gayden T, Mikael LG, Bechet D, Karamboulas C, et al. Impaired H3K36 methylation defines a subset of head and neck squamous cell carcinomas. Nat Genet. 2017;49:180–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Peri S, Izumchenko E, Schubert AD, Slifker MJ, Ruth K, Serebriiskii IG, et al. NSD1- and NSD2-damaging mutations define a subset of laryngeal tumors with favorable prognosis. Nat Commun. 2017;8:1772.

    Article  Google Scholar 

  18. Lucio-Eterovic AK, Carpenter PB. An open and shut case for the role of NSD proteins as oncogenes. Transcription. 2011;2:158–61.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Cellosaurus HaCaT (CVCL_0038). Accessed 28 Dec 2018.

  20. Bui N, Huang JK, Bojorquez-Gomez A, Licon K, Sanchez KS, Tang SN, et al. Disruption of NSD1 in head and neck Cancer promotes favorable chemotherapeutic responses linked to Hypomethylation. Mol Cancer Ther. 2018;17:1585–94.

    Article  CAS  Google Scholar 

  21. Salgado D, Desvignes J-P, Rai G, Blanchard A, Miltgen M, Pinard A, et al. UMD-predictor: a high-throughput sequencing compliant system for pathogenicity prediction of any human cDNA substitution. Hum Mutat. 2016;37:439–46.

    Article  CAS  Google Scholar 

  22. Genomics of Drug Sensitivity in Cancer. Accessed 1 Dec 2018.

  23. DG Pfister, S Spencer. NCCN clinical practice guidelines in oncology: head and neck cancers (version 1. 2018). Accessed 11 Feb 2018.

  24. Povirk LF. DNA damage and mutagenesis by radiomimetic DNA-cleaving agents: bleomycin, neocarzinostatin and other enediynes. Mutat Res. 1996;355:71–89.

    Article  Google Scholar 

  25. Montecucco A, Zanetta F, Biamonti G. Molecular mechanisms of etoposide. EXCLI J. 2015;14:95–108.

    PubMed  PubMed Central  Google Scholar 

  26. Chen SH, Chan N-L, Hsieh T. New mechanistic and functional insights into DNA topoisomerases. Annu Rev Biochem. 2013;82:139–70.

    Article  CAS  Google Scholar 

  27. Sun X, Li F, Sun N, Shukui Q, Baoan C, Jifeng F, et al. Polymorphisms in XRCC1 and XPG and response to platinum-based chemotherapy in advanced non-small cell lung cancer patients. Lung Cancer. 2009;65:230–6.

    Article  Google Scholar 

  28. Ceccaldi R, O’Connor KW, Mouw KW, Li AY, Matulonis UA, D’Andrea AD, et al. A unique subset of epithelial ovarian cancers with platinum sensitivity and PARP inhibitor resistance. Cancer Res. 2015.

    Article  CAS  Google Scholar 

  29. Bambury RM, Jordan E, Zabor EC, Bouvier N, Al-Ahmadie H, Boyd ME, et al. Association of somatic mutations in DNA damage repair (DDR) genes with efficacy of platinum-based chemotherapy in advanced urothelial carcinoma. JCO 2015;33 15_suppl:4532–4532. doi:

    Article  Google Scholar 

  30. Wagner EJ, Carpenter PB. Understanding the language of Lys36 methylation at histone H3. Nat Rev Mol Cell Biol. 2012;13:115–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Huang N, vom Baur E, Garnier JM, Lerouge T, Vonesch JL, Lutz Y, et al. Two distinct nuclear receptor interaction domains in NSD1, a novel SET protein that exhibits characteristics of both corepressors and coactivators. EMBO J. 1998;17:3398–412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kurotaki N, Imaizumi K, Harada N, Masuno M, Kondoh T, Nagai T, et al. Haploinsufficiency of NSD1 causes Sotos syndrome. Nat Genet. 2002;30:365–6.

    Article  CAS  Google Scholar 

  33. Imaizumi K, Kimura J, Matsuo M, Kurosawa K, Masuno M, Niikawa N, et al. Sotos syndrome associated with a de novo balanced reciprocal translocation t (5;8)(q35;q24.1). Am J Med Genet. 2002;107:58–60.

    Article  Google Scholar 

  34. Douglas J, Hanks S, Temple IK, Davies S, Murray A, Upadhyaya M, et al. NSD1 mutations are the major cause of Sotos syndrome and occur in some cases of weaver syndrome but are rare in other overgrowth phenotypes. Am J Hum Genet. 2003;72:132–43.

    Article  CAS  Google Scholar 

  35. Tonon G, Wong K-K, Maulik G, Brennan C, Feng B, Zhang Y, et al. High-resolution genomic profiles of human lung cancer. Proc Natl Acad Sci U S A. 2005;102:9625–30.

    Article  CAS  Google Scholar 

  36. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013;499:43–9.

    Article  Google Scholar 

  37. Su X, Zhang J, Mouawad R, Compérat E, Rouprêt M, Allanic F, et al. NSD1 inactivation and SETD2 mutation drive a convergence toward loss of function of H3K36 writers in clear cell renal cell carcinomas. Cancer Res. 2017;77:4835–45.

    Article  CAS  Google Scholar 

  38. Quintana RM, Dupuy AJ, Bravo A, Casanova ML, Alameda JP, Page A, et al. A transposon-based analysis of gene mutations related to skin cancer development. J Invest Dermatol. 2013;133:239–48.

    Article  CAS  Google Scholar 

  39. Wang GG, Cai L, Pasillas MP, Kamps MP. NUP98-NSD1 links H3K36 methylation to Hox-a gene activation and leukaemogenesis. Nat Cell Biol. 2007;9:804–12.

    Article  CAS  Google Scholar 

  40. Berdasco M, Ropero S, Setien F, Fraga MF, Lapunzina P, Losson R, et al. Epigenetic inactivation of the Sotos overgrowth syndrome gene histone methyltransferase NSD1 in human neuroblastoma and glioma. PNAS. 2009;106:21830–5.

    Article  PubMed  Google Scholar 

  41. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.

    Article  Google Scholar 

  42. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:pl1.

    Article  Google Scholar 

  43. Ghasemi F, Prokopec SD, MacNeil D, Mundi N, Gameiro SF, Howlett C, et al. Mutational analysis of head and neck squamous cell carcinoma stratified by smoking status. JCI Insight. 2019;4(1):e123443.

    Article  Google Scholar 

Download references


The authors thank the Department of Otolaryngology/Head and Neck Surgery at the University of North Carolina at Chapel Hill and the Yale Department of Surgery, Otolaryngology Division for the support.


Financial support was provided by the Yale Department of Surgery Otolaryngology Division and the Department of Otolaryngology/Head and Neck Surgery University of North Carolina at Chapel Hill. Funding was also provided by the NIH-NHLBI Medical Student Research Fellowship, the William U. Gardner Memorial Student Research Fellowship, and the Rickard K. Gershon, MD, Student Research Fellowship.

Author information

Authors and Affiliations



CP: Study design, acquisition of data, analysis and interpretation of data, and article preparation. SI: Acquisition of data. WGY: Study concept and design, analysis and interpretation of data, and article preparation. NI: Study concept and design, acquisition of data; analysis and interpretation of data, and article preparation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Natalia Issaeva.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pan, C., Izreig, S., Yarbrough, W.G. et al. NSD1 mutations by HPV status in head and neck cancer: differences in survival and response to DNA-damaging agents. Cancers Head Neck 4, 3 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: