Next Article in Journal
Disruption of CCL2 in Mesenchymal Stem Cells as an Anti-Tumor Approach against Prostate Cancer
Previous Article in Journal
Dexamethasone-Induced Fatty Acid Oxidation and Autophagy/Mitophagy Are Essential for T-ALL Glucocorticoid Resistance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma

by
Alexandra Ruth Glathar
1,
Akinsola Oyelakin
1,
Kasturi Bala Nayak
1,
Jennifer Sosa
1,
Rose-Anne Romano
2 and
Satrajit Sinha
1,*
1
Department of Biochemistry, Jacobs School of Medicine and Biomedical Sciences, State University of New York at Buffalo, Buffalo, NY 14203, USA
2
Department of Oral Biology, School of Dental Medicine, State University of New York at Buffalo, Buffalo, NY 14214, USA
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(2), 446; https://doi.org/10.3390/cancers15020446
Submission received: 9 December 2022 / Revised: 4 January 2023 / Accepted: 7 January 2023 / Published: 10 January 2023

Abstract

:

Simple Summary

Oral squamous cell carcinomas (OSCC) are the most common malignancies affecting the oral cavity and account for 40% of all head and neck squamous cell carcinoma cases. Unfortunately, patient outcomes are generally unfavorable, as diagnosis normally occurs at a late stage of the disease, as well as a lack of effective targeted treatments. Our identification of potential prognostic biomarkers is of particular importance to the field and could lead to the generation of effective targeted therapies which may improve patient outcomes and survival. Here, we utilize a mouse model of OSCC to explore the role of p63 as an important oncogenic transcription factor in the control of OSCC proliferation and migration. We also generate a p63-driven gene expression signature for mouse OSCC which identifies both novel and conserved genes and pathways, which may be relevant in human disease and which may serve as potential biomarkers and targets for future therapeutics.

Abstract

Oral squamous cell carcinoma (OSCC) is the most common malignancy of the oral cavity and is linked to tobacco exposure, alcohol consumption, and human papillomavirus infection. Despite therapeutic advances, a lack of molecular understanding of disease etiology, and delayed diagnoses continue to negatively affect survival. The identification of oncogenic drivers and prognostic biomarkers by leveraging bulk and single-cell RNA-sequencing datasets of OSCC can lead to more targeted therapies and improved patient outcomes. However, the generation, analysis, and continued utilization of additional genetic and genomic tools are warranted. Tobacco-induced OSCC can be modeled in mice via 4-nitroquinoline 1-oxide (4NQO), which generates a spectrum of neoplastic lesions mimicking human OSCC and upregulates the oncogenic master transcription factor p63. Here, we molecularly characterized established mouse 4NQO treatment-derived OSCC cell lines and utilized RNA and chromatin immunoprecipitation-sequencing to uncover the global p63 gene regulatory and signaling network. We integrated our p63 datasets with published bulk and single-cell RNA-sequencing of mouse 4NQO-treated tongue and esophageal tumors, respectively, to generate a p63-driven gene signature that sheds new light on the role of p63 in murine OSCC. Our analyses reveal known and novel players, such as COTL1, that are regulated by p63 and influence various oncogenic processes, including metastasis. The identification of new sets of potential biomarkers and pathways, some of which are functionally conserved in human OSCC and can prognosticate patient survival, offers new avenues for future mechanistic studies.

Graphical Abstract

1. Introduction

Oral squamous cell carcinomas (OSCCs) are the most common malignancies arising from the mucous membrane of the oral cavity and typically affect the buccal and labial mucosae, gingiva, floor of the mouth, hard palate, lips, tongue, and cheeks [1,2,3]. OSCC is linked to tobacco exposure, alcohol consumption, and human papillomavirus infection and accounts for 40% of all head and neck squamous cell carcinoma (HNSCC) cases [1,2,3]. The outcomes of patients with OSCC are generally unfavorable because diagnosis typically occurs at a late stage of the disease and there are no effective treatments [1,2,3]. To this end, recent genomic and transcriptomic characterizations of human OSCC/HNSCC tumors have offered crucial molecular insights into the disease’s etiology and progression. [4,5,6]. Harnessing the power of such data-rich resources to identify biomarkers for targeted therapy, however, remains a challenge and necessitates parallel studies of OSCC, particularly in well-defined genetic models
In mice, 4-nitroquinoline 1-oxide (4NQO)-induced carcinogenesis produces a spectrum of preneoplastic and neoplastic lesions that mimic the progression, molecular, and histological changes of OSCC in humans [7,8]. Administration of 4NQO in drinking water or by topical application induces sequential changes in murine oral epithelium, from hyperplasia and dysplasia to the development of in situ carcinomas [7,8]. These lesions exhibit up to 94% similarity with the mutational landscape of human carcinogen-induced HNSCC [9], highlighting the potential of this model for identifying human-relevant biomarkers. Several studies have used 4NQO-induced carcinogenesis to uncover the molecular drivers of OSCC [9,10,11,12,13]. Deletion of p16INK4a, a gene in one of the most commonly mutated genomic loci in human OSCC, has been found to upregulate the expression of the oncogenic transcription factor Trp63 (p63) and enhance the development of 4NQO-induced tumors [12]. p63 is a master regulator of epithelial development, shown in many mouse and human studies to drive oncogenesis [12,13,14,15,16,17,18,19]. The human gene TP63 encodes two major isoforms, ΔNp63 and TAp63, with ΔNp63α being the predominant isoform expressed in most cells of epithelial origin and epithelial cancers [20,21,22,23]. TP63 is amplified and overexpressed in most OSCC tumors and promotes tumor growth by multiple mechanisms [12,13,14,15,16,17,18,19,20,21,22,23]. Interestingly, the downregulation of p63 in late-stage and metastatic tumors and a 4NQO model of HNSCC facilitates the metastatic migration of cancer cells to secondary sites [24,25]. Despite the overwhelming evidence of a strong link between p63 and various facets of epithelial cancers, the underlying molecular mechanisms and the full repertoire of the downstream target genes, by which p63 promotes and contributes to the development of OSCC, are not fully known.
Here, we have utilized a comprehensive repertoire of transcriptomic and genomic tools to uncover novel p63-driven oncogenic pathways and targets to better understand its role in murine OSCC (mOSCC). We perform the first comprehensive transcriptomic characterization using bulk RNA sequencing (RNA-seq) of 4NQO-induced mOSCC cell lines, highlighting the diversity of murine OSCC, as well as providing a robust dataset for future studies employing these cell lines. To investigate the molecular function of p63 in mOSCC and identify its transcriptional targets, we utilized both loss- and gain-of-function models of p63 to identify a well-defined p63-driven gene signature. We also show that our identified signature coincides well with signatures identified in recent publicly available bulk and single-cell RNA-seq of 4NQO-induced mOSCC tumors and mouse esophageal squamous cell carcinoma (ESCC) tumors, respectively. Results from our RNA-seq, Chromatin immunoprecipitation-sequencing (ChIP-seq), and phenotypic experiments with mOSCC cell lines demonstrate that p63 regulates a broad range of pathways and targets, which affects downstream pathways including cell adhesion, migration, and metastasis, as well as the novel target, Cotl1. Our findings reaffirm the notion that dynamic expression of p63, specifically the ΔNp63 isoform, regulates an extensive swath of transcriptional targets, as well as maintains broad oncogenic signaling. These gene regulatory networks are essential for a wide range of oncogenic processes in mOSCC and are also likely to be relevant in the human disease.

2. Materials and Methods

2.1. Cell Culture Studies

Cell lines B7E3, B4B8, and B7E11 were obtained from Dr. Carter Van Waes [26]. All cell lines were grown and maintained in high-glutamine DMEM with the following supplements: 10% FBS, 1% nonessential amino acids, and antibiotics. All cell lines were tested by the eMycoPlus Mycoplasma PCR detection kit (BulldogBio) to ensure they were bereft of any mycoplasma infection.

2.2. Generation of Knockdown and Overexpression of Cell Lines

Lentivirus-mediated depletion of p63 in B7E11 cells was performed using the pGIPZ system. GIPZ lentiviruses containing short hairpin RNAs (shRNAs) (clone IDs: V3LMM_508694 (sh1), V3LMM_523845 (sh2), and V3LMM_523849 (sh3)) targeting Trp63 were generated with the help of the Gene Modulation Services Shared Core at Roswell Park Comprehensive Cancer Center. The pINDUCER20 vector comprising a neomycin resistance gene was used as previously described [27] for doxycycline-inducible expression of FLAG-tagged ΔNp63α in B7E3 cells. The FLAG-tagged ΔNp63α cDNA (a gift from Caterina Missero) was transferred to pInducer20 using recombination Gateway LR clonase Enzyme mix kit (Thermofisher, Waltham, MA, USA) and the resulting plasmid was confirmed by DNA sequencing. Viral infection and selection with either puromycin or neomycin were performed as described before [28].

2.3. Western Blot Analysis

Western blots were generated according to a previously described protocol [28]. Briefly, 5 µL protein lysates were loaded onto SDS-polyacrylamide gels and transferred to Immun-Blot polyvinylidene difluoride membranes (Bio-Rad Laboratories, Hercules, CA, USA). After blocking in 5% milk, the membranes were incubated first in primary antibodies against p63 (4A4, 1:20,000), COTL1 (Proteintech, 1:10,000), K14 (a gift from Dr. Rose-Anne Romano) [29], Vimentin (CST, 1:5000), MMP9 (Proteintech, 1:10,000), Fibronectin (SinoBiological, 1:5000), ITGB4 (Proteintech, 1:10,000), E-cadherin (CST, 1:5000), and K6 (a gift from Dr. Julie Segre), then with horseradish peroxidase-conjugated secondary antibodies corresponding to the host of the primary antibody, and then washed in Tris-buffered saline with 0.05% Tween-20. Protein expression was detected with the LumiGLO peroxidase chemiluminescent substrate kit (SeraCare, Milford, MA, USA), and membranes were imaged using a Bio-Rad ChemiDoc imaging system. Uncropped Western blot images can be found in Figures S7–S9.

2.4. ChIP-seq of p63/ΔNp63

The iDeal ChIP-seq kit for transcription factors (C01010055; Diagenode, Denville, NJ, USA) and the associated protocols were used to perform ChIP-seq. B7E11 cells were grown to ~90% confluency and cross-linked in the supplied fixation buffer supplemented with 0.5% formaldehyde for 10 min. Lysates from the fixed cells were subsequently sonicated with a Diagenode Bioruptor to obtain sheared chromatin with an approximate fragment length of 150–400 bp. The ChIP experiments for p63 were carried out using 2 µg of pan-p63 4A4 antibody (Santa Cruz Biotechnology) and 2 µg of ΔNp63 antibody (E6Q3O; Cell Signaling Technology, Danvers, MA, USA). Libraries were prepared as described before [30]. ChIP DNA and input controls were then subjected to 50 bp single-end sequencing on an Illumina HiSeq 2500 instrument (Illumina Inc., San Diego, CA, USA), which resulted in 15–25 million reads per sample.

