DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

DNA binding and mitotic phosphorylation protect polyglutamine proteins from assembly formation

Polyglutamine (polyQ) expansion is associated with pathogenic protein aggregation in neurodegenerative disorders. However, long polyQ tracts are also found in many transcription factors (TFs), such as FOXP2, a TF implicated in human speech. Here, we explore how FOXP2 and other glutamine-rich TFs avoid unscheduled assembly. Throughout interphase, DNA binding, irrespective of sequence specificity, has a solubilizing effect. During mitosis, multiple phosphorylation events promote FOXP2's eviction from chromatin and supplant the solubilizing function of DNA. Further, human-specific amino acid substitutions linked to the evolution of speech map to a mitotic phospho-patch, the "EVO patch," and reduce the propensity of the human FOXP2 to assemble. Fusing the pathogenic form of Huntingtin to either a DNA-binding domain, a phosphomimetic variant of this EVO patch, or a negatively charged peptide is sufficient to diminish assembly formation, suggesting that hijacking mechanisms governing solubility of glutamine-rich TFs may offer new strategies for treatment of polyQ expansion diseases.

The ALX4 dimer structure provides insight into how disease alleles impact function

How homeodomain proteins gain sufficient DNA binding specificity to regulate diverse processes is a long-standing question. Here, we determine how the ALX4 Paired-like protein achieves DNA binding specificity for a TAAT-NNN-ATTA dimer site. We first show that ALX4 binds this motif independently of its co-factor, TWIST1, in cranial neural crest cells. Structural analysis identifies seven ALX4 residues that participate in dimer binding, many of which are conserved across the Paired-like family, but not other homeodomain proteins. Unexpectedly, the two ALX4 proteins within the dimer use distinct residues to form asymmetric protein-protein and protein-DNA interactions and mediate cooperativity. Moreover, we find that ALX4 cooperativity is required for transcriptional activation and that ALX4 disease variants cause distinct molecular defects that include loss of cooperativity. These findings provide insights into how Paired-like factors gain DNA specificity and show how disease variants can be stratified based on their molecular defects.

Auricular malformations are driven by copy number variations in a hierarchical enhancer cluster and a dominant enhancer recapitulates human pathogenesis

Enhancers, through the combinatorial action of transcription factors (TFs), dictate both the spatial specificity and the levels of gene expression, and their aberrations can result in diseases. While a HMX1 downstream enhancer is associated with ear malformations, the mechanisms underlying bilateral constricted ear (BCE) remain unclear. Here, we identify a copy number variation (CNV) containing three enhancers-collectively termed the positional identity hierarchical enhancer cluster (PI-HEC)-that drives BCE by coordinately regulating HMX1 expression. Each enhancer exhibits distinct activity-location-structure features, and the dominant enhancer with high mobility group (HMG)-box combined with Coordinator and homeodomain TF motifs modulating its activity and specificity, respectively. Mouse models demonstrate that neural crest-derived fibroblasts with aberrant Hmx1 expression in the basal pinna, along with ectopic distal pinna expression, disrupt outer ear development, affecting cartilage, muscle, and epidermis. Our findings elucidate mammalian ear morphogenesis and underscore the complexity of synergistic regulation among enhancers and between enhancers and transcription factors.

Transcription factor cooperativity at a GATA3 tandem DNA sequence determines oncogenic enhancer-mediated activation

The TAL1 oncogene driving T cell lymphoblastic leukemia is frequently activated through mutated cis-regulatory elements, whereby small insertions or deletions (indels) create a binding site for the transcription factor MYB. Unraveling how non-coding mutations create oncogenic enhancers is key to understanding cancer biology and can provide important insights into fundamental mechanisms of gene regulation. Utilizing a CRISPR-Cas9 screening approach, we identify GATA3 as the key transcriptional regulator of enhancer-mediated TAL1 overexpression. CRISPR-Cas9 engineering of the mutant enhancer reveals a tandem GATA3 site that is required for binding of GATA3, chromatin accessibility, and MYB recruitment. Reciprocally, MYB binding to its motif is required for GATA3 recruitment, consistent with a transcription factor cooperativity model. Importantly, we show that GATA3 stabilizes a TAL1-MYB interaction and that complex formation requires GATA3 binding to DNA. Our work sheds light on the mechanisms of enhancer-mediated oncogene activation, where key transcription factors cooperate to achieve maximal transcriptional output, thereby supporting leukemogenesis.

