The transcription factor Grainy head primes epithelial enhancers for spatiotemporal activation by displacing nucleosomes

Exploring the reciprocity between pioneer factors and development

Development is regulated by coordinated changes in gene expression. Control of these changes in expression is largely governed by the binding of transcription factors to specific regulatory elements. However, the packaging of DNA into chromatin prevents the binding of many transcription factors. Pioneer factors overcome this barrier owing to unique properties that enable them to bind closed chromatin, promote accessibility and, in so doing, mediate binding of additional factors that activate gene expression. Because of these properties, pioneer factors act at the top of gene-regulatory networks and drive developmental transitions. Despite the ability to bind target motifs in closed chromatin, pioneer factors have cell type-specific chromatin occupancy and activity. Thus, developmental context clearly shapes pioneer-factor function. Here, we discuss this reciprocal interplay between pioneer factors and development: how pioneer factors control changes in cell fate and how cellular environment influences pioneer-factor binding and activity.

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.

Pioneer factors: nature or nurture?

Chromatin is densely packed with nucleosomes, which limits the accessibility of many chromatin-associated proteins. Pioneer factors (PFs) are usually viewed as a special group of sequence-specific transcription factors (TFs) that can recognize nucleosome-embedded motifs, invade compact chromatin, and generate open chromatin regions. Through this process, PFs initiate a cascade of events that play key roles in gene regulation and cell differentiation. A current debate in the field is if PFs belong to a unique subset of TFs with intrinsic "pioneering activity", or if all TFs have the potential to function as PFs within certain cellular contexts. There are also different views regarding the key feature(s) that define pioneering activity. In this review, we present evidence from the literature related to these alternative views and discuss how to potentially reconcile them. It is possible that both intrinsic properties, like tight nucleosome binding and structural compatibility, and cellular conditions, like concentration and co-factor availability, are important for PF function.

Divergence of Grainy head affects chromatin accessibility, gene expression, and embryonic viability in Drosophila melanogaster

Pioneer factors are critical for gene regulation and development because they bind chromatin and make DNA more accessible for binding by other transcription factors. The pioneer factor Grainy head (Grh) is present across metazoans and has been shown to retain a role in epithelium development in fruit flies, nematodes, and mice despite extensive divergence in both amino acid sequence and length. Here, we investigate the evolution of Grh function by comparing the effects of the fly ( Drosophila melanogaster ) and worm ( Caenorhabditis elegans ) Grh orthologs on chromatin accessibility, gene expression, embryonic development, and viability in transgenic D. melanogaster . We found that the Caenorhabditis elegans ortholog rescued cuticle development but not full embryonic viability in Drosophila melanogaster grh null mutants. At the molecular level, the C. elegans ortholog only partially rescued chromatin accessibility and gene expression. Divergence in the disordered N-terminus of the Grh protein contributes to these differences in embryonic viability and molecular phenotypes. These data show how pioneer factors can diverge in sequence and function at the molecular level while retaining conserved developmental functions at the organismal level.

SUMMARY STATEMENT

Despite divergence in a disordered region that affects function at both molecular and organismal levels, the Caenorhabditis elegans Grainy head (Grh) protein rescued cuticle morphology in D. melanogaster embryos.

Protein-intrinsic properties and context-dependent effects regulate pioneer factor binding and function

Chromatin is a barrier to the binding of many transcription factors. By contrast, pioneer factors access nucleosomal targets and promote chromatin opening. Despite binding to target motifs in closed chromatin, many pioneer factors display cell-type-specific binding and activity. The mechanisms governing pioneer factor occupancy and the relationship between chromatin occupancy and opening remain unclear. We studied three Drosophila transcription factors with distinct DNA-binding domains and biological functions: Zelda, Grainy head and Twist. We demonstrated that the level of chromatin occupancy is a key determinant of pioneering activity. Multiple factors regulate occupancy, including motif content, local chromatin and protein concentration. Regions outside the DNA-binding domain are required for binding and chromatin opening. Our results show that pioneering activity is not a binary feature intrinsic to a protein but occurs on a spectrum and is regulated by a variety of protein-intrinsic and cell-type-specific features.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

Targeted design of synthetic enhancers for selected tissues in the Drosophila embryo

Enhancers control gene expression and have crucial roles in development and homeostasis1-3. However, the targeted de novo design of enhancers with tissue-specific activities has remained challenging. Here we combine deep learning and transfer learning to design tissue-specific enhancers for five tissues in the Drosophila melanogaster embryo: the central nervous system, epidermis, gut, muscle and brain. We first train convolutional neural networks using genome-wide single-cell assay for transposase-accessible chromatin with sequencing (ATAC-seq) datasets and then fine-tune the convolutional neural networks with smaller-scale data from in vivo enhancer activity assays, yielding models with 13% to 76% positive predictive value according to cross-validation. We designed and experimentally assessed 40 synthetic enhancers (8 per tissue) in vivo, of which 31 (78%) were active and 27 (68%) functioned in the target tissue (100% for central nervous system and muscle). The strategy of combining genome-wide and small-scale functional datasets by transfer learning is generally applicable and should enable the design of tissue-, cell type- and cell state-specific enhancers in any system.

Pioneer factors: roles and their regulation in development

