Identification of CT-based non-invasive radiomic biomarkers for overall survival prediction in oral cavity squamous cell carcinoma.

This study addresses the limited non-invasive tools for Oral Cavity Squamous Cell Carcinoma (OSCC) survival prediction by identifying Computed Tomography (CT)-based biomarkers to improve prognosis prediction. A retrospective analysis was conducted on data from 149 OSCC patients, including CT radiomics and clinical information. An ensemble approach involving correlation analysis, score screening, and the Sparse-L1 algorithm was used to select functional features, which were then used to build Cox Proportional Hazards models (CPH). Our CPH achieved a 0.70 concordance index in testing. The model identified two CT-based radiomics features, Gradient-Neighboring-Gray-Tone-Difference-Matrix-Strength (GNS) and normalized-Wavelet-LLL-Gray-Level-Dependence-Matrix-Large-Dependence-High-Gray-Level-Emphasis (HLE), as well as stage and alcohol usage, as survival biomarkers. The GNS group with values above 14 showed a hazard ratio of 0.12 and a 3-year survival rate of about 90%. Conversely, the GNS group with values less than or equal to 14 had a 49% survival rate. For normalized HLE, the high-end group (HLE > - 0.415) had a hazard ratio of 2.41, resulting in a 3-year survival rate of 70%, while the low-end group (HLE ≤ - 0.415) had a 36% survival rate. These findings contribute to our knowledge of how radiomics can be used to predict the outcome so that treatment plans can be tailored for patients people with OSCC to improve their survival.

Identification of CT-based non-invasive Radiographic Biomarkers for Overall Survival Stratification in Oral Cavity Squamous Cell Carcinoma.

This study addresses the limited non-invasive tools for Oral Cavity Squamous Cell Carcinoma OSCC survival prediction by identifying Computed Tomography (CT)-based biomarkers for improved prognosis. A retrospective analysis was conducted on data from 149 OSCC patients, including radiomics and clinical. An ensemble approach involving correlation analysis, score screening, and the Sparse-L1 algorithm was used to select functional features, which were then used to build Cox Proportional Hazards models (CPH). Our CPH achieved a 0.70 concordance index in testing. The model identified two CT-based radiomics features, Gradient-Neighboring-Gray-Tone-Difference-Matrix-Strength (GNS) and normalized-Wavelet-LLL-Gray-Level-Dependence-Matrix-Large-Dependence-High-Gray-Level-Emphasis (HLE), as well as smoking and alcohol usage, as survival biomarkers. The GNS group with values above 14 showed a hazard ratio of 0.12 and a 3-year survival rate of about 90%. Conversely, the GNS group with values less than or equal to 14 had a 49% survival rate. For normalized HLE, the high-end group (HLE > -0.415) had a hazard ratio of 2.41, resulting in a 3-year survival rate of 70%, while the low-end group (HLE <= -0.415) had a 36% survival rate. These findings contribute to our knowledge of how radiomics can be used to anticipate the outcome and tailor treatment plans from people with OSCC.

Simultaneously targeting ErbB family kinases and PI3K in HPV-positive head and neck squamous cell carcinoma.

OBJECTIVES: To identify the most effective PI3K and EGFR inhibitors in HPV-positive head and neck squamous cell carcinoma (HNSCC) and investigate the efficacy of a combination of an ErbB family kinase inhibitor and a PI3K inhibitor to inhibit cell proliferation of HPV-positive HNSCC. MATERIALS AND METHOD: HPV-positive HNSCC cell lines were treated with the FDA approved ErbB kinase inhibitor, Afatinib or FDA-approved PI3K inhibitor, Copanlisib, alone or in combination, and phosphorylation and total protein levels of cells were assessed by Western blot analysis.Cell proliferation and apoptosis were examined by MTS assay, flow cytometry, and Western blots, respectively. RESULTS: Copanlisib more effectively inhibited cell proliferation in comparison to other PI3K inhibitors tested. HPV-positive HNSCC cells differentially responded to cisplatin, Afatinib, or Copanlisib. The combination of Afatinib and Copanlisib more effectively suppressed cell proliferation and induced apoptosis compared to either treatment alone. Mechanistically, the combination of Afatinib and Copanlisib completely blocked phosphorylation of EGFR, HER2, HER3, and Akt as well as significantly decreased the HPV E7 expression compared to either treatment alone. CONCLUSION: Afatinib and Copanlisib more effectively suppress cell proliferation and survival of HPV-positive HNSCC in comparison to either treatment alone.