Thalidomide-induced limb malformations: an update and reevaluation

Historically, thalidomide-induced congenital malformations have served as an important example of the enhanced susceptibility of developing embryos to chemical perturbation. The compound produced a wide variety of congenital malformations in humans, which were initially detected by an association with a relatively rare limb defect labeled phocomelia. Although true phocomelia in the most severe form is a transverse defect with intercalary absence of limb regions, it is proposed that thalidomide produces a longitudinal limb phenotype in humans under usual circumstances that can become transverse in severe cases with a preferential sensitivity of forelimb over hindlimb, preaxial over postaxial, and left more impacted than the corresponding non-autopod limb bones on the right. The thalidomide-induced limb phenotype in humans is described and followed by a hierarchical comparison with various laboratory animal species. Mechanistic studies have been hampered by the fact that only non-human primates and rabbits have malformations that are anatomically similar to humans. Included in this review are unpublished data on limb malformations produced by thalidomide in rhesus monkeys from experiments performed more than 50 years ago. The critical period in gestation for the induction of phocomelia may initiate prior to the development of the embryonic limb bud, which contrasts with other chemical and physical agents that are known to produce this phenotype. The importance of toxicokinetic parameters is reviewed including dose, enantiomers, absorption, distribution, and both non-enzymatic and enzymatic biotransformations. The limb embryopathy mechanism that provides a partial explanation of the limb phenotype is that cereblon binds to thalidomide creating a protein complex that ubiquitinates protein substrates (CRL4CRBN) that are not targets for the complex in the absence of the thalidomide. One of these neosubstrates is SALL4 which when mutated causes a syndrome that phenocopies aspects of thalidomide embryopathy. Other candidate neosubstrates for the complex that have been found in non-human species may contribute to an understanding of the limb defect including PLZF, p63, and various zinc finger transcription factors. It is proposed that it is important to consider the species-specificity of the compound when considering potential mechanistic pathways and that some of the more traditional mechanisms for explaining the embryopathy, such as anti-angiogenesis and redox perturbation, may contribute to a full understanding of this teratogen.

DNA-guided transcription factor interactions extend human gene regulatory code

In the same way that the mRNA-binding specificities of transfer RNAs define the genetic code, the DNA-binding specificities of transcription factors (TFs) form the molecular basis of the gene regulatory code1,2. The human gene regulatory code is much more complex than the genetic code, in particular because there are more than 1,600 TFs that commonly interact with each other. TF-TF interactions are required for specifying cell fate and executing cell-type-specific transcriptional programs. Despite this, the landscape of interactions between DNA-bound TFs is poorly defined. Here we map the biochemical interactions between DNA-bound TFs using CAP-SELEX, a method that can simultaneously identify individual TF binding preferences, TF-TF interactions and the DNA sequences that are bound by the interacting complexes. A screen of more than 58,000 TF-TF pairs identified 2,198 interacting TF pairs, 1,329 of which preferentially bound to their motifs arranged in a distinct spacing and/or orientation. We also discovered 1,131 TF-TF composite motifs that were markedly different from the motifs of the individual TFs. In total, we estimate that the screen identified between 18% and 47% of all human TF-TF motifs. The novel composite motifs we found were enriched in cell-type-specific elements, active in vivo and more likely to be formed between developmentally co-expressed TFs. Furthermore, TFs that define embryonic axes commonly interacted with different TFs and bound to distinct motifs, explaining how TFs with a similar specificity can define distinct cell types along developmental axes.

Ginsenoside Rb1 Ameliorates Heart Failure Ventricular Remodeling by Regulating the Twist1/PGC-1α/PPARα Signaling Pathway

