Developmental factor IRF6 exhibits tumor suppressor activity in squamous cell carcinomas

The multiple roles of interferon regulatory factor family in health and disease

Interferon Regulatory Factors (IRFs), a family of transcription factors, profoundly influence the immune system, impacting both physiological and pathological processes. This review explores the diverse functions of nine mammalian IRF members, each featuring conserved domains essential for interactions with other transcription factors and cofactors. These interactions allow IRFs to modulate a broad spectrum of physiological processes, encompassing host defense, immune response, and cell development. Conversely, their pivotal role in immune regulation implicates them in the pathophysiology of various diseases, such as infectious diseases, autoimmune disorders, metabolic diseases, and cancers. In this context, IRFs display a dichotomous nature, functioning as both tumor suppressors and promoters, contingent upon the specific disease milieu. Post-translational modifications of IRFs, including phosphorylation and ubiquitination, play a crucial role in modulating their function, stability, and activation. As prospective biomarkers and therapeutic targets, IRFs present promising opportunities for disease intervention. Further research is needed to elucidate the precise mechanisms governing IRF regulation, potentially pioneering innovative therapeutic strategies, particularly in cancer treatment, where the equilibrium of IRF activities is of paramount importance.

Breast Cancer Stem Cells Upregulate IRF6 in Stromal Fibroblasts to Induce Stromagenesis

The microenvironment of a cancer stem cell (CSC) niche is often found in coexistence with cancer-associated fibroblasts (CAFs). Here, we show the first in-depth analysis of the interaction between primary triple-negative breast cancer stem cells (BCSCs) with fibroblasts. Using 2D co-culture models with specific seeding ratios, we identified stromal fibroblast aggregation at the BCSC cluster periphery, and, on closer observation, the aggregated fibroblasts was found to encircle BCSC clusters in nematic organization. In addition, collagen type I and fibronectin accumulation were also found at the BCSC-stromal periphery. MACE-Seq analysis of BCSC-encapsulating fibroblasts displayed the transformation of stromal fibroblasts to CAFs and the upregulation of fibrosis regulating genes of which the Interferon Regulatory Factor 6 (IRF6) gene was identified. Loss of function experiments with the IRF6 gene decreased fibroblast encapsulation around BCSC clusters in 2D co-cultures. In BCSC xenografts, fibroblast IRF6 expression led to an increase in the stromal area and fibroblast density in tumors, in addition to a reduction in necrotic growth. Based on our findings, we propose that fibroblast IRF6 function is an important factor in the development of the stromal microenvironment and in sustaining the BCSC tumor niche.

Metabolic landscape of the healthy pancreas and pancreatic tumor microenvironment

Pancreatic cancer, one of the deadliest human malignancies, is characterized by a fibro-inflammatory tumor microenvironment and wide array of metabolic alterations. To comprehensively map metabolism in a cell type-specific manner, we harnessed a unique single-cell RNA-sequencing dataset of normal human pancreata. This was compared with human pancreatic cancer samples using a computational pipeline optimized for this study. In the cancer cells we observed enhanced biosynthetic programs. We identified downregulation of mitochondrial programs in several immune populations, relative to their normal counterparts in healthy pancreas. Although granulocytes, B cells, and CD8+ T cells all downregulated oxidative phosphorylation, the mechanisms by which this occurred were cell type specific. In fact, the expression pattern of the electron transport chain complexes was sufficient to identify immune cell types without the use of lineage markers. We also observed changes in tumor-associated macrophage (TAM) lipid metabolism, with increased expression of enzymes mediating unsaturated fatty acid synthesis and upregulation in cholesterol export. Concurrently, cancer cells exhibited upregulation of lipid/cholesterol receptor import. We thus identified a potential crosstalk whereby TAMs provide cholesterol to cancer cells. We suggest that this may be a new mechanism boosting cancer cell growth and a therapeutic target in the future.

Sevoflurane inhibits lung cancer development by promoting FUS1 transcription via downregulating IRF6

Lung cancer is a major contributor to cancer deaths worldwide and is on the rise. Although surgical resection has been widely used as a standard therapy for lung cancer patients, the relapse rate after surgery is high. It is still unclear whether there is a potential drug that can reduce the probability of postsurgical recurrence in lung cancer patients. We used 5 typical lung cancer cell lines as well as 41 lung cancer tissue samples and paracancer tissue samples to investigate the expression levels of interferon regulatory factor 6 (IRF6) and tumor suppressor candidate 2 (TUSC2, also known as FUS1). We also treated lung cancer cells (H322 and A549) with different concentrations of sevoflurane to study its influence on lung cancer cell tumorigenesis. Lentivirus-mediated gain-of-function studies of IRF6 and FUS1 were applied to validate the role of IRF6 and FUS1 in lung cancer. Next, we used short hairpin RNA-mediated loss of function of IRF6 and luciferase, chromatin immunoprecipitation assays to validate the regulatory role of IRF6 on FUS1. Our findings reported that IRF6 was upregulated in lung cancer tissues, while FUS1 was downregulated. Functional assays revealed that sevoflurane inhibits lung cancer development by downregulating IRF6 expression. Luciferase and chromatin immunoprecipitation-quantitative real-time PCR assays uncovered that IRF6 represses FUS1 transcriptional expression in lung cancer cells. We have shown that sevoflurane prevents lung cancer development by downregulating IRF6 to stimulate FUS1 transcription, indicating that sevoflurane can be used as the potential anesthetic drug in surgical resection to reduce postoperative tumor relapse in lung cancer patients.

ABCC1 Is a ΔNp63 Target Gene Overexpressed in Squamous Cell Carcinoma

The transcription factor ΔNp63 plays a pivotal role in maintaining the integrity of stratified epithelial tissues by regulating the expression of distinct target genes involved in lineage specification, cell stemness, cell proliferation and differentiation. Here, we identified the ABC transporter subfamily member ABCC1 as a novel ΔNp63 target gene. We found that in immortalized human keratinocytes and in squamous cell carcinoma (SCC) cells, ∆Np63 induces the expression of ABCC1 by physically occupying a p63-binding site (p63 BS) located in the first intron of the ABCC1 gene locus. In cutaneous SCC and during the activation of the keratinocyte differentiation program, ∆Np63 and ABCC1 levels are positively correlated raising the possibility that ABCC1 might be involved in the regulation of the proliferative/differentiative capabilities of squamous tissue. However, we did not find any gross alteration in the structure and morphology of the epidermis in humanized hABCC1 knock-out mice. Conversely, we found that the genetic ablation of ABCC1 led to a marked reduction in inflammation-mediated proliferation of keratinocytes, suggesting that ABCC1 might be involved in the regulation of keratinocyte proliferation upon inflammatory/proliferative signals. In line with these observations, we found a significant increase in ABCC1 expression in squamous cell carcinomas (SCCs), a tumor type characterized by keratinocyte hyper-proliferation and a pro-inflammatory tumor microenvironment. Collectively, these data uncover ABCC1 as an additional ∆Np63 target gene potentially involved in those skin diseases characterized by dysregulation of proliferation/differentiation balance.

DNA Methylation Effects on Van der Woude Syndrome Phenotypic Variability

OBJECTIVE

Van der Woude Syndrome (VWS) presents with combinations of lip pits (LP) and cleft lip and/or cleft palate (CL/P, CPO). VWS phenotypic heterogeneity even amongst relatives, suggests that epigenetic factors may act as modifiers. IRF6, causal for 70% of VWS cases, and TP63 interact in a regulatory loop coordinating epithelial proliferation and differentiation in palatogenesis. We hypothesize that differential DNA methylation within IRF6 and TP63 regulatory regions underlie VWS phenotypic discordance.

METHODS

DNA methylation of CpG sites in IRF6 and TP63 promoters and in an IRF6 enhancer element was compared amongst blood or saliva DNA samples of 78 unrelated cases. Analyses were done separately for blood and saliva, within each sex and in combination, and to address cleft type (CL/P ± LP vs. CPO ± LP) and phenotypic severity (any cleft + LP vs. any cleft only).

RESULTS

For cleft type, blood samples showed higher IRF6 and TP63 promoter methylation on males with CPO ± LP compared to CL/P ± LP and on individuals with CPO ± LP compared to those with CL/P ± LP, respectively. Saliva samples showed higher IRF6 enhancer methylation on individuals with CPO ± LP compared to CL/P ± LP and contrary to above, lower TP63 promoter methylation on CPO ± LP compared to CL/P ± LP. For phenotypic severity, blood samples showed no differences; however, saliva samples showed higher IRF6 promoter methylation in individuals with any cleft + LP compared to those without lip pits.

CONCLUSION

We observed differential methylation in IRF6 and TP63 regulatory regions associated with cleft type and phenotypic severity, indicating that epigenetic changes in IRF6 and TP63 can contribute to phenotypic heterogeneity in VWS.

Interferon regulatory factor 6 (IRF6) determines intestinal epithelial cell development and immunity