Epigenetic Regulations of Perineural Invasion in Head and Neck Squamous Cell Carcinoma.

Carcinomas of the oral cavity and oropharynx belong among the ten most common malignancies in the human population. The prognosis of head and neck squamous cell carcinoma (HNSCC) is determined by the degree of invasiveness of the primary tumor and by the extent of metastatic spread into regional and distant lymph nodes. Moreover, the level of the perineural invasion itself associates with tumor localization, invasion's extent, and the presence of nodal metastases. Here, we summarize the current knowledge about different aspects of epigenetic changes, which can be associated with HNSCC while focusing on perineural invasion (PNI). We review epigenetic modifications of the genes involved in the PNI process in HNSCC from the omics perspective and specific epigenetic modifications in OSCC or other neurotropic cancers associated with perineural invasion. Moreover, we summarize DNA methylation status of tumor-suppressor genes, methylation and demethylation enzymes and histone post-translational modifications associated with PNI. The influence of other epigenetic factors on the HNSCC incidence and perineural invasion such as tobacco, alcohol and oral microbiome is overviewed and HPV infection is discussed as an epigenetic factor associated with OSCC and related perineural invasion. Understanding epigenetic regulations of axon growth that lead to tumorous spread or uncovering the molecular control of axon interaction with cancer tissue can help to discover new therapeutic targets for these tumors.

Entinostat Decreases Immune Suppression to Promote Antitumor Responses in a HER2+ Breast Tumor Microenvironment.

Therapeutic combinations to alter immunosuppressive, solid tumor microenvironments (TME), such as in breast cancer, are essential to improve responses to immune checkpoint inhibitors (ICI). Entinostat, an oral histone deacetylase inhibitor, has been shown to improve responses to ICIs in various tumor models with immunosuppressive TMEs. The precise and comprehensive alterations to the TME induced by entinostat remain unknown. Here, we employed single-cell RNA sequencing on HER2-overexpressing breast tumors from mice treated with entinostat and ICIs to fully characterize changes across multiple cell types within the TME. This analysis demonstrates that treatment with entinostat induced a shift from a protumor to an antitumor TME signature, characterized predominantly by changes in myeloid cells. We confirmed myeloid-derived suppressor cells (MDSC) within entinostat-treated tumors associated with a less suppressive granulocytic (G)-MDSC phenotype and exhibited altered suppressive signaling that involved the NFκB and STAT3 pathways. In addition to MDSCs, tumor-associated macrophages were epigenetically reprogrammed from a protumor M2-like phenotype toward an antitumor M1-like phenotype, which may be contributing to a more sensitized TME. Overall, our in-depth analysis suggests that entinostat-induced changes on multiple myeloid cell types reduce immunosuppression and increase antitumor responses, which, in turn, improve sensitivity to ICIs. Sensitization of the TME by entinostat could ultimately broaden the population of patients with breast cancer who could benefit from ICIs.

Extracellular Vesicles Released by Tumor Endothelial Cells Spread Immunosuppressive and Transforming Signals Through Various Recipient Cells.

