TFAP2 paralogs facilitate chromatin access for MITF at pigmentation and cell proliferation genes

GET: a foundation model of transcription across human cell types

Transcriptional regulation, involving the complex interplay between regulatory sequences and proteins, directs all biological processes. Computational models of transcription lack generalizability to accurately extrapolate in unseen cell types and conditions. Here, we introduce GET, an interpretable foundation model designed to uncover regulatory grammars across 213 human fetal and adult cell types. Relying exclusively on chromatin accessibility data and sequence information, GET achieves experimental-level accuracy in predicting gene expression even in previously unseen cell types. GET showcases remarkable adaptability across new sequencing platforms and assays, enabling regulatory inference across a broad range of cell types and conditions, and uncovering universal and cell type specific transcription factor interaction networks. We evaluated its performance on prediction of regulatory activity, inference of regulatory elements and regulators, and identification of physical interactions between transcription factors. Specifically, we show GET outperforms current models in predicting lentivirus-based massive parallel reporter assay readout with reduced input data. In fetal erythroblasts, we identify distal (>1Mbp) regulatory regions that were missed by previous models. In B cells, we identified a lymphocyte-specific transcription factor-transcription factor interaction that explains the functional significance of a leukemia-risk predisposing germline mutation. In sum, we provide a generalizable and accurate model for transcription together with catalogs of gene regulation and transcription factor interactions, all with cell type specificity.

Transcription factor activating enhancer-binding protein 2ε (AP2ε) modulates phenotypic plasticity and progression of malignant melanoma

Malignant melanoma, the most aggressive form of skin cancer, is often incurable once metastatic dissemination of cancer cells to distant organs has occurred. We investigated the role of Transcription Factor Activating Enhancer-Binding Protein 2ε (AP2ε) in the progression of metastatic melanoma. Here, we observed that AP2ε is a potent activator of metastasis and newly revealed AP2ε to be an important player in melanoma plasticity. High levels of AP2ε lead to worsened prognosis of melanoma patients. Using a transgenic melanoma mouse model with a specific loss of AP2ε expression, we confirmed the impact of AP2ε to modulate the dynamic switch from a migratory to a proliferative phenotype. AP2ε deficient melanoma cells show a severely reduced migratory potential in vitro and reduced metastatic behavior in vivo. Consistently, we revealed increased activity of AP2ε in quiescent and migratory cells compared to heterogeneously proliferating cells in bioprinted 3D models. In conclusion, these findings disclose a yet-unknown role of AP2ε in maintaining plasticity and migration in malignant melanoma cells.

Regulation of Enhancers by SUMOylation Through TFAP2C Binding and Recruitment of HDAC Complex to the Chromatin

Enhancers are fundamental to gene regulation. Post-translational modifications by the small ubiquitin-like modifiers (SUMO) modify chromatin regulation enzymes, including histone acetylases and deacetylases. However, it remains unclear whether SUMOylation regulates enhancer marks, acetylation at the 27th lysine residue of the histone H3 protein (H3K27Ac). To investigate whether SUMOylation regulates H3K27Ac, we performed genome-wide ChIP-seq analyses and discovered that knockdown (KD) of the SUMO activating enzyme catalytic subunit UBA2 reduced H3K27Ac at most enhancers. Bioinformatic analysis revealed that TFAP2C-binding sites are enriched in enhancers whose H3K27Ac was reduced by UBA2 KD. ChIP-seq analysis in combination with molecular biological methods showed that TFAP2C binding to enhancers increased upon UBA2 KD or inhibition of SUMOylation by a small molecule SUMOylation inhibitor. However, this is not due to the SUMOylation of TFAP2C itself. Proteomics analysis of TFAP2C interactome on the chromatin identified histone deacetylation (HDAC) and RNA splicing machineries that contain many SUMOylation targets. TFAP2C KD reduced HDAC1 binding to chromatin and increased H3K27Ac marks at enhancer regions, suggesting that TFAP2C is important in recruiting HDAC machinery. Taken together, our findings provide insights into the regulation of enhancer marks by SUMOylation and TFAP2C and suggest that SUMOylation of proteins in the HDAC machinery regulates their recruitments to enhancers.

DNA damage remodels the MITF interactome to increase melanoma genomic instability

