Grainyhead-like (Grhl) Target Genes in Development and Cancer

Expression Patterns of Grainyhead-Like 2 and Ovo-Like 2 in Mouse Mammary Gland Development During Pregnancy, Lactation, and Weaning

Grainyhead-like 2 (Grhl2) is a transcription factor that regulates cell adhesion genes in mammary ductal development and serves as a repressor of the epithelial-mesenchymal transition. Conversely, Ovo-like2 (Ovol2) is a target gene of Grhl2 but functions as a substitute in Grhl2-deficient mice, facilitating successful epithelial barrier formation and lumen expansion in kidney-collecting ductal epithelial cells. Our objective was to examine the expression patterns of Grhl2, Ovol2, and their associated genes during the intricate phases of mouse mammary gland development. The mRNA expression of Grhl2 and Ovol2 increased after pregnancy. We observed Grhl2 protein presence in the epithelial cell's region, coinciding with acini formation, and its signal significantly correlated with E-cadherin (Cdh1) expression. However, Ovol2 was present in the epithelial region without a correlation with Cdh1. Similarly, Zeb1, a mesenchymal transcription factor, showed Cdh1-independent expression. Subsequently, we explored the interaction between Rab25, a small G protein, and Grhl2/Ovol2. The expressions of Grhl2 and Ovol2 exhibited a strong correlation with Rab25 and claudin-4, a tight junction protein. These findings suggest that Grhl2 and Ovol2 may collaborate to regulate genes associated with cell adhesion and are crucial for maintaining epithelial integrity during the different phases of mammary gland development.

Grainyhead-like-2, an epithelial master programmer, promotes interferon induction and suppresses breast cancer recurrence

Type-I and -III interferons play a central role in immune rejection of pathogens and tumors, thus promoting immunogenicity and suppressing tumor recurrence. Double strand RNA is an important ligand that stimulates tumor immunity via interferon responses. Differentiation of embryonic stem cells to pluripotent epithelial cells activates the interferon response during development, raising the question of whether epithelial vs. mesenchymal gene signatures in cancer potentially regulate the interferon pathway as well. Here, using genomics and signaling approaches, we show that Grainyhead-like-2 (GRHL2), a master programmer of epithelial cell identity, promotes type-I and -III interferon responses to double-strand RNA. GRHL2 enhanced the activation of IRF3 and relA/NF-kB and the expression of IRF1; a functional GRHL2 binding site in the IFNL1 promoter was also identified. Moreover, time to recurrence in breast cancer correlated positively with GRHL2 protein expression, indicating that GRHL2 is a tumor recurrence suppressor, consistent with its enhancement of interferon responses. These observations demonstrate that epithelial cell identity supports interferon responses in the context of cancer.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Exploration of the Correlation Between GRHL1 Expression and Tumor Microenvironment in Endometrial Cancer and Immunotherapy

Introduction

GRHL1 belongs to the family of Grainyhead-like (GRHL). Previous studies have shown that dysregulation of growth and survival pathways is associated with the GRHL family of gene cancers. Immunotherapy with checkpoint inhibitors has changed the treatment paradigm for many tumors, including endometrial cancer (EC). However, the effect of GRHL1 on immunotherapy in EC and its relationship with immune cell infiltration are poorly understood.

Methods

Differential expression of GRHL1 between EC and normal EC tissues was analyzed by searching the TCGA database, and the results were verified utilizing immunohistochemistry analyses. Next, the relationship between GRHL1, CD8+ T cells and tumor microenvironment (TME) was also investigated, and the effect of GRHL1 expression on immunotherapy in EC was evaluated.

Results

According to the findings, EC tissues had elevated expression levels of GRHL1 relative to normal tissues. Patients with EC who expressed GRHL1 at high levels experienced worse overall survival (OS) and Progression-free survival (PFS) than those whose expression was lower. In addition, GRHL1 expression was negatively correlated with CD8+ T cells, and patients with high GRHL1 expression were less effective in receiving immunotherapy.

Conclusion

The expression of GRHL1 was high in EC patients, and high expression of GRHL1 inhibits the proliferation of CD8+ T cells in the tumor microenvironment of EC and affect the efficacy of immunotherapy.

Transcriptome and open chromatin analysis reveals the process of myocardial cell development and key pathogenic target proteins in Long QT syndrome type 7

OBJECTIVE