Intestinal epithelial cell (IEC) responses to interferon (IFN) favor antiviral defense with minimal cytotoxicity, but IEC-specific factors that regulate these responses remain poorly understood. Interferon regulatory factors (IRFs) are a family of nine related transcription factors, and IRF6 is preferentially expressed by epithelial cells, but its roles in IEC immunity are unknown. In this study, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) screens found that Irf6 deficiency enhanced IFN-stimulated antiviral responses in transformed mouse IECs but not macrophages. Furthermore, knockout (KO) of Irf6 in IEC organoids resulted in profound changes to homeostasis and immunity gene expression. Irf6 KO organoids grew more slowly, and single-cell ribonucleic acid sequencing indicated reduced expression of genes in epithelial differentiation and immunity pathways. IFN-stimulated gene expression was also significantly different in Irf6 KO organoids, with increased expression of stress and apoptosis-associated genes. Functionally, the transcriptional changes in Irf6 KO organoids were associated with increased cytotoxicity upon IFN treatment or inflammasome activation. These data indicate a previously unappreciated role for IRF6 in IEC biology, including regulation of epithelial development and moderation of innate immune responses to minimize cytotoxicity and maintain barrier function.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

A novel noncanonical function for IRF6 in the recycling of E-cadherin

Interferon Regulatory Factor 6 (IRF6) is a transcription factor essential for keratinocyte cell-cell adhesions. Previously, we found that recycling of E-cadherin was defective in the absence of IRF6, yet total E-cadherin levels were not altered, suggesting a previously unknown, nontranscriptional function for IRF6. IRF6 protein contains a DNA binding domain (DBD) and a protein binding domain (PBD). The transcriptional function of IRF6 depends on its DBD and PBD, however, whether the PBD is necessary for the interaction with cytoplasmic proteins has yet to be demonstrated. Here, we show that an intact PBD is required for recruitment of cell-cell adhesion proteins at the plasma membrane, including the recycling of E-cadherin. Colocalizations and coimmunoprecipitations reveal that IRF6 forms a complex in recycling endosomes with Rab11, Myosin Vb, and E-cadherin, and that the PBD is required for this interaction. These data indicate that IRF6 is a novel effector of the endosomal recycling of E-cadherin and demonstrate a non-transcriptional function for IRF6 in regulating cell-cell adhesions.

Neural crest and periderm-specific requirements of Irf6 during neural tube and craniofacial development

IRF6 is a key genetic determinant of syndromic and non-syndromic cleft lip and palate. The ability to interrogate post-embryonic requirements of Irf6 has been hindered, as global Irf6 ablation in the mouse causes neonatal lethality. Prior work analyzing Irf6 in mouse models defined its role in the embryonic surface epithelium and periderm where it is required to regulate cell proliferation and differentiation. Several reports have also described Irf6 gene expression in other cell types, such as muscle, and neuroectoderm. However, analysis of a functional role in non-epithelial cell lineages has been incomplete due to the severity and lethality of the Irf6 knockout model and the paucity of work with a conditional Irf6 allele. Here we describe the generation and characterization of a new Irf6 floxed mouse model and analysis of Irf6 ablation in periderm and neural crest lineages. This work found that loss of Irf6 in periderm recapitulates a mild Irf6 null phenotype, suggesting that Irf6-mediated signaling in periderm plays a crucial role in regulating embryonic development. Further, conditional ablation of Irf6 in neural crest cells resulted in an anterior neural tube defect of variable penetrance. The generation of this conditional Irf6 allele allows for new insights into craniofacial development and new exploration into the post-natal role of Irf6.

The involvement of RIPK4 in TNF-α-stimulated IL-6 and IL-8 production by melanoma cells

PURPOSE

The receptor-interacting protein kinase (RIPK4) has an oncogenic function in melanoma, regulates NF-κB and Wnt/β-catenin pathways, and is sensitive to the BRAF inhibitors: vemurafenib and dabrafenib which lead to its decreased level. As its role in melanoma remains not fully understood, we examined the effects of its downregulation on the transcriptomic profile of melanoma.

METHODS

Applying RNA-seq, we revealed global alterations in the transcriptome of WM266.4 cells with RIPK4 silencing. Functional partners of RIPK4 were evaluated using STRING and GeneMANIA databases. Cells with transient knockdown (via siRNA) and stable knockout (via CRISPR/Cas9) of RIPK4 were stimulated with TNF-α. The expression levels of selected proteins were assessed using Western blot, ELISA, and qPCR.

RESULTS

Global analysis of gene expression changes indicates a complex role for RIPK4 in regulating adhesion, migration, proliferation, and inflammatory processes in melanoma cells. Our study highlights potential functional partners of RIPK4 such as BIRC3, TNF-α receptors, and MAP2K6. Data from RIPK4 knockout cells suggest a putative role for RIPK4 in modulating TNF-α-induced production of IL-8 and IL-6 through two distinct signaling pathways-BIRC3/NF-κB and p38/MAPK. Furthermore, increased serum TNF-α levels and the correlation of RIPK4 with NF-κB were revealed in melanoma patients.

CONCLUSION

These data reveal a complex role for RIPK4 in regulating the immune signaling network in melanoma cells and suggest that this kinase may represent an alternative target for melanoma-targeted adjuvant therapy.

Biophysical and functional characterization of N-terminal domain of Human Interferon Regulatory Factor 6

BACKGROUND

Interferon regulatory factor 6 (IRF6) has a key function in palate fusion during palatogenesis during embryonic development, and mutations in IRF6 cause orofacial clefting disorders.

METHODS AND RESULTS

The in silico analysis of IRF6 is done to obtain leads for the domain boundaries and subsequently the sub-cloning of the N-terminal domain of IRF6 into the pGEX-2TK expression vector and successfully optimized the overexpression and purification of recombinant glutathione S-transferase-fused NTD-IRF6 protein under native conditions. After cleavage of the GST tag, NTD-IRF6 was subjected to protein folding studies employing Circular Dichroism and Intrinsic fluorescence spectroscopy at variable pH, temperature, and denaturant. CD studies showed predominantly alpha-helical content and the highest stability of NTD-IRF6 at pH 9.0. A comparison of native and renatured protein depicts loss in the secondary structural content. Intrinsic fluorescence and quenching studies have identified that tryptophan residues are majorly present in the buried areas of the protein and a small fraction was on or near the protein surface. Upon the protein unfolding with a higher concentration of denaturant urea, the peak of fluorescence intensity decreased and red shifted, confirming that tryptophan residues are majorly present in a more polar environment. While regulating IFNβ gene expression during viral infection, the N-terminal domain binds to the promoter region of Virus Response Element-Interferon beta (VRE-IFNβ). Along with the protein folding analysis, this study also aimed to identify the DNA-binding activity and determine the binding affinities of NTD-IRF6 with the VRE-IFNβ promoter region. The protein-DNA interaction is specific as demonstrated by gel retardation assay and the kinetics of molecular interactions as quantified by Biolayer Interferometry showed a strong affinity with an affinity constant (KD) value of 7.96 × 10-10 M.

CONCLUSION

NTD-IRF6 consists of a mix of α-helix and β-sheets that show temperature-dependent cooperative unfolding between 40 °C and 55 °C. Urea-induced unfolding shows moderate tolerance to urea as the mid-transition concentration of urea (Cm) is 3.2 M. The tryptophan residues are majorly buried as depicted by fluorescence quenching studies. NTD-IRF6 has a specific and high affinity toward the promoter region of VRE-IFNβ.

Plasticity-induced repression of Irf6 underlies acquired resistance to cancer immunotherapy in pancreatic ductal adenocarcinoma

Acquired resistance to immunotherapy remains a critical yet incompletely understood biological mechanism. Here, using a mouse model of pancreatic ductal adenocarcinoma (PDAC) to study tumor relapse following immunotherapy-induced responses, we find that resistance is reproducibly associated with an epithelial-to-mesenchymal transition (EMT), with EMT-transcription factors ZEB1 and SNAIL functioning as master genetic and epigenetic regulators of this effect. Acquired resistance in this model is not due to immunosuppression in the tumor immune microenvironment, disruptions in the antigen presentation machinery, or altered expression of immune checkpoints. Rather, resistance is due to a tumor cell-intrinsic defect in T-cell killing. Molecularly, EMT leads to the epigenetic and transcriptional silencing of interferon regulatory factor 6 (Irf6), rendering tumor cells less sensitive to the pro-apoptotic effects of TNF-α. These findings indicate that acquired resistance to immunotherapy may be mediated by programs distinct from those governing primary resistance, including plasticity programs that render tumor cells impervious to T-cell killing.

DNA methylation effects on Van der Woude Syndrome phenotypic variability

OBJECTIVE