2.5. ChIP-seq Analysis

The raw ChIP-seq reads from all experiments were mapped to the Mus musculus genome (mm10) using Bowtie2 v2.3.4.1. Peak calling was then performed using MACS2 v2.1.0, with a minimum false-discovery rate cutoff of 0.05. Sequenced genomic input was used as the control for each experiment, and resultant peaks were matched to the nearest gene using GREAT analysis with default settings [31,32]. To visualize ChIP peaks, bam files were processed with deepTools v3.3.2 to generate bigwig files, which were then uploaded to IGV [33]. The R tool ChIPseeker was used to annotate p63 ChIP-seq peaks to the nearest genomic feature of the mm10 genome assembly [34]. Adobe Illustrator was used for final image processing.

2.6. RNA Isolation and Library Preparation for RNA-seq

Total RNA from cell lines was extracted using a Direct-zol RNA miniprep kit (Zymo Research, Irvine, CA, USA). The extracted RNA was snap-frozen on dry ice and stored at −80 °C until library preparation. For each RNA sample, cDNA libraries were prepared using the TrueSeq RNA sample preparation kit (Illumina). The libraries underwent 50 bp single-end sequencing on an Illumina HiSeq 2500 instrument. Quality control metrics were performed on raw sequencing reads using the FASTQC v0.11.9 application.

2.7. RNA-seq Analysis

Reads were mapped to the appropriate reference genome, GRCm38/mm10 build, with HISAT v2.1.0 [35]. Reads aligning to the reference genome were quantified with featureCounts v1.5.3 to generate a matrix of raw counts, which was then processed in R, a free program for statistical computing and graphics, to generate normalized expression values in transcripts per million, according to the method proposed by Wagner et al. [36]. Differential gene expression analysis comparing control to p63 knockdown was carried out using DESeq2 v1.24.0 [37]. Differentially expressed genes with a false-discovery rate of ≤0.1 were considered statistically significant.

2.8. qRT-PCR Analysis

Total RNA from B7E11 knockdown and B7E3 overexpression cell lines were extracted using the Direct-zol RNA miniprep kit. The RNA was reverse transcribed with the Bio-Rad iScript cDNA synthesis kit, according to the manufacturer’s instructions. The resulting cDNA was used for qPCR with Bio-Rad iQ SYBR green Supermix. A list of the qRT-PCR primers can be found in Table S9.

2.9. HNSCC and mOSCC Dataset Analysis

RNA-seq data from patient samples were obtained from The Cancer Genome Atlas (TCGA). TCGA RNA-seq expression and survival datasets were downloaded from cBioPortal and the UCSC Xena Browser [38,39,40]. RNA-seq fastq files from mouse 4NQO-generated oral tumors were downloaded from GSE54246 [11] and realigned with HISAT v2.1.0 to the mm10 genome, and then quantified using featureCounts v1.5.3. Differential gene expression analysis was performed using DeSeq2 between the control and 4NQO-treated mouse tongues. Processed single-cell RNA-sequencing data for mouse esophageal squamous cell carcinoma generated by Wang et al. [41] were obtained from the GSA (CRA002118). Cells with less than 200 genes expressed and more than 10% mitochondrial genes expressed were filtered out to ensure that only high-quality cells were utilized in downstream analyses. Genes expressed in less than 3 cells were also filtered out of the analysis. After filtering, data were normalized and variable features were identified utilizing the selection method: “vst”. Cell clustering was performed using nearest neighbor construction and then clustered. Identified clusters were then annotated to identify known cell populations. For the analyses performed here, cells were filtered to retain only tumor cells, as identified through our annotation, and the distribution of p63 expression was evaluated to identify p63low and p63high clusters, and differential gene expression analysis was performed between the identified cell populations. Heatmaps for the scRNA-seq data were generated using the Seurat package [42].

2.10. Invasion Assay

B7E3 pIND-Trp63 cells were treated with 0, 50, or 200 ng of doxycycline (dox) for 24 h to induce p63 expression. Cells were plated at a density of 200,000 cells in the top of a Boyden chamber (Corning, product No. 354480) with a medium containing no serum and were allowed to invade and migrate for 12 h through a Matrigel matrix to the bottom chamber, which contained fetal bovine serum as an attractant. Cells were then fixed in formaldehyde and permeabilized with methanol according to a standard protocol. Cells on the bottom of the chamber were then stained with a mixture of Evan’s blue and methylene blue, and imaged with a Cytation 1 imaging plate reader (BioTek of Agilent, Santa Clara, CA); the stained area was calculated using ImageJ.

2.11. Spheroid Assay

B7E3 cells were grown in 96-well culture plates coated with 50% Matrigel (Corning, product No. 356237) for 6–8 days. To induce p63 expression, the cells were treated with 100 μg/mL doxycycline. The resulting spheroids were imaged using a 4× lens objective on a Cytation 1 imaging device; images were analyzed using ImageJ software. Graphing and statistical analysis were performed using Microsoft Excel (Microsoft Corp., Redmond, WA, USA). Student’s t-tests for samples with equal variance were used to determine statistical significance.

2.12. Motif Enrichment Analysis of Enhancers

The MEME Suite tool CentriMo was used for local motif enrichment analysis to determine which DNA-binding transcription factor motifs from the HOCOMOCO mouse (v11 FULL) database [43] were enriched within regions obtained from the p63 ChIP-seq in B7E11 cells. Motifs were ranked according to p-value.

2.13. Gene Ontology/Pathway Enrichment Analysis

The GREAT tool was used to annotate binding loci to the nearest gene [32]. Identified genes were then subjected to either the KEGG pathway analysis or WikiPathway analysis utilizing the DAVID functional annotation tool or the Broad Institute Gene Set Enrichment Analysis (GSEA) web tool [44,45,46,47,48].

2.14. Immunofluorescence and Immunohistochemistry

De-identified SCC patient samples were obtained from the archives of the Department of Oral and Maxillofacial Pathology, the University of New York at the University at Buffalo. An oral cavity tumor with a normal tissue microarray (US Biomax, Inc.—Tissue Array, T271b, Derwood, MD, USA) and SCC tissue slides were deparaffinized and sequentially rehydrated in decreasing concentrations of ethanol in water. After heat-induced antigen retrieval in sodium citrate, tissues were blocked and then incubated overnight with primary antibodies specific to p63 (4A4, 1:300), COTL1 (Proteintech, 1:200), or K14 (1:300), according to standard protocols. For immunofluorescence, the slides were incubated with appropriate fluorescent secondary antibodies, and coverslips were mounted with a mounting medium containing DAPI (Vector Laboratories, Newark, CA, USA) to stain nuclei, and then they were sealed for imaging. For immunohistochemistry, antibody labeling was visualized with an Impact DAB substrate kit (Vector Laboratories). Counterstaining was conducted using hematoxylin (Vector labs) before the slides were rinsed in tap water, air dried, and coverslipped with Permount mounting medium (Thermo Fisher Scientific, Waltham, MA, USA).

2.15. Statistics

Statistical analyses were performed using R. A Shapiro–Wilk test was performed to check the normality of data, and then either Student’s t-test or a Wilcoxon signed-rank test was performed according to whether the data were normally distributed. A p-value lower than 0.05 was considered significant.

3. Results

3.1. The Global Transcriptome of 4NQO-Generated mOSCC Cell Lines

To perform molecular studies of p63, we chose three tumorigenic murine oral SCC cell lines derived from BALB/c oral keratinocytes exposed to 4NQO [26]. These cell lines, B7E3, B4B8, and B7E11, showed low (undetectable), medium, and high p63 expression, respectively (Figure 1A), thus presenting excellent models to examine the molecular function of p63 in mOSCC. These cell lines are known to differ in their invasive abilities and the associated expression of epithelial–mesenchymal transition (EMT)-associated genes [49]. Thus, we first chose to characterize these cell lines by exploring the protein expression of EMT-associated factors through Western blotting. B7E3 and B7E11 cells had the highest and lowest levels of vimentin, respectively (Figure 1A), consistent with previous immunofluorescence findings [49]. B7E11 cells also had higher keratin levels than B4B8 and B7E3 cells (Figure 1A). B7E11 cells displayed high expression of epithelial EMT markers, whereas B7E3 cells showed high levels of mesenchymal-associated markers (Figure 1A). Interestingly, B4B8 cells appeared to display a hybrid EMT phenotype, with high levels of both epithelial- and mesenchymal-associated markers (Figure 1A). This partial EMT state is supported by previous findings showing that B4B8 cells migrate more than both B7E3 and B7E11 cells and also produce larger tumors than B7E11 cells [49].
First, to explore the variability in gene expression between these cell lines, we performed dimension reduction using principal component analysis (PCA) using the top 2500 most variable genes. PCA revealed that the most variation in gene expression was between B7E3 and B7E11 cells, with B4B8 cells lying in between, matching our Western blot results (Figure S1). Next, we performed RNA-seq of each cell line in order to explore their transcriptional landscape. Our differential gene expression analysis identified 9171 differentially expressed genes (DEGs) between B7E3 and B7E11 cells, and 8693 DEGs between B7E3 and B4B8 cells (Table S1). To integrate the gene expression differences across the three cell lines and to simplify the downstream analysis, we identified the genes showing consistent upregulation and downregulation in the low p63-expressing (p63low) and high p63-expressing (p63high) cell lines. This combined analysis resulted in 2588 genes enriched in B7E3 cells and 2376 genes enriched in both B4B8 and B7E11 cells (Figure 1B, Table S1). GSEA of the top 500 enriched genes from this analysis found that p63low B7E3 cells were enriched in genes associated with inflammatory response, whereas the p63high cells showed enrichment in EMT, hypoxia, and xenobiotic metabolism genes (Figure 1C,D). Together, these results revealed that the B4B8 and B7E11 cells share a similar gene expression profile that is distinct from that of the B7E3 mOSCC cells. Interestingly, this shared B4B8/B7E11 gene expression profile may be associated with the transcription factor p63, as evidenced by the enrichment of pathways shown to be regulated by p63, including the p53 pathway, Wnt beta-catenin signaling, and EMT [50].

3.2. Generating a p63 Cistrome in the 4NQO-Transformed B7E11 mOSCC Cell Line