Since genome instability can drive cancer initiation and progression, cells have evolved highly effective and ubiquitous DNA damage response (DDR) programs. However, some cells (for example, in skin) are normally exposed to high levels of DNA-damaging agents. Whether such high-risk cells possess lineage-specific mechanisms that tailor DNA repair to the tissue remains largely unknown. Using melanoma as a model, we show here that the microphthalmia-associated transcription factor MITF, a lineage addition oncogene that coordinates many aspects of melanocyte and melanoma biology, plays a nontranscriptional role in shaping the DDR. On exposure to DNA-damaging agents, MITF is phosphorylated at S325, and its interactome is dramatically remodeled; most transcription cofactors dissociate, and instead MITF interacts with the MRE11-RAD50-NBS1 (MRN) complex. Consequently, cells with high MITF levels accumulate stalled replication forks and display defects in homologous recombination-mediated repair associated with impaired MRN recruitment to DNA damage. In agreement with this, high MITF levels are associated with increased single-nucleotide and copy number variant burdens in melanoma. Significantly, the SUMOylation-defective MITF-E318K melanoma predisposition mutation recapitulates the effects of DNA-PKcs-phosphorylated MITF. Our data suggest that a nontranscriptional function of a lineage-restricted transcription factor contributes to a tissue-specialized modulation of the DDR that can impact cancer initiation.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Dynamic regulation of chromatin accessibility during melanocyte stem cell activation

Melanocyte stem cells (McSCs) of the hair follicle are necessary for hair pigmentation and can serve as melanoma cells of origin when harboring cancer-driving mutations. McSCs can be released from quiescence, activated, and undergo differentiation into pigment-producing melanocytes during the hair cycle or due to environmental stimuli, such as ultraviolet-B (UVB) exposure. However, our current understanding of the mechanisms regulating McSC stemness, activation, and differentiation remains limited. Here, to capture the differing possible states in which murine McSCs can exist, we sorted melanocyte nuclei from quiescent (telogen) skin, skin actively producing hair shafts (anagen), and skin exposed to UVB. With these sorted nuclei, we then utilized single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing (snATAC-seq) and characterized three melanocyte lineages: quiescent McSCs (qMcSCs), activated McSCs (aMcSCs), and differentiated melanocytes (dMCs) that co-exist in all three skin conditions. Furthermore, we successfully identified differentially accessible genes and enriched transcription factor binding motifs for each melanocyte lineage. Our findings reveal potential gene regulators that determine these melanocyte cell states and provide new insights into how aMcSC chromatin states are regulated differently under divergent intrinsic and extrinsic cues. We also provide a publicly available online tool with a user-friendly interface to explore this comprehensive dataset, which will provide a resource for further studies on McSC regulation upon natural or UVB-mediated stem cell activation.

Comprehensive Analysis of Prognostic Value and Immune Infiltration of TFAP2 Family Members in Bladder Cancer from Database and FFPE Sample

Background: Bladder cancer (BLCA) is one of the common malignant tumors worldwide. Recent studies have shown that Transcription factor activating protein-2(TFAP2) family proteins plays a bidirectional regulatory role in the process of tumorigenesis versus evolution by regulating the expression of tumor associated genes. However, little is known about the function of distinct TFAP2s proteins in patient with BLCA. Methods: Formalin-fixed paraffin-embedded (FFPE) sample tissues and clinical data of 240 patients with bladder cancer were collected for immunohistochemical analysis. The Human Protein Atlas, Gene Expression Profiling Interactive Analysis (GEPIA), Shiny Methylation Analysis Resource Tool (SMART), Kaplan-Meier plotter, cBioPortal, Metascape, LinkedOmics, TIMER and CIBERSORT were utilized to analyze differential expression, prognostic value, genetic alteration and immune cell infiltration of TFAP2 family in patients with BLCA. Results: Our study found that TFAP2 family proteins are generally expressed higher in BLCA tissues than in normal tissues. However, they show different trends in the growth, metastasis and survival prognosis of BLCA. TFAP2A and TFAP2C was associated with worse clinical stage and prognosis in BLCA patients, while TFAP2B, TFAP2D and TFAP2E showed the opposite trend. Importantly, the functions of the differentially expressed TFAP2s were primarily related to the developmental process, reproductive process, response to stimulus and immune system process, etc. Moreover, TFAP2 family was significantly correlated with the infiltration of six immune cell types and might regulate TAM polarization. Conclusion: TFAP2 family might be an important regulator of immune cell infiltration and a valuable prognostic biomarker in patients with BLCA.

Structural basis for specific DNA sequence motif recognition by the TFAP2 transcription factors