Background: Heart failure (HF), the terminal stage of cardiovascular disease with high morbidity and mortality, remains poorly managed by current therapies. Ventricular remodeling in HF is fundamentally characterized by myocardial fibrosis. While ginsenoside Rb1 has demonstrated anti-fibrotic effects in HF, the underlying mechanism remains unclear. Twist1, an upstream regulator of energy metabolism factors PGC-1α and PPARα, may attenuate fibrosis by preserving systemic energy homeostasis, suggesting its pivotal role in HF pathogenesis. This study explores ginsenoside Rb1's anti-HF mechanisms through the regulation of ginsenoside Rb1 on these metabolic regulators. Methods: Sprague Dawley rats were subjected to a ligation of the left anterior descending coronary artery to induce an HF model, followed by ginsenoside Rb1 treatment for 6 weeks. Therapeutic effects were evaluated through cardiac function assessment, myocardial histopathological staining (HE, Masson, immunofluorescence, immunohistochemistry), mitochondrial morphology observation (transmission electron microscopy), energy metabolism analysis (electron transport chain efficiency, mitochondrial membrane potential, ATP content), and protein expression profiling (Twist1, PGC-1α, PPARα, GLUT4, PPARγ). Additionally, H9c2 cells induced with endothelin-1 to model HF were employed as an in vitro model to further investigate ginsenoside Rb1's regulatory effects on the Twist1/PGC-1α/PPARα signaling pathway. Results: Ginsenoside Rb1 can restore cardiac function in HF rats, improve mitochondrial function, alleviate energy metabolism disorders, and inhibit ventricular remodeling. By modulating the Twist1/PGC-1α/PPARα signaling pathway, ginsenoside Rb1 suppressed the abnormal overexpression of Twist1 and maintained normal expression of downstream PGC-1α and PPARα. In vitro experiments further demonstrated that ginsenoside Rb1 significantly inhibited Twist1 expression in H9c2 cardiomyocytes with HF while promoting PGC-1α and PPARα expression, thereby restoring myocardial energy metabolism and mitigating ventricular remodeling in HF. Conclusions: Ginsenoside Rb1 can inhibit the upregulation of Twist1 and activate the expression of its downstream PGC-1α and PPARα expression, by modulating the Twist1/PGC-1α/PPARα signaling pathway, alleviating ventricular remodeling in HF patients and improving myocardial energy metabolism dysfunction. Twist1 may be a key target for the treatment of HF. This study not only elucidates the mechanism by which ginsenoside Rb1 alleviates HF, but also provides new insights into the clinical treatment of HF.

The molecular properties of the bHLH TCF4 protein as an intrinsically disordered hub transcription factor

BACKGROUND

Transcription factor 4 (TCF4) is a member of the basic helix-loop-helix (bHLH) family of transcription factors that guides proper embryogenesis, particularly neurogenesis, myogenesis, heart development and hematopoiesis. The interaction of TCF4 with DNA is dependent on the presence of a conserved bHLH domain, particularly the presence of a basic (b) motif. Most mutations in the Tcf4 gene are either associated with the development of serious nervous system disorders, such as Pitt-Hopkins syndrome or schizophrenia, or are lethal. Although TCF4 is essential for the proper development and function of the human body, there is a lack of fundamental knowledge about the structure of TCF4 since structural studies were previously limited exclusively to its bHLH.

METHODS

Recombinant full-length TCF4 was expressed in bacterial cells and purified using chromatographic techniques. To compare the properties of TCF4 in its apo and holo form, we determined the dissociation constant (KD) of the TCF4:DNA complex using independent methods, including fluorescence polarization (FP), electrophoretic mobility shift assay (EMSA), and fluorescence correlation spectroscopy (FCS). Then we compared the properties of TCF4 in its apo and holo form in relation to the changes of the conformation of the polypeptide chain (hydrogen/deuterium exchange mass spectrometry; HDX-MS), hydrodynamic properties (e.g., sedimentation-velocity analytical ultracentrifugation; SV-AUC), and stability (thermal shift, circular dichroism; CD).

RESULTS

We demonstrate the molecular characteristics of TCF4, the dimer of which is one of the largest intrinsically disordered proteins (IDPs) described to date. According to our findings, the structure of TCF4 is extensively disordered. Only the bHLH domain exhibits a stable fold. Strikingly, Ephrussi-box (E-box) binding via the bHLH domain has no significant effect on the disordered nature of TCF4, but it does influence the dynamic of bHLH and stability of the protein.

CONCLUSIONS

We suggest that bHLH plays the role of an anchor localizing TCF4 to specific gene sequences. The dual nature of the TCF4 structure and the fact that the intrinsically disordered regions (IDRs) represent most of the protein sequence, suggest that TCF4 may act as a hub transcription factor regulating the expression of specific genes through the interaction of IDRs with gene-specific partners.

Enhanced insights into the genetic architecture of 3D cranial vault shape using pleiotropy-informed GWAS