To further characterize the p63–regulated oncogenic program in mOSCC cells, we next sought to identify target genes, whose regulatory regions are bound by p63. To this end, we performed independent ChIP-seq experiments in the p63high B7E11 cells utilizing two p63-specific antibodies, 4A4 and ΔNp63 E6Q30 (Figure 1A). The 4A4 antibody recognizes all p63 isoforms and identified 29,001 genomic sites, whereas the ΔNp63-specific E6Q3O antibody identified 6414 sites (Figure 2A, Table S2). Next, we used the intersectBed function in bedtools to identify the 6396 high-confidence p63-bound sites common to both antibodies (Table S2). As expected, a motif analysis performed on these high-confidence peaks revealed that the consensus p63 motif (p = 1.4 × 10−3985) was the most enriched, followed closely by the p73 (p = 2.5 × 10−3815) and the p53 (p = 4.8 × 10−2832) motifs (Figure 2B). The distribution of p63 peaks relative to the transcriptional start site indicates that p63 preferentially targets distal intergenic and intronic regions, which are likely to act as enhancers (Figure 2C).
To establish a p63 gene regulatory network, we integrated the p63 cistrome with the cell line DEG datasets. This resulted in the identification of 1897 direct p63 targets and 3067 indirect targets (Figure 2D and Table S2). For the direct targets, we found an enrichment of several oncogenic programs, including vascular endothelial growth factor (VEGF) signaling, focal adhesion, and TGF-β signaling in the p63high B7E11 cells (Figure 2E). Conversely, both p63low B7E3 and B4B8 cells were enriched for biological processes regulated by p63, such as cytodifferentiation, phosphoinositide-3 kinase pathway, EMT, and Notch1 signaling (Figure 2F).

3.3. Generating a p63-Driven mOSCC Transcriptomic Signature

To examine molecular processes regulated by p63 in B7E11 cells, we next performed knockdown studies. For these experiments, we used three stable, independent lentiviral shRNAs targeting Trp63. Western blotting confirmed that all three shRNAs (sh1, sh2, and sh3) reduced p63 expression, albeit at varying levels (Figure 3A). We performed RNA-seq on sh2 and sh3 cells, which affected the most robust p63 depletion. sh2 resulted in 100 DEGs, while sh3 resulted in 587 DEGs (Table S3). We focused our analyses on sh3 as it produced the strongest knockdown of p63 at both the mRNA and protein levels.
Among the DEGs resulting from sh3-mediated p63 knockdown, 287 genes were upregulated and 300 were downregulated (Figure 3B,C). An analysis of enriched KEGG pathways using the DAVID bioinformatics database of these DEGs revealed cell cycle, epithelial–mesenchymal transition, apical junction, and the p53 signaling pathway among genes that were downregulated (Figure 3D). Conversely, upregulated DEGs were enriched in immune system-related pathways, including inflammatory responses and TNF signaling via NF-κB (Figure 3D). The upregulation of many cytokines and chemokines upon the loss of p63 is interesting because dysregulation of p63 has been associated with altered immune responses in SCC and other cancers [14,30,51,52,53].
To complement our p63 knockdown studies, we performed gain-of-function studies in the p63low B7E3 cell line. For this purpose, we established a lentiviral-based doxycycline-inducible system to drive the expression of a FLAG epitope-tagged version of ΔNp63α. First, we confirmed that doxycycline (dox) induces ΔNp63α overexpression in a dose-dependent manner (Figure 4A). Next, we carried out RNA-seq experiments on cells in the presence and absence of dox. Differential gene expression analysis identified 147 DEGs (Figure 4B,C, Table S4). To explore the pathways enriched in the 64 downregulated DEGs from the B7E3 p63-overexpressing cells, we performed GSEA using the hallmark gene set. The top enriched pathways in the downregulated genes were involved in IL2 signaling, hypoxia, and myogenesis (Figure 4D). Conversely, upregulated genes were enriched in p53 signaling, inflammatory responses, and EMT (Figure 4E). To validate the results from the RNA-seq results, we performed qRT-PCR for several of these genes, which confirmed the trends observed in both knockdown and overexpression experiments, indicating that p63 regulates the expression of these genes (Figures S2 and S3).
To identify a consensus set of genes that are most reliant on p63 expression, we combined the DEGs from both the loss- and gain-of-function experiments, revealing a set of 37 genes common to both datasets (Table S5). To further identify genes that were directly regulated by p63, we incorporated our p63 ChIP-seq data, which revealed 22 genes directly bound by p63 in B7E11 cells (Table S5). These analyses provided a p63-driven gene expression signature for 4NQO-induced mOSCC.

3.4. Evaluating the p63-Driven Gene Signature in an Independent 4NQO-Induced mOSCC Model

We explored if our p63-driven signature is also enriched in an independent RNA-seq dataset generated by Tang et al. [11] from tongue tumors in 6-week-old C57BL/6 mice treated with 100 μg/mL of 4NQO for 10 weeks. We reprocessed this data to perform DEG analysis between data from the 4NQO-treated and control tongues (mice treated with the vehicle) and identified 8171 DEGs (4014 upregulated and 4157 downregulated) at an adjusted p-value of ≤0.1 between the two conditions (Figure 5A and Table S6). We focused on DEGs showing a log2 fold change of >1 to identify enriched pathways (Figure 5B). As expected, we found that Trp63 was upregulated, in agreement with observations in other 4NQO studies (Table S6). Overall, 2588 DEGs were downregulated upon tumor induction by 4NQO, and these were enriched for metabolic pathways, salivary secretion, and focal adhesion (Figure 5C). The enrichment of these pathways among downregulated genes matches that observed upon induction of p63 expression, suggesting that p63 plays a key role in this downregulation and that these pathways are key to 4NQO-induced tumorigenesis. By contrast, 2274 genes were upregulated in 4NQO-induced tumors (Table S6). These upregulated DEGs displayed enrichment of pathways involved in cytokine–cytokine receptor interaction, chemokine signaling, cell cycle, and p53 signaling (Figure 5D). These pathways are similar to those identified for the p63-induced DEGs, highlighting the potentially important role of p63-signaling in 4NQO-induced tumors, particularly in regulating the tumor immune microenvironment.
We compared our p63-based signature with the list of DEGs from the tongue tumor dataset by Tang et al. [11] and identified 22 genes that were common to both (Figure 6A, Table S7). These genes are involved in the p53 signaling pathway, cell adhesion pathways, and metabolic pathways (Figure 6A). Of the 22 genes, 14 are direct targets of p63 according to our high-confidence p63 target list (Figure 6A and Table S7). Having ascertained a core signature for p63 signaling in murine OSCC, we turned our attention to human OSCC to determine if any of our identified targets may be conserved in the human disease context. In addition to human orthologs previously identified as playing key roles in OSCC, p63 signaling, and other cancers, namely, FAT2, K14, and PERP, we identified several novel regulators, including Cotl1, Bcam, Adipor2, and Wnt7b, (Figure 6A) [54,55,56,57,58]. Finally, to determine if the gene expression changes translate to differences in protein levels, we performed Western blotting on a subset of proteins encoded by DEGs. Overexpression of p63 in B7E3 cells resulted in an increase in the levels of Cotl1, Krt14, and Krt6 (Figure 6B), similar to what was seen in the RNA-seq data and demonstrating that the observed gene expression changes result in alterations to protein expression (Figure 6B).

3.5. scRNA-seq Analysis of 4NQO-Induced Mouse ESCC

In the absence of publicly available scRNA-seq data for 4NQO-generated OSCC, we turned our attention to esophageal squamous cell carcinoma (ESCC) [59]. Treatment of mouse esophageal tissues with 4NQO induces SCC in a manner that mimics the tumorigenic processes of ESCC in humans, similar to a 4NQO-induced mouse OSCC [59]. Therefore, we also employed a recently published dataset from scRNA-seq of a 4NQO-induced ESCC generated by Yao et al. [41] to explore the possible enrichment of our identified p63-based gene expression signature. First, we reanalyzed the ESCC scRNA-seq data generated by Yao et al. [41] and performed clustering analysis to identify nine overarching cell-type clusters (Figure 7A). We then assessed p63 expression across the tumor cell cluster and found a gradient of p63 expression across this cluster, with a clear population of cells with high expression and a population of cells with low expression (Figure 7B). Next, we performed DEG analysis on the identified p63low cell population and p63high population and found 480 genes, five of which overlapped our p63-based gene expression signature (Figure 7B, Table S8). We also performed GSEA using the hallmark set of genes enriched in either the p63low or p63high cell population (Figure 7C). GSEA data from our p63 knockdown study markedly overlapped that from the p63low population, sharing enrichment of KRAS signaling and interferon gamma response pathways (Figure 3E and Figure 7C). Likewise, in the p63high cell population, the hallmark pathways overlap those in the cells with induced p63 overexpression, such as EMT, TNF signaling via NF-κB, and p53 signaling (Figure 4E and Figure 7D).

3.6. p63-Driven mOSCC Gene Signature Enriched in the Human TCGA Dataset

The identification of a p63 gene expression signature in mOSCC prompted us to examine its relevance in human tumors. To this end, we next explored the extent of cross-species conservation of the p63-driven signature in the TCGA-HNSCC datasets, keeping in mind the underlying differences in disease etiology and genetic complexity between these two cohorts. We compared the expression of 22 of the signature genes in the mouse dataset to that of the TCGA-HNSCC data. This analysis identified 16 genes that showed significant differences in expression between normal and tumor tissue (Figure S4). We expected six of these to be downregulated in tumors compared to expression in normal tissues; however, only two genes, CES2 and EMP1, matched the pattern of expression identified in the mOSCC datasets utilized. (Figure S4). Similarly, 10 of the 15 genes upregulated in our signature, including FAT2, COTL1, and KRT14, were also upregulated in human tissues (Figure S4).
To determine the prognostic value of our signature, we performed Kaplan–Meier analyses for each of the 21 genes in our p63-driven signature (Figure S5). Only two genes were associated with differences in overall patient survival: higher expression of BCAM and WNT7B were associated with worse survival (Figure S5). A few other genes in our signature, such as AMOTL2 and ANXA3, displayed trends toward a worse prognosis at higher levels, but these were not statistically significant (Figure S5). Overall, our analysis identified known and novel players in the p63 network that may be further investigated as drivers of the OSCC oncogenic program in both mice and humans.

3.7. p63 Affects mOSCC Cell Line Migration and Proliferation

p63 has been reported to regulate the proliferation and metastasis of human HNSCC cells, but its effect on mouse OSCC has not been explored [53,60,61,62]. Our previous pathway analysis revealed that B7E3 cells, which have undetectable levels of p63, are enriched for genes associated with processes such as cell adhesion and the cell cycle (compared to that of cells with high levels of p63). Thus, we examined if induced p63 expression influences the migratory and proliferative potential of B7E3 cells. We performed cell migration assays using a Transwell system and found that dox-induced p63 expression impaired the migration of cells compared to that of control cells (no dox) (Figure 8A), suggesting that p63 blocks migration. This finding is consistent with the downregulation of p63 expression in SCCs prior to cell migration and metastasis [25,63].
Next, to assess the effect of p63 expression on the proliferation of B7E3 cells in an anchorage-independent environment, we performed spheroid growth assays. B7E3 cells were grown either under control conditions with no dox or with 100 ng of dox to induce p63 overexpression. After 9 days, spheroids formed by p63-expressing B7E3 cells were significantly larger on average than spheroids grown from control cells (Figure 8B), confirming that p63 expression increases the proliferative ability of cells, including mOSCC cell lines.