Van der Woude Syndrome (VWS) classically presents with combinations of lip pits (LP) and orofacial clefts, with marked phenotypic discordance even amongst individuals carrying the same mutation. Such discordance suggests a possible role for epigenetic factors as phenotypic modifiers. Both IRF6 , causal for 70% of VWS cases, and TP63 interact in a regulatory loop to coordinate epithelial proliferation and differentiation for palatogenesis. We hypothesize that differential DNA methylation (DNAm) in CpG sites within regulatory regions of IRF6 and TP63 are associated with VWS phenotypic discordance.

METHODS

We measured DNAm levels of CpG sites located in the promoter regions of IRF6 and TP63 and in an IRF6 enhancer element (MCS9.7) in 83 individuals with VWS grouped within 5 phenotypes for primary analysis: 1=CL+/-P+LP, 2=CL+/-P, 3=CP+LP, 4=CP, 5=LP and 2 phenotypes for secondary analysis: 1=any cleft and LP, 2= any cleft without LP. DNA samples were bisulfite converted and pyrosequenced with target-specific primers. Methylation levels were compared amongst phenotypes.

RESULTS

CpG sites in the IRF6 promoter showed statistically significant differences in methylation among phenotypic groups in both analyses (P<0.05). Individuals with any form of cleft (Groups 1-4) had significantly higher methylation levels than individuals with lip pits only (Group 5). In the secondary analysis, individuals in Group 1 (cleft+LP) had significantly higher methylation than Group 2 (cleft only).

CONCLUSION

Results indicated that hypermethylation of the IRF6 promoter is associated with more severe phenotypes (any cleft +/- lip pits); thus, possibly impacting an already genetically weakened IRF6 protein and leading to a more severe phenotype.

Plasticity-induced repression of Irf6 underlies acquired resistance to cancer immunotherapy

Acquired resistance to immune checkpoint immunotherapy remains a critical yet incompletely understood biological mechanism. Here, using a mouse model of pancreatic ductal adenocarcinoma (PDAC) to study tumor relapse following immunotherapy-induced responses, we found that tumors underwent an epithelial-to-mesenchymal transition (EMT) that resulted in reduced sensitivity to T cell-mediated killing. EMT-transcription factors (EMT-TFs) ZEB1 and SNAIL function as master genetic and epigenetic regulators of this tumor-intrinsic effect. Acquired resistance was not due to immunosuppression in the tumor immune microenvironment, disruptions in the antigen presentation machinery, or altered expression of immune checkpoints. Rather, EMT was associated with epigenetic and transcriptional silencing of interferon regulatory factor 6 (Irf6), which renders tumor cells less sensitive to the pro-apoptotic effects of TNF-α. These findings show how resistance to immunotherapy in PDAC can be acquired through plasticity programs that render tumor cells impervious to T cell killing.

PRMT1 Inhibition Selectively Targets BNC1-Dependent Proliferation, but not Migration in Squamous Cell Carcinoma

Squamous Cell Carcinoma (SCC) develops in stratified epithelial tissues and demonstrates frequent alterations in transcriptional regulators. We sought to discover SCC-specific transcriptional programs and identified the transcription factor Basonuclin 1 (BNC1) as highly expressed in SCC compared to other tumor types. RNA-seq and ChIP-seq analysis identified pro-proliferative genes activated by BNC1 in SCC cells and keratinocytes. Inhibition of BNC1 in SCC cells suppressed proliferation and increased migration via FRA1. In contrast, BNC1 reduction in keratinocytes caused differentiation, which was abrogated by IRF6 knockdown, leading to increased migration. Protein interactome analysis identified PRMT1 as a co-activator of BNC1-dependent proliferative genes. Inhibition of PRMT1 resulted in a dose-dependent reduction in SCC cell proliferation without increasing migration. Importantly, therapeutic inhibition of PRMT1 in SCC xenografts significantly reduced tumor size, resembling functional effects of BNC1 knockdown. Together, we identify BNC1-PRMT1 as an SCC-lineage specific transcriptional axis that promotes cancer growth, which can be therapeutically targeted to inhibit SCC tumorigenesis.

DNA methylation differences in monozygotic twins with Van der Woude syndrome

Introduction

Van der Woude Syndrome (VWS) is an autosomal dominant disorder responsible for 2% of all syndromic orofacial clefts (OFCs) with IRF6 being the primary causal gene (70%). Cases may present with lip pits and either cleft lip, cleft lip with cleft palate, or cleft palate, with marked phenotypic discordance even among individuals carrying the same mutation. This suggests that genetic or epigenetic modifiers may play additional roles in the syndrome's etiology and variability in expression. We report the first DNA methylation profiling of 2 pairs of monozygotic twins with VWS. Our goal is to explore epigenetic contributions to VWS etiology and variable phenotypic expressivity by comparing DNAm profiles in both twin pairs. While the mutations that cause VWS in these twins are known, the additional mechanism behind their phenotypic risk and variability in expression remains unclear.

Methods

We generated whole genome DNAm data for both twin pairs. Differentially methylated positions (DMPs) were selected based on: (1) a coefficient of variation in DNAm levels in unaffected individuals < 20%, and (2) intra-twin pair absolute difference in DNAm levels >5% (delta beta > | 0.05|). We then divided the DMPs in two subgroups for each twin pair for further analysis: (1) higher methylation levels in twin A (Twin A > Twin B); and (2) higher methylation levels in twin B (Twin B >Twin A).

Results and Discussion

Gene ontology analysis revealed a list of enriched genes that showed significant differential DNAm, including clef-associated genes. Among the cleft-associated genes, TP63 was the most significant hit (p=7.82E-12). Both twin pairs presented differential DNAm levels in CpG sites in/near TP63 (Twin 1A > Twin 1B and Twin 2A < Twin 2B). The genes TP63 and IRF6 function in a biological regulatory loop to coordinate epithelial proliferation and differentiation in a process that is critical for palatal fusion. The effects of the causal mutations in IRF6 can be further impacted by epigenetic dysregulation of IRF6 itself, or genes in its pathway. Our data shows evidence that changes in DNAm is a plausible mechanism that can lead to markedly distinct phenotypes, even among individuals carrying the same mutation.

ARHGAP29 is required for keratinocyte proliferation and migration

BACKGROUND

RhoA GTPase plays critical roles in actin cytoskeletal remodeling required for controlling a diverse range of cellular functions including cell proliferation, cell adhesions, migration and changes in cell shape. RhoA cycles between an active GTP-bound and an inactive GDP-bound form, a process that is regulated by guanine nucleotide exchange factors (GEFs), and GTPase-activating proteins (GAPs). ARHGAP29 is a GAP expressed in keratinocytes of the skin and is decreased in the absence of Interferon Regulator Factor 6, a critical regulator of cell proliferation and migration. However, the role for ARHGAP29 in keratinocyte biology is unknown.

RESULTS

Novel ARHGAP29 knockdown keratinocyte cell lines were generated using both CRISPR/Cas9 and shRNA technologies. Knockdown cells exhibited significant reduction of ARHGAP29 protein (50-80%) and displayed increased filamentous actin (stress fibers), phospho-myosin light chain (contractility), cell area and population doubling time. Furthermore, we found that ARHGAP29 knockdown keratinocytes displayed significant delays in scratch wound closure in both single cell and collective cell migration conditions. Particularly, our results show a reduction in path lengths, speed, directionality and persistence in keratinocytes with reduced ARHGAP29. The delay in scratch closure was rescued by both adding back ARHGAP29 or adding a ROCK inhibitor to ARHGAP29 knockdown cells.

CONCLUSIONS

These data demonstrate that ARHGAP29 is required for keratinocyte morphology, proliferation and migration mediated through the RhoA pathway.

The NOTCH-RIPK4-IRF6-ELOVL4 Axis Suppresses Squamous Cell Carcinoma

Receptor-interacting serine/threonine protein kinase 4 (RIPK4) and its kinase substrate the transcription factor interferon regulatory factor 6 (IRF6) play critical roles in the development and maintenance of the epidermis. In addition, ourselves and others have previously shown that RIPK4 is a NOTCH target gene that suppresses the development of cutaneous and head and neck squamous cell carcinomas (HNSCCs). In this study, we used autochthonous mouse models, where the expression of Pik3ca^H1047R oncogene predisposes the skin and oral cavity to tumor development, and show that not only loss of Ripk4, but also loss of its kinase substrate Irf6, triggers rapid SCC development. In vivo rescue experiments using Ripk4 or a kinase-dead Ripk4 mutant showed that the tumor suppressive function of Ripk4 is dependent on its kinase activity. To elucidate critical mediators of this tumor suppressive pathway, we performed transcriptional profiling of Ripk4-deficient epidermal cells followed by multiplexed in vivo CRISPR screening to identify genes with tumor suppressive capabilities. We show that Elovl4 is a critical Notch-Ripk4-Irf6 downstream target gene, and that Elovl4 loss itself triggers SCC development. Importantly, overexpression of Elovl4 suppressed tumor growth of Ripk4-deficient keratinocytes. Altogether, our work identifies a potent Notch1-Ripk4-Irf6-Elovl4 tumor suppressor axis.