Large-scale GWAS studies have uncovered hundreds of genomic loci linked to facial and brain shape variation, but only tens associated with cranial vault shape, a largely overlooked aspect of the craniofacial complex. Surrounding the neocortex, the cranial vault plays a central role during craniofacial development and understanding its genetics are pivotal for understanding craniofacial conditions. Experimental biology and prior genetic studies have generated a wealth of knowledge that presents opportunities to aid further genetic discovery efforts. Here, we use the conditional FDR method to leverage GWAS data of facial shape, brain shape, and bone mineral density to enhance SNP discovery for cranial vault shape. This approach identified 120 independent genomic loci at 1% FDR, nearly tripling the number discovered through unconditioned analysis and implicating crucial craniofacial transcription factors and signaling pathways. These results significantly advance our genetic understanding of cranial vault shape and craniofacial development more broadly.

Transfer learning reveals sequence determinants of the quantitative response to transcription factor dosage

Deep learning models have advanced our ability to predict cell-type-specific chromatin patterns from transcription factor (TF) binding motifs, but their application to perturbed contexts remains limited. We applied transfer learning to predict how concentrations of the dosage-sensitive TFs TWIST1 and SOX9 affect regulatory element (RE) chromatin accessibility in facial progenitor cells, achieving near-experimental accuracy. High-affinity motifs that allow for heterotypic TF co-binding and are concentrated at the center of REs buffer against quantitative changes in TF dosage and predict unperturbed accessibility. Conversely, low-affinity or homotypic binding motifs distributed throughout REs drive sensitive responses with minimal impact on unperturbed accessibility. Both buffering and sensitizing features display purifying selection signatures. We validated these sequence features through reporter assays and demonstrated that TF-nucleosome competition can explain low-affinity motifs' sensitizing effects. This combination of transfer learning and quantitative chromatin response measurements provides a novel approach for uncovering additional layers of the cis-regulatory code.

Distinct gene regulatory dynamics drive skeletogenic cell fate convergence during vertebrate embryogenesis

Cell type repertoires have expanded extensively in metazoan animals, with some clade-specific cells being crucial to evolutionary success. A prime example are the skeletogenic cells of vertebrates. Depending on anatomical location, these cells originate from three different precursor lineages, yet they converge developmentally towards similar cellular phenotypes. Furthermore, their 'skeletogenic competency' arose at distinct evolutionary timepoints, thus questioning to what extent different skeletal body parts rely on truly homologous cell types. Here, we investigate how lineage-specific molecular properties are integrated at the gene regulatory level, to allow for skeletogenic cell fate convergence. Using single-cell functional genomics, we find that distinct transcription factor profiles are inherited from the three precursor states and incorporated at lineage-specific enhancer elements. This lineage-specific regulatory logic suggests that these regionalized skeletogenic cells are distinct cell types, rendering them amenable to individualized selection, to define adaptive morphologies and biomaterial properties in different parts of the vertebrate skeleton.

Twist1 Acts Upstream of the Dlx5-Hand2 Pathway to Pattern the Mammalian Jaw

Both the upper and lower jaws develop from cranial neural crest cells (CNCCs) populating the first pharyngeal arch in all gnathostomes. Previous studies showed that the Edn1/Ednra-Dlx5/Dlx6-Hand2 signaling pathway is necessary for lower jaw formation and that ectopic expression of Edn1 or Hand2 throughout the CNCCs partly transformed the upper jaw to lower jaw structures, but the molecular mechanisms regulating upper jaw development remain unclear. Here we show that the basic helix-loop-helix transcription factor Twist1 is required for upper jaw development. Whereas the Twist1fl/fl;Wnt1-Cre mouse embryos, with tissue-specific inactivation of Twist1 in premigratory CNCCs, exhibited aberrantly persistent expression of the key neuroglial lineage regulator Sox10 in the postmigratory CNCCs populating the facial primordia, we found that genetic inactivation of Sox10 did not rescue the defects in CNCC survival and patterning in Twist1fl/fl;Wnt1-Cre embryos. However, analysis of Sox10fl/+;Twist1fl/fl;Wnt1-Cre mice revealed duplicated mandibular structures, including ectopic Meckel's cartilage, in place of the maxilla. Both Sox10fl/+;Twist1fl/fl;Wnt1-Cre and Sox10fl/fl;Twist1fl/fl;Wnt1-Cre embryos exhibited ectopic expression of Dlx5 and Hand2 in the developing maxillary processes at E10.5. Furthermore, we found that Twist1fl/fl;Wnt1-Cre embryos also expressed Dlx5 and Hand2 ectopically in the maxillary domain at E10.5 and subsequently developed Meckel's cartilage-like cartilage rods bilaterally at the maxillary region. However, the expression of Edn1 was unaltered in the developing Twist1fl/fl;Wnt1-Cre embryos, indicating that Twist1 functions in the CNCC-derived facial mesenchyme to regulate the Dlx5-Hand2 pathway without affecting Edn1 expression in the epithelial and mesodermal compartments. We further show that Twist1 represses reporter gene activation driven by the Dlx5/Dlx6 intergenic enhancer known to drive Dlx5/Dlx6 expression in the developing mandibular arch. Together, these data identify a new role of Twist1 in patterning the regional identities of the CNCC-derived facial mesenchyme and provide novel insight into the pathogenic mechanisms underlying TWIST1-related craniofacial developmental disorders.