Pioneer factors are a subclass of transcription factors that can bind and initiate opening of silent chromatin regions. Pioneer factors subsequently regulate lineage-specific genes and enhancers and, thus, activate the zygotic genome after fertilization, guide cell fate transitions during development, and promote various forms of human cancers. As such, pioneer factors are useful in directed cell reprogramming. In this review, we define the structural and functional characteristics of pioneer factors, how they bind and initiate opening of closed chromatin regions, and the consequences for chromatin dynamics and gene expression during cell differentiation. We also discuss emerging mechanisms that modulate pioneer factors during development.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types. We found that zonation affects gene expression and chromatin accessibility in hepatocytes, among other cell types. These states are driven by the repressors TCF7L1 and TBX3, alongside other core hepatocyte transcription factors, such as HNF4A, CEBPA, FOXA1 and ONECUT1. To examine the architecture of the enhancers driving these cell states, we trained a hierarchical deep learning model called DeepLiver. Our study provides a multimodal understanding of the regulatory code underlying hepatocyte identity and their zonation state that can be used to engineer enhancers with specific activity levels and zonation patterns.

Genomic and epigenomic integrative subtypes of renal cell carcinoma in a Japanese cohort

Renal cell carcinoma (RCC) comprises several histological types characterised by different genomic and epigenomic aberrations; however, the molecular pathogenesis of each type still requires further exploration. We perform whole-genome sequencing of 128 Japanese RCC cases of different histology to elucidate the significant somatic alterations and mutagenesis processes. We also perform transcriptomic and epigenomic sequencing to identify distinguishing features, including assay for transposase-accessible chromatin sequencing (ATAC-seq) and methyl sequencing. Genomic analysis reveals that the mutational signature differs among the histological types, suggesting that different carcinogenic factors drive each histology. From the ATAC-seq results, master transcription factors are identified for each histology. Furthermore, clear cell RCC is classified into three epi-subtypes, one of which expresses highly immune checkpoint molecules with frequent loss of chromosome 14q. These genomic and epigenomic features may lead to the development of effective therapeutic strategies for RCC.

GCparagon: evaluating and correcting GC biases in cell-free DNA at the fragment level

Analyses of cell-free DNA (cfDNA) are increasingly being employed for various diagnostic and research applications. Many technologies aim to increase resolution, e.g. for detecting early-stage cancer or minimal residual disease. However, these efforts may be confounded by inherent base composition biases of cfDNA, specifically the over - and underrepresentation of guanine (G) and cytosine (C) sequences. Currently, there is no universally applicable tool to correct these effects on sequencing read-level data. Here, we present GCparagon, a two-stage algorithm for computing and correcting GC biases in cfDNA samples. In the initial step, length and GC base count parameters are determined. Here, our algorithm minimizes the inclusion of known problematic genomic regions, such as low-mappability regions, in its calculations. In the second step, GCparagon computes weights counterbalancing the distortion of cfDNA attributes (correction matrix). These fragment weights are added to a binary alignment map (BAM) file as alignment tags for individual reads. The GC correction matrix or the tagged BAM file can be used for downstream analyses. Parallel computing allows for a GC bias estimation below 1 min. We demonstrate that GCparagon vastly improves the analysis of regulatory regions, which frequently show specific GC composition patterns and will contribute to standardized cfDNA applications.

Single-cell transcriptomics illuminates regulatory steps driving anterior-posterior patterning of Drosophila embryonic mesoderm

Single-cell technologies promise to uncover how transcriptional programs orchestrate complex processes during embryogenesis. Here, we apply a combination of single-cell technology and genetic analysis to investigate the dynamic transcriptional changes associated with Drosophila embryo morphogenesis at gastrulation. Our dataset encompassing the blastoderm-to-gastrula transition provides a comprehensive single-cell map of gene expression across cell lineages validated by genetic analysis. Subclustering and trajectory analyses revealed a surprising stepwise progression in patterning to transition zygotic gene expression and specify germ layers as well as uncovered an early role for ecdysone signaling in epithelial-to-mesenchymal transition in the mesoderm. We also show multipotent progenitors arise prior to gastrulation by analyzing the transcription trajectory of caudal mesoderm cells, including a derivative that ultimately incorporates into visceral muscles of the midgut and hindgut. This study provides a rich resource of gastrulation and elucidates spatially regulated temporal transitions of transcription states during the process.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

Widespread regulatory specificities between transcriptional co-repressors and enhancers in Drosophila

Gene expression is controlled by the precise activation and repression of transcription. Repression is mediated by specialized transcription factors (TFs) that recruit co-repressors (CoRs) to silence transcription, even in the presence of activating cues. However, whether CoRs can dominantly silence all enhancers or display distinct specificities is unclear. In this work, we report that most enhancers in Drosophila can be repressed by only a subset of CoRs, and enhancers classified by CoR sensitivity show distinct chromatin features, function, TF motifs, and binding. Distinct TF motifs render enhancers more resistant or sensitive to specific CoRs, as we demonstrate by motif mutagenesis and addition. These CoR-enhancer compatibilities constitute an additional layer of regulatory specificity that allows differential regulation at close genomic distances and is indicative of distinct mechanisms of transcriptional repression.

Resistance to mesenchymal reprogramming sustains clonal propagation in metastatic breast cancer