Long QT syndrome type 7 (Andersen-Tawil syndrome, ATS), which is caused by KCNJ2 gene mutation, often leads to ventricular arrhythmia, periodic paralysis and skeletal malformations. The development, differentiation and electrophysiological maturation of cardiomyocytes (CMs) changes promote the pathophysiology of Long QT syndrome type 7(LQT7). We aimed to specifically reproduce the ATS disease phenotype and study the pathogenic mechanism.

METHODS AND RESULTS

We established a cardiac cell model derived from human induced pluripotent stem cells (hiPSCs) to the phenotypes and electrophysiological function, and the establishment of a human myocardial cell model that specifically reproduces the symptoms of ATS provides a reliable platform for exploring the mechanism of this disease or potential drugs. The spontaneous pulsation rate of myocardial cells in the mutation group was significantly lower than that in the repair CRISPR group, the action potential duration was prolonged, and the Kir2.1 current of the inward rectifier potassium ion channel was decreased, which is consistent with the clinical symptoms of ATS patients. Only ZNF528, a chromatin-accessible TF related to pathogenicity, was continuously regulated beginning from the cardiac mesodermal precursor cell stage (day 4), and continued to be expressed at low levels, which was identified by WGCNA method and verified with ATAC-seq data in the mutation group. Subsequently, it indicated that seven pathways were downregulated (all p < 0.05) by used single sample Gene Set Enrichment Analysis to evaluate the overall regulation of potassium-related pathways enriched in the transcriptome and proteome of late mature CMs. Among them, the three pathways (GO: 0008076, GO: 1990573 and GO: 0030007) containing the mutated gene KCNJ2 is involved that are related to the whole process by which a potassium ion enters the cell via the inward rectifier potassium channel to exert its effect were inhibited. The other four pathways are related to regulation of the potassium transmembrane pathway and sodium:potassium exchange ATPase (p < 0.05). ZNF528 small interfering (si)-RNA was applied to hiPSC-derived cardiomyocytes for CRISPR group to explore changes in potassium ion currents and growth and development related target protein levels that affect disease phenotype. Three consistently downregulated proteins (KCNJ2, CTTN and ATP1B1) associated with pathogenicity were verificated through correlation and intersection analysis.

CONCLUSION

This study uncovers TFs and target proteins related to electrophysiology and developmental pathogenicity in ATS myocardial cells, obtaining novel targets for potential therapeutic candidate development that does not rely on gene editing.

Human Pluripotent Stem Cell Derived Organoids Reveal a Role for WNT Signaling in Dorsal-Ventral Patterning of the Hindgut

The cloaca is a transient structure that forms in the terminal hindgut giving rise to the rectum dorsally and the urogenital sinus ventrally. Similarly, human hindgut cultures derived from human pluripotent stem cells generate human colonic organoids (HCOs) which also contain co-developing urothelial tissue. In this study, our goal was to identify pathways involved in cloacal patterning and apply this to human hindgut cultures. RNA-seq data comparing dorsal versus ventral cloaca in e10.5 mice revealed that WNT signaling was elevated in the ventral versus dorsal cloaca. Inhibition of WNT signaling in hindgut cultures biased their differentiation towards a colorectal fate. WNT activation biased differentiation towards a urothelial fate, giving rise to human urothelial organoids (HUOs). HUOs contained cell types present in human urothelial tissue. Based on our results, we propose a mechanism whereby WNT signaling patterns the ventral cloaca, prior to cloacal septation, to give rise to the urogenital sinus.

Grainyhead-like 2 interacts with noggin to regulate tissue fusion in mouse

Defective tissue fusion during mammalian embryogenesis results in congenital anomalies, such as exencephaly, spina bifida and cleft lip and/or palate. The highly conserved transcription factor grainyhead-like 2 (Grhl2) is a crucial regulator of tissue fusion, with mouse models lacking GRHL2 function presenting with a fully penetrant open cranial neural tube, facial and abdominal clefting (abdominoschisis), and an open posterior neuropore. Here, we show that GRHL2 interacts with the soluble morphogen protein and bone morphogenetic protein (BMP) inhibitor noggin (NOG) to impact tissue fusion during development. The maxillary prominence epithelium in embryos lacking Grhl2 shows substantial morphological abnormalities and significant upregulation of NOG expression, together with aberrantly distributed pSMAD5-positive cells within the neural crest cell-derived maxillary prominence mesenchyme, indicative of disrupted BMP signalling. Reducing this elevated NOG expression (by generating Grhl2-/-;Nog+/- embryos) results in delayed embryonic lethality, partial tissue fusion rescue, and restoration of tissue form within the craniofacial epithelia. These data suggest that aberrant epithelial maintenance, partially regulated by noggin-mediated regulation of BMP-SMAD pathways, may underpin tissue fusion defects in Grhl2-/- mice.