Mechanisms of enhanced or impaired DNA target selectivity driven by protein dimerization

Successful DNA transcription demands coordination between proteins that bind DNA while simultaneously binding to one another into dimers or higher-order complexes. Measurements that report on the lifetime or occupancy of an individual protein on DNA thus represent a convolution over the protein interactions with specific DNA, nonspecific DNA, or protein partners on DNA. For some DNA-binding proteins, dimerization is considered an essential step for stable DNA association, but here we show that protein dimerization can also reduce dwell times on specific DNA targets, enhance or impair occupancy on target sequences, and spatially redistribute proteins on DNA. We use mass-action kinetic models of pairwise association reactions between proteins and DNA (specific and nonspecific) and protein dimers, along with theory and spatial stochastic simulations to isolate the role of dimerization on observed dwell times and occupancy. For proteins binding a spatially localized cluster of targets, dimerization can drive up dwell time by 1000-fold and produce high selectivity for clustered over isolated targets. However, this effect can become negligible when proteins outnumber target sequences. In contrast, for isolated DNA targets, dimerization often reduces dwell times by sequestering proteins from their target sites, in some cases thus reducing overall occupancy. The ability of these proteins to bind DNA nonspecifically and diffuse in 1D to exploit dimensional reduction is a key determinant controlling degree of enhancement, despite the presence of nucleosome barriers to 1D diffusion. By comparison with ChipSeq data, our model explains how the distribution of the GAF pioneer proteins throughout the genome is highly selective for clustered targets due to protein interactions and provides a framework to predict how even weak dimerization can redistribute or stabilize proteins on DNA.

Liquid condensates: a new barrier to loop extrusion?

Liquid-liquid phase separation (LLPS), driven by dynamic, low-affinity multivalent interactions of proteins and RNA, results in the formation of macromolecular condensates on chromatin. These structures are likely to provide high local concentrations of effector factors responsible for various processes including transcriptional regulation and DNA repair. In particular, enhancers, super-enhancers, and promoters serve as platforms for condensate assembly. In the current paradigm, enhancer-promoter (EP) interaction could be interpreted as a result of enhancer- and promoter-based condensate contact/fusion. There is increasing evidence that the spatial juxtaposition of enhancers and promoters could be provided by loop extrusion (LE) by SMC complexes. Here, we propose that condensates may act as barriers to LE, thereby contributing to various nuclear processes including spatial contacts between regulatory genomic elements.

Combinatorial transcription factor binding encodes cis-regulatory wiring of mouse forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to cis-regulatory elements, orchestrating transcriptional programs. Although studies of chromatin state and chromosomal interactions have demonstrated dynamic neurodevelopmental cis-regulatory landscapes, parallel understanding of TF interactions lags. To elucidate combinatorial TF binding driving mouse basal ganglia development, we integrated chromatin immunoprecipitation sequencing (ChIP-seq) for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and gene expression data, and functional enhancer assays. We identified sets of putative regulatory elements with shared TF binding (TF-pRE modules) that orchestrate distinct processes of GABAergic neurogenesis and suppress other cell fates. The majority of pREs were bound by one or two TFs; however, a small proportion were extensively bound. These sequences had exceptional evolutionary conservation and motif density, complex chromosomal interactions, and activity as in vivo enhancers. Our results provide insights into the combinatorial TF-pRE interactions that activate and repress expression programs during telencephalon neurogenesis and demonstrate the value of TF binding toward modeling developmental transcriptional wiring.

Tissue-specific chromatin accessibility and transcriptional regulation in maize cold stress response