3.8. COTL1 Is a Novel Target of p63 and Is Driven by p63 Expression

The enrichment of genes for several cell adhesion molecules in our combined p63 signature prompted us to look more closely at a novel p63 target, coactosin-like protein (Cotl1). COTL1 belongs to the ADF/cofilin family and likely binds to actin filaments to regulate the cytoskeleton [64]. COTL1 has been found to regulate the migration of mouse neocortical neurons and potentially plays a role in T-cell activation [65]. COTL1 also plays a role in the progression of cancer, promoting the proliferation of lung cancer and glioblastoma cells in vitro and in vivo [66,67]. COTL1 mutations are one of the most prevalent subclonal mutations in human OSCC, with increased copy numbers in 10% of OSCC cases [68]. We assessed COTL1 expression in B7E3 cells by immunofluorescence staining using antibodies specific to p63 and COTL1 in B7E3 cells with and without dox-induced p63 expression. Although previous studies showed COTL1 expression to be exclusively cytoplasmic, we observed strong expression in nuclei (Figure 9). The nuclear staining for COTL1 overlapped that of DAPI in areas with low DNA density and remained in discreet aggregates within nuclei under all conditions, except in the case of actively dividing cells, in which it was more diffuse and surrounded the condensed DNA (Figure S6). The intensity of cytoplasmic staining of COTL1 matched that of p63, with stronger staining of cytoplasmic COTL1 with greater p63 induction (Figure 9). These findings are consistent with both our RNA-seq and Western blot results, showing that COTL1 expression is regulated by p63 expression.
To determine if these patterns of COTL1 and p63 expression in mOSCC cells occur in human OSCC tissues, we performed immunofluorescence staining of SCC tissues representing well-differentiated, moderately differentiated, and poorly differentiated SCC, because previous work has shown that p63 expression increases and expands as SCC becomes less differentiated [20]. Accordingly, we observed that p63 showed increased expression as well as expanded coverage of expression across the tissue from well-differentiated tissues to poorly differentiated tissues (Figure 10). We saw a similar pattern of expression for COTL1. Interestingly, COTL1 expression was tightly defined around endothelial tissue, suggesting the cytoplasmic COTL1 may play a role in maintaining basement membrane integrity (Figure 10).
Finally, we performed immunohistochemistry on a tissue microarray of normal and cancer human oral tissues, staining for p63, COTL1, and K14 (Figure 11). Similar to the results shown in Figure 10, we observed overall increased p63 staining in tumors compared to that in normal tissues. Notably, we found a stage-dependent increase in p63 expression, with advanced-stage tumors expressing the highest levels of p63 regardless of the anatomical site of the tumor (Figure 11). Similarly, the expression of K14, a known target of p63, was increased in cancer compared to that in normal tissue. Consistent with the results obtained from the previous immunofluorescence experiments, the pattern of COTL1 staining matched that of p63, with increased expression in cancer tissue (Figure 11). Taken together, these staining results confirm the gene expression-based findings of 4NQO-derived cell lines and tumors and further support the notion of an oncogenic network of p63 and its targets that operate in oral tumors.

4. Discussion

Tobacco smoking claims the lives of more than 6 million people every year worldwide and is one of the leading causes of cancer deaths in the United States [69,70]. The compound 4NQO is a precursor carcinogen that mimics some effects of tobacco, thus 4NQO-mediated chemical carcinogenesis models serve as valuable tools for the mechanistic exploration of tobacco-associated cancers, such as OSCC [71]. The molecular characterization of three representative 4NQO-induced OSCC cell lines through RNA-seq reported here provides a valuable resource and complements similar studies of the well-established MOC1 cell line [72,73]. One interesting observation from our studies is the difference in epithelial and mesenchymal gene expression patterns across the three mOSCC cell lines, in agreement with a prominent role for EMT in OSCC [49]. Specifically, we note that the B4B8 cell line displayed characteristics of a partial, or hybrid EMT state, as markers associated with both epithelial and mesenchymal states were highly expressed at the mRNA and protein levels. The hybrid EMT state of B4B8 cells may also predispose them to be particularly invasive, as shown previously [49], and might also be relevant in other contexts; for example, where it has been used as a murine HNSCC model to explore EGFR/ERBB-dependent growth. [74]
The p63-null B7E3 and the p63high B7E11 cells offered us an excellent toolkit to study the role of the oncogenic ΔNp63 isoform in mOSCC and identify its transcriptional targets. Although previous studies have suggested p63 as a potential therapeutic target that is upregulated upon the deletion of the commonly mutated p16INK4a locus in 4NQO-induced OSCC, the molecular and cellular processes that are regulated by p63 in this context have not been identified [12,24]. Our RNA-seq and ChIP-seq experiments in the p63high B7E11 led us to the identification of genes regulated directly and indirectly by p63 and to a consensus p63-dependent core gene signature. By incorporating results from other independent and complementary datasets, we also examined the relevance of this mouse 4NQO-derived OSCC signature to human tumors [11,41]. The resultant signature includes genes well-known in OSCC and p63-dependent signaling networks, as well as several novel factors that shed light on the specific role of p63 in OSCC and even more broadly in SCC.
It is important to highlight that, unlike the genetically well-defined 4NQO models, human tumors result from multiple carcinogenic insults. Thus, one caveat is that cross-species comparisons of molecular signatures of OSCC are likely to reveal both similarities and differences. This is exemplified by Bcam, which encodes a cell adhesion protein that acts as a receptor for LAMA5, laminin that is a major component of the basement membrane [75]. BCAM is highly expressed in KRAS-mutant hepatic metastases from colorectal cancer, and inhibition of BCAM/LAMA5 interferes with the adhesion of colorectal cells to vascular endothelial cells, thereby reducing metastatic growth [76]. Intriguingly, while BCAM expression is increased in human HNSCC tumors, we found that Bcam expression is reduced in 4NQO-induced tumors, highlighting the differences between human and murine tumors. We observed higher Bcam expression in the p63high cell population compared to the p63low cells in the scRNA-seq generated by Yao et al. [41], further supporting the notion that this gene is regulated by p63. Since survival analyses suggest that BCAM expression does affect the outcome of human HNSCC, this might be a p63 target worthy of follow-up investigations. In addition, we also found significant enrichment of genes involved in the p53 signaling pathway, some of which have known interactions with p63, such as Perp and Sfn [77,78]. Interestingly, our analysis showed an upregulation of Perp in 4NQO-induced tongue tumors even though Perp is proapoptotic, suggesting a complicated balance of proliferative and apoptotic signaling [78]. Similarly, Sfn (Stratifin, also known as 14-3-3 protein sigma) was also increased in 4NQO-induced tumors, matching results from prior studies, which showed Sfn to be over-expressed in human HNSCC tumor samples in relation to non-cancerous head and neck tissues [79,80].
One of our notable findings was that the expression of several promigratory genes, including Fat2 and Cotl1, paralleled that of p63 in 4NQO-induced tumors. p63 is known to induce the expression of human FAT2, as well as the mesenchymal gene Slug, to promote tumor invasion in breast cancer, whereas COTL1 increases the migratory ability of both breast and non-small cell lung cancer cells [56,81,82]. Interestingly, the p63 gene signature we identified also includes genes associated with invasion, such as AMOTL2 that are paradoxically suppressed by p63. In humans, AMOTL2 has conflicting roles in migration that seemingly depend on the cellular context, such that knockdown inhibits the migration of human umbilical vein endothelial cells but promotes EMT of mammary epithelial cells [83]. The role of Amotl2 in mOSCC may more closely align with that in mammary epithelial cells, but that remains an area for further exploration. Similarly, we also observed that high p63 expression was associated with the downregulation of Msln, which encodes mesothelin, a protein that promotes EMT and invasion in breast, lung, and ovarian cancers [84,85,86]. A similar dichotomy of p63 function in cell migration and invasion has been reported before and has been attributed to an oscillatory expression of the ΔNp63 isoform that results in its tumor suppressive activities [25].
Finally, our exploration of the novel p63 target, COTL1 in the mouse cell lines and human HNSCC, suggests that it might have an oncogenic role that is conserved in both mouse and human disease. Notably, COTL1 and p63 have similar expression profiles in OSCC tissues, and higher expression of both was observed in more advanced disease stages, suggesting they play a role in the clinical outcome of OSCC patients. In support of this notion, COTL1 levels were amplified in up to 10% of patients with OSCC tumors in one study [68]. Previous immunofluorescence staining for COTL1 has shown diffuse localization in the cytoplasm; however, we observed both nuclear localization, as well as cytoplasmic staining of COTL1, which suggests an intriguing nuclear function for COTL1 [87]. The discrete aggregates of COTL1 in the nucleus, which seemingly disappear during cell division, hint at a possible role for COTL1 in cell division. Although the specific molecular function of COTL1 in OSCC requires additional studies, our results described here, combined with those previously reported, suggest COTL1 is involved in OSCC metastasis and patient survival. Additional support for this notion comes from a recent study identifying COTL1 as one of the 52 differentially expressed proteins in a proteomic-based comparison of metastatic and non-metastatic lymph nodes from HNSCC patients [88].
Overall, our studies reaffirm the previously identified role for p63 in modulating cellular processes that underlie tumor progression and metastasis, such as cell adhesion in OSCC, similar to what has been reported in other SCC and in epithelial development [16,89]. The p63-driven mOSCC signature and data generated from the 4NQO mOSCC cell lines set the stage for future studies into p63 function in OSCC and for the identification of actionable targets for future therapeutics. One major current limitation in this field is the lack of robust genomics datasets, in particular for 4NQO-derived tumors at different stages of tumor development and progression. Future transcriptomic and epigenomic studies addressing this shortcoming, particularly at single-cell resolutions, will offer valuable insights into the p63 gene regulatory and signaling networks and their role in the etiology of OSCC.

5. Conclusions