Head and neck squamous cell carcinoma (HNSCC) has a high recurrence and metastatic rate with an unknown mechanism of cancer spread. Tumor inflammation is the most critical processes of cancer onset, growth, and metastasis. We hypothesize that the release of extracellular vesicles (EVs) by tumor endothelial cells (TECs) induce reprogramming of immune cells as well as stromal cells to create an immunosuppressive microenvironment that favor tumor spread. We call this mechanism as non-metastatic contagious carcinogenesis. Extracellular vesicles were collected from primary HNSCC-derived endothelial cells (TEC-EV) and were used for stimulation of peripheral blood mononuclear cells (PBMCs) and primary adipose mesenchymal stem cells (ASCs). Regulation of ASC gene expression was investigated by RNA sequencing and protein array. PBMC, stimulated with TEC-EV, were analyzed by enzyme-linked immunosorbent assay and fluorescence-activated cell sorting. We validated in vitro the effects of TEC-EV on ASCs or PBMC by measuring invasion, adhesion, and proliferation. We found and confirmed that TEC-EV were able to change ASC inflammatory gene expression signature within 24-48 h. TEC-EV were also able to enhance the secretion of TGF-ß1 and IL-10 by PBMC and to increase T regulatory cell (Treg) expansion. TEC-EV carry specific proteins and RNAs that are responsible for Treg differentiation and immune suppression. ASCs and PBMC, treated with TEC-EV, enhanced proliferation, adhesion of tumor cells, and their invasion. These data indicate that TEC-EV exhibit a mechanism of non-metastatic contagious carcinogenesis that regulates tumor microenvironment and reprograms immune cells to sustain tumor growth and progression.

HPV E2, E4, E5 drive alternative carcinogenic pathways in HPV positive cancers.

The dominant paradigm for HPV carcinogenesis includes integration into the host genome followed by expression of E6 and E7 (E6/E7). We explored an alternative carcinogenic pathway characterized by episomal E2, E4, and E5 (E2/E4/E5) expression. Half of HPV positive cervical and pharyngeal cancers comprised a subtype with increase in expression of E2/E4/E5, as well as association with lack of integration into the host genome. Models of the E2/E4/E5 carcinogenesis show p53 dependent enhanced proliferation in vitro, as well as increased susceptibility to induction of cancer in vivo. Whole genomic expression analysis of the E2/E4/E5 pharyngeal cancer subtype is defined by activation of the fibroblast growth factor receptor (FGFR) pathway and this subtype is susceptible to combination FGFR and mTOR inhibition, with implications for targeted therapy.

Correction: Integrated single-cell and bulk gene expression and ATAC-seq reveals heterogeneity and early changes in pathways associated with resistance to cetuximab in HNSCC-sensitive cell lines.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Integrated single-cell and bulk gene expression and ATAC-seq reveals heterogeneity and early changes in pathways associated with resistance to cetuximab in HNSCC-sensitive cell lines.

BACKGROUND: Identifying potential resistance mechanisms while tumour cells still respond to therapy is critical to delay acquired resistance. METHODS: We generated the first comprehensive multi-omics, bulk and single-cell data in sensitive head and neck squamous cell carcinoma (HNSCC) cells to identify immediate responses to cetuximab. Two pathways potentially associated with resistance were focus of the study: regulation of receptor tyrosine kinases by TFAP2A transcription factor, and epithelial-to-mesenchymal transition (EMT). RESULTS: Single-cell RNA-seq demonstrates heterogeneity, with cell-specific TFAP2A and VIM expression profiles in response to treatment and also with global changes to various signalling pathways. RNA-seq and ATAC-seq reveal global changes within 5 days of therapy, suggesting early onset of mechanisms of resistance; and corroborates cell line heterogeneity, with different TFAP2A targets or EMT markers affected by therapy. Lack of TFAP2A expression is associated with HNSCC decreased growth, with cetuximab and JQ1 increasing the inhibitory effect. Regarding the EMT process, short-term cetuximab therapy has the strongest effect on inhibiting migration. TFAP2A silencing does not affect cell migration, supporting an independent role for both mechanisms in resistance. CONCLUSION: Overall, we show that immediate adaptive transcriptional and epigenetic changes induced by cetuximab are heterogeneous and cell type dependent; and independent mechanisms of resistance arise while tumour cells are still sensitive to therapy.