Maize, a vital crop globally, faces significant yield losses due to its sensitivity to cold stress, especially in temperate regions. Understanding the molecular mechanisms governing maize response to cold stress is crucial for developing strategies to enhance cold tolerance. However, the precise chromatin-level regulatory mechanisms involved remain largely unknown. In this study, we employed DNase-seq and RNA-seq techniques to investigate chromatin accessibility and gene expression changes in maize root, stem, and leaf tissues subjected to cold treatment. We discovered widespread changes in chromatin accessibility and gene expression across these tissues, with strong tissue specificity. Cold stress-induced DNase I hypersensitive sites (coiDHSs) were associated with differentially expressed genes, suggesting a direct link between chromatin accessibility and gene regulation under cold stress. Motif enrichment analysis identified ERF transcription factors (TFs) as central regulators conserved across tissues, with ERF5 emerging as pivotal in the cold response regulatory network. Additionally, TF co-localization analysis highlighted six TF pairs (ERF115-SHN3, ERF9-LEP, ERF7-SHN3, LEP-SHN3, LOB-SHN3, and AS2-LOB) conserved across tissues but showing tissue-specific binding preferences. These findings indicate intricate regulatory networks in maize cold response. Overall, our study provides insights into the chromatin-level regulatory mechanisms underpinning maize adaptive response to cold stress, offering potential targets for enhancing cold tolerance in agricultural contexts.

Target gene responses differ when transcription factor levels are acutely decreased by nuclear export versus degradation

Defining the time of action for morphogens requires tools capable of temporally controlled perturbations. To study how the transcription factor Dorsal affects patterning of the Drosophila embryonic dorsal-ventral axis, we used two light-inducible tags that trigger either nuclear export or degradation of Dorsal under blue light. Nuclear export of Dorsal leads to loss of the high-threshold, ventrally expressed target gene snail (sna), while the low-threshold, laterally expressed target gene short-gastrulation (sog) is retained. In contrast, degradation of Dorsal results in retention of sna, loss of sog, and lower nuclear levels compared to when Dorsal is exported from the nucleus. To understand why nuclear export causes loss of sna but degradation does not, we investigated Dorsal kinetics using photobleaching and found that it rapidly re-enters the nucleus even under blue-light conditions favoring export. The associated kinetics of Dorsal being rapidly imported and exported continuously are likely responsible for loss of sna but, alternatively, can support sog. Collectively, our results indicate that this dynamic patterning process is influenced by both Dorsal concentration and nuclear retention.

A DNA base-specific sequence interposed between CRX and NRL contributes to RHODOPSIN expression

Gene expression emerges from DNA sequences through the interaction of transcription factors (TFs) with DNA cis-regulatory sequences. In eukaryotes, TFs bind to transcription factor binding sites (TFBSs) with differential affinities, enabling cell-specific gene expression. In this view, DNA enables TF binding along a continuum ranging from low to high affinity depending on its sequence composition; however, it is not known whether evolution has entailed a further level of entanglement between DNA-protein interaction. Here we found that the composition and length (22 bp) of the DNA sequence interposed between the CRX and NRL retinal TFs in the proximal promoter of RHODOPSIN (RHO) largely controls the expression levels of RHO. Mutagenesis of CRX-NRL DNA linking sequences (here termed "DNA-linker") results in uncorrelated gene expression variation. In contrast, mutual exchange of naturally occurring divergent human and mouse Rho cis-regulatory elements conferred similar yet species-specific Rho expression levels. Two orthogonal DNA-binding proteins targeted to the DNA-linker either activate or repress the expression of Rho depending on the DNA-linker orientation relative to the CRX and NRL binding sites. These results argue that, in this instance, DNA itself contributes to CRX and NRL activities through a code based on specific base sequences of a defined length, ultimately determining optimal RHO expression levels.

Detection of regions of homozygosity in an unusual case of frontonasal dysplasia

We present the case of a 7-year-old Ecuadorian mestizo girl with multiple orofacial malformations. The patient is the product of a first-degree relationship (father-daughter). A cytogenetic study revealed a normal karyotype. The genetic mapping array study identified 0.73 Gb of alterations, 727,087,295 bp involved in regions of homozygosity (ROH) in all chromosomes (25.2% of the genome) and 764,028 bp in gains in chromosomes 9 and 14. Genes from the TGFB, BMP, FGF, SHH and WNT families, among others, were identified in the ROH. They are related to craniofacial development and their protein products showed a strong association in the interactome analysis.

ARID1A-BAF coordinates ZIC2 genomic occupancy for epithelial-to-mesenchymal transition in cranial neural crest specification