Consistent downregulation of the cleft lip/palate-associated genes IRF6 and GRHL3 in carcinomas

Interferon Regulatory Factor 6 (IRF6) and Grainyhead Like Transcription Factor 3 (GRHL3) are transcription factors that orchestrate gene regulatory networks required for the balance between keratinocyte differentiation and proliferation. Absence of either protein results in the lack of a normal stratified epidermis with keratinocytes failing to stop proliferating and to terminally differentiate. Numerous pathological variants within IRF6 and GRHL3 have been identified in orofacial cleft-affected individuals and expression of the two transcription factors has been found to be often dysregulated in cancers. However, whether orofacial cleft-associated IRF6 and GRHL3 variants in patients might also affect their cancer risk later in life, is not clear yet. The fact that the role of IRF6 and GRHL3 in cancer remains controversial makes this question even more challenging. Some studies identified IRF6 and GRHL3 as oncogenes, while others could attribute tumor suppressive functions to them. Trying to solve this apparent conundrum, we herein aimed to characterize IRF6 and GRHL3 function in various types of carcinomas. We screened multiple cancer and normal cell lines for their expression, and subsequently proceeded with functional assays in cancer cell lines. Our data uncovered consistent downregulation of IRF6 and GRHL3 in all types of carcinomas analyzed. Reduced levels of IRF6 and GRHL3 were found to be associated with several tumorigenic properties, such as enhanced cell proliferation, epithelial mesenchymal transition, migration and reduced differentiation capacity. Based on our findings, IRF6 and GRHL3 can be considered as tumor suppressor genes in various carcinomas, which makes them potential common etiological factors for cancer and CLP in a fraction of CLP-affected patients.

SOX9 maintains human foetal lung tip progenitor state by enhancing WNT and RTK signalling

The balance between self-renewal and differentiation in human foetal lung epithelial progenitors controls the size and function of the adult organ. Moreover, progenitor cell gene regulation networks are employed by both regenerating and malignant lung cells, where modulators of their effects could potentially be of therapeutic value. Details of the molecular networks controlling human lung progenitor self-renewal remain unknown. We performed the first CRISPRi screen in primary human lung organoids to identify transcription factors controlling progenitor self-renewal. We show that SOX9 promotes proliferation of lung progenitors and inhibits precocious airway differentiation. Moreover, by identifying direct transcriptional targets using Targeted DamID, we place SOX9 at the centre of a transcriptional network, which amplifies WNT and RTK signalling to stabilise the progenitor cell state. In addition, the proof-of-principle CRISPRi screen and Targeted DamID tools establish a new workflow for using primary human organoids to elucidate detailed functional mechanisms underlying normal development and disease.

ErbB2/Her2-dependent downregulation of a cell death-promoting protein BLNK in breast cancer cells is required for 3D breast tumor growth

A significant proportion of breast cancers are driven by ErbB2/Her2 oncoprotein that they overexpress. These malignancies are typically treated with various ErbB2-targeted drugs, but many such cancers develop resistance to these agents and become incurable. Conceivably, treatment of ErbB2-positive cancers could be facilitated by use of agents blocking oncogenic signaling mechanisms downstream of ErbB2. However, current understanding of these mechanisms is limited. The ability of solid tumor cells to resist anoikis, cell death triggered by cell detachment from the extracellular matrix (ECM), is thought to be critical for 3D tumor growth. In an effort to understand the mechanisms of ErbB2-driven breast cancer cell anoikis resistance we found that detachment of non-malignant breast epithelial cells from the ECM upregulates a cell death-promoting tumor suppressor adapter protein BLNK and that ErbB2 blocks this upregulation by reducing tumor cell levels of transcription factor IRF6. We further observed that trastuzumab, a therapeutic anti-ErbB2 antibody, upregulates BLNK in human trastuzumab-sensitive but not trastuzumab-resistant ErbB2-positive breast cancer cells. Moreover, we established that BLNK promotes anoikis by activating p38 MAP kinase and that ErbB2-dependent BLNK downregulation blocks breast cancer cell anoikis. In search for pharmacological approaches allowing to upregulate BLNK in tumor cells we found that clinically approved proteasome inhibitor bortezomib upregulates IRF6 and BLNK in human breast cancer cells and inhibits their 3D growth in a BLNK-dependent manner. In addition, we found that BLNK upregulation in human ErbB2-positive breast cancer cells blocks their ability to form tumors in mice. Furthermore, we used publicly available data on mRNA levels in multiple breast cancers to demonstrate that increased BLNK mRNA levels correlate with increased relapse-free survival in a cohort of approximately 400 patients with ErbB2-positive breast cancer. In summary, we discovered a novel mechanism of ErbB2-driven 3D breast tumor growth mediated by ErbB2-dependent BLNK downregulation.

Identifying the critical states and dynamic network biomarkers of cancers based on network entropy

Background

There are sudden deterioration phenomena during the progression of many complex diseases, including most cancers; that is, the biological system may go through a critical transition from one stable state (the normal state) to another (the disease state). It is of great importance to predict this critical transition or the so-called pre-disease state so that patients can receive appropriate and timely medical care. In practice, however, this critical transition is usually difficult to identify due to the high nonlinearity and complexity of biological systems.

Methods

In this study, we employed a model-free computational method, local network entropy (LNE), to identify the critical transition/pre-disease states of complex diseases. From a network perspective, this method effectively explores the key associations among biomolecules and captures their dynamic abnormalities.

Results

Based on LNE, the pre-disease states of ten cancers were successfully detected. Two types of new prognostic biomarkers, optimistic LNE (O-LNE) and pessimistic LNE (P-LNE) biomarkers, were identified, enabling identification of the pre-disease state and evaluation of prognosis. In addition, LNE helps to find “dark genes” with nondifferential gene expression but differential LNE values.

Conclusions

The proposed method effectively identified the critical transition states of complex diseases at the single-sample level. Our study not only identified the critical transition states of ten cancers but also provides two types of new prognostic biomarkers, O-LNE and P-LNE biomarkers, for further practical application. The method in this study therefore has great potential in personalized disease diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03445-0.

High-throughput functional evaluation of human cancer-associated mutations using base editors

Comprehensive phenotypic characterization of the many mutations found in cancer tissues is one of the biggest challenges in cancer genomics. In this study, we evaluated the functional effects of 29,060 cancer-related transition mutations that result in protein variants on the survival and proliferation of non-tumorigenic lung cells using cytosine and adenine base editors and single guide RNA (sgRNA) libraries. By monitoring base editing efficiencies and outcomes using surrogate target sequences paired with sgRNA-encoding sequences on the lentiviral delivery construct, we identified sgRNAs that induced a single primary protein variant per sgRNA, enabling linking those mutations to the cellular phenotypes caused by base editing. The functions of the vast majority of the protein variants (28,458 variants, 98%) were classified as neutral or likely neutral; only 18 (0.06%) and 157 (0.5%) variants caused outgrowing and likely outgrowing phenotypes, respectively. We expect that our approach can be extended to more variants of unknown significance and other tumor types.

IRF7 inhibits the Warburg effect via transcriptional suppression of PKM2 in osteosarcoma

Osteosarcoma (OS) is a malignant bone tumor among adolescents and young adults. IRF7 belongs to the transcription factor family of interferon regulatory factors (IRFs) and has previously been described to function as a tumor suppressor in multiple cancer types. However, the biological functions and cellular mechanism of IRF7 in OS remain elusive. In this study, by quantitative real-time PCR (qRT-PCR) and western blotting, we found that IRF7 was downregulated in OS, and the higher expression of IRF7 was correlated with a better survival prognosis. Moreover, loss-of-function and gain-of-function studies have proved the critical functions of IRF7 in suppressing aerobic glycolysis of osteosarcoma cells as evidenced by glucose uptake, lactate production, extracellular acidification rate, and oxygen consumption rate. Mechanistically, IRF7 inhibited the expression of key glycolytic gene PKM2 via direct transcriptional regulation. Moreover, the in vitro and in vivo tumor-suppressive roles of IRF7 were uncovered in OS and these effects were largely glycolysis-dependent. Therefore, our study unveils a previous unprecedented role of IRF7 in glucose metabolism reprogram and suggests that IRF7 might serve as a potential therapeutic target for patients with OS.

RIPK4 regulates cell–cell adhesion in epidermal development and homeostasis