The dynamic regulation of p63 and its vast network of targets in both mouse 4NQO-derived cell lines and human tissues highlights its importance to OSCC. Our data re-emphasizes the crucial role of p63, specifically the ΔNp63 isoform in regulating networks of known and yet-to-be-discovered genes and a myriad of oncogenic processes, such as cell adhesion, migration, and metastasis that are important in OSCC. This is best illustrated by the identification of COTL1 as a novel factor regulated by p63 that can affect the clinical outcome of OSCC patients. How COTL1 functions in the cellular context and how its biological actions affect tumor cell behavior remains largely unknown. This remains the focus of mechanistic follow-up studies that are likely to be very informative. One limitation of our genomic studies is that they are based on cells grown in culture and thus likely do not represent an in vivo tumor. In this regard, syngeneic mouse models bearing tumors generated from the 4NQO-derived parental and modified cell lines can be extremely valuable to further validate our findings, especially the p63 gene signature. Taken together, the results from this study add to our understanding of the underlying molecular alterations associated with OSCC and provide potential targets for better diagnostics and therapy.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cancers15020446/s1. Table S1: Combined Table showing results of mOSCC cell line DEG analyses; Table S2: B7E11 p63 ChIP-seq Combined Sites; Table S3: B7E11 p63 knockdown DEG list; Table S4: B7E11 p63 inducible DEG list; Table S5: Combined loss- and gain-of-function p63-based mOSCC Signature; Table S6: 4NQO-induced mOSCC DEG analysis; Table S7: Combined overall p63-based gene expression signature; Table S8: scRNA-seq DEG analysis based on p63high versus p63low cell Populations; Table S9: List of primers utilized in qRT-PCR analysis; Figure S1: Principal component analysis (PCA) of mouse 4NQO-induced oral squamous cell carcinoma cell lines. Top 2500 variable genes were used to generate PCA; Figure S2: qRT-PCR analysis of a selection of DEGs identified in B7E11 p63 knockdown cell line; Figure S3: qRT-PCR analysis of a selection of DEGs identified in B7E3 p63-inducible cell line; Figure S4: Boxplot of expression levels of the genes that constitute the combined p63-based mOSCC gene expression signature in the human TCGA HNSCC cohort; Figure S5: Kaplan-Meier survival analysis of the genes that constitute the combined p63-based mOSCC gene expression signature in the human TCGA HNSCC cohort; Figure S6: Immunofluorescent staining of COTL1 shows novel nuclear staining in areas of low DNA density, except for actively dividing cells in which diffused COTL1 expression is observed surrounding the condensed DNA. Yellow arrows highlight the punctual staining of COTL1, while white arrows highlight the diffuse staining of COTL1; Figure S7: Full Western blots of parental mouse 4NQO-induced cell lines B7E3, B4B8, and B7E11 as seen in Figure 1; Figure S8: Full Western blots of B7E11 p63 knockdown cell line as seen in Figure 3; Figure S9: Full Western blots of B7E3 inducible p63 cell lines as seen in Figure 4 and Figure 6.

Author Contributions

Conceptualization, A.R.G. and S.S.; Investigation, A.R.G., A.O., K.B.N. and J.S.; Visualization, A.R.G. and A.O.; Formal Analysis, A.R.G. and A.O.; Writing—original draft, A.R.G.; Writing—review and editing, A.R.G., R.-A.R. and S.S.; Supervision, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

A.O. acknowledges support from National Institute of Health, grants DE028480 and DE023526.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets presented in this study are available at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE217670 (accessed on 1 December 2022).

Acknowledgments