Correction: Novel Insight into Mutational Landscape of Head and Neck Squamous Cell Carcinoma.

[This corrects the article DOI: 10.1371/journal.pone.0093102.].

Chromatin structure regulates cancer-specific alternative splicing events in primary HPV-related oropharyngeal squamous cell carcinoma.

Human papillomavirus-related oropharyngeal squamous cell carcinoma (HPV+ OPSCC) represents a unique disease entity within head and neck cancer with rising incidence. Previous work has shown that alternative splicing events (ASEs) are prevalent in HPV+ OPSCC, but further validation is needed to understand the regulation of this process and its role in these tumours. In this study, eleven ASEs (GIT2, CTNNB1, MKNK2, MRPL33, SIPA1L3, SNHG6, SYCP2, TPRG1, ZHX2, ZNF331, and ELOVL1) were selected for validation from 109 previously published candidate ASEs to elucidate the post-transcriptional mechanisms of oncogenesis in HPV+ disease. In vitro qRT-PCR confirmed differential expression of 9 of 11 ASE candidates, and in silico analysis within the TCGA cohort confirmed 8 of 11 candidates. Six ASEs (MRPL33, SIPA1L3, SNHG6, TPRG1, ZHX2, and ELOVL1) showed significant differential expression across both methods. Further evaluation of chromatin modification revealed that ASEs strongly correlated with cancer-specific distribution of acetylated lysine 27 of histone 3 (H3K27ac). Subsequent epigenetic treatment of HPV+ HNSCC cell lines (UM-SCC-047 and UPCI-SCC-090) with JQ1 not only induced downregulation of cancer-specific ASE isoforms, but also growth inhibition in both cell lines. The UPCI-SCC-090 cell line, with greater ASE expression, also showed more significant growth inhibition after JQ1 treatment. This study confirms several novel cancer-specific ASEs in HPV+OPSCC and provides evidence for the role of chromatin modifications in regulation of alternative splicing in HPV+OPSCC. This highlights the role of epigenetic changes in the oncogenesis of HPV+OPSCC, which represents a unique, unexplored target for therapeutics that can alter the global post-transcriptional landscape.

Expression of Programmed Cell Death Ligand 1 and Associated Lymphocyte Infiltration in Olfactory Neuroblastoma.

BACKGROUND: Programmed cell death ligand 1 (PD-L1) is a transmembrane glycoprotein that interacts with the receptor programmed cell death 1 (PD-1) to suppress T-cell activation, reduce adjacent tissue damage, and promote tolerance to self-antigens. Tumors may express PD-L1 as a mechanism to evade immune detection. Recent clinical trials have demonstrated the efficacy of PD-L1/PD-1 antagonists through activation of tumor-infiltrated CD8+ T cells. The aim of this study was to determine the expression pattern of PD-L1 and PD-1 in olfactory neuroblastoma (ONB) tumor cells and to determine the presence of PD-1+ and CD8+ lymphocytes in the ONB immune microenvironment. METHODS: Immunohistochemistry for expression of PD-L1, PD-1, and CD8 was performed on paraffin-embedded ONB tissue. RESULTS: Of the 10 primary site ONB samples, 4 demonstrated positive PD-L1 expression. Of PD-L1+ tumors, the 2 highest expressing samples were found to contain PD-1+ tumor cells. Of the 4 available metastatic samples, all of which arose from PD-L1- primary site ONB, 3 were positive for PD-L1 and contained PD-1+ tumor cells. PD-L1+ primary and metastatic tumors also demonstrated increased PD-1+ infiltrating lymphocytes in the tumor and stroma (11.6- and 4.62-fold increase) compared with PD-L1- samples (P < 0.05 and P = 0.068 respectively). PD-L1+ specimens demonstrated increased CD8+ lymphocytes in the tumor and stroma (7.46- and 2.14-fold increase) compared with PD-L1- tumors (P < 0.05 for both). CONCLUSIONS: These data demonstrate that a proportion of ONB primary and metastatic tumors express PD-L1 and possess an associated tumor and stromal infiltrate of PD-1+ and CD8+ lymphocytes.