The acquisition of mesenchymal traits is considered a hallmark of breast cancer progression. However, the functional relevance of epithelial-to-mesenchymal transition (EMT) remains controversial and context dependent. Here, we isolate epithelial and mesenchymal populations from human breast cancer metastatic biopsies and assess their functional potential in vivo. Strikingly, progressively decreasing epithelial cell adhesion molecule (EPCAM) levels correlate with declining disease propagation. Mechanistically, we find that persistent EPCAM expression marks epithelial clones that resist EMT induction and propagate competitively. In contrast, loss of EPCAM defines clones arrested in a mesenchymal state, with concomitant suppression of tumorigenicity and metastatic potential. This dichotomy results from distinct clonal trajectories impacting global epigenetic programs that are determined by the interplay between human ZEB1 and its target GRHL2. Collectively, our results indicate that susceptibility to irreversible EMT restrains clonal propagation, whereas resistance to mesenchymal reprogramming sustains disease spread in multiple models of human metastatic breast cancer, including patient-derived cells in vivo.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

Protein-intrinsic properties and context-dependent effects regulate pioneer-factor binding and function

Chromatin is a barrier to the binding of many transcription factors. By contrast, pioneer factors access nucleosomal targets and promote chromatin opening. Despite binding to target motifs in closed chromatin, many pioneer factors display cell-type specific binding and activity. The mechanisms governing pioneer-factor occupancy and the relationship between chromatin occupancy and opening remain unclear. We studied three Drosophila transcription factors with distinct DNA-binding domains and biological functions: Zelda, Grainy head, and Twist. We demonstrated that the level of chromatin occupancy is a key determinant of pioneering activity. Multiple factors regulate occupancy, including motif content, local chromatin, and protein concentration. Regions outside the DNA-binding domain are required for binding and chromatin opening. Our results show that pioneering activity is not a binary feature intrinsic to a protein but occurs on a spectrum and is regulated by a variety of protein-intrinsic and cell-type-specific features.

GRHL2 and AP2a coordinate early surface ectoderm lineage commitment during development

Ectodermal dysplasias including skin abnormalities and cleft lip/palate result from improper surface ectoderm (SE) patterning. However, the connection between SE gene regulatory networks and disease remains poorly understood. Here, we dissect human SE differentiation with multiomics and establish GRHL2 as a key mediator of early SE commitment, which acts by skewing cell fate away from the neural lineage. GRHL2 and master SE regulator AP2a balance early cell fate output, with GRHL2 facilitating AP2a binding to SE loci. In turn, AP2a restricts GRHL2 DNA binding away from de novo chromatin contacts. Integration of these regulatory sites with ectodermal dysplasia-associated genomic variants annotated within the Biomedical Data Commons identifies 55 loci previously implicated in craniofacial disorders. These include ABCA4/ARHGAP29 and NOG regulatory regions where disease-linked variants directly affect GRHL2/AP2a binding and gene transcription. These studies elucidate the logic underlying SE commitment and deepen our understanding of human oligogenic disease pathogenesis.

The transcription factor optomotor-blind restricts apterous expression through TrxG and PcG genes

The establishment of body pattern is a fundamental process in developmental biology. In Drosophila, the wing disc is subdivided into dorsal (D) and ventral (V) compartments by the D/V boundary. The dorsal fate is adopted by expressing the selector gene apterous (ap). ap expression is regulated by three combinational cis-regulatory modules which are activated by EGFR pathway, Ap-Vg auto-regulatory and epigenetic mechanisms. Here, we found that the Tbx family transcription factor Optomotor-blind (Omb) restricted ap expression in the ventral compartment. Loss of omb induced autonomous initiation of ap expression in the middle third instar larvae in the ventral compartment. Oppositely, over-activation of omb inhibited ap in the medial pouch. All three enhancers apE, apDV and apP were upregulated in omb null mutants, indicating a combinational regulation of ap modulators. However, Omb affected ap expression neither by directly regulating EGFR signaling, nor via Vg regulation. Therefore, a genetic screen of epigenetic regulators, including the Trithorax group (TrxG) and Polycomb group (PcG) genes was performed. We found that knocking down the TrxG gene kohtalo (kto), domino (dom) or expressing the PcG gene grainy head (grh), the ectopic ap in omb mutants was repressed. The inhibition of apDV by kto knockdown and grh activation could contribute to ap repression. Moreover, Omb and the EGFR pathway are genetically parallel in ap regulation in the ventral compartment. Collectively, Omb is a repressive signal for ap expression in the ventral compartment, which requires TrxG and PcG genes.

Increased chromatin accessibility promotes the evolution of a transcriptional silencer in Drosophila

The loss of discrete morphological traits, the most common evolutionary transition, is typically driven by changes in developmental gene expression. Mutations accumulating in regulatory elements of these genes can disrupt DNA binding sites for transcription factors patterning their spatial expression, or delete entire enhancers. Regulatory elements, however, may be silenced through changes in chromatin accessibility or the emergence of repressive elements. Here, we show that increased chromatin accessibility at the gene yellow, combined with the gain of a repressor site, underlies the loss of a wing spot pigmentation pattern in a Drosophila species. The gain of accessibility of this repressive element is regulated by E93, a transcription factor governing the progress of metamorphosis. This convoluted evolutionary scenario contrasts with the parsimonious mutational paths generally envisioned and often documented for morphological losses. It illustrates how evolutionary changes in chromatin accessibility may directly contribute to morphological diversification.

