An integrated genomic approach identifies follistatin as a target of the p63-epidermal growth factor receptor oncogenic network in head and neck squamous cell carcinoma.

Although numerous putative oncogenes have been associated with the etiology of head and neck squamous cell carcinoma (HNSCC), the mechanisms by which these oncogenes and their downstream targets mediate tumor progression have not been fully elucidated. We performed an integrative analysis to identify a crucial set of targets of the oncogenic transcription factor p63 that are common across multiple transcriptomic datasets obtained from HNSCC patients, and representative cell line models. Notably, our analysis revealed FST which encodes follistatin, a secreted glycoprotein that inhibits the transforming growth factor TGFß/activin signaling pathways, to be a direct transcriptional target of p63. In addition, we found that FST expression is also driven by epidermal growth factor receptor EGFR signaling, thus mediating a functional link between the TGF-ß and EGFR pathways. We show through loss- and gain-of-function studies that FST predominantly imparts a tumor-growth and migratory phenotype in HNSCC cells. Furthermore, analysis of single-cell RNA sequencing data from HNSCC patients unveiled cancer cells as the dominant source of FST within the tumor microenvironment and exposed a correlation between the expression of FST and its regulators with immune infiltrates. We propose FST as a prognostic biomarker for patient survival and a compelling candidate mediating the broad effects of p63 on the tumor and its associated microenvironment.

EHF is a novel regulator of cellular redox metabolism and predicts patient prognosis in HNSCC.

Head and Neck Squamous Cell Carcinoma (HNSCC) is a heterogeneous disease with relatively high morbidity and mortality rates. The lack of effective therapies, high recurrence rates and drug resistance driven in part, by tumor heterogeneity, contribute to the poor prognosis for patients diagnosed with this cancer. This problem is further exacerbated by the fact that key regulatory factors contributing to the disease diversity remains largely elusive. Here, we have identified EHF as an important member of the ETS family of transcription factors that is highly expressed in normal oral tissues, but lost during HNSCC progression. Interestingly, HNSCC tumors and cell lines exhibited a dichotomy of high and low EHF expression, and patients whose tumors retained EHF expression showed significantly better prognosis, suggesting a potential tumor suppressive role for EHF. To address this, we have performed gain and loss of function studies and leveraged bulk and single-cell cancer genomic datasets to identify global EHF targets by RNA-sequencing (RNA-seq) and Chromatin Immunoprecipitation and next generation sequencing (ChIP-seq) experiments of HNSCC cell lines. These mechanistic studies have revealed that EHF, acts as a regulator of a broad spectrum of metabolic processes, specifically targeting regulators of redox homeostasis such as NRF2 and SOX2. Our immunostaining results confirm the mutually exclusive expression patterns of EHF and SOX2 in HNSCC tumors and suggest a possible role for these two factors in establishing discrete metabolic states within the tumor microenvironment. Taken together, EHF may serve as a novel prognostic marker for classifying HNSCC patients for actionable and targeted therapeutic intervention.

p63 Directs Subtype-Specific Gene Expression in HPV+ Head and Neck Squamous Cell Carcinoma.

The complex heterogeneity of head and neck squamous cell carcinoma (HNSCC) reflects a diverse underlying etiology. This heterogeneity is also apparent within Human Papillomavirus-positive (HPV+) HNSCC subtypes, which have distinct gene expression profiles and patient outcomes. One aggressive HPV+ HNSCC subtype is characterized by elevated expression of genes involved in keratinization, a process regulated by the oncogenic transcription factor ΔNp63. Furthermore, the human TP63 gene locus is a frequent HPV integration site and HPV oncoproteins drive ΔNp63 expression, suggesting an unexplored functional link between ΔNp63 and HPV+ HNSCC. Here we show that HPV+ HNSCCs can be molecularly stratified according to ΔNp63 expression levels and derive a ΔNp63-associated gene signature profile for such tumors. We leveraged RNA-seq data from p63 knockdown cells and ChIP-seq data for p63 and histone marks from two ΔNp63high HPV+ HNSCC cell lines to identify an epigenetically refined ΔNp63 cistrome. Our integrated analyses reveal crucial ΔNp63-bound super-enhancers likely to mediate HPV+ HNSCC subtype-specific gene expression that is anchored, in part, by the PI3K-mTOR pathway. These findings implicate ΔNp63 as a key regulator of essential oncogenic pathways in a subtype of HPV+ HNSCC that can be exploited as a biomarker for patient stratification and treatment choices.