Differentially Methylated Super-Enhancers Regulate Target Gene Expression in Human Cancer.

Current literature suggests that epigenetically regulated super-enhancers (SEs) are drivers of aberrant gene expression in cancers. Many tumor types are still missing chromatin data to define cancer-specific SEs and their role in carcinogenesis. In this work, we develop a simple pipeline, which can utilize chromatin data from etiologically similar tumors to discover tissue-specific SEs and their target genes using gene expression and DNA methylation data. As an example, we applied our pipeline to human papillomavirus-related oropharyngeal squamous cell carcinoma (HPV + OPSCC). This tumor type is characterized by abundant gene expression changes, which cannot be explained by genetic alterations alone. Chromatin data are still limited for this disease, so we used 3627 SE elements from public domain data for closely related tissues, including normal and tumor lung, and cervical cancer cell lines. We integrated the available DNA methylation and gene expression data for HPV + OPSCC samples to filter the candidate SEs to identify functional SEs and their affected targets, which are essential for cancer development. Overall, we found 159 differentially methylated SEs, including 87 SEs that actively regulate expression of 150 nearby genes (211 SE-gene pairs) in HPV + OPSCC. Of these, 132 SE-gene pairs were validated in a related TCGA cohort. Pathway analysis revealed that the SE-regulated genes were associated with pathways known to regulate nasopharyngeal, breast, melanoma, and bladder carcinogenesis and are regulated by the epigenetic landscape in those cancers. Thus, we propose that gene expression in HPV + OPSCC may be controlled by epigenetic alterations in SE elements, which are common between related tissues. Our pipeline can utilize a diversity of data inputs and can be further adapted to SE analysis of diseased and non-diseased tissues from different organisms.

Chromatin dysregulation and DNA methylation at transcription start sites associated with transcriptional repression in cancers.

Although promoter-associated CpG islands have been established as targets of DNA methylation changes in cancer, previous studies suggest that epigenetic dysregulation outside the promoter region may be more closely associated with transcriptional changes. Here we examine DNA methylation, chromatin marks, and transcriptional alterations to define the relationship between transcriptional modulation and spatial changes in chromatin structure. Using human papillomavirus-related oropharyngeal carcinoma as a model, we show aberrant enrichment of repressive H3K9me3 at the transcriptional start site (TSS) with methylation-associated, tumor-specific gene silencing. Further analysis identifies a hypermethylated subtype which shows a functional convergence on MYC targets and association with CREBBP/EP300 mutation. The tumor-specific shift to transcriptional repression associated with DNA methylation at TSSs was confirmed in multiple tumor types. Our data may show a common underlying epigenetic dysregulation in cancer associated with broad enrichment of repressive chromatin marks and aberrant DNA hypermethylation at TSSs in combination with MYC network activation.

Publisher Correction: Chromatin dysregulation and DNA methylation at transcription start sites associated with transcriptional repression in cancers.

The original version of this Article contained an error in the author affiliations. Trey Ideker was incorrectly associated with 'Department of Medicine (Oncology), Stanford University School of Medicine, 875 Blake Wilbur Dr, Palo Alto, CA 94304, USA.' This has now been corrected in both the PDF and HTML versions of the Article.

PDGF enhances the protective effect of adipose stem cell-derived extracellular vesicles in a model of acute hindlimb ischemia.