Enhancer architecture and chromatin accessibility constrain phenotypic space during Drosophila development

Developmental enhancers bind transcription factors and dictate patterns of gene expression during development. Their molecular evolution can underlie phenotypical evolution, but the contributions of the evolutionary pathways involved remain little understood. Here, using mutation libraries in Drosophila melanogaster embryos, we observed that most point mutations in developmental enhancers led to changes in gene expression levels but rarely resulted in novel expression outside of the native pattern. In contrast, random sequences, often acting as developmental enhancers, drove expression across a range of cell types; random sequences including motifs for transcription factors with pioneer activity acted as enhancers even more frequently. Our findings suggest that the phenotypic landscapes of developmental enhancers are constrained by enhancer architecture and chromatin accessibility. We propose that the evolution of existing enhancers is limited in its capacity to generate novel phenotypes, whereas the activity of de novo elements is a primary source of phenotypic novelty.

Differential Grainy head binding correlates with variation in chromatin structure and gene expression in Drosophila melanogaster

Phenotypic evolution is often caused by variation in gene expression resulting from altered gene regulatory mechanisms. Genetic variation affecting chromatin remodeling has been identified as a potential source of variable gene expression; however, the roles of specific chromatin remodeling factors remain unclear. Here, we address this knowledge gap by examining the relationship between variation in gene expression, variation in chromatin structure, and variation in binding of the pioneer factor Grainy head between imaginal wing discs of two divergent strains of Drosophila melanogaster and their F₁ hybrid. We find that (1) variation in Grainy head binding is mostly due to sequence changes that act in cis but are located outside of the canonical Grainy head binding motif, (2) variation in Grainy head binding correlates with changes in chromatin accessibility, and (3) this variation in chromatin accessibility, coupled with variation in Grainy head binding, correlates with variation in gene expression in some cases but not others. Interactions among these three molecular layers is complex, but these results suggest that genetic variation affecting the binding of pioneer factors contributes to variation in chromatin remodeling and the evolution of gene expression.

Hox function and specificity – A tissue centric view

Since their discovery, the Hox genes, with their incredible power to reprogram the identity of complete body regions, a phenomenon called homeosis, have captured the fascination of many biologists. Recent research has provided new insights into the function of Hox proteins in different germ layers and the mechanisms they employ to control tissue morphogenesis. We focus in this review on the ectoderm and mesoderm to highlight new findings and discuss them with regards to established concepts of Hox target gene regulation. Furthermore, we highlight the molecular mechanisms involved the transcriptional repression of specific groups of Hox target genes, and summarize the role of Hox mediated gene silencing in tissue development. Finally, we reflect on recent findings identifying a large number of tissue-specific Hox interactor partners, which open up new avenues and directions towards a better understanding of Hox function and specificity in different tissues.

GRHL2 Enhances Phosphorylated Estrogen Receptor (ER) Chromatin Binding and Regulates ER-Mediated Transcriptional Activation and Repression

Phosphorylation of estrogen receptor α (ER) at serine 118 (pS118-ER) is induced by estrogen and is the most abundant posttranslational mark associated with a transcriptionally active receptor. Cistromic analysis of pS118-ER from our group revealed enrichment of the GRHL2 motif near pS118-ER binding sites. In this study, we used cistromic and transcriptomic analyses to interrogate the relationship between GRHL2 and pS118-ER. We found that GRHL2 is bound to chromatin at pS118-ER/GRHL2 co-occupancy sites prior to ligand treatment, and GRHL2 binding is required for maximal pS118-ER recruitment. pS118-ER/GRHL2 co-occupancy sites were enriched at active enhancers marked by H3K27ac and H3K4me1, along with FOXA1 and p300, compared to sites where each factor binds independently. Transcriptomic analysis yielded four subsets of ER/GRHL2-coregulated genes revealing that GRHL2 can both enhance and antagonize E2-mediated ER transcriptional activity. Gene ontology analysis indicated that coregulated genes are involved in cell migration. Accordingly, knockdown of GRHL2, combined with estrogen treatment, resulted in increased cell migration but no change in proliferation. These results support a model in which GRHL2 binds to selected enhancers and facilitates pS118-ER recruitment to chromatin, which then results in differential activation and repression of genes that control estrogen-regulated ER-positive breast cancer cell migration.

The regulatory role of pioneer factors during cardiovascular lineage specification – A mini review

Cardiovascular disease (CVD) remains the number one cause of death worldwide. Ischemic heart disease contributes to heart failure and has considerable morbidity and mortality. Therefore, alternative therapeutic strategies are urgently needed. One class of epigenetic regulators known as pioneer factors has emerged as an important tool for the development of regenerative therapies for the treatment of CVD. Pioneer factors bind closed chromatin and remodel it to drive lineage specification. Here, we review pioneer factors within the cardiovascular lineage, particularly during development and reprogramming and highlight the implications this field of research has for the future development of cardiac specific regenerative therapies.

Divergent expression of paralogous genes by modification of shared enhancer activity through a promoter-proximal silencer