Epidermal development and maintenance are finely regulated events requiring a strict balance between proliferation and differentiation. Alterations in these processes give rise to human disorders such as cancer or syndromes with skin and annexes defects, known as ectodermal dysplasias (EDs). Here, we studied the functional effects of two novel receptor-interacting protein kinase 4 (RIPK4) missense mutations identified in siblings with an autosomal recessive ED with cutaneous syndactyly, palmoplantar hyperkeratosis and orofacial synechiae. Clinical overlap with distinct EDs caused by mutations in transcription factors (i.e. p63 and interferon regulatory factor 6, IRF6) or nectin adhesion molecules was noticed. Impaired activity of the RIPK4 kinase resulted both in altered epithelial differentiation and defective cell adhesion. We showed that mutant RIPK4 resulted in loss of PVRL4/nectin-4 expression in patient epidermis and primary keratinocytes, and demonstrated that PVRL4 is transcriptionally regulated by IRF6, a RIPK4 phosphorylation target. In addition, defective RIPK4 altered desmosome morphology through modulation of plakophilin-1 and desmoplakin. In conclusion, this work implicates RIPK4 kinase function in the p63-IRF6 regulatory loop that controls the proliferation/differentiation switch and cell adhesion, with implications in ectodermal development and cancer.

IRF6 Regulates the Delivery of E-Cadherin to the Plasma Membrane

IRF6 is a transcription factor that is required for craniofacial development and epidermal morphogenesis. Specifically, Irf6-deficient mice lack the terminally differentiated epidermal layers, leading to an absence of barrier function. This phenotype also includes intraoral adhesions due to the absence of the oral periderm, leading to the mislocalization of E-cadherin and other cell‒cell adhesion proteins of the oral epithelium. However, the mechanisms by which IRF6 controls the localization of cell adhesion proteins are not understood. In this study, we show that in human and murine keratinocytes, loss of IRF6 leads to a breakdown of epidermal sheets after mechanical stress. This defect is due to a reduction of adhesion proteins at the plasma membrane. Dynamin inhibitors rescued the IRF6-dependent resistance of epidermal sheets to mechanical stress, but only inhibition of clathrin-mediated endocytosis rescued the localization of junctional proteins at the membrane. Our data show that E-cadherin recycling but not its endocytosis is affected by loss of IRF6. Overall, we demonstrate a role for IRF6 in the delivery of adhesion proteins to the cell membrane.

Update on Management Recommendations for Advanced Cutaneous Squamous Cell Carcinoma

Cutaneous squamous cell carcinoma (cSCC) is the second most common form of skin cancer, the incidence of which has risen over the last years. Although cSCC rarely metastasizes, early detection and treatment of primary tumours are critical to limit progression and local invasion. Several prognostic factors related to patients' clinicopathologic profile and tumour features have been identified as high-risk markers and included in the stratification scales, but their association with regional control or survival is uncertain. Therefore, decision-making on the diagnosis and management of cSCC should be made based on each individual patient's characteristics. Recent advances in non-invasive imaging techniques and molecular testing have enhanced clinical diagnostic accuracy. Surgical excision is the mainstay of local treatment, whereas radiotherapy (RT) is recommended for patients with inoperable disease or in specific circumstances. Novel systemic treatments including immunotherapies and targeted therapies have changed the therapeutic landscape for cSCC. The anti-PD-1 agent cemiplimab is currently the only FDA/EMA-approved first-line therapy for patients with locally advanced or metastatic cSCC who are not candidates for curative surgery or RT. Given the likelihood of recurrence and the increased risk of developing multiple cSCC, close follow-up should be performed during the first years of treatment and continued long-term surveillance is warranted.

Comprehensive analysis of expression profile and prognostic significance of interferon regulatory factors in pancreatic cancer

BACKGROUND

Pancreatic cancer (PC) is a highly lethal disease and an increasing cause of cancer-associated mortality worldwide. Interferon regulatory factors (IRFs) play vital roles in immune response and tumor cellular biological processes. However, the specific functions of IRFs in PC and tumor immune response are far from systematically clarified. This study aimed to explorer the expression profile, prognostic significance, and biological function of IRFs in PC.

RESULTS

We observed that the levels of IRF2, 6, 7, 8, and 9 were elevated in tumor compared to normal tissues in PC. IRF7 expression was significantly associated with patients' pathology stage in PC. PC patients with high IRF2, low IRF3, and high IRF6 levels had significantly poorer overall survival. High mRNA expression, amplification and, deep deletion were the three most common types of genetic alterations of IRFs in PC. Low expression of IRF2, 4, 5, and 8 was resistant to most of the drugs or small molecules from Genomics of Drug Sensitivity in Cancer. Moreover, IRFs were positively correlated with the abundance of tumor infiltrating immune cells in PC, including B cells, CD8+ T cells, CD4+ T cells, macrophages, Neutrophil, and Dendritic cells. Functional analysis indicated that IRFs were involved in T cell receptor signaling pathway, immune response, and Toll-like receptor signaling pathway.

CONCLUSIONS

Our results indicated that certain IRFs could serve as potential therapeutic targets and prognostic biomarkers for PC patients. Further basic and clinical studies are needed to validate our findings and generalize the clinical application of IRFs in PC.

Decreased interferon regulatory factor 6 expression due to DNA hypermethylation predicts an unfavorable prognosis in clear cell renal cell carcinoma

Background: Emerging evidences have indicated that IRF6, as a member of the Interferon regulatory factors (IRFs) family, plays important roles in a variety of tumors. However, the expression status of IRF6 and its prognostic value in clear cell renal cell carcinoma (ccRCC) remain unclear. Methods: In this study, we used TCGA-KIRC, GEO and TIP databases and immunohistochemistry staining to determine the expression profile, clinico-pathological features and prognostic value of IRF6 in ccRCC. MSP and demethylation analysis were utilized to verify the regulatory effect of DNA methylation on IRF6 expression. Results: Our results found that IRF6 expression was downregulated in ccRCC tissues and cell lines, and decreased IRF6 expression was associated with worse clinicopathological features and poorer prognosis. Besides, the results of multivariate Cox regression analysis also confirmed that decreased IRF6 expression was an independently risk factor predictor of shorter Overall Survival (OS) (HR: 0.8524, 95%CI: 0.7614-0.9543, P=0.0056) and Disease Free Survival (DFS) (HR: 0.7024, 95%CI: 0.6087-0.8104, P<0.0001) in ccRCC patients. Moreover, the results of MSP and demethylation analysis validated that decreased IRF6 expression was caused by DNA hypermethylation. Furthermore, our results showed that IRF6 expression was associated with the infiltration levels of multiple immune cells in ccRCC. Conclusions: These findings demonstrated that IRF6 expression was significantly reduced in ccRCC and DNA hypermethylation played an important role in decreased IRF6 expression. In addition, the decrease of IRF6 was related to the unfavorable prognosis of ccRCC patients and the alterations of tumor immune cells infiltration.

Lack of IRF6 Disrupts Human Epithelial Homeostasis by Altering Colony Morphology, Migration Pattern, and Differentiation Potential of Keratinocytes

Variants within the gene encoding for the transcription factor Interferon Regulatory Factor 6 (IRF6) are associated with syndromic and non-syndromic Cleft Lip/Palate (CLP) cases. IRF6 plays a vital role in the regulation of the proliferation/differentiation balance in keratinocytes and is involved in wound healing and migration. Since a fraction of CLP patients undergoing corrective cleft surgery experience wound healing complications, IRF6 represents an interesting candidate gene linking the two processes. However, Irf6 function has been mainly studied in mice and knowledge on IRF6 in human cells remains sparse. Here, we aimed to elucidate the role of IRF6 in human postnatal skin- and oral mucosa-derived keratinocytes. To do so, we applied CRISPR/Cas9 to ablate IRF6 in two TERT-immortalized keratinocyte cultures, which we used as model cell lines. We show that IRF6 controls the appearance of single cells and colonies, with the latter being less cohesive in its absence. Consequently, IRF6 knockout keratinocytes often moved as single cells instead of a collective epithelial sheet migration but maintained their epithelial character. Lack of IRF6 triggered severe keratinocyte differentiation defects, which were already apparent in the stratum spinosum and extended to the stratum corneum in 3D organotypic skin cultures, while it did not alter their growth rate. Finally, proteomics revealed that most of the differentially expressed proteins in the absence of IRF6 could be associated with differentiation, cell-cell adhesion as well as immune response. Our data expand the knowledge on IRF6 in human postnatal keratinocytes, which will help to better understand IRF6-related pathologies.

Association analysis of SNPs in GRHL3, FAF1, and KCNJ2 with NSCPO sub-phenotypes in Han Chinese

OBJECTIVES

Non-syndromic cleft palate only (NSCPO) is a common congenital deformity with complex etiologies. GRHL3, FAF1, and KCNJ2 have been reported to be involved in the pathogenesis of NSCPO. Up till now, there have been no replication studies based on large Han Chinese. Therefore, this study aimed to investigate associations between GRHL3, FAF1, KCNJ2, and NSCPO sub-phenotypes patients in Han Chinese.

MATERIALS AND METHODS