We thank Kirsten Smalley for excellent technical help. We also thank the UB Genomics and Bioinformatics Core for their next-generation sequencing services.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leemans, C.R.; Braakhuis, B.J.M.; Brakenhoff, R.H. The Molecular Biology of Head and Neck Cancer. Nat. Rev. Cancer 2011, 11, 9–22. [Google Scholar] [CrossRef] [PubMed]
  2. Johnson, D.E.; Burtness, B.; Leemans, C.R.; Lui, V.W.Y.; Bauman, J.E.; Grandis, J.R. Head and Neck Squamous Cell Carcinoma. Nat. Rev. Dis. Prim. 2020, 6, 92. [Google Scholar] [CrossRef] [PubMed]
  3. Howard, A.; Agrawal, N.; Gooi, Z. Lip and Oral Cavity Squamous Cell Carcinoma. Hematol. Oncol. Clin. North Am. 2021, 35, 895–911. [Google Scholar] [CrossRef]
  4. Lawrence, M.S.; Sougnez, C.; Lichtenstein, L.; Cibulskis, K.; Lander, E.; Gabriel, S.B.; Getz, G.; Ally, A.; Balasundaram, M.; Birol, I.; et al. Comprehensive Genomic Characterization of Head and Neck Squamous Cell Carcinomas. Nature 2015, 517, 576–582. [Google Scholar] [CrossRef] [Green Version]
  5. Puram, S.V.; Tirosh, I.; Parikh, A.S.; Patel, A.P.; Yizhak, K.; Gillespie, S.; Rodman, C.; Luo, C.L.; Mroz, E.A.; Emerick, K.S.; et al. Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer. Cell 2017, 171, 1611–1624.e24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kürten, C.H.L.; Kulkarni, A.; Cillo, A.R.; Santos, P.M.; Roble, A.K.; Onkar, S.; Reeder, C.; Lang, S.; Chen, X.; Duvvuri, U.; et al. Investigating Immune and Non-Immune Cell Interactions in Head and Neck Tumors by Single-Cell RNA Sequencing. Nat. Commun. 2021, 12, 7338. [Google Scholar] [CrossRef]
  7. Vered, M.; Yarom, N.; Dayan, D. 4NQO Oral Carcinogenesis: Animal Models, Molecular Markers and Future Expectations. Oral Oncol. 2005, 41, 337–339. [Google Scholar] [CrossRef]
  8. Kanojia, D.; Vaidya, M.M. 4-Nitroquinoline-1-Oxide Induced Experimental Oral Carcinogenesis. Oral Oncol. 2006, 42, 655–667. [Google Scholar] [CrossRef]
  9. Wang, Z.; Wu, V.H.; Allevato, M.M.; Gilardi, M.; He, Y.; Luis Callejas-Valera, J.; Vitale-Cross, L.; Martin, D.; Amornphimoltham, P.; Mcdermott, J.; et al. Syngeneic Animal Models of Tobacco-Associated Oral Cancer Reveal the Activity of in Situ Anti-CTLA-4. Nat. Commun. 2019, 10, 5546. [Google Scholar] [CrossRef] [Green Version]
  10. Kondo, S. A Test for Mutation Theory of Cancer: Carcinogenesis by Misrepair of Dna Damaged by 4-Nitroquinoline 1-Oxide. Br. J. Cancer 1977, 35, 595–601. [Google Scholar] [CrossRef]
  11. Tang, X.H.; Osei-Sarfo, K.; Urvalek, A.M.; Zhang, T.; Scognamiglio, T.; Gudas, L.J. Combination of Bexarotene and the Retinoid CD1530 Reduces Murine Oral-Cavity Carcinogenesis Induced by the Carcinogen 4-Nitroquinoline 1-Oxide. Proc. Natl. Acad. Sci. USA 2014, 111, 8907–8912. [Google Scholar] [CrossRef] [Green Version]
  12. Ishida, K.; Tomita, H.; Kanayama, T.; Noguchi, K.; Niwa, A.; Kawaguchi, M.; Miyai, M.; Matsuo, M.; Imaizumi, Y.; Kato, K.; et al. Specific Deletion of P16INK4a with Retention of P19ARF Enhances the Development of Invasive Oral Squamous Cell Carcinoma. Am. J. Pathol. 2020, 190, 1332–1342. [Google Scholar] [CrossRef]
  13. Chen, W.; Kang, K.L.; Alshaikh, A.; Varma, S.; Lin, Y.; Shin, K.; Kim, R.; Wang, C.; Park, N.; Walentin, K.; et al. Grainyhead-like 2 (GRHL2) Knockout Abolishes Oral Cancer Development through Reciprocal Regulation of the MAP Kinase and TGF-β Signaling Pathways. Oncogenesis 2018, 7, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Du, J.; Romano, R.A.; Si, H.; Mattox, A.; Bian, Y.; Yang, X.; Sinha, S.; Van Waes, C.; Chen, Z. Epidermal Overexpression of Transgenic Δnp63 Promotes Type 2 Immune and Myeloid Inflammatory Responses and Hyperplasia via NF-ΚB Activation. J. Pathol. 2014, 232, 356–368. [Google Scholar] [CrossRef] [PubMed]
  15. Sethi, I.; Gluck, C.; Zhou, H.; Buck, M.J.; Sinha, S. Evolutionary Re-Wiring of P63 and the Epigenomic Regulatory Landscape in Keratinocytes and Its Potential Implications on Species-Specific Gene Expression and Phenotypes. Nucleic Acids Res. 2017, 45, 8208–8224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Romano, R.; Smalley, K.; Magraw, C.; Serna, V.A.; Kurita, T.; Raghavan, S.; Sinha, S. DNp63 Knockout Mice Reveal Its Indispensable Role as a Master Regulator of Epithelial Development and Differentiation. Development 2012, 139, 772–782. [Google Scholar] [CrossRef] [Green Version]
  17. Keyes, W.M.; Pecoraro, M.; Aranda, V.; Vernersson-Lindahl, E.; Li, W.; Vogel, H.; Guo, X.; Garcia, E.L.; Michurina, T.V.; Enikolopov, G.; et al. ΔNp63α Is an Oncogene That Targets Chromatin Remodeler Lsh to Drive Skin Stem Cell Proliferation and Tumorigenesis. Cell Stem Cell 2011, 8, 164–176. [Google Scholar] [CrossRef] [Green Version]
  18. Botchkarev, V.; Flores, E. P53/P63/P73 in the Epidermis in Health and Disease. Cold Spring Harb. Perspect. Med. 2014, 4, 1260–1266. [Google Scholar] [CrossRef]
  19. Romano, R.; Ortt, K.; Birkaya, B.; Smalley, K.; Sinha, S. An Active Role of the D N Isoform of P63 in Regulating Basal Keratin Genes K5 and K14 and Directing Epidermal Cell Fate. PLoS ONE 2009, 4, e5623. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Lo Muzio, L.; Santarelli, A.; Caltabiano, R.; Rubini, C.; Pieramici, T.; Trevisiol, L.; Carinci, F.; Leonardi, R.; De Lillo, A.; Lanzafame, S.; et al. P63 Overexpression Associates With Poor Prognosis in Head and Neck Squamous Cell Carcinoma. Hum. Pathol. 2005, 36, 187–194. [Google Scholar] [CrossRef]
  21. Somerville, T.D.D.; Xu, Y.; Miyabayashi, K.; Tiriac, H.; Cleary, C.R.; Maia-Silva, D.; Milazzo, J.P.; Tuveson, D.A.; Vakoc, C.R. TP63-Mediated Enhancer Reprogramming Drives the Squamous Subtype of Pancreatic Ductal Adenocarcinoma. Cell Rep. 2018, 25, 1741–1755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Hsieh, M.-H.; Choe, J.H.; Gadhvi, J.; Kim, Y.J.; Arguez, M.A.; Palmer, M.; Gerold, H.; Nowak, C.; Do, H.; Mazambani, S.; et al. P63 and SOX2 Dictate Glucose Reliance and Metabolic Vulnerabilities in Squamous Cell Carcinomas. Cell Rep. 2019, 28, 1860–1878.e9. [Google Scholar] [CrossRef] [Green Version]
  23. Sethi, I.; Romano, R.A.; Gluck, C.; Smalley, K.; Vojtesek, B.; Buck, M.J.; Sinha, S. A Global Analysis of the Complex Landscape of Isoforms and Regulatory Networks of P63 in Human Cells and Tissues. BMC Genomics 2015, 16, 584. [Google Scholar] [CrossRef] [Green Version]
  24. Lakshmanachetty, S.; Balaiya, V.; High, W.A.; Koster, M.I. Loss of TP63 Promotes the Metastasis of Head and Neck Squamous Cell Carcinoma by Activating MAPK and STAT3 Signaling. Mol. Cancer Res. 2019, 17, 1279–1293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bui, N.H.B.; Napoli, M.; Davis, A.J.; Abbas, H.A.; Rajapakshe, K.; Coarfa, C.; Flores, E.R. Spatiotemporal Regulation of ΔNp63 by TGFβ-Regulated MiRNAs Is Essential for Cancer Metastasis. Cancer Res. 2020, 80, 2833–2847. [Google Scholar] [CrossRef] [Green Version]
  26. Thomas, G.R.; Chen, Z.; Oechli, M.N.; Hendler, F.J.; Van Waes, C. Decreased Expression of CD80 Is a Marker for Increased Tumorigenicity in a New Murine Model of Oral Squamous-Cell Carcinoma. Int. J. Cancer 1999, 82, 377–384. [Google Scholar] [CrossRef]
  27. Oyelakin, A.; Nayak, K.B.; Glathar, A.R.; Gluck, C.; Wrynn, T.; Tugores, A.; Romano, R.; Sinha, S. EHF Is a Novel Regulator of Cellular Redox Metabolism and Predicts Patient Prognosis in HNSCC. NAR cancer 2022, 4, zcac017. [Google Scholar] [CrossRef]
  28. Gluck, C.; Glathar, A.; Tsompana, M.; Nowak, N.; Garrett-Sinha, L.A.; Buck, M.J.; Sinha, S. Molecular Dissection of the Oncogenic Role of ETS1 in the Mesenchymal Subtypes of Head and Neck Squamous Cell Carcinoma. PLoS Genet. 2019, 15, e1008250. [Google Scholar] [CrossRef] [Green Version]
  29. Rizzo, J.M.; Oyelakin, A.; Min, S.; Smalley, K.; Bard, J.; Luo, W.; Nyquist, J.; Guttman-Yassky, E.; Yoshida, T.; De Benedetto, A.; et al. ΔNp63 Regulates IL-33 and IL-31 Signaling in Atopic Dermatitis. Cell Death Differ. 2016, 23, 1073–1085. [Google Scholar] [CrossRef]
  30. Glathar, A.R.; Oyelakin, A.; Gluck, C.; Bard, J.; Sinha, S. P63 Directs Subtype-Specific Gene Expression in HPV+ Head and Neck Squamous Cell Carcinoma. Front. Oncol. 2022, 12, 879054. [Google Scholar] [CrossRef]
  31. Zhang, Y.; Liu, T.; Meyer, C.A.; Eeckhoute, J.; Johnson, D.S.; Bernstein, B.E.; Nusbaum, C.; Myers, R.M.; Brown, M.; Li, W.; et al. Model-Based Analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9, R137. [Google Scholar] [CrossRef] [Green Version]
  32. Mclean, C.Y.; Bristor, D.; Hiller, M.; Clarke, S.L.; Schaar, B.T.; Wenger, A.M.; Bejerano, G. GREAT Improves Functional Interpretation of Cis-Regulatory Regions. Nat Biotechnol. 2010, 28, 495–501. [Google Scholar] [CrossRef] [Green Version]
  33. Ramírez, F.; Ryan, D.P.; Grüning, B.; Bhardwaj, V.; Kilpert, F.; Richter, A.S.; Heyne, S.; Dündar, F.; Manke, T. DeepTools2: A next Generation Web Server for Deep-Sequencing Data Analysis. Nucleic Acids Res. 2016, 44, W160–W165. [Google Scholar] [CrossRef] [PubMed]
  34. Yu, G.; Wang, L.G.; He, Q.Y. ChIPseeker: An R/Bioconductor Package for ChIP Peak Annotation, Comparison and Visualization. Bioinformatics 2015, 31, 2382–2383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-Based Genome Alignment and Genotyping with HISAT2 and HISAT-Genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef]
  36. Wagner, G.P.; Kin, K.; Lynch, V.J. Measurement of MRNA Abundance Using RNA-Seq Data: RPKM Measure Is Inconsistent among Samples. Theory Biosci. 2012, 131, 281–285. [Google Scholar] [CrossRef]
  37. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E. Yevgeniy Antipin; Boris Reva; Arthur P. Goldberg; Chris Sander; Nikolaus Schultz. The CBio Cancer Genomics. Cancer Discov. 2017, 32, 736–740. [Google Scholar] [CrossRef]
  39. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the CBioPortal. Sci. Signal. 2014, 6, pl1. [Google Scholar] [CrossRef] [Green Version]
  40. Goldman, M.J.; Craft, B.; Hastie, M.; Repecka, K.; McDade, F.; Kamath, A.; Banerjee, A.; Luo, Y.; Rogers, D.; Brooks, A.N.; et al. Overcoming Barriers to Early Disease Intervention. Nat. Biotechnol. 2020, 38, 675–678. [Google Scholar] [CrossRef] [PubMed]
  41. Yao, J.; Cui, Q.; Fan, W.; Ma, Y.; Chen, Y.; Liu, T.; Zhang, X.; Xi, Y.; Wang, C.; Peng, L.; et al. Single-Cell Transcriptomic Analysis in a Mouse Model Deciphers Cell Transition States in the Multistep Development of Esophageal Cancer. Nat. Commun. 2020, 11, 3715. [Google Scholar] [CrossRef] [PubMed]
  42. Hao, Y.; Hao, S.; Andersen-Nissen, E.; Mauck, W.M.; Zheng, S.; Butler, A.; Lee, M.J.; Wilk, A.J.; Darby, C.; Zager, M.; et al. Integrated Analysis of Multimodal Single-Cell Data. Cell 2021, 184, 3573–3587.e29. [Google Scholar] [CrossRef]
  43. Bailey, T.L.; MacHanick, P. Inferring Direct DNA Binding from ChIP-Seq. Nucleic Acids Res. 2012, 40, e128. [Google Scholar] [CrossRef] [Green Version]
  44. Kanehisa, M.; Furumichi, M.; Sato, Y.; Ishiguro-Watanabe, M.; Tanabe, M. KEGG: Integrating Viruses and Cellular Organisms. Nucleic Acids Res. 2021, 49, D545–D551. [Google Scholar] [CrossRef] [PubMed]
  45. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  46. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics Enrichment Tools: Paths toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [Green Version]
  47. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  48. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.-F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstrale, M.; Laurila, E.; et al. PGC-1α-Responsive Genes Involved in Oxidative Phosphorylation Are Coordinately Downregulated in Human Diabetes. Nat. Genet. 2003, 34, 267–273. [Google Scholar] [CrossRef]
  49. Mogi, S.; Dang, D.; Waes, C.V.A.N.; Ellis, D.; Atakilit, A.; Ramos, D.M. The Expression of Integrin AvB6 Promotes the Epithelial Cell Morphology and Suppresses Invasive Behavior in Transformed Oral Keratinocyte. Anticancer Res. 2005, 756, 751–755. [Google Scholar]
  50. Fisher, M.L.; Balinth, S.; Mills, A.A. ΔNp63α in Cancer: Importance and Therapeutic Opportunities. Trends Cell Biol. 2022. [Google Scholar] [CrossRef]
  51. Kumar, S.; Gabrilovich, D.; Chakrabarti, R.; Kumar, S.; Wilkes, D.W.; Samuel, N.; Blanco, M.A.; Nayak, A.; Alicea-torres, K.; Gluck, C.; et al. ΔNp63-Driven Recruitment of Myeloid-Derived Suppressor Cells Promotes Metastasis in Triple- Negative Breast Cancer. J. Clin. Invest. 2018, 128, 5095–5109. [Google Scholar] [CrossRef]
  52. Kubo, T.; Ichimiya, S.; Tonooka, A.; Nagashima, T.; Kikuchi, T.; Sato, N. P63 Induces CD4+ T-Cell Chemoattractant TARC/CCL17 in Human Epithelial Cells. J. Interf. Cytokine Res. 2008, 28, 725–732. [Google Scholar] [CrossRef]
  53. Mitchell, M.J.; King, M.R. ΔNp63 Versatilely Regulates a Broad NF-ΚB Gene Program and Promotes Squamous Epithelial Proliferation, Migration and Inflammation. Cancer Res. 2014, 71, 3688–3700. [Google Scholar] [CrossRef] [Green Version]
  54. Soares, E.; Zhou, H. Master Regulatory Role of P63 in Epidermal Development and Disease. Cell. Mol. Life Sci. 2018, 75, 1179–1190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Romano, R.; Birkaya, B.; Sinha, S. A Functional Enhancer of Keratin14 Is a Direct Transcriptional Target of D Np63. J. Invest. Dermatol. 2007, 127, 1175–1186. [Google Scholar] [CrossRef] [Green Version]
  56. Dang, T.T.; Westcott, J.M.; Maine, E.A.; Kanchwala, M.; Xing, C.; Pearson, G.W. ΔNp63α Induces the Expression of FAT2 and Slug to Promote Tumor Invasion. Oncotarget 2016, 7, 28592–28611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yi, M.; Tan, Y.; Wang, L.; Cai, J.; Li, X.; Zeng, Z.; Xiong, W.; Li, G.; Li, X.; Tan, P.; et al. TP63 Links Chromatin Remodeling and Enhancer Reprogramming to Epidermal Differentiation and Squamous Cell Carcinoma Development. Cell. Mol. Life Sci. 2020, 77, 4325–4346. [Google Scholar] [CrossRef]
  58. Hamdan, F.H.; Johnsen, S.A. DeltaNp63-Dependent Super Enhancers Define Molecular Identity in Pancreatic Cancer by an Interconnected Transcription Factor Network. Proc. Natl. Acad. Sci. USA 2018, 115, 201812915. [Google Scholar] [CrossRef] [Green Version]
  59. Tang, X.H.; Knudsen, B.; Bemis, D.; Tickoo, S.; Gudas, L.J. Oral Cavity and Esophageal Carcinogenesis Modeled in Carcinogen-Treated Mice. Clin. Cancer Res. 2004, 10, 301–313. [Google Scholar] [CrossRef] [Green Version]
  60. Devos, M.; Declercq, W. ΔNp63α Acts as a Lineage-Survival Oncogene in Squamous Cell Carcinoma. J. Med. Surg. Pathol. 2017, 2, 7–9. [Google Scholar] [CrossRef]
  61. Nekulova, M.; Holcakova, J.; Coates, P.; Vojtesek, B. The Role of P63 in Cancer, Stem Cells and Cancer Stem Cells. Cell. Mol. Biol. Lett. 2011, 16, 296–327. [Google Scholar] [CrossRef] [PubMed]
  62. Lambert, A.W.; Fiore, C.; Chutake, Y.; Verhaar, E.R.; Strasser, P.C.; Chen, M.W.; Farouq, D.; Das, S.; Li, X.; Eaton, E.N.; et al. ΔNp63/P73 Drive Metastatic Colonization by Controlling a Regenerative Epithelial Stem Cell Program in Quasi-Mesenchymal Cancer Stem Cells. Dev. Cell 2022, 57, 2714–2730.e8. [Google Scholar] [CrossRef]
  63. Barbieri, C.E.; Tang, L.J.; Brown, K.A.; Pietenpol, J.A. Loss of P63 Leads to Increased Cell Migration and Up-Regulation of Genes Involved in Invasion and Metastasis. Cancer Res. 2006, 66, 7589–7597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Provost, P.; Doucet, J.; Stock, A.; Gerisch, G.; Samuelsson, B.; Rådmark, O. Coactosin-like Protein, a Human F-Actin-Binding Protein: Critical Role of Lysine-75. Biochem. J. 2001, 359, 255–263. [Google Scholar] [CrossRef] [PubMed]
  65. Li, G.; Yin, Y.; Chen, J.; Fan, Y.; Ma, J.; Huang, Y.; Chen, C.; Dai, P.; Chen, S.; Zhao, S. Coactosin-like Protein 1 Inhibits Neuronal Migration during Mouse Corticogenesis. J. Vet. Sci. 2018, 19, 21–26. [Google Scholar] [CrossRef] [PubMed]
  66. Hye-Cheol, J.; Gwang-Il, K.; Sang-Ho, C.; Kwang-Hyung, L.; Jung-Jae, K.; Jeong-Hee, K.; Kwang-Hoe, C. Proteomic Analysis of Human Small Cell Lung Cancer Tissues: Up-Regulation of Coactosin-like Protein-1. J. Proteome Res. 2011, 10, 269–276. [Google Scholar] [CrossRef]
  67. Shao, S.; Fan, Y.; Zhong, C.; Zhu, X.; Zhu, J. Coactosin-like Protein (Cotl1) Promotes Glioblastoma (Gbm) Growth in Vitro and in Vivo. Cancer Manag. Res. 2020, 12, 10909–10917. [Google Scholar] [CrossRef]
  68. Wang, L.; Wang, X.; Jia, Y.; Guo, F.; Zhengjun, S.; Shao, Z. Intratumoural Heterogeneity and Clone Evolution of Oral Squamous Cell Carcinoma. Mol. Carcinog. 2021, 60, 758–768. [Google Scholar] [CrossRef]
  69. Lim, S.S.; Vos, T.; Flaxman, A.D.; Danaei, G.; Shibuya, K.; Adair-Rohani, H.; Amann, M.; Anderson, H.R.; Andrews, K.G.; Aryee, M.; et al. A Comparative Risk Assessment of Burden of Disease and Injury Attributable to 67 Risk Factors and Risk Factor Clusters in 21 Regions, 1990-2010: A Systematic Analysis for the Global Burden of Disease Study 2010. Lancet 2012, 380, 2224–2260. [Google Scholar] [CrossRef] [Green Version]
  70. Agudo, A.; Bonet, C.; Travier, N.; González, C.A.; Vineis, P.; Bueno-de-Mesquita, H.B.; Trichopoulos, D.; Boffetta, P.; Clavel-Chapelon, F.; Boutron-Ruault, M.C.; et al. Impact of Cigarette Smoking on Cancer Risk in the European Prospective Investigation into Cancer and Nutrition Study. J. Clin. Oncol. 2012, 30, 4550–4557. [Google Scholar] [CrossRef]
  71. Li, Q.; Dong, H.; Yang, G.; Song, Y.; Mou, Y.; Ni, Y. Mouse Tumor-Bearing Models as Preclinical Study Platforms for Oral Squamous Cell Carcinoma. Front. Oncol. 2020, 10, 212. [Google Scholar] [CrossRef] [PubMed]
  72. Judd, N.P.; Winkler, A.E.; Murillo-Sauca, O.; Brotman, J.J.; Law, J.H.; Lewis, J.S., Jr.; Dunn, G.P.; Bui, J.D.; Sunwoo, J.B.; Uppaluri, R. ERK1/2 Regulation of CD44 Modulates Oral Cancer Aggressiveness. Cancer Res. 2012, 72, 365–374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Zhou, L.; Zeng, Z.; Egloff, A.M.; Zhang, F.; Guo, F.; Campbell, K.M.; Du, P.; Fu, J.; Zolkind, P.; Ma, X.; et al. Checkpoint Blockade-Induced CD8+ T Cell Differentiation in Head and Neck Cancer Responders. J. Immunother. Cancer 2022, 10, e004034. [Google Scholar] [CrossRef] [PubMed]
  74. Korpela, S.P.; Hinz, T.K.; Oweida, A.; Kim, J.; Calhoun, J.; Ferris, R.; Nemenoff, R.A.; Karam, S.D.; Clambey, E.T.; Heasley, L.E. Role of Epidermal Growth Factor Receptor Inhibitor-Induced Interferon Pathway Signaling in the Head and Neck Squamous Cell Carcinoma Therapeutic Response. J. Transl. Med. 2021, 19, 43. [Google Scholar] [CrossRef] [PubMed]
  75. Udani, M.; Zen, Q.; Cottman, M.; Leonard, N.; Jefferson, S.; Daymont, C.; Truskey, G.; Telen, M.J. Basal Cell Adhesion Molecule/Lutheran Protein: The Receptor Critical for Sickle Cell Adhesion to Laminin. J. Clin. Invest. 1998, 101, 2550–2558. [Google Scholar] [CrossRef] [Green Version]
  76. Bartolini, A.; Cardaci, S.; Lamba, S.; Oddo, D.; Marchiò, C.; Cassoni, P.; Amoreo, C.A.; Corti, G.; Testori, A.; Bussolino, F.; et al. BCAM and LAMA5 Mediate the Recognition between Tumor Cells and the Endothelium in the Metastatic Spreading of KRAS-Mutant Colorectal Cancer. Clin. Cancer Res. 2016, 22, 4923–4933. [Google Scholar] [CrossRef] [Green Version]
  77. Awais, R.; Spiller, D.G.; White, M.R.H.; Paraoan, L. P63 Is Required beside P53 for PERP-Mediated Apoptosis in Uveal Melanoma. Br. J. Cancer 2016, 115, 983–992. [Google Scholar] [CrossRef] [Green Version]
  78. Beaudry, V.G.; Jiang, D.; Dusek, R.L.; Park, E.J.; Knezevich, S.; Ridd, K.; Vogel, H.; Bastian, B.C.; Attardi, L.D. Loss of the P53/P63 Regulated Desmosomal Protein Perp Promotes Tumorigenesis. PLoS Genet. 2010, 6, e1001168. [Google Scholar] [CrossRef] [Green Version]
  79. Ralhan, R.; DeSouza, L.V.; Matta, A.; Tripathi, S.C.; Ghanny, S.; Datta Gupta, S.; Bahadur, S.; Siu, K.W.M. Discovery and Verification of Head-and-Neck Cancer Biomarkers by Differential Protein Expression Analysis Using ITRAQ Labeling, Multidimensional Liquid Chromatography, and Tandem Mass Spectrometry. Mol. Cell. Proteomics 2008, 7, 1162–1173. [Google Scholar] [CrossRef] [Green Version]
  80. Danilov, A.V.; Neupane, D.; Nagaraja, A.S.; Feofanova, E.V.; Humphries, L.A.; DiRenzo, J.; Korc, M. DeltaNp63alpha-Mediated Induction of Epidermal Growth Factor Receptor Promotes Pancreatic Cancer Cell Growth and Chemoresistance. PLoS ONE 2011, 6, e26815. [Google Scholar] [CrossRef]
  81. Guo, S.; Yang, P.; Jiang, X.; Li, X.; Wang, Y.; Zhang, X.; Sun, B.; Zhang, Y.; Jia, Y. Genetic and Epigenetic Silencing of MircoRNA-506-3p Enhances COTL1 Oncogene Expression to Foster Non-Small Lung Cancer Progression. Oncotarget 2017, 8, 644–657. [Google Scholar] [CrossRef] [PubMed]
  82. Pei, B.; Li, T.; Qian, Q.; Fan, W.; He, X.; Zhu, Y.; Xu, L. Downregulation of MicroRNA-30c-5p Was Responsible for Cell Migration and Tumor Metastasis via COTL1-Mediated Microfilament Arrangement in Breast Cancer. Gland Surg. 2020, 9, 747–758. [Google Scholar] [CrossRef]
  83. Lv, M.; Shen, Y.; Yang, J.; Li, S.; Wang, B.; Chen, Z.; Li, P.; Liu, P.; Yang, J. Angiomotin Family Members: Oncogenes or Tumor Suppressors? Int. J. Biol. Sci. 2017, 13, 772–781. [Google Scholar] [CrossRef] [PubMed]
  84. He, X.; Wang, L.; Riedel, H.; Wang, K.; Yang, Y.; Dinu, C.Z.; Rojanasakul, Y. Mesothelin Promotes Epithelial-to-Mesenchymal Transition and Tumorigenicity of Human Lung Cancer and Mesothelioma Cells. Mol. Cancer 2017, 16, 63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Chang, M.C.; Chen, C.A.; Chen, P.J.; Chiang, Y.C.; Chen, Y.L.; Mao, T.L.; Lin, H.W.; Lin Chiang, W.H.; Cheng, W.F. Mesothelin Enhances Invasion of Ovarian Cancer by Inducing MMP-7 through MAPK/ERK and JNK Pathways. Biochem. J. 2012, 442, 293–302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Wang, Y.; Wang, L.; Li, D.; Wang, H.B.; Chen, Q.F. Mesothelin Promotes Invasion and Metastasis in Breast Cancer Cells. J. Int. Med. Res. 2012, 40, 2109–2116. [Google Scholar] [CrossRef]
  87. Kim, J.; Shapiro, M.J.; Bamidele, A.O.; Gurel, P.; Thapa, P.; Higgs, H.N.; Hedin, K.E.; Shapiro, V.S.; Billadeau, D.D. Coactosin-like 1 Antagonizes Cofilin to Promote Lamellipodial Protrusion at the Immune Synapse. PLoS ONE 2014, 9, e85090. [Google Scholar] [CrossRef]
  88. Vidotto, A.; Polachini, G.M.; De Paula-Silva, M.; Oliani, S.M.; Henrique, T.; López, R.V.M.; Cury, P.M.; Nunes, F.D.; Góis-Filho, J.F.; De Carvalho, M.B.; et al. Differentially Expressed Proteins in Positive versus Negative HNSCC Lymph Nodes. BMC Med. Genomics 2018, 11, 73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Nekulova, M.; Holcakova, J.; Gu, X.; Hrabal, V.; Galtsidis, S.; Orzol, P.; Liu, Y.; Logotheti, S.; Zoumpourlis, V.; Nylander, K.; et al. ΔNp63α Expression Induces Loss of Cell Adhesion in Triple-Negative Breast Cancer Cells. BMC Cancer 2016, 16, 782. [Google Scholar] [CrossRef]