The duplication of genes and their associated cis-regulatory elements, or enhancers, is a key contributor to genome evolution and biological complexity. Moreover, many paralogs, particularly tandem duplicates, are fixed for long periods of time under the control of shared enhancers. However, in most cases, the mechanism by which gene expression and function diverge following duplication is not known. Here, we dissect the regulation and function of the paralogous nubbin/pdm2 genes during wing development in Drosophila melanogaster. We show that these paralogs play a redundant role in the wing and that their expression relies on a single shared wing enhancer. However, the two genes differ in their ability to respond to this enhancer, with nub responding in all wing progenitor cells and pdm2 only in a small subset. This divergence is a result of a pdm2-specific silencer element at the pdm2 promoter that receives repressive input from the transcription factor Rotund. Repression through this silencer also depends on nub, allowing pdm2 to fully respond to the wing enhancer when nub expression is perturbed and functional compensation to occur. Thus, expression divergence downstream of a shared enhancer arises as a consequence of silencing the promoter of one paralog.

Cis-regulatory chromatin loops analysis identifies GRHL3 as a master regulator of surface epithelium commitment

Understanding the regulatory network of cell fate acquisition remains a major challenge. Using the induction of surface epithelium (SE) from human embryonic stem cells as a paradigm, we show that the dynamic changes in morphology-related genes (MRGs) closely correspond to SE fate transitions. The marked remodeling of cytoskeleton indicates the initiation of SE differentiation. By integrating promoter interactions, epigenomic features, and transcriptome, we delineate an SE-specific cis-regulatory network and identify grainyhead-like 3 (GRHL3) as an initiation factor sufficient to drive SE commitment. Mechanically, GRHL3 primes the SE chromatin accessibility landscape and activates SE-initiating gene expression. In addition, the evaluation of GRHL3-mediated promoter interactions unveils a positive feedback loop of GRHL3 and bone morphogenetic protein 4 on SE fate decisions. Our work proposes a concept that MRGs could be used to identify cell fate transitions and provides insights into regulatory principles of SE lineage development and stem cell-based regenerative medicine.

Hox genes collaborate with helix-loop-helix factor Grainyhead to promote neuroblast apoptosis along the anterior-posterior axis of the Drosophila larval central nervous system

Hox genes code for a family of a homeodomain (HD) containing transcription factors that use TALE-HD containing factors Pbx/Exd and Meis/Hth to specify the development of the anterior-posterior (AP) axis of an organism. However, the absence of TALE-HD containing factors from specific tissues emphasizes the need to identify and validate new Hox cofactors. In Drosophila central nervous system (CNS), Hox execute segment-specific apoptosis of neural stem cells (neuroblasts-NBs) and neurons. In abdominal segments of larval CNS, Hox gene Abdominal-A (AbdA) mediates NB apoptosis with the help of Exd and bHLH factor Grainyhead (Grh) using a 717 bp apoptotic enhancer. In this study, we show that this enhancer is critical for abdominal NB apoptosis and relies on two separable set of DNA binding motifs responsible for its initiation and maintenance. Our results also show that AbdA and Grh interact through their highly conserved DNA binding domains, and the DNA binding specificity of AbdA-HD is important for it to interact with Grh and essential for it to execute NB apoptosis in CNS. We also establish that Grh is required for Hox-dependent NB apoptosis in Labial and Sex Combs Reduced (Scr) expressing regions of the CNS, and it can physically interact with all the Hox proteins in vitro. Our biochemical and functional data collectively support the idea that Grh can function as a Hox cofactor and help them carry out their in vivo roles during development.

Epigenetic Mechanisms Influencing Therapeutic Response in Breast Cancer

The majority of breast cancers are estrogen receptor (ER)+ and agents targeting the ER signaling pathway have markedly increased survival for women with breast cancer for decades. However, therapeutic resistance eventually emerges, especially in the metastatic setting. In the past decade disrupted epigenetic regulatory processes have emerged as major contributors to carcinogenesis in many cancer types. Aberrations in chromatin modifiers and transcription factors have also been recognized as mediators of breast cancer development and therapeutic outcome, and new epigenetic-based therapies in combination with targeted therapies have been proposed. Here we will discuss recent progress in our understanding of the chromatin-based mechanisms of breast tumorigenesis, how these mechanisms affect therapeutic response to standard of care treatment, and discuss new strategies towards therapeutic intervention to overcome resistance.

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Evolutionarily conserved transcription factors as regulators of longevity and targets for geroprotection

Aging is the single largest risk factor for many debilitating conditions, including heart diseases, stroke, cancer, diabetes, and neurodegenerative disorders. While far from understood in its full complexity, it is scientifically well-established that aging is influenced by genetic and environmental factors, and can be modulated by various interventions. One of aging's early hallmarks are aberrations in transcriptional networks, controlling for example metabolic homeostasis or the response to stress. Evidence in different model organisms abounds that a number of evolutionarily conserved transcription factors, which control such networks, can affect lifespan and healthspan across species. These transcription factors thus potentially represent conserved regulators of longevity and are emerging as important targets in the challenging quest to develop treatments to mitigate age-related diseases, and possibly even to slow aging itself. This review provides an overview of evolutionarily conserved transcription factors that impact longevity or age-related diseases in at least one multicellular model organism (nematodes, flies, or mice), and/or are tentatively linked to human aging. Discussed is the general evidence for transcriptional regulation of aging and disease, followed by a more detailed look at selected transcription factor families, the common metabolic pathways involved, and the targeting of transcription factors as a strategy for geroprotective interventions.