Firstly, we selected 2 SNPs based on previous literatures: FAF1 (rs3827730) and GRHL3 (rs41268753). Also, we selected 8 tagSNPs in GRHL3 (rs557811, rs609352, rs10903078, rs6659209, rs12401714, rs12568599, rs3887581, rs12024148) and 2 tagSNPs in KCNJ2 (rs75855040 and rs236514). Afterward, we evaluated these SNPs among 1668 NSCPO patients and 1811 normal controls from Han Chinese. Following data were analyzed by PLINK and Haploview program.

RESULTS

Association analysis under additive model showed that allele A at rs12568599 in GRHL3 gene is significantly associated with NSCPO (p = 0.0034, OR = 1.38 and 95%CI: 1.11-1.72) and its sub-phenotype incomplete cleft palate (ICP) (p = 0.0039, OR = 1.4 and 95%CI: 1.11-1.75), and it could increase the risk of both NSCPO and ICP.

CONCLUSIONS

This study firstly found that rs12568599 in GRHL3 is associated with NSCPO and ICP in Han Chinese, indicating that sub-phenotypes of NSCPO have different etiologies.

Inhibition of KIF20A by transcription factor IRF6 affects the progression of renal clear cell carcinoma

Background

Renal clear cell carcinoma (ccRCC) is one of the most common malignant tumors, whose incidence is increasing year by year. IRF6 plays an important role in the occurrence of tumors, although there is yet no report on its expression in ccRCC.

Methods

The expression of IRF6 and KIF20A in ccRCC was predicted by GEPIA and HAP databases. In addition, GEPIA database predicted the relationship between IRF6 and KIF20A expressions and the pathological staging, overall survival, and disease-free survival of ccRCC. The possible binding sites of IRF6 and KIF20A promoters were predicted by JASPAR database and verified by luciferase and ChIP assays. The specific effects of IRF6 on ccRCC cell proliferation, invasion and apoptosis were subsequently examined at both cellular level and animal level.

Results

The database predicted down-regulated IRF6 expression in renal carcinoma tissues and its correlation with poor prognosis. IRF6 overexpression inhibited cRCC cell proliferation, invasion and migration. In addition, up-regulated KIF20A expression in renal carcinoma tissues and its association with prognosis were also predicted. Interference with KIF20A inhibited the proliferation, invasion, and migration of ccRCC cells. Finally, we confirmed that KIF20A is a functional target of IRF6 and can partially reverse the effects of IRF6 on the proliferation, invasion and migration of ccRCC cells. Conclusion: Inhibition of KIF20A by transcription factor IRF6 affects cell proliferation, invasion and migration in renal clear cell carcinoma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01879-y.

Cyproheptadine, an epigenetic modifier, exhibits anti-tumor activity by reversing the epigenetic silencing of IRF6 in urothelial carcinoma

BACKGROUND

Urothelial carcinoma (UC) is the second most common malignancy of the urinary system with high rate of recurrence, UC patients therefore needed to be treated with surgery followed by chemotherapy. Development of novel therapeutics with minimal side-effect is an urgent issue. Our previous study showed that cyproheptadine (CPH), an anti-histamine, exhibited antitumor activity in UC in vitro and in an xenograft model. However, the molecular mechanism of how CPH inhibits tumor progression is not fully understood.

METHODS

Genes that were upregulated after treatment with CPH in UC cells, were examined by RNA-Seq. Real-time quantitative PCR (RT-qPCR) was employed to detect IRF6 expression while COBRA assay and bisulphite pyrosequencing were used to examine promoter methylation of IRF6. Enrichment of total H3K27 acetylation and H3K4 mono-methylation were detected by western blotting. Colony formation and flow cytometry were used to examine proliferation and apoptosis in UC cells overexpressed or depleted with IRF6. Nude mice xenograft model was used to examine the effect of IRF6 in UC.

RESULTS

Our result showed that several genes, including IRF6 were upregulated after treatment with CPH in BFTC905 UC cells. Further experiments found that treatment of CPH could restore the expression of IRF6 in several other UC cell lines, probably due to promoter hypomethylation and enrichment of H3K27 acetylation and H3K4 mono-methylation. These results may be due to the fact that CPH could alter the activity, but not the expression of epigenetic modifiers. Finally, re-expression of IRF6 in UC inhibited tumor growth in vitro and in an xenograft mouse model, by inducing apoptosis.

CONCLUSION

In conclusion, our results suggested that CPH may be an epigenetic modifier, modulating the expression of the potential tumor suppressor IRF6, in inhibiting tumor growth in UC.

Tumor levels of the mediators of ErbB2-driven anoikis resistance correlate with breast cancer relapse in patients receiving trastuzumab-based therapies

PURPOSE

Patients with ErbB2/Her2 oncoprotein-positive breast cancers often receive neoadjuvant therapies (NATs) containing the anti-ErbB2 antibody trastuzumab. Tumors that are still present after NATs are resected, and patients continue receiving trastuzumab. These cancers are associated with high relapse risk. Whether relapse will occur cannot be presently reliably predicted. The ability to make such predictions could improve disease management. We found previously that ErbB2 blocks breast tumor cell anoikis, apoptosis induced by cell detachment from the extracellular matrix, by downregulating the pro-apoptotic protein Irf6 and upregulating the anti-apoptotic protein Epidermal Growth Factor Receptor (EGFR) in the cells and, thus, promotes their three-dimensional growth. We now tested whether tumor levels of these proteins before and after NATs correlate with patients' relapse-free survival (RFS) and overall survival (OS).

METHODS

We selected archival breast tumor samples collected from 37 women with ErbB2-positive stages II and III breast cancer before and after NATs. We used immunohistochemistry to test whether levels of the indicated proteins in respective tumors correlate with RFS and OS.

RESULTS

We observed that the presence of high Irf6 levels in the tumors following NATs correlated with reduced RFS and OS. Perhaps not by coincidence, we noticed that trastuzumab-sensitive ErbB2-positive breast cancer cells selected for the ability to overproduce exogenous Irf6 in culture acquired trastuzumab resistance. Finally, EGFR presence in patients' tumors before or after NATs was associated with decreased RFS and OS.

CONCLUSIONS

This study could help identify patients with ErbB2-positive tumors that are at increased risk of disease relapse following NATs.

Interferon regulatory factor 6 correlates with the progression of non-small cell lung cancer and can be regulated by miR-320

OBJECTIVES

The expression of interferon regulatory factor 6 (IRF6) has been reported in several cancer types, but its roles underlying the progression of lung cancer have not been detailedly investigated.

METHODS

The pairs of lung cancer tissues and para-carcinoma tissues and The Cancer Genome Atlas database were collected to detect IRF6 expression. Cell counting kit-8, transwell and terminal-deoxynucleoitidyl transferase-mediated nick end labelling assays were used to evaluate cell proliferation, migration and apoptosis.

KEY FINDINGS

A significant up-regulation of IRF6 in both lung adenocarcinoma and lung squamous cell carcinoma tissues compared with normal non-tumor tissues. Subsequently, Immunostaining also revealed that canceration of lung tissues predisposed to evoke IRF6 expression. In vitro experiments revealed the antitumour effects, including growth and migration inhibition, of IRF6 siRNA transfection. Considering miR-320 as an endogenous inhibitor to IRF6, miR-320 mimics transfection led to the inhibition of proliferation and migration of lung cancer cells. However, overexpressed IRF6 neutralized the antineoplastic activities of miR-320 in lung cancer cells.

CONCLUSIONS

The miR-320/IRF6 signalling axis was implicated in pulmonary canceration. miR-320 as an endogenous inhibitor of IRF6 provided a novel therapeutic strategy for the treatment of lung cancer.

Delineating the roles of Grhl2 in craniofacial development through tissue-specific conditional deletion and epistasis approaches in mouse

BACKGROUND

The highly conserved Grainyhead-like (Grhl) family of transcription factors play critical roles in the development of the neural tube and craniofacial skeleton. In particular, deletion of family member Grainyhead-like 2 (Grhl2) leads to mid-gestational embryonic lethality, maxillary clefting, abdominoschisis, and both cranial and caudal neural tube closure defects. These highly pleiotropic and systemic defects suggest that Grhl2 plays numerous critical developmental roles to ensure correct morphogenesis and patterning.

RESULTS

Here, using four separate Cre-lox conditional deletion models, as well as one genetic epistasis approach (Grhl2^+/- ;Edn1^+/- double heterozygous mice) we have investigated tissue-specific roles of Grhl2 in embryonic development, with a particular focus on the craniofacial skeleton. We find that loss of Grhl2 in the pharyngeal epithelium (using the Shh^Cre driver) leads to low-penetrance micrognathia, whereas deletion of Grhl2 within the ectoderm of the pharynx (Nestin^Cre ) leads to small, albeit significant, differences in the proximal-distal elongation of both the maxilla and mandible. Loss of Grhl2 in endoderm (Sox17-2a^iCre ) resulted in noticeable lung defects and a single instance of secondary palatal clefting, although formation of other endoderm-derived organs such as the stomach, bladder and intestines was not affected. Lastly, deletion of Grhl2 in cells of the neural crest (Wnt1^Cre ) did not lead to any discernible defects in craniofacial development, and similarly, our epistasis approach did not detect any phenotypic consequences of loss of a single allele of both Grhl2 and Edn1.