In Situ Hybridization to Characterize Neurulation and Midbrain-Hindbrain Boundary Formation in Zebrafish

Whole-mount in situ hybridization is cable to harness the inherent advantages of zebrafish as a model organism for developmental biology, particularly when visualizing the formation of the neural tube, specifically at the level of the midbrain-hindbrain boundary. The size and transparency of developing zebrafish embryos allow for the visualization of neural markers in vivo along the length of the developing zebrafish central nervous system. In practice, this technique is useful for examining defects in neurulation and midbrain-hindbrain boundary formation that may arise following gene manipulation, for example, CRISPR mutagenesis. This method describes the process of embryo collection and preparation, RNA probe transcription, probe hybridization in vivo, as well as the process of probe detection and visualization.

Discovery of a non-canonical GRHL1 binding site using deep convolutional and recurrent neural networks

BACKGROUND

Transcription factors regulate gene expression by binding to transcription factor binding sites (TFBSs). Most models for predicting TFBSs are based on position weight matrices (PWMs), which require a specific motif to be present in the DNA sequence and do not consider interdependencies of nucleotides. Novel approaches such as Transcription Factor Flexible Models or recurrent neural networks consequently provide higher accuracies. However, it is unclear whether such approaches can uncover novel non-canonical, hitherto unexpected TFBSs relevant to human transcriptional regulation.

RESULTS

In this study, we trained a convolutional recurrent neural network with HT-SELEX data for GRHL1 binding and applied it to a set of GRHL1 binding sites obtained from ChIP-Seq experiments from human cells. We identified 46 non-canonical GRHL1 binding sites, which were not found by a conventional PWM approach. Unexpectedly, some of the newly predicted binding sequences lacked the CNNG core motif, so far considered obligatory for GRHL1 binding. Using isothermal titration calorimetry, we experimentally confirmed binding between the GRHL1-DNA binding domain and predicted GRHL1 binding sites, including a non-canonical GRHL1 binding site. Mutagenesis of individual nucleotides revealed a correlation between predicted binding strength and experimentally validated binding affinity across representative sequences. This correlation was neither observed with a PWM-based nor another deep learning approach.

CONCLUSIONS

Our results show that convolutional recurrent neural networks may uncover unanticipated binding sites and facilitate quantitative transcription factor binding predictions.

HSPA12A was identified as a key driver in colorectal cancer GWAS loci 10q26.12 and modulated by an enhancer-promoter interaction

Although genome-wide association studies (GWASs) have identified over 100 colorectal cancer (CRC) risk loci, an understanding of causal genes or risk variants and their biological functions in these loci remain unclear. Recently, genomic loci 10q26.12 with lead SNP rs1665650 was identified as an essential CRC risk loci of Asian populations. However, the functional mechanism of this region has not been fully clarified. Here, we applied an RNA interfering-based on-chip approach to screen for the genes essential for cell proliferation in the CRC risk loci 10q26.12. Notably, HSPA12A had the most significant effect among the identified genes and functioned as a crucial oncogene facilitating cell proliferation. Moreover, we conducted an integrative fine-mapping analysis to identify putative casual variants and further explored their association with CRC risk in a large-scale Chinese population consisting of 4054 cases and 4054 controls and also independently validated in 5208 cases and 20,832 controls from the UK biobank cohort. We identified a risk SNP rs7093835 in the intron of HSPA12A that was significantly associated with an increased risk of CRC (OR 1.23, 95% CI 1.08-1.41, P = 1.92 × 10-3). Mechanistically, the risk variant could facilitate an enhancer-promoter interaction mediated by the transcriptional factor (TF) GRHL1 and ultimately upregulate HSPA12A expression, which provides functional evidence to support our population findings. Collectively, our study reveals the important role of HSPA12A in CRC development and illustrates a novel enhancer-promoter interaction module between HSPA12A and its regulatory elements rs7093835, providing new insights into the etiology of CRC.

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.