Reactivation of super-enhancers by KLF4 in human Head and Neck Squamous Cell Carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is a disease of significant morbidity and mortality and rarely diagnosed in early stages. Despite extensive genetic and genomic characterization, targeted therapeutics and diagnostic markers of HNSCC are lacking due to the inherent heterogeneity and complexity of the disease. Herein, we have generated the global histone mark based epigenomic and transcriptomic cartogram of SCC25, a representative cell type of mesenchymal HNSCC and its normal oral keratinocyte counterpart. Examination of genomic regions marked by differential chromatin states and associated with misregulated gene expression led us to identify SCC25 enriched regulatory sequences and transcription factors (TF) motifs. These findings were further strengthened by ATAC-seq based open chromatin and TF footprint analysis which unearthed Krüppel-like Factor 4 (KLF4) as a potential key regulator of the SCC25 cistrome. We reaffirm the results obtained from in silico and chromatin studies in SCC25 by ChIP-seq of KLF4 and identify ΔNp63 as a co-oncogenic driver of the cancer-specific gene expression milieu. Taken together, our results lead us to propose a model where elevated KLF4 levels sustains the oncogenic state of HNSCC by reactivating repressed chromatin domains at key downstream genes, often by targeting super-enhancers.

Molecular dissection of the oncogenic role of ETS1 in the mesenchymal subtypes of head and neck squamous cell carcinoma.

Head and Neck Squamous Cell Carcinoma (HNSCC) is a heterogeneous disease of significant mortality and with limited treatment options. Recent genomic analysis of HNSCC tumors has identified several distinct molecular classes, of which the mesenchymal subtype is associated with Epithelial to Mesenchymal Transition (EMT) and shown to correlate with poor survival and drug resistance. Here, we utilize an integrated approach to characterize the molecular function of ETS1, an oncogenic transcription factor specifically enriched in Mesenchymal tumors. To identify the global ETS1 cistrome, we have performed integrated analysis of RNA-Seq, ChIP-Seq and epigenomic datasets in SCC25, a representative ETS1high mesenchymal HNSCC cell line. Our studies implicate ETS1 as a crucial regulator of broader oncogenic processes and specifically Mesenchymal phenotypes, such as EMT and cellular invasion. We found that ETS1 preferentially binds cancer specific regulator elements, in particular Super-Enhancers of key EMT genes, highlighting its role as a master regulator. Finally, we show evidence that ETS1 plays a crucial role in regulating the TGF-ß pathway in Mesenchymal cell lines and in leading-edge cells in primary HNSCC tumors that are endowed with partial-EMT features. Collectively our study highlights ETS1 as a key regulator of TGF-ß associated EMT and reveals new avenues for sub-type specific therapeutic intervention.

Targeting CD39 Toward Activated Platelets Reduces Systemic Inflammation and Improves Survival in Sepsis: A Preclinical Pilot Study.