CONCLUSION

Taken together, our study identifies a pharyngeal-epithelium intrinsic, non-cell-autonomous role for Grhl2 in the patterning and formation of the craniofacial skeleton, as well as an endoderm-specific role for Grhl2 in the formation and establishment of the mammalian lung.

Identification of a dynamic gene regulatory network required for pluripotency factor-induced reprogramming of mouse fibroblasts and hepatocytes

The generation of induced pluripotent stem cells (iPSCs) from somatic cells provides an excellent model to study mechanisms of transcription factor-induced global alterations of the epigenome and genome function. Here, we have investigated the early transcriptional events of cellular reprogramming triggered by the co-expression of Oct4, Sox2, Klf4, and c-Myc (OSKM) in mouse embryonic fibroblasts (MEFs) and mouse hepatocytes (mHeps). In this analysis, we identified a gene regulatory network composed of nine transcriptional regulators (9TR; Cbfa2t3, Gli2, Irf6, Nanog, Ovol1, Rcan1, Taf1c, Tead4, and Tfap4), which are directly targeted by OSKM, in vivo. Functional studies using single and double shRNA knockdowns of any of these factors caused disruption of the network and dramatic reductions in reprogramming efficiency, indicating that this network is essential for the induction and establishment of pluripotency. We demonstrate that the stochastic co-expression of 9TR network components occurs in a remarkably small number of cells, approximating the percentage of terminally reprogrammed cells as a result of dynamic molecular events. Thus, the early DNA-binding patterns of OSKM and the subsequent probabilistic co-expression of essential 9TR components in subpopulations of cells undergoing reprogramming steer the reconstruction of a gene regulatory network marking the transition to pluripotency.

An Irf6-Esrp1/2 regulatory axis controls midface morphogenesis in vertebrates

Irf6 and Esrp1 are important for palate development across vertebrates. In zebrafish, we found that irf6 regulates the expression of esrp1 We detailed overlapping Irf6 and Esrp1/2 expression in mouse orofacial epithelium. In zebrafish, irf6 and esrp1/2 share expression in periderm, frontonasal ectoderm and oral epithelium. Genetic disruption of irf6 and esrp1/2 in zebrafish resulted in cleft of the anterior neurocranium. The esrp1/2 mutant also developed cleft of the mouth opening. Lineage tracing of cranial neural crest cells revealed that the cleft resulted not from migration defect, but from impaired chondrogenesis. Analysis of aberrant cells within the cleft revealed expression of sox10, col1a1 and irf6, and these cells were adjacent to krt4+ and krt5+ cells. Breeding of mouse Irf6; Esrp1; Esrp2 compound mutants suggested genetic interaction, as the triple homozygote and the Irf6; Esrp1 double homozygote were not observed. Further, Irf6 heterozygosity reduced Esrp1/2 cleft severity. These studies highlight the complementary analysis of Irf6 and Esrp1/2 in mouse and zebrafish, and identify a unique aberrant cell population in zebrafish expressing sox10, col1a1 and irf6 Future work characterizing this cell population will yield additional insight into cleft pathogenesis.

miR-587 promotes cervical cancer by repressing interferon regulatory factor 6

BACKGROUND

Interferon regulatory factor 6 (IRF6) exhibits tumor-suppressive functions in several cancer types. In the present study, the antitumor properties and related pathway mechanism of IRF6 were investigated in cervical cancer.

METHODS

Forty-one pairs of cervical cancer specimens and para-carcinoma tissues were collected to evaluate IRF6 expression using immunohistochemical staining and miR-587. The effects of miR-587 and IRF6 on cervical cancer cell growth were explored by MTT assays and in a HeLa tumor xenograft mouse model. The migration and invasion of cervical cancer cells were monitored using transwell assays.

RESULTS

IRF6 expression in cervical cancer specimens and cell lines was significantly reduced compared to that in the corresponding control group. In addition, IRF6 expression was negatively correlated with miR-587 in cervical cancer tissues. Bioinformatics algorithms and luciferase assays revealed that IRF6 is a potential target of miR-587, and miR-587 mimic transfection led to a significant repression of IRF6 protein levels in cervical cancer cells. We also discovered that the antineoplastic properties of IRF6 could be reversed by overexpressing miR-587 in cervical cancer cells. The up-regulation of miR-587 was correlated with poor overall survival in cervical cancer. In an in vivo experiment, miR-587 silencing induced HeLa tumor growth inhibition, which was associated with the up-regulation of IRF6 protein in the tumor.

CONCLUSIONS

miR-587 post-transcriptionally represses IRF6 protein expression to abrogate the antineoplastic activity of IRF6. The miR-587/IRF6 signaling pathway plays a crucial role in the progression of cervical cancer and serves as a potential therapeutic target for the treatment of cervical cancer.

A Novel Van der Woude Syndrome-Causing IRF6 Variant Is Subject to Incomplete Non-sense-Mediated mRNA Decay Affecting the Phenotype of Keratinocytes

Van der Woude syndrome (VWS) is a genetic syndrome that leads to typical phenotypic traits, including lower lip pits and cleft lip/palate (CLP). The majority of VWS-affected patients harbor a pathogenic variant in the gene encoding for the transcription factor interferon regulatory factor 6 (IRF6), a crucial regulator of orofacial development, epidermal differentiation and tissue repair. However, most of the underlying mechanisms leading from pathogenic IRF6 gene variants to phenotypes observed in VWS remain poorly understood and elusive. The availability of one VWS individual within our cohort of CLP patients allowed us to identify a novel VWS-causing IRF6 variant and to functionally characterize it. Using VWS patient-derived keratinocytes, we reveal that most of the mutated IRF6_VWS transcripts are subject to a non-sense-mediated mRNA decay mechanism, resulting in IRF6 haploinsufficiency. While moderate levels of IRF6_VWS remain detectable in the VWS keratinocytes, our data illustrate that the IRF6_VWS protein, which lacks part of its protein-binding domain and its whole C-terminus, is noticeably less stable than its wild-type counterpart. Still, it maintains transcription factor function. As we report and characterize a so far undescribed VWS-causing IRF6 variant, our results shed light on the physiological as well as pathological role of IRF6 in keratinocytes. This acquired knowledge is essential for a better understanding of the molecular mechanisms leading to VWS and CLP.

Insight Into the Function of RIPK4 in Keratinocyte Differentiation and Carcinogenesis

The receptor-interacting protein kinase 4 (RIPK4), a member of the RIPK family, was originally described as an interaction partner of protein kinase C (PKC) β and PKCδ. RIPK4 is identified as a key regulator of keratinocyte differentiation, cutaneous inflammation, and cutaneous wound repair. The mechanism by which RIPK4 integrates upstream signals to initiate specific responses remains elusive. Previous studies have indicated that RIPK4 can regulate several signaling pathways, including the NF-κB, Wnt/β-catenin, and RAF/MEK/ERK pathways. Furthermore, RIPK4-related biological signaling pathways interact with each other to form a complex network. Mounting evidence suggests that RIPK4 is aberrantly expressed in various kinds of cancers. In several types of squamous cell carcinoma (SCC), the mutations that drive aggressive SCC have been found in RIPK4. In addition, the function of RIPK4 in carcinogenesis is probably tissue-specific, since RIPK4 can play a dual role as both a tumor promoter and a tumor suppressor in different tumor types. Therefore, RIPK4 may represent as an independent prognostic factor and a promising novel therapeutic target, which can be used to identify the risks of patients and guide personalized treatments. In future, RIPK4-interacting pathways and precise molecular targets need to be investigated in order to further elucidate the mechanisms underlying epidermal differentiation and carcinogenesis.

Non-random distribution of deleterious mutations in the DNA and protein-binding domains of IRF6 are associated with Van Der Woude syndrome

BACKGROUND

The development of the face occurs during the early days of intrauterine life by the formation of facial processes from the first Pharyngeal arch. Derangement in these well-organized fusion events results in Orofacial clefts (OFC). Van der Woude syndrome (VWS) is one of the most common causes of syndromic cleft lip and/or palate accounting for 2% of all cases. Mutations in the IRF6 gene account for 70% of cases with the majority of these mutations located in the DNA-binding (exon 3, 4) or protein-binding domains (exon 7-9). The current study was designed to update the list of IRF6 variants reported for VWS by compiling all the published mutations from 2013 to date as well as including the previously unreported VWS cases from Africa and Puerto Rico.

METHODS

We used PubMed with the search terms; "Van der Woude syndrome," "Popliteal pterygium syndrome," "IRF6," and "Orofacial cleft" to identify eligible studies. We compiled the CADD score for all the mutations to determine the percentage of deleterious variants.

RESULTS

