YEATS2 is upregulated in head and neck cancer.

(A) Computational strategy used to find epigenetic factors with poor patient prognosis overexpressed in HNC. (B) Heatmap depicting the TCGA mRNA expression of top 10 epigenetic factors shortlisted using strategy mentioned in (A). (C) RT-qPCR result showing mRNA expression of YEATS2 in HNC samples (n=8). (D) Kaplan-Meier curve of YEATS2 showing association with overall head and neck cancer patient survival. (E-F) Gene expression profile of YEATS2 in publicly available HNC microarray datasets, GSE30784 (E) and GSE9844 (F). (G) Immunoblot showing protein levels of YEATS2 in nuclear protein lysates extracted from HNC tissues (n=23). (H) Volcano plot of differentially expressed genes in RNA-seq analysis of shControl vs. shYEATS2 in BICR10 cells. (I) Results of overrepresentation analysis of genes significantly downregulated in shControl vs. shYEATS2 RNA-seq data. Error bars, min to max; two-tailed t test, ∗∗∗p < 0.001, N-Normal, T-Tumor.

YEATS2 drives EMT in head and neck cancer cells.

(A) Immunoblot showing the expression levels of various EMT factors upon YEATS2 knockdown in BICR10. (B-C) Results of invasion assay with quantification (below), after knockdown (B) or overexpression (C) of YEATS2 in BICR10 (Scale bar, 200 μm). (D-E) Wound healing assay performed after knockdown (D) or overexpression (E) of YEATS2 in BICR10 (Scale bar, 275 μm). (F-G) Results of 3D invasion assay showing change in invasive potential of BICR10 cells in collagen matrix after silencing (F) or overexpression (G) of YEATS2 (quantification shown below). (Scale bar: 4X, 200 μm; 10X, 50 μm) Error bars, mean ± SEM; two-tailed t test, ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001, n= 3 biological replicates.

Regulation of EMT by YEATS2 is SP1-dependent in HNC

(A) Luciferase assay results showing the difference in relative luciferase activity of the two YEATS2 promoter deletion constructs in BICR10. (B) Plot showing decrease in mRNA expression of YEATS2 on SP1-knockdown in BICR10 cells. (C) Immunoblot showing the reduced expression level of YEATS2 upon SP1-knockdown in BICR10. (D) Plot depicting SP1 binding on YEATS2 promoter in SP1-ChIP-qPCR assay in BICR10. (E) Plot showing difference in relative luciferase activity of YEATS2 Luc-508 in shControl vs. shSP1 BICR10 cells. (F) Schematic showing mutation of SP1-binding site sequence in YEATS2 Luc-508 construct (above), and relative luciferase activity of wild-type (WT) vs. mutant (Mut) YEATS2 Luc-508 in BICR10 (below). (G) Immunoblot depicting the decreased Twist1 levels on SP1 knockdown and its subsequent rescue of expression upon YEATS2 overexpression in BICR10 (* indicates endogenous YEATS2 band). (H) Invasion assay images (with quantification on right) showing decrease and rescue of the percentage of invaded cells in shSP1 BICR10 cells, and shSP1 cells with YEATS2 overexpression, respectively. Scale bar, 200 μm. Error bars, mean ± SEM; two-tailed t test, ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001, n = 3 biological replicates.

YEATS2 and GCDH regulate histone crotonylation in HNC.

(A) GSEA plot showing enrichment of KEGG LYSINE_DEGRADATION_PATHWAY in TCGA samples stratified as YEATS2_high as compared to YEATS2_low samples (NES-net enrichment score, NOM p-nominal p-value). (B) Schematic depicting the canonical lysine degradation pathway highlighting the step that leads to crotonyl-CoA production via GCDH. (C) RT-qPCR result showing enhanced mRNA expression of GCDH in tumor vs. normal samples (n=8). (D) Immunoblot showing enhanced levels of H3K27cr (quantification on right) in nuclear lysates extracted from HNC tumor vs. normal samples (n=23). (E) Representative IHC images (quantification below) showing the levels of ECHS1, GCDH, YEATS2, H3K27cr and Twist1 in HNC normal vs. tumor tissue samples (n=8) (Scale bar, 100 μm). (F and H) Immunoblot depicting the decrease in H3K27cr levels on (F) YEATS2- and (H) GCDH-knockdown in BICR10 cells. (G and I) Immunoflorescence images depicting the decrease in H3K27cr levels on (G) YEATS2- and (I) GCDH-knockdown in BICR10 cells (Scale bar, 10 μm). Error bars, mean ± SEM; two-tailed t test, ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001. n = 3 biological replicates.

YEATS2 regulates expression of EMT-related SPARC in HNC.