The BAF chromatin remodeler regulates lineage commitment including cranial neural crest cell (CNCC) specification. Variants in BAF subunits cause Coffin-Siris syndrome (CSS), a congenital disorder characterized by coarse craniofacial features and intellectual disability. Approximately 50% of individuals with CSS harbor variants in one of the mutually exclusive BAF subunits, ARID1A/ARID1B. While Arid1a deletion in mouse neural crest causes severe craniofacial phenotypes, little is known about the role of ARID1A in CNCC specification. Using CSS-patient-derived ARID1A+/- induced pluripotent stem cells to model CNCC specification, we discovered that ARID1A-haploinsufficiency impairs epithelial-to-mesenchymal transition (EMT), a process necessary for CNCC delamination and migration from the neural tube. Furthermore, wild-type ARID1A-BAF regulates enhancers associated with EMT genes. ARID1A-BAF binding at these enhancers is impaired in heterozygotes while binding at promoters is unaffected. At the sequence level, these EMT enhancers contain binding motifs for ZIC2, and ZIC2 binding at these sites is ARID1A-dependent. When excluded from EMT enhancers, ZIC2 relocates to neuronal enhancers, triggering aberrant neuronal gene activation. In mice, deletion of Zic2 impairs NCC delamination, while ZIC2 overexpression in chick embryos at post-migratory neural crest stages elicits ectopic delamination from the neural tube. These findings reveal an essential ARID1A-ZIC2 axis essential for EMT and CNCC delamination.

Pioneer factors: Emerging rules of engagement for transcription factors on chromatinized DNA

Pioneering transcription factors (TFs) can drive cell fate changes by binding their DNA motifs in a repressive chromatin environment. Recent structures illustrate emerging rules for nucleosome engagement: TFs distort the nucleosomal DNA to gain access or employ alternative DNA-binding modes with smaller footprints, they preferentially access solvent-exposed motifs near the entry/exit sites, and frequently interact with histones. The extent of TF-histone interactions, in turn, depends on the motif location on the nucleosome, the type of DNA-binding fold, and adjacent domains present. TF-histone interactions can phase TF motifs relative to nucleosomes, and we discuss how these complex and surprisingly diverse interactions between nucleosomes and TFs contribute to function.

Prediction and design of transcriptional repressor domains with large-scale mutational scans and deep learning

Regulatory proteins have evolved diverse repressor domains (RDs) to enable precise context-specific repression of transcription. However, our understanding of how sequence variation impacts the functional activity of RDs is limited. To address this gap, we generated a high-throughput mutational scanning dataset measuring the repressor activity of 115,000 variant sequences spanning more than 50 RDs in human cells. We identified thousands of clinical variants with loss or gain of repressor function, including TWIST1 HLH variants associated with Saethre-Chotzen syndrome and MECP2 domain variants associated with Rett syndrome. We also leveraged these data to annotate short linear interacting motifs (SLiMs) that are critical for repression in disordered RDs. Then, we designed a deep learning model called TENet ( T ranscriptional E ffector Net work) that integrates sequence, structure and biochemical representations of sequence variants to accurately predict repressor activity. We systematically tested generalization within and across domains with varying homology using the mutational scanning dataset. Finally, we employed TENet within a directed evolution sequence editing framework to tune the activity of both structured and disordered RDs and experimentally test thousands of designs. Our work highlights critical considerations for future dataset design and model training strategies to improve functional variant prioritization and precision design of synthetic regulatory proteins.

Mutagenesis Sensitivity Mapping of Human Enhancers In Vivo

Distant-acting enhancers are central to human development. However, our limited understanding of their functional sequence features prevents the interpretation of enhancer mutations in disease. Here, we determined the functional sensitivity to mutagenesis of human developmental enhancers in vivo. Focusing on seven enhancers active in the developing brain, heart, limb and face, we created over 1700 transgenic mice for over 260 mutagenized enhancer alleles. Systematic mutation of 12-basepair blocks collectively altered each sequence feature in each enhancer at least once. We show that 69% of all blocks are required for normal in vivo activity, with mutations more commonly resulting in loss (60%) than in gain (9%) of function. Using predictive modeling, we annotated critical nucleotides at base-pair resolution. The vast majority of motifs predicted by these machine learning models (88%) coincided with changes to in vivo function, and the models showed considerable sensitivity, identifying 59% of all functional blocks. Taken together, our results reveal that human enhancers contain a high density of sequence features required for their normal in vivo function and provide a rich resource for further exploration of human enhancer logic.