We previously have shown that platelet-derived growth factor (PDGF) modulates the biological activity of extracellular vesicles released by adipose-derived mesenchymal stem cells (ASC-EVs). ASC-EVs may interact with blood and vessel cells by transferring proteins and nucleic acids and regulate their functions. In this study, we investigated immunomodulatory activity and protection from acute hindlimb ischemia of EVs released by PDGF-stimulated ASC (PDGF-EVs). PDGF treatment of ASC changed protein and RNA composition of released EVs by enhancing the expression of anti-inflammatory and immunomodulatory factors. In vitro, control EVs (cEVs) derived from non-stimulated ASC increased the secretion of both the IL-1b, IL-17, IFNγ, TNFα pro-inflammatory factors and the IL-10 anti-inflammatory factor, and enhanced the in vitro peripheral blood mononuclear cell (PBMC) adhesion on endothelium. In contrast, PDGF-EVs enhanced IL-10 secretion and induced TGF-ß1 secretion by PBMC. Moreover, PDGF-EVs stimulated the formation of T regulatory cells. In vivo, PDGF-EVs protected muscle tissue from acute ischemia, reduced infiltration of inflammatory cells and increased T regulatory cell infiltration in respect to cEVs. Our results suggest that PDGF-EVs are enriched in anti-inflammatory and immunomodulatory factors and induced in PBMC an enhanced production of IL-10 and TGF-ß1 resulting in protection of muscle from acute ischemia in vivo.

Functional characterization of alternatively spliced GSN in head and neck squamous cell carcinoma.

We have recently performed the characterization of alternative splicing events (ASEs) in head and neck squamous cell carcinoma, which allows dysregulation of protein expression common for cancer cells. Such analysis demonstrated a high ASE prevalence among tumor samples, including tumor-specific alternative splicing in the GSN gene.In vitro studies confirmed that overall expression of either ASE-GSN or wild-type GSN (WT-GSN) isoform inversely correlated with cell proliferation, whereas the high ratio of ASE-GSN to WT-GSN correlated with increased cellular invasion. Additionally, a change in expression of either isoform caused compensatory changes in expression of the other isoform. Our results suggest that the overall expression and the balance between GSN isoforms are mediating factors in proliferation, while increased overall expression of ASE-GSN is specific to cancer tissues. As a result, we propose ASE-GSN can serve not only as a biomarker of disease and disease progression, but also as a neoantigen for head and neck squamous cell carcinoma treatment, for which only a limited number of disease-specific targeted therapies currently exist.

HEY1 is expressed independent of NOTCH1 and is associated with poor prognosis in head and neck squamous cell carcinoma.

OBJECTIVES: Notch signaling is frequently altered in head and neck squamous cell carcinoma (HNSCC). However, the nature and clinical implications of this dysregulation are not well understood. We previously described an association of transcriptionally active NOTCH1 Intracellular Domain (NICD1) immunohistochemical (IHC) expression pattern with high-risk pathologic characteristics. Here we further characterize Notch signaling in HNSCC. MATERIALS AND METHODS: IHC expression patterns and clinicopathologic associations of Notch pathway molecules were evaluated among 78 tumors with known NOTCH1 mutation status. IHC was performed for JAG1, a NOTCH1 activating ligand, and HEY1, an NICD1 transcriptional target and Notch pathway activation marker. IHC pattern and H-score (% staining × intensity) were recorded and compared to clinicopathologic characteristics and survival. Survival was analyzed using Kaplan Meier method and Cox proportional hazards models (HR). RESULTS: JAG1 and NICD1 expression patterns were highly concordant among tumors without truncating NOTCH1 mutations (p < 0.001), but were dissimilar among tumors with truncating NOTCH1 mutations (p = 0.24). There was evidence for JAG1-independent NOTCH1 activation among seven tumors, all with wild-type NOTCH1. HEY1 expression was associated with neither JAG1 nor NICD1 expression, but was associated with NOTCH1 mutation status (p = 0.03). Twelve (16%) tumors expressed HEY1 but not NICD1. Higher HEY1 H-score was significantly associated with worse overall (adjusted hazard ratio [aHR] 2.0, 95% CI = 1.0-4.2) and disease-specific (aHR = 3.3, 95% CI = 1.4-7.9) survival, whereas JAG1 and NICD1 expression were not associated with survival. CONCLUSIONS: These findings suggest both NOTCH1-dependent and -independent HEY1 regulation, and imply a previously unrecognized prognostic role for HEY1 in HNSCC.