Twenty-one new mutations were identified from nine papers. The majority of these mutations were in exon 4. Mutations in exon 3 and 4 had CADD scores between 20 and 30 and mutations in exon 7-9 had CADD scores between 30 and 40. The presence of higher CADD scores in the protein-binding domain (exon 7-9) further confirms the crucial role played by this domain in the function of IRF6. In the new cases, we identified five IRF6 mutations, three novel missense mutations (p.Phe36Tyr, p.Lys109Thr, and p.Gln438Leu), and two previously reported nonsense mutations (p.Ser424*and p.Arg250*).

CONCLUSION

Mutations in the protein and DNA-binding domains of IRF6 ranked among the top 0.1% and 1% most deleterious genetic mutations, respectively. Overall, these findings expand the range of VWS mutations and are important for diagnostic and counseling purposes.

Association Analysis of Driver Gene-Related Genetic Variants Identified Novel Lung Cancer Susceptibility Loci with 20,871 Lung Cancer Cases and 15,971 Controls

BACKGROUND

A substantial proportion of cancer driver genes (CDG) are also cancer predisposition genes. However, the associations between genetic variants in lung CDGs and the susceptibility to lung cancer have rarely been investigated.

METHODS

We selected expression-related single-nucleotide polymorphisms (eSNP) and nonsynonymous variants of lung CDGs, and tested their associations with lung cancer risk in two large-scale genome-wide association studies (20,871 cases and 15,971 controls of European descent). Conditional and joint association analysis was performed to identify independent risk variants. The associations of independent risk variants with somatic alterations in lung CDGs or recurrently altered pathways were investigated using data from The Cancer Genome Atlas (TCGA) project.

RESULTS

We identified seven independent SNPs in five lung CDGs that were consistently associated with lung cancer risk in discovery (P < 0.001) and validation (P < 0.05) stages. Among these loci, rs78062588 in TPM3 (1q21.3) was a new lung cancer susceptibility locus (OR = 0.86, P = 1.65 × 10-6). Subgroup analysis by histologic types further identified nine lung CDGs. Analysis of somatic alterations found that in lung adenocarcinomas, rs78062588[C] allele (TPM3 in 1q21.3) was associated with elevated somatic copy number of TPM3 (OR = 1.16, P = 0.02). In lung adenocarcinomas, rs1611182 (HLA-A in 6p22.1) was associated with truncation mutations of the transcriptional misregulation in cancer pathway (OR = 0.66, P = 1.76 × 10-3).

CONCLUSIONS

Genetic variants can regulate functions of lung CDGs and influence lung cancer susceptibility.

IMPACT

Our findings might help unravel biological mechanisms underlying lung cancer susceptibility.

Lin28A promotes IRF6-regulated aerobic glycolysis in glioma cells by stabilizing SNHG14

Warburg effect is a hallmark of cancer cells, wherein glycolysis is preferred over oxidative phosphorylation even in aerobic conditions. Reprogramming of glycometabolism is especially crucial for malignancy in glioma. RNA-binding proteins and long noncoding RNAs are important for aerobic glycolysis during malignant transformation. Thus, we determined the expression and function of RNA-binding protein Lin28A, long noncoding RNA SNHG14, and transcription factor IRF6 in human glioma cells to elucidate the mechanism(s) underlying their role in glycolysis. Quantitative real-time polymerase chain reaction and western blotting showed that Lin28A and SNHG14 were overexpressed and IRF6 was downregulated in glioma. Depleting Lin28A from cells decreased the stability and expression of SNHG14. Furthermore, depleting SNHG14 reduced IRF6 mRNA degradation by targeting its 3' untranslated region and inhibiting STAU1-mediated degradation, thereby increasing the expression of IRF6. PKM2 is an important enzyme in aerobic glycolysis, and GLUT1 is the primary transporter that facilitates glucose uptake. IRF6 inhibited the transcription of PKM2 and GLUT1, thereby impairing glycolysis and cell proliferation and inducing apoptosis in glioma. Notably, depleting Lin28A and SNHG14 and overexpressing IRF6 reduced the growth of xenograft tumors in vivo and prolonged the survival of nude mice. Taken together, our data revealed that the Lin28A/SNHG14/IRF6 axis is crucial for reprogramming glucose metabolism and stimulating tumorigenesis in glioma cells. Thus, targeting this axis might help in the development of a novel therapeutic strategy for glioma metabolism.

SPECC1L regulates palate development downstream of IRF6

SPECC1L mutations have been identified in patients with rare atypical orofacial clefts and with syndromic cleft lip and/or palate (CL/P). These mutations cluster in the second coiled-coil and calponin homology domains of SPECC1L and severely affect the ability of SPECC1L to associate with microtubules. We previously showed that gene-trap knockout of Specc1l in mouse results in early embryonic lethality. We now present a truncation mutant mouse allele, Specc1lΔC510, that results in perinatal lethality. Specc1lΔC510/ΔC510 homozygotes showed abnormal palate rugae but did not show cleft palate. However, when crossed with a gene-trap allele, Specc1lcGT/ΔC510 compound heterozygotes showed a palate elevation delay with incompletely penetrant cleft palate. Specc1lcGT/ΔC510 embryos exhibit transient oral epithelial adhesions at E13.5, which may delay shelf elevation. Consistent with oral adhesions, we show periderm layer abnormalities, including ectopic apical expression of adherens junction markers, similar to Irf6 hypomorphic mutants and Arhgap29 heterozygotes. Indeed, SPECC1L expression is drastically reduced in Irf6 mutant palatal shelves. Finally, we wanted to determine if SPECC1L deficiency also contributed to non-syndromic (ns) CL/P. We sequenced 62 Caucasian, 89 Filipino, 90 Ethiopian, 90 Nigerian and 95 Japanese patients with nsCL/P and identified three rare coding variants (p.Ala86Thr, p.Met91Iso and p.Arg546Gln) in six individuals. These variants reside outside of SPECC1L coiled-coil domains and result in milder functional defects than variants associated with syndromic clefting. Together, our data indicate that palate elevation is sensitive to deficiency of SPECC1L dosage and function and that SPECC1L cytoskeletal protein functions downstream of IRF6 in palatogenesis.

Loss of S100A14 expression at the tumor-invading front correlates with poor differentiation and worse prognosis in oral squamous cell carcinoma

BACKGROUND

We previously showed a tumor-suppressive function of S100A14 in oral squamous cell carcinoma (OSCC). This study aimed to examine the prognostic significance and differentiation-related function of S100A14 in OSCC.

METHODS

S100A14 expression was examined in 170 OSCCs from Norwegian and Nepalese populations using immunohistochemistry. Pro-differentiation function was investigated by overexpressing and silencing S100A14 expression in OSCC-derived cells. External transcriptomic datasets were used to validate association between S100A14 and differentiation markers in OSCC.

RESULT

Loss of S100A14 expression at the invading tumor fronts significantly correlated with poor differentiation and reduced 10-years survival of OSCC-patients. Multivariate Cox analysis identified S100A14 to be an independent prognostic factor. Modulation of S100A14 expression in OSCC-derived cells positively correlated with the expression of differentiation markers. Analysis of external datasets supported the pro-differentiation function of S100A14.

CONCLUSION

These results indicate that S100A14 is a pro-differentiation protein and its expression might be useful as a prognostic marker in OSCC.

Interrogating the Grainyhead-like 2 (Grhl2) genomic locus identifies an enhancer element that regulates palatogenesis in mouse

The highly-conserved Grainyhead-like (Grhl) transcription factors are critical regulators of embryogenesis that regulate cellular survival, proliferation, migration and epithelial integrity, especially during the formation of the craniofacial skeleton. Family member Grhl2 is expressed throughout epithelial tissues during development, and loss of Grhl2 function leads to significant defects in neurulation, abdominal wall closure, formation of the face and fusion of the maxilla/palate. Whereas numerous downstream target genes of Grhl2 have been identified, very little is known about how this crucial developmental transcription factor itself is regulated. Here, using in silico and in utero expression analyses and functional deletion in mice, we have identified a novel 2.4 ​kb enhancer element (mm1286) that drives reporter gene expression in a pattern that strongly recapitulates endogenous Grhl2 in the craniofacial primordia, modulates Grhl2 expression in these tissues, and augments Grhl2-mediated closure of the secondary palate. Deletion of this genomic element, in the context of inactivation of one allele of Grhl2 (through generation of double heterozygous Grhl2^+/-;mm1286^+/- mice), results in a significant predisposition to palatal clefting at birth. Moreover, we found that a highly conserved 325 bp region of mm1286 is both necessary and sufficient for mediating the craniofacial-specific enhancer activity of this region, and that an extremely well-conserved 12-bp sequence within this element (CTGTCAAACAGGT) substantially determines full enhancer function. Together, these data provide valuable new insights into the upstream genomic regulatory landscape responsible for transcriptional control of Grhl2 during palatal closure.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.