Synergistic activation by Glass and Pointed promotes neuronal identity in the Drosophila eye disc

The integration of extrinsic signaling with cell-intrinsic transcription factors can direct progenitor cells to differentiate into distinct cell fates. In the developing Drosophila eye, differentiation of photoreceptors R1-R7 requires EGFR signaling mediated by the transcription factor Pointed, and our single-cell RNA-Seq analysis shows that the same photoreceptors require the eye-specific transcription factor Glass. We find that ectopic expression of Glass and activation of EGFR signaling synergistically induce neuronal gene expression in the wing disc in a Pointed-dependent manner. Targeted DamID reveals that Glass and Pointed share many binding sites in the genome of developing photoreceptors. Comparison with transcriptomic data shows that Pointed and Glass induce photoreceptor differentiation through intermediate transcription factors, including the redundant homologs Scratch and Scrape, as well as directly activating neuronal effector genes. Our data reveal synergistic activation of a multi-layered transcriptional network as the mechanism by which EGFR signaling induces neuronal identity in Glass-expressing cells.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Transfer learning reveals sequence determinants of the quantitative response to transcription factor dosage

Deep learning approaches have made significant advances in predicting cell type-specific chromatin patterns from the identity and arrangement of transcription factor (TF) binding motifs. However, most models have been applied in unperturbed contexts, precluding a predictive understanding of how chromatin state responds to TF perturbation. Here, we used transfer learning to train and interpret deep learning models that use DNA sequence to predict, with accuracy approaching experimental reproducibility, how the concentration of two dosage-sensitive TFs (TWIST1, SOX9) affects regulatory element (RE) chromatin accessibility in facial progenitor cells. High-affinity motifs that allow for heterotypic TF co-binding and are concentrated at the center of REs buffer against quantitative changes in TF dosage and strongly predict unperturbed accessibility. In contrast, motifs with low-affinity or homotypic binding distributed throughout REs lead to sensitive responses with minimal contributions to unperturbed accessibility. Both buffering and sensitizing features show signatures of purifying selection. We validated these predictive sequence features using reporter assays and showed that a biophysical model of TF-nucleosome competition can explain the sensitizing effect of low-affinity motifs. Our approach of combining transfer learning and quantitative measurements of the chromatin response to TF dosage therefore represents a powerful method to reveal additional layers of the cis-regulatory code.

Target gene responses differ when transcription factor levels are acutely decreased by nuclear export versus degradation

Defining the time of action for morphogens requires tools capable of temporally controlled perturbations. To study how the transcription factor Dorsal affects patterning of the Drosophila embryonic dorsal-ventral axis, we used two light-inducible tags that result in either nuclear export or degradation of Dorsal when exposed to blue light. Nuclear export of Dorsal results in loss of expression for the high threshold, ventrally-expressed target gene snail (sna) but retention of the low threshold, laterally-expressed target gene short-gastrulation (sog). In contrast, degradation of Dorsal results in retention of sna, loss of sog, and lower nuclear levels than when Dorsal is exported from the nucleus. To elucidate how nuclear export results in loss of sna but degradation does not, we investigated Dorsal kinetics using photobleaching and found it reenters the nucleus even under conditions of blue-light when export is favored. The associated kinetics of being imported and exported continuously are likely responsible for loss of sna but, alternatively, can support sog. Collectively, our results show that this dynamic patterning process is influenced by both Dorsal concentration and nuclear retention.

Spatial enhancer activation influences inhibitory neuron identity during mouse embryonic development

The mammalian telencephalon contains distinct GABAergic projection neuron and interneuron types, originating in the germinal zone of the embryonic basal ganglia. How genetic information in the germinal zone determines cell types is unclear. Here we use a combination of in vivo CRISPR perturbation, lineage tracing and ChIP-sequencing analyses and show that the transcription factor MEIS2 favors the development of projection neurons by binding enhancer regions in projection-neuron-specific genes during mouse embryonic development. MEIS2 requires the presence of the homeodomain transcription factor DLX5 to direct its functional activity toward the appropriate binding sites. In interneuron precursors, the transcription factor LHX6 represses the MEIS2-DLX5-dependent activation of projection-neuron-specific enhancers. Mutations of Meis2 result in decreased activation of regulatory enhancers, affecting GABAergic differentiation. We propose a differential binding model where the binding of transcription factors at cis-regulatory elements determines differential gene expression programs regulating cell fate specification in the mouse ganglionic eminence.

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.