OBJECTIVES: Sepsis is associated with a systemic inflammatory reaction, which can result in a life-endangering organ dysfunction. Pro-inflammatory responses during sepsis are characterized by increased activation of leukocytes and platelets, formation of platelet-neutrophil aggregates, and cytokine production. Sequestration of platelet-neutrophil aggregates in the microvasculature contributes to tissue damage during sepsis. At present no effective therapeutic strategy to ameliorate these events is available. In this preclinical pilot study, a novel anti-inflammatory approach was evaluated, which targets nucleoside triphosphate hydrolase activity toward activated platelets via a recombinant fusion protein combining a single-chain antibody against activated glycoprotein IIb/IIIa and the extracellular domain of CD39 (targ-CD39). DESIGN: Experimental animal study and cell culture study. SETTING: University-based experimental laboratory. SUBJECTS: Human dermal microvascular endothelial cells 1, human platelets and neutrophils, and C57BL/6NCrl mice. INTERVENTIONS: Platelet-leukocyte-endothelium interactions were evaluated under inflammatory conditions in vitro and in a murine lipopolysaccharide-induced sepsis model in vivo. The outcome of polymicrobial sepsis was evaluated in a murine cecal ligation and puncture model. To evaluate the anti-inflammatory potential of activated platelet targeted nucleoside triphosphate hydrolase activity, we employed a potato apyrase in vitro and in vivo, as well as targ-CD39 and as a control, nontarg-CD39 in vivo. MEASUREMENTS AND MAIN RESULTS: Under conditions of sepsis, agents with nucleoside triphosphate hydrolase activity decreased platelet-leukocyte-endothelium interaction, transcription of pro-inflammatory cytokines, microvascular platelet-neutrophil aggregate sequestration, activation marker expression on platelets and neutrophils contained in these aggregates, leukocyte extravasation, and organ damage. Targ-CD39 had the strongest effect on these variables and retained hemostasis in contrast to nontarg-CD39 and potato apyrase. Most importantly, targ-CD39 improved survival in the cecal ligation and puncture model to a stronger extent then nontarg-CD39 and potato apyrase. CONCLUSIONS: Targeting nucleoside triphosphate hydrolase activity (CD39) toward activated platelets is a promising new treatment concept to decrease systemic inflammation and mortality of sepsis. This innovative therapeutic approach warrants further development toward clinical application.

ΔNp63-driven recruitment of myeloid-derived suppressor cells promotes metastasis in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is particularly aggressive, with enhanced incidence of tumor relapse, resistance to chemotherapy, and metastases. As the mechanistic basis for this aggressive phenotype is unclear, treatment options are limited. Here, we showed an increased population of myeloid-derived immunosuppressor cells (MDSCs) in TNBC patients compared with non-TNBC patients. We found that high levels of the transcription factor ΔNp63 correlate with an increased number of MDSCs in basal TNBC patients, and that ΔNp63 promotes tumor growth, progression, and metastasis in human and mouse TNBC cells. Furthermore, we showed that MDSC recruitment to the primary tumor and metastatic sites occurs via direct ΔNp63-dependent activation of the chemokines CXCL2 and CCL22. CXCR2/CCR4 inhibitors reduced MDSC recruitment, angiogenesis, and metastasis, highlighting a novel treatment option for this subset of TNBC patients. Finally, we found that MDSCs secrete prometastatic factors such as MMP9 and chitinase 3-like 1 to promote TNBC cancer stem cell function, thereby identifying a nonimmunologic role for MDSCs in promoting TNBC progression. These findings identify a unique crosstalk between ΔNp63+ TNBC cells and MDSCs that promotes tumor progression and metastasis, which could be exploited in future combined immunotherapy/chemotherapy strategies for TNBC patients.

Functional characterization and genomic studies of a novel murine submandibular gland epithelial cell line.

A better understanding of the normal and diseased biology of salivary glands (SG) has been hampered, in part, due to difficulties in cultivating and maintaining salivary epithelial cells. Towards this end, we have generated a mouse salivary gland epithelial cell (mSGc) culture system that is well-suited for the molecular characterization of SG cells and their differentiation program. We demonstrate that mSGc can be maintained for multiple passages without a loss of proliferation potential, readily form 3D-spheroids and importantly express a panel of well-established salivary gland epithelial cell markers. Moreover, mSGc 3D-spheroids also exhibit functional maturation as evident by robust agonist-induced intracellular calcium signaling. Finally, transcriptomic characterization of mSGc by RNA-seq and hierarchical clustering analysis with adult organ RNA-seq datasets reveal that mSGc retain most of the molecular attributes of adult mouse salivary gland. This well-characterized mouse salivary gland cell line will fill a critical void in the field by offering a valuable resource to examine various mechanistic aspects of mouse salivary gland biology.

Evolutionary re-wiring of p63 and the epigenomic regulatory landscape in keratinocytes and its potential implications on species-specific gene expression and phenotypes.