(A) Venn diagram showing SPARC gene obtained after integration of RNA-seq data (genes downregulated in shControl vs. shYEATS2), YEATS2 ChIP-seq data and, hallmark EPITHELIAL_MESENCHYMAL_TRANSITION gene signature. (B) Integrative genome viewer (IGV) plot showing decrease in SPARC expression in shControl vs. shYEATS2 RNA-seq data. (C) YEATS2-ChIP-qPCR results showing decreased binding of YEATS2 on SPARC promoter in shYEATS2 BICR10 cells. (D) RT-qPCR results showing decreased expression of SPARC on YEATS2-knockdown in BICR10. (E) Immunoblot showing decreased expression of SPARC in conditioned media derived from shControl and shYEATS2 BICR10 cells. (F) p300-ChIP-qPCR results showing decreased binding of p300 on SPARC promoter in shYEATS2 BICR10 cells. (G) Immunoblot showing decreased SPARC levels on p300 knockdown in conditioned media from shControl and shEP300 BICR10 cells. (H) Immunoblot showing co-immunoprecipitation of YEATS2 by p300 in YEATS2-overexpressed BICR10 cells (endogenous YEATS2 band is highlighted by *). (I-J) Immunoblot depicting the overexpression of Flag-tagged SPARC in shYEATS2 BICR10 cells (I), and Invasion assay images (J) (with quantification below) showing decrease and rescue of the percentage of invaded cells in shControl vs. shSP1 BICR10 cells, and shSP1 cells with YEATS2 overexpression, respectively. Scale bar, 200 μm. ∗∗p < 0.01, ∗∗∗p < 0.001, n = 3 biological replicates.

Maintenance of H3K27cr marks is dependent on YEATS2-mediated recruitment of p300 on SPARC promoter.

(A) H3K27cr-ChIP-qPCR results showing decrease in H3K27cr enrichment on SPARC promoter on YEATS2 knockdown. (B) H3K27ac-ChIP-qPCR results showing non-significant change in H3K27ac enrichment on SPARC promoter on YEATS2 knockdown. (C) RNA Pol II ChIP-qPCR showing decrease in Pol II occupancy on SPARC after YEATS2 downregulation in BICR10. (D) H3K27cr-ChIP-qPCR results showing increase in H3K27cr enrichment on SPARC promoter on treating BICR10 cells with 2.5 mM sodium crotonate (NaCr). (E-F) RT-qPCR (E) and immunoblot (F) showing enhanced SPARC expression in untreated BICR10 cells vs. BICR10 cells treated with 2.5 mM NaCr. (G) Invasion assay images (with quantification below) showing increased invasion on treating BICR10 cells with 2.5 mM NaCr (Scale bar, 200 μm). (H) H3K27cr-ChIP-qPCR data depicting the lack of significant difference in H3K27cr levels between shYEATS2 vs. shYEATS2+ 2.5 mM NaCr cells. (I) Immunoblot showing inability of NaCr treatment to rescue SPARC expression in YEATS2-knockdown BICR10 cells. (J) H3K27cr ChIP-qPCR showing decrease in H3K27cr levels in BICR10 shControl vs. shEP300 cells. (K-L) p300 (K) and H3K27cr ChIP-qPCR (L) showing decrease and subsequent rescue in p300 binding and H3K27cr enrichment on SP1 knockdown and YEATS2 overexpression, respectively. Error bars, mean ± SEM; two-tailed t test, ns-non-significant, ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001.

GCDH expression is SP1-dependent and regulates H3K27cr-mediated SPARC expression with YEATS2 synergistically.

(A) RT-qPCR results showing decrease in SPARC expression on GCDH-knockdown in BICR10 cells. (B) Immunoblot showing the reduced expression of SPARC upon GCDH-knockdown in SCC9 cells. (C) Decrease in H3K27cr levels on SPARC promoter in shGCDH BICR10 cells. (D) Immunoblot showing the reduced expression of GCDH on SP1-knockdown in BICR10 cells. (E) Plot showing SP1 binding on GCDH promoter in SP1-ChIP assay in BICR10. (F) Invasion assay images (quantification shown below) showing decrease in invasion of BICR10 cells on GCDH knockdown (Scale bar, 200 μm). (G) IF images (quantification shown below) showing reduced nuclear localization of GCDH in shYEATS2 BICR10 cells. Scale bar, 5 μm. (H) Immunoblot depicting decreased SPARC expression on SP1 knockdown and its subsequent rescue upon dual YEATS2 and GCDH overexpression in shSP1 BICR10 cells. Error bars, mean ± SEM; two-tailed t test, ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001, n = 3 biological replicates.

YEATS2 regulates gene expression by maintaining H3K27cr levels globally.

(A) Differential ChIP-seq profile of H3K27cr in shControl vs. shYEATS2. (B) Overrepresentation analysis showing pathways enriched among genes with decreased H3K27cr enrichment in H3K27cr ChIP-seq data. (C) Overlap of genes with reduced H3K27cr on their promoter and downregulated genes in RNA-seq. (D-E) IGV plot showing representative examples of genes (CREB3L2 and ERG) common in (C).

Schematic showing YEATS2- and GCDH-mediated regulation of EMT through H3K27cr-dependent SPARC upregulation in head and neck cancer.