Pioneer factors as master regulators of the epigenome and cell fate

Pioneer factors are transcription factors with the unique ability to initiate opening of closed chromatin. The stability of cell identity relies on robust mechanisms that maintain the epigenome and chromatin accessibility to transcription factors. Pioneer factors counter these mechanisms to implement new cell fates through binding of DNA target sites in closed chromatin and introduction of active-chromatin histone modifications, primarily at enhancers. As master regulators of enhancer activation, pioneers are thus crucial for the implementation of correct cell fate decisions in development, and as such, they hold tremendous potential for therapy through cellular reprogramming. The power of pioneer factors to reshape the epigenome also presents an Achilles heel, as their misexpression has major pathological consequences, such as in cancer. In this Review, we discuss the emerging mechanisms of pioneer factor functions and their roles in cell fate specification, cellular reprogramming and cancer.

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

Grainyhead-like (GRHL) factors are essential, highly conserved transcription factors (TFs) that regulate processes common to both natural cellular behaviours during embryogenesis, and de-regulation of growth and survival pathways in cancer. Serving to drive the transcription, and therefore activation of multiple co-ordinating pathways, the three GRHL family members (GRHL1-3) are a critical conduit for modulating the molecular landscape that guides cellular decision-making processes during proliferation, epithelial-mesenchymal transition (EMT) and migration. Animal models and in vitro approaches harbouring GRHL loss or gain-of-function are key research tools to understanding gene function, which gives confidence that resultant phenotypes and cellular behaviours may be translatable to humans. Critically, identifying and characterising the target genes to which these factors bind is also essential, as they allow us to discover and understand novel genetic pathways that could ultimately be used as targets for disease diagnosis, drug discovery and therapeutic strategies. GRHL1-3 and their transcriptional targets have been shown to drive comparable cellular processes in Drosophila, C. elegans, zebrafish and mice, and have recently also been implicated in the aetiology and/or progression of a number of human congenital disorders and cancers of epithelial origin. In this review, we will summarise the state of knowledge pertaining to the role of the GRHL family target genes in both development and cancer, primarily through understanding the genetic pathways transcriptionally regulated by these factors across disparate disease contexts.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences.

Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Grainyhead 1 acts as a drug-inducible conserved transcriptional regulator linked to insulin signaling and lifespan

Aging is impacted by interventions across species, often converging on metabolic pathways. Transcription factors regulate longevity yet approaches for their pharmacological modulation to exert geroprotection remain sparse. We show that increased expression of the transcription factor Grainyhead 1 (GRH-1) promotes lifespan and pathogen resistance in Caenorhabditis elegans. A compound screen identifies FDA-approved drugs able to activate human GRHL1 and promote nematodal GRH-1-dependent longevity. GRHL1 activity is regulated by post-translational lysine methylation and the phosphoinositide (PI) 3-kinase C2A. Consistently, nematodal longevity following impairment of the PI 3-kinase or insulin/IGF-1 receptor requires grh-1. In BXD mice, Grhl1 expression is positively correlated with lifespan and insulin sensitivity. In humans, GRHL1 expression positively correlates with insulin receptor signaling and also with lifespan. Fasting blood glucose levels, including in individuals with type 2 diabetes, are negatively correlated with GRHL1 expression. Thereby, GRH-1/GRHL1 is identified as a pharmacologically malleable transcription factor impacting insulin signaling and lifespan.

Life- and healthspan of organisms can be modulated by dietary, genetic, or pharmacological interventions, which often affect metabolic pathways. Here the authors report that Grainyhead 1 is an evolutionarily conserved, drug-inducible transcription factor that promotes longevity in C. elegans, and thus a potential target for the development of geroprotective drugs.

Decoding gene regulation in the fly brain

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis^1-6, three-dimensional morphological classification⁷ and electron microscopy mapping of the connectome^8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation.

Quantitative-enhancer-FACS-seq (QeFS) reveals epistatic interactions among motifs within transcriptional enhancers in developing Drosophila tissue

Understanding the contributions of transcription factor DNA binding sites to transcriptional enhancers is a significant challenge. We developed Quantitative enhancer-FACS-Seq for highly parallel quantification of enhancer activities from a genomically integrated reporter in Drosophila melanogaster embryos. We investigate the contributions of the DNA binding motifs of four poorly characterized TFs to the activities of twelve embryonic mesodermal enhancers. We measure quantitative changes in enhancer activity and discover a range of epistatic interactions among the motifs, both synergistic and alleviating. We find that understanding the regulatory consequences of TF binding motifs requires that they be investigated in combination across enhancer contexts.

Cell-type-specific Hox regulatory strategies orchestrate tissue identity

Hox proteins are homeodomain transcription factors that diversify serially homologous segments along the animal body axis, as revealed by the classic bithorax phenotype of Drosophila melanogaster, in which mutations in Ultrabithorax (Ubx) transform the third thoracic segment into the likeness of the second thoracic segment. To specify segment identity, we show that Ubx both increases and decreases chromatin accessibility, coinciding with its dual role as both an activator and repressor of transcription. However, the choice of transcriptional activity executed by Ubx is spatially regulated and depends on the availability of cofactors, with Ubx acting as a repressor in some populations and as an activator in others. Ubx-mediated changes to chromatin accessibility positively and negatively affect the binding of Scalloped (Sd), a transcription factor that is required for appendage development in both segments. These findings illustrate how a single Hox protein can modify complex gene regulatory networks to transform the identity of an entire tissue.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