Although epidermal keratinocyte development and differentiation proceeds in similar fashion between humans and mice, evolutionary pressures have also wrought significant species-specific physiological differences. These differences between species could arise in part, by the rewiring of regulatory network due to changes in the global targets of lineage-specific transcriptional master regulators such as p63. Here we have performed a systematic and comparative analysis of the p63 target gene network within the integrated framework of the transcriptomic and epigenomic landscape of mouse and human keratinocytes. We determined that there exists a core set of ∼1600 genomic regions distributed among enhancers and super-enhancers, which are conserved and occupied by p63 in keratinocytes from both species. Notably, these DNA segments are typified by consensus p63 binding motifs under purifying selection and are associated with genes involved in key keratinocyte and skin-centric biological processes. However, the majority of the p63-bound mouse target regions consist of either murine-specific DNA elements that are not alignable to the human genome or exhibit no p63 binding in the orthologous syntenic regions, typifying an occupancy lost subset. Our results suggest that these evolutionarily divergent regions have undergone significant turnover of p63 binding sites and are associated with an underlying inactive and inaccessible chromatin state, indicative of their selective functional activity in the transcriptional regulatory network in mouse but not human. Furthermore, we demonstrate that this selective targeting of genes by p63 correlates with subtle, but measurable transcriptional differences in mouse and human keratinocytes that converges on major metabolic processes, which often exhibit species-specific trends. Collectively our study offers possible molecular explanation for the observable phenotypic differences between the mouse and human skin and broadly informs on the prevailing principles that govern the tug-of-war between evolutionary forces of rigidity and plasticity over transcriptional regulatory programs.

A global analysis of the complex landscape of isoforms and regulatory networks of p63 in human cells and tissues.

BACKGROUND: The transcription factor p63 belongs to the p53/p63/p73 family and plays key functional roles during normal epithelial development and differentiation and in pathological states such as squamous cell carcinomas. The human TP63 gene, located on chromosome 3q28 is driven by two promoters that generate the full-length transactivating (TA) and N-terminal truncated (ΔN) isoforms. Furthermore alternative splicing at the C-terminus gives rise to additional α, ß, γ and likely several other minor variants. Teasing out the expression and biological function of each p63 variant has been both the focus of, and a cause for contention in the p63 field. RESULTS: Here we have taken advantage of a burgeoning RNA-Seq based genomic data-sets to examine the global expression profiles of p63 isoforms across commonly utilized human cell-lines and major tissues and organs. Consistent with earlier studies, we find ΔNp63 transcripts, primarily that of the ΔNp63α isoforms, to be expressed in most cells of epithelial origin such as those of skin and oral tissues, mammary glands and squamous cell carcinomas. In contrast, TAp63 is not expressed in the majority of normal cell-types and tissues; rather it is selectively expressed at moderate to high levels in a subset of Burkitt's and diffuse large B-cell lymphoma cell lines. We verify this differential expression pattern of p63 isoforms by Western blot analysis, using newly developed ΔN and TA specific antibodies. Furthermore using unsupervised clustering of human cell lines, tissues and organs, we show that ΔNp63 and TAp63 driven transcriptional networks involve very distinct sets of molecular players, which may underlie their different biological functions. CONCLUSIONS: In this study we report comprehensive and global expression profiles of p63 isoforms and their relationship to p53/p73 and other potential transcriptional co-regulators. We curate publicly available data generated in part by consortiums such as ENCODE, FANTOM and Human Protein Atlas to delineate the vastly different transcriptomic landscapes of ΔNp63 and TAp63. Our studies help not only in dispelling prevailing myths and controversies on p63 expression in commonly used human cell lines but also augur new isoform- and cell type-specific activities of p63.

Suberoylanilide hydroxamic acid limits migration and invasion of glioma cells in two and three dimensional culture.

High grade gliomas are aggressive cancers that are not well addressed by current chemotherapies, in large measure because these drugs do not curtail the diffuse invasion of glioma cells into brain tissue surrounding the tumor. Here, we investigate the effects of suberoylanilide hydroxamic acid (SAHA) on glioma cells in 2D and 3D in vitro assays, as SAHA has previously been shown to significantly increase apoptosis, decrease proliferation, and interfere with migration in other cell lines. We find that SAHA has significant independent effects on proliferation, migration, and invasion. These effects are seen in both 2D and 3D culture. In 3D culture, with glioma spheroids embedded in collagen I matrices, SAHA independently limits both glioma invasion and the reorganization of the tumor surroundings that usually proceeds such invasion. The decreased matrix reorganization and invasion is not accompanied by decreased production or activity of matrix-metalloproteases but instead may be related to increased cell-cell adhesion.