The TFAP2 family regulates gene expression during differentiation, development, and organogenesis, and includes five homologs in humans. They all possess a highly conserved DNA binding domain (DBD) followed by a helix-span-helix (HSH) domain. The DBD-HSH tandem domain specifically binds to a GCC(N3)GGC consensus sequence, but the precise recognition mechanisms remain unclear. Here, we found that TFAP2 preferred binding to the GCC(N3)GGC sequence, and the pseudo-palindromic GCC and GGC motifs and the length of the central spacer between the two motifs determined their binding specificity. Structural studies revealed that the two flat amphipathic α-helical HSH domains of TFAP2A stacked with each other to form a dimer via hydrophobic interactions, while the stabilized loops from both DBD domains inserted into two neighboring major grooves of the DNA duplex to form base-specific interactions. This specific DNA binding mechanism controlled the length of the central spacer and determined the DNA sequence specificity of TFAP2. Mutations of the TFAP2 proteins are implicated in various diseases. We illustrated that reduction or disruption of the DNA binding ability of the TFAP2 proteins is the primary cause of TFAP2 mutation-associated diseases. Thus, our findings also offer valuable insights into the pathogenesis of disease-associated mutations in TFAP2 proteins.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

DNA damage-induced interaction between a lineage addiction oncogenic transcription factor and the MRN complex shapes a tissue-specific DNA Damage Response and cancer predisposition

Since genome instability can drive cancer initiation and progression, cells have evolved highly effective and ubiquitous DNA Damage Response (DDR) programs. However, some cells, in skin for example, are normally exposed to high levels of DNA damaging agents. Whether such high-risk cells possess lineage-specific mechanisms that tailor DNA repair to the tissue remains largely unknown. Here we show, using melanoma as a model, that the microphthalmia-associated transcription factor MITF, a lineage addition oncogene that coordinates many aspects of melanocyte and melanoma biology, plays a non-transcriptional role in shaping the DDR. On exposure to DNA damaging agents, MITF is phosphorylated by ATM/DNA-PKcs, and unexpectedly its interactome is dramatically remodelled; most transcription (co)factors dissociate, and instead MITF interacts with the MRE11-RAD50-NBS1 (MRN) complex. Consequently, cells with high MITF levels accumulate stalled replication forks, and display defects in homologous recombination-mediated repair associated with impaired MRN recruitment to DNA damage. In agreement, high MITF levels are associated with increased SNV burden in melanoma. Significantly, the SUMOylation-defective MITF-E318K melanoma predisposition mutation recapitulates the effects of ATM/DNA-PKcs-phosphorylated MITF. Our data suggest that a non-transcriptional function of a lineage-restricted transcription factor contributes to a tissue-specialised modulation of the DDR that can impact cancer initiation.

Zebrafish pigment cells develop directly from persistent highly multipotent progenitors

Neural crest cells are highly multipotent stem cells, but it remains unclear how their fate restriction to specific fates occurs. The direct fate restriction model hypothesises that migrating cells maintain full multipotency, whilst progressive fate restriction envisages fully multipotent cells transitioning to partially-restricted intermediates before committing to individual fates. Using zebrafish pigment cell development as a model, we show applying NanoString hybridization single cell transcriptional profiling and RNAscope in situ hybridization that neural crest cells retain broad multipotency throughout migration and even in post-migratory cells in vivo, with no evidence for partially-restricted intermediates. We find that leukocyte tyrosine kinase early expression marks a multipotent stage, with signalling driving iridophore differentiation through repression of fate-specific transcription factors for other fates. We reconcile the direct and progressive fate restriction models by proposing that pigment cell development occurs directly, but dynamically, from a highly multipotent state, consistent with our recently-proposed Cyclical Fate Restriction model.

Epigenetic and pharmacological control of pigmentation via Bromodomain Protein 9 (BRD9)

Lineage-specific differentiation programs are activated by epigenetic changes in chromatin structure. Melanin-producing melanocytes maintain a gene expression program ensuring appropriate enzymatic conversion of metabolites into the pigment, melanin, and transfer to surrounding cells. During neuroectodermal development, SMARCA4 (BRG1), the catalytic subunit of SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodeling complexes, is essential for lineage specification. SMARCA4 is also required for development of multipotent neural crest precursors into melanoblasts, which differentiate into pigment-producing melanocytes. In addition to the catalytic domain, SMARCA4 and several SWI/SNF subunits contain bromodomains which are amenable to pharmacological inhibition. We investigated the effects of pharmacological inhibitors of SWI/SNF bromodomains on melanocyte differentiation. Strikingly, treatment of murine melanoblasts and human neonatal epidermal melanocytes with selected bromodomain inhibitors abrogated melanin synthesis and visible pigmentation. Using functional genomics, iBRD9, a small molecule selective for the bromodomain of BRD9 was found to repress pigmentation-specific gene expression. Depletion of BRD9 confirmed a requirement for expression of pigmentation genes in the differentiation program from melanoblasts into pigmented melanocytes and in melanoma cells. Chromatin immunoprecipitation assays showed that iBRD9 disrupts the occupancy of BRD9 and the catalytic subunit SMARCA4 at melanocyte-specific loci. These data indicate that BRD9 promotes melanocyte pigmentation whereas pharmacological inhibition of BRD9 is repressive.