GRHL3 activates FSCN1 to relax cell-cell adhesions between migrating keratinocytes during wound reepithelialization

The migrating keratinocyte wound front is required for skin wound closure. Despite significant advances in wound healing research, we do not fully understand the molecular mechanisms that orchestrate collective keratinocyte migration. Here, we show that, in the wound front, the epidermal transcription factor Grainyhead like-3 (GRHL3) mediates decreased expression of the adherens junction protein E-cadherin; this results in relaxed adhesions between suprabasal keratinocytes, thus promoting collective cell migration and wound closure. Wound fronts from mice lacking GRHL3 in epithelial cells (Grhl3-cKO) have lower expression of Fascin-1 (FSCN1), a known negative regulator of E-cadherin. Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) on wounded keratinocytes shows decreased wound-induced chromatin accessibility near the Fscn1 gene in Grhl3-cKO mice, a region enriched for GRHL3 motifs. These data reveal a wound-induced GRHL3/FSCN1/E-cadherin pathway that regulates keratinocyte-keratinocyte adhesion during wound-front migration; this pathway is activated in acute human wounds and is altered in diabetic wounds in mice, suggesting translational relevance.

Transcription Control of Liver Development

During liver organogenesis, cellular transcriptional profiles are constantly reshaped by the action of hepatic transcriptional regulators, including FoxA1-3, GATA4/6, HNF1α/β, HNF4α, HNF6, OC-2, C/EBPα/β, Hex, and Prox1. These factors are crucial for the activation of hepatic genes that, in the context of compact chromatin, cannot access their targets. The initial opening of highly condensed chromatin is executed by a special class of transcription factors known as pioneer factors. They bind and destabilize highly condensed chromatin and facilitate access to other "non-pioneer" factors. The association of target genes with pioneer and non-pioneer transcription factors takes place long before gene activation. In this way, the underlying gene regulatory regions are marked for future activation. The process is called "bookmarking", which confers transcriptional competence on target genes. Developmental bookmarking is accompanied by a dynamic maturation process, which prepares the genomic loci for stable and efficient transcription. Stable hepatic expression profiles are maintained during development and adulthood by the constant availability of the main regulators. This is achieved by a self-sustaining regulatory network that is established by complex cross-regulatory interactions between the major regulators. This network gradually grows during liver development and provides an epigenetic memory mechanism for safeguarding the optimal expression of the regulators.

EGFR-ERK induced activation of GRHL1 promotes cell cycle progression by up-regulating cell cycle related genes in lung cancer

Grainyhead-like 1 (GRHL1) is a transcription factor involved in embryonic development. However, little is known about the biological functions of GRHL1 in cancer. In this study, we found that GRHL1 was upregulated in non-small cell lung cancer (NSCLC) and correlated with poor survival of patients. GRHL1 overexpression promoted the proliferation of NSCLC cells and knocking down GRHL1 inhibited the proliferation. RNA sequencing showed that a series of cell cycle-related genes were altered when knocking down GRHL1. We further demonstrated that GRHL1 could regulate the expression of cell cycle-related genes by binding to the promoter regions and increasing the transcription of the target genes. Besides, we also found that EGF stimulation could activate GRHL1 and promoted its nuclear translocation. We identified the key phosphorylation site at Ser76 on GRHL1 that is regulated by the EGFR-ERK axis. Taken together, these findings elucidate a new function of GRHL1 on regulating the cell cycle progression and point out the potential role of GRHL1 as a drug target in NSCLC.

Pioneering the developmental frontier

Coordinated changes in gene expression allow a single fertilized oocyte to develop into a complex multi-cellular organism. These changes in expression are controlled by transcription factors that gain access to discrete cis-regulatory elements in the genome, allowing them to activate gene expression. Although nucleosomes present barriers to transcription factor occupancy, pioneer transcription factors have unique properties that allow them to bind DNA in the context of nucleosomes, define cis-regulatory elements, and facilitate the subsequent binding of additional factors that determine gene expression. In this capacity, pioneer factors act at the top of gene-regulatory networks to control developmental transitions. Developmental context also influences pioneer factor binding and activity. Here we discuss the interplay between pioneer factors and development, their role in driving developmental transitions, and the influence of the cellular environment on pioneer factor binding and activity.

Grainyhead-like transcription factors: guardians of the skin barrier

There has been selective pressure to maintain a skin barrier since terrestrial animals evolved 360 million years ago. These animals acquired an unique integumentary system with a keratinized, stratified, squamous epithelium surface barrier. The barrier protects against dehydration and entry of microbes and toxins. The skin barrier centres on the stratum corneum layer of the epidermis and consists of cornified envelopes cemented by the intercorneocyte lipid matrix. Multiple components of the barrier undergo cross-linking by transglutaminase (TGM) enzymes, while keratins provide additional mechanical strength. Cellular tight junctions also are crucial for barrier integrity. The grainyhead-like (GRHL) transcription factors regulate the formation and maintenance of the integument in diverse species. GRHL3 is essential for formation of the skin barrier during embryonic development, whereas GRHL1 maintains the skin barrier postnatally. This is achieved by transactivation of Tgm1 and Tgm5, respectively. In addition to its barrier function, GRHL3 plays key roles in wound repair and as an epidermal tumour suppressor. In its former role, GRHL3 activates the planar cell polarity signalling pathway to mediate wound healing by providing directional migration cues. In squamous epithelium, GRHL3 regulates the balance between proliferation and differentiation, and its loss induces squamous cell carcinoma (SCC). In the skin, this is mediated through increased expression of MIR21, which reduces the expression levels of GRHL3 and its direct target, PTEN, leading to activation of the PI3K-AKT signalling pathway. These data position the GRHL family as master regulators of epidermal homeostasis across a vast gulf of evolutionary history.

Interpretation of allele-specific chromatin accessibility using cell state-aware deep learning

Genomic sequence variation within enhancers and promoters can have a significant impact on the cellular state and phenotype. However, sifting through the millions of candidate variants in a personal genome or a cancer genome, to identify those that impact cis-regulatory function, remains a major challenge. Interpretation of noncoding genome variation benefits from explainable artificial intelligence to predict and interpret the impact of a mutation on gene regulation. Here we generate phased whole genomes with matched chromatin accessibility, histone modifications, and gene expression for 10 melanoma cell lines. We find that training a specialized deep learning model, called DeepMEL2, on melanoma chromatin accessibility data can capture the various regulatory programs of the melanocytic and mesenchymal-like melanoma cell states. This model outperforms motif-based variant scoring, as well as more generic deep learning models. We detect hundreds to thousands of allele-specific chromatin accessibility variants (ASCAVs) in each melanoma genome, of which 15%-20% can be explained by gains or losses of transcription factor binding sites. A considerable fraction of ASCAVs are caused by changes in AP-1 binding, as confirmed by matched ChIP-seq data to identify allele-specific binding of JUN and FOSL1. Finally, by augmenting the DeepMEL2 model with ChIP-seq data for GABPA, the TERT promoter mutation, as well as additional ETS motif gains, can be identified with high confidence. In conclusion, we present a new integrative genomics approach and a deep learning model to identify and interpret functional enhancer mutations with allelic imbalance of chromatin accessibility and gene expression.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Core transcription regulatory circuitry orchestrates corneal epithelial homeostasis

Adult stem cell identity, plasticity, and homeostasis are precisely orchestrated by lineage-restricted epigenetic and transcriptional regulatory networks. Here, by integrating super-enhancer and chromatin accessibility landscapes, we delineate core transcription regulatory circuitries (CRCs) of limbal stem/progenitor cells (LSCs) and find that RUNX1 and SMAD3 are required for maintenance of corneal epithelial identity and homeostasis. RUNX1 or SMAD3 depletion inhibits PAX6 and induces LSCs to differentiate into epidermal-like epithelial cells. RUNX1, PAX6, and SMAD3 (RPS) interact with each other and synergistically establish a CRC to govern the lineage-specific cis-regulatory atlas. Moreover, RUNX1 shapes LSC chromatin architecture via modulating H3K27ac deposition. Disturbance of RPS cooperation results in cell identity switching and dysfunction of the corneal epithelium, which is strongly linked to various human corneal diseases. Our work highlights CRC TF cooperativity for establishment of stem cell identity and lineage commitment, and provides comprehensive regulatory principles for human stratified epithelial homeostasis and pathogenesis.

Corneal epithelium shares similar molecular signatures to other stratified epithelia. Here, the authors map super-enhancers and accessible chromatin in corneal epithelium, identifying a transcription regulatory circuit, including RUNX1, PAX6, and SMAD3, required for corneal epithelial identity and homeostasis.

Loss of FOXC1 contributes to the corneal epithelial fate switch and pathogenesis

Forkhead box C1 (FOXC1) is required for neural crest and ocular development, and mutations in FOXC1 lead to inherited Axenfeld-Rieger syndrome. Here, we find that FOXC1 and paired box 6 (PAX6) are co-expressed in the human limbus and central corneal epithelium. Deficiency of FOXC1 and alternation in epithelial features occur in patients with corneal ulcers. FOXC1 governs the fate of the corneal epithelium by directly binding to lineage-specific open promoters or enhancers marked by H3K4me2. FOXC1 depletion not only activates the keratinization pathway and reprograms corneal epithelial cells into skin-like epithelial cells, but also disrupts the collagen metabolic process and interferon signaling pathways. Loss of interferon regulatory factor 1 and PAX6 induced by FOXC1 dysfunction is linked to the corneal ulcer. Collectively, our results reveal a FOXC1-mediated regulatory network responsible for corneal epithelial homeostasis and provide a potential therapeutic target for corneal ulcer.

Mesenchymal-to-Epithelial Transitions in Development and Cancer

The evolutionary emergence of the mesenchymal phenotype greatly increased the complexity of tissue architecture and composition in early Metazoan species. At the molecular level, an epithelial-to-mesenchymal transition (EMT) was permitted by the innovation of specific transcription factors whose expression is sufficient to repress the epithelial transcriptional program. The reverse process, mesenchymal-to-epithelial transition (MET), involves direct inhibition of EMT transcription factors by numerous mechanisms including tissue-specific MET-inducing transcription factors (MET-TFs), micro-RNAs, and changes to cell and tissue architecture, thus providing an elegant solution to the need for tight temporal and spatial control over EMT and MET events during development and adult tissue homeostasis.