Integrated time course omics analysis distinguishes immediate therapeutic response from acquired resistance.

BACKGROUND: Targeted therapies specifically act by blocking the activity of proteins that are encoded by genes critical for tumorigenesis. However, most cancers acquire resistance and long-term disease remission is rarely observed. Understanding the time course of molecular changes responsible for the development of acquired resistance could enable optimization of patients' treatment options. Clinically, acquired therapeutic resistance can only be studied at a single time point in resistant tumors. METHODS: To determine the dynamics of these molecular changes, we obtained high throughput omics data (RNA-sequencing and DNA methylation) weekly during the development of cetuximab resistance in a head and neck cancer in vitro model. The CoGAPS unsupervised algorithm was used to determine the dynamics of the molecular changes associated with resistance during the time course of resistance development. RESULTS: CoGAPS was used to quantify the evolving transcriptional and epigenetic changes. Applying a PatternMarker statistic to the results from CoGAPS enabled novel heatmap-based visualization of the dynamics in these time course omics data. We demonstrate that transcriptional changes result from immediate therapeutic response or resistance, whereas epigenetic alterations only occur with resistance. Integrated analysis demonstrates delayed onset of changes in DNA methylation relative to transcription, suggesting that resistance is stabilized epigenetically. CONCLUSIONS: Genes with epigenetic alterations associated with resistance that have concordant expression changes are hypothesized to stabilize the resistant phenotype. These genes include FGFR1, which was associated with EGFR inhibitors resistance previously. Thus, integrated omics analysis distinguishes the timing of molecular drivers of resistance. This understanding of the time course progression of molecular changes in acquired resistance is important for the development of alternative treatment strategies that would introduce appropriate selection of new drugs to treat cancer before the resistant phenotype develops.

Splice Expression Variation Analysis (SEVA) for inter-tumor heterogeneity of gene isoform usage in cancer.

Motivation: Current bioinformatics methods to detect changes in gene isoform usage in distinct phenotypes compare the relative expected isoform usage in phenotypes. These statistics model differences in isoform usage in normal tissues, which have stable regulation of gene splicing. Pathological conditions, such as cancer, can have broken regulation of splicing that increases the heterogeneity of the expression of splice variants. Inferring events with such differential heterogeneity in gene isoform usage requires new statistical approaches. Results: We introduce Splice Expression Variability Analysis (SEVA) to model increased heterogeneity of splice variant usage between conditions (e.g. tumor and normal samples). SEVA uses a rank-based multivariate statistic that compares the variability of junction expression profiles within one condition to the variability within another. Simulated data show that SEVA is unique in modeling heterogeneity of gene isoform usage, and benchmark SEVA's performance against EBSeq, DiffSplice and rMATS that model differential isoform usage instead of heterogeneity. We confirm the accuracy of SEVA in identifying known splice variants in head and neck cancer and perform cross-study validation of novel splice variants. A novel comparison of splice variant heterogeneity between subtypes of head and neck cancer demonstrated unanticipated similarity between the heterogeneity of gene isoform usage in HPV-positive and HPV-negative subtypes and anticipated increased heterogeneity among HPV-negative samples with mutations in genes that regulate the splice variant machinery. These results show that SEVA accurately models differential heterogeneity of gene isoform usage from RNA-seq data. Availability and implementation: SEVA is implemented in the R/Bioconductor package GSReg. Contact: bahman@jhu.edu or favorov@sensi.org or ejfertig@jhmi.edu. Supplementary information: Supplementary data are available at Bioinformatics online.