Figure 1. Profile of 4NQO-generated mOSCC cell lines. (A) Western blots for levels of p63 and EMT-associated factors in B7E3, B4B8, and B7E11 cell lines. (B) Heat map of the expression of genes consistently enriched in p63high or p63low cells. Enriched hallmark terms found in the p63low B7E3 DEGs (C), the p63high B4B8, and B7E11 DEGs (D).
Figure 1. Profile of 4NQO-generated mOSCC cell lines. (A) Western blots for levels of p63 and EMT-associated factors in B7E3, B4B8, and B7E11 cell lines. (B) Heat map of the expression of genes consistently enriched in p63high or p63low cells. Enriched hallmark terms found in the p63low B7E3 DEGs (C), the p63high B4B8, and B7E11 DEGs (D).
Cancers 15 00446 g001
Figure 2. Identification of the p63 global targets in B7E11 cells by ChIP-seq. (A) Heat map of the average ChIP-seq signals from p63 binding sites for two p63 antibodies across the genome. (B) Top transcription factor motifs derived from CentriMo motif analysis on B7E11 consensus p63 ChIP peaks. (C) Distribution pattern of genomic features associated with p63 binding sites across the genome. (D) Integration of the identified p63 cistrome with the previously identified cell line-based DEG dataset. Enriched WikiPathways identified for p63high (E) and p63low (F) target genes.
Figure 2. Identification of the p63 global targets in B7E11 cells by ChIP-seq. (A) Heat map of the average ChIP-seq signals from p63 binding sites for two p63 antibodies across the genome. (B) Top transcription factor motifs derived from CentriMo motif analysis on B7E11 consensus p63 ChIP peaks. (C) Distribution pattern of genomic features associated with p63 binding sites across the genome. (D) Integration of the identified p63 cistrome with the previously identified cell line-based DEG dataset. Enriched WikiPathways identified for p63high (E) and p63low (F) target genes.
Cancers 15 00446 g002
Figure 3. p63 knockdown in p63high B7E11 cells. (A) Western blot analysis of p63 expression in B7E11 cells expressing either p63-targeting shRNAs or a nontargeting shRNA (shCTL). GAPDH was used as a loading control. A pie chart (B) and a volcano plot (C) of DEGs resulting from sh3 in B7E11 cells. Enriched hallmark terms in downregulated (D) and upregulated (E) DEGs in B7E11 cells with p63 knockdown.
Figure 3. p63 knockdown in p63high B7E11 cells. (A) Western blot analysis of p63 expression in B7E11 cells expressing either p63-targeting shRNAs or a nontargeting shRNA (shCTL). GAPDH was used as a loading control. A pie chart (B) and a volcano plot (C) of DEGs resulting from sh3 in B7E11 cells. Enriched hallmark terms in downregulated (D) and upregulated (E) DEGs in B7E11 cells with p63 knockdown.
Cancers 15 00446 g003
Figure 4. p63 overexpression in p63low B7E3 cells. (A) Western blot analysis of p63 expression in B7E3 cells with dox-inducible p63 expression. GAPDH was used as a loading control. A pie chart (B) and a volcano plot (C) of DEGs resulting from induced p63 overexpression in B7E3 cells. Enriched hallmark terms in downregulated (D) and upregulated (E) DEGs in p63-overexpressing B7E3 cells.
Figure 4. p63 overexpression in p63low B7E3 cells. (A) Western blot analysis of p63 expression in B7E3 cells with dox-inducible p63 expression. GAPDH was used as a loading control. A pie chart (B) and a volcano plot (C) of DEGs resulting from induced p63 overexpression in B7E3 cells. Enriched hallmark terms in downregulated (D) and upregulated (E) DEGs in p63-overexpressing B7E3 cells.
Cancers 15 00446 g004
Figure 5. Investigating the p63-driven gene signature in an independent 4NQO-induced mOSCC model. A pie chart (A) and a volcano plot (B) of DEGs in the dataset by Tang et al. [11]. Enriched KEGG pathways in downregulated (C) and upregulated (D) DEGs from the dataset by Tang et al.
Figure 5. Investigating the p63-driven gene signature in an independent 4NQO-induced mOSCC model. A pie chart (A) and a volcano plot (B) of DEGs in the dataset by Tang et al. [11]. Enriched KEGG pathways in downregulated (C) and upregulated (D) DEGs from the dataset by Tang et al.
Cancers 15 00446 g005
Figure 6. Combined mOSCC p63-based signature. (A) Heat map of the combined 22-gene signature in the B7E11 cells with p63 knockdown, B7E3 cells with induced p63 expression, and the Tang et al. [11] 4NQO mOSCC dataset. Genes in boldface are direct targets of p63, according to our high-confidence p63 target list. (B) Western blot analysis of the levels of proteins encoded by DEGs identified in our combined p63-based signature in B7E3 p63 overexpressing cells. GAPDH was used as a loading control.
Figure 6. Combined mOSCC p63-based signature. (A) Heat map of the combined 22-gene signature in the B7E11 cells with p63 knockdown, B7E3 cells with induced p63 expression, and the Tang et al. [11] 4NQO mOSCC dataset. Genes in boldface are direct targets of p63, according to our high-confidence p63 target list. (B) Western blot analysis of the levels of proteins encoded by DEGs identified in our combined p63-based signature in B7E3 p63 overexpressing cells. GAPDH was used as a loading control.
Cancers 15 00446 g006
Figure 7. scRNA-seq analysis of 4NQO-induced mouse ESCC shows enrichment of our p63-based OSCC gene signature. (A) UMAP clustering of different cell types in mouse ESCC tumors. (B) UMAPs shows a zoomed-in cancer cell cluster highlighted from overall UMAP clustering. UMAP distribution of p63 expression across tumor cell clusters highlights p63high and p63low populations. DEG analysis of p63high versus p63low cells shows enrichment of several DEGs previously identified in the p63-based gene expression signature. UMAP projections display the expression of these factors across the tumor cell population. Enriched hallmark terms in DEGs enriched in p63low (C) and p63high (D) cells.
Figure 7. scRNA-seq analysis of 4NQO-induced mouse ESCC shows enrichment of our p63-based OSCC gene signature. (A) UMAP clustering of different cell types in mouse ESCC tumors. (B) UMAPs shows a zoomed-in cancer cell cluster highlighted from overall UMAP clustering. UMAP distribution of p63 expression across tumor cell clusters highlights p63high and p63low populations. DEG analysis of p63high versus p63low cells shows enrichment of several DEGs previously identified in the p63-based gene expression signature. UMAP projections display the expression of these factors across the tumor cell population. Enriched hallmark terms in DEGs enriched in p63low (C) and p63high (D) cells.
Cancers 15 00446 g007
Figure 8. p63 regulates mOSCC cell line migration and proliferation. (A) Transwell migration assay with B7E3 cells with p63 expression induced by 50 or 200 ng dox. Cells that migrated through the Matrigel matrix by 12 h were then imaged and quantified by using ImageJ. (B) Spheroid assay with B7E3 cells with inducible p63 expression. Cells were treated with no dox or 100 ng of dox and allowed to form spheroids for 9 days. The sizes of the top 500 largest spheroids in both groups were then quantified using ImageJ. Values are presented as mean ± SD. Significance level was determined using two-tailed Student’s t-test for samples of equal variance.  p < 0.05, ★★ p < 0.01, and ★★★★ p < 0.001.
Figure 8. p63 regulates mOSCC cell line migration and proliferation. (A) Transwell migration assay with B7E3 cells with p63 expression induced by 50 or 200 ng dox. Cells that migrated through the Matrigel matrix by 12 h were then imaged and quantified by using ImageJ. (B) Spheroid assay with B7E3 cells with inducible p63 expression. Cells were treated with no dox or 100 ng of dox and allowed to form spheroids for 9 days. The sizes of the top 500 largest spheroids in both groups were then quantified using ImageJ. Values are presented as mean ± SD. Significance level was determined using two-tailed Student’s t-test for samples of equal variance.  p < 0.05, ★★ p < 0.01, and ★★★★ p < 0.001.
Cancers 15 00446 g008
Figure 9. Co-staining of p63 and COTL1 in B7E3 cells with induced p63 expression. Immunofluorescence images of p63 (green), COTL1 (red), and DAPI (blue) in B7E3 cells treated with no dox (A), 50 ng dox (B), or 200 ng dox (C).
Figure 9. Co-staining of p63 and COTL1 in B7E3 cells with induced p63 expression. Immunofluorescence images of p63 (green), COTL1 (red), and DAPI (blue) in B7E3 cells treated with no dox (A), 50 ng dox (B), or 200 ng dox (C).
Cancers 15 00446 g009
Figure 10. Co-immunofluorescence of p63 and COTL1 in human OSCC tissue. Two independent samples of well-differentiated, moderately differentiated, and poorly differentiated human OSCC tissues were co-stained for p63 (green) and COTL1 (red), and imaged at 20× magnification. DAPI (blue) was used as a nuclear stain.
Figure 10. Co-immunofluorescence of p63 and COTL1 in human OSCC tissue. Two independent samples of well-differentiated, moderately differentiated, and poorly differentiated human OSCC tissues were co-stained for p63 (green) and COTL1 (red), and imaged at 20× magnification. DAPI (blue) was used as a nuclear stain.
Cancers 15 00446 g010
Figure 11. Immunohistochemical staining in HNSCC tumor microarray tissues. Staining of p63 (A), COTL1 (B), and K14 (C) across normal, malignant tumor stage II, and malignant tumor stage III tissues at 10× magnification.
Figure 11. Immunohistochemical staining in HNSCC tumor microarray tissues. Staining of p63 (A), COTL1 (B), and K14 (C) across normal, malignant tumor stage II, and malignant tumor stage III tissues at 10× magnification.
Cancers 15 00446 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Glathar, A.R.; Oyelakin, A.; Nayak, K.B.; Sosa, J.; Romano, R.-A.; Sinha, S. A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma. Cancers 2023, 15, 446. https://doi.org/10.3390/cancers15020446

AMA Style

Glathar AR, Oyelakin A, Nayak KB, Sosa J, Romano R-A, Sinha S. A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma. Cancers. 2023; 15(2):446. https://doi.org/10.3390/cancers15020446

Chicago/Turabian Style

Glathar, Alexandra Ruth, Akinsola Oyelakin, Kasturi Bala Nayak, Jennifer Sosa, Rose-Anne Romano, and Satrajit Sinha. 2023. "A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma" Cancers 15, no. 2: 446. https://doi.org/10.3390/cancers15020446

APA Style

Glathar, A. R., Oyelakin, A., Nayak, K. B., Sosa, J., Romano, R. -A., & Sinha, S. (2023). A Systemic and Integrated Analysis of p63-Driven Regulatory Networks in Mouse Oral Squamous Cell Carcinoma. Cancers, 15(2), 446. https://doi.org/10.3390/cancers15020446

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop