DNA-mediated cooperativity facilitates the co-selection of cryptic enhancer sequences by SOX2 and PAX6 transcription factors

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

Transcriptional kinetic synergy: A complex landscape revealed by integrating modeling and synthetic biology

Transcription factors (TFs) control gene expression, often acting synergistically. Classical thermodynamic models offer a biophysical explanation for synergy based on binding cooperativity and regulated recruitment of RNA polymerase. Because transcription requires polymerase to transition through multiple states, recent work suggests that "kinetic synergy" can arise through TFs acting on distinct steps of the transcription cycle. These types of synergy are not mutually exclusive and are difficult to disentangle conceptually and experimentally. Here, we model and build a synthetic circuit in which TFs bind to a single shared site on DNA, such that TFs cannot synergize by simultaneous binding. We model mRNA production as a function of both TF binding and regulation of the transcription cycle, revealing a complex landscape dependent on TF concentration, DNA binding affinity, and regulatory activity. We use synthetic TFs to confirm that the transcription cycle must be integrated with recruitment for a quantitative understanding of gene regulation.

Improving cooperativity of transcription activators by oligomerization domains in mammalian cells

Cooperative activation is critical for the applications of synthetic biology in mammalian cells. In this study, we have developed cooperative transcription factor by fusing oligomerization domain in mammalian cells. Firstly, we demonstrated that two oligomerized domains (CI434 and CI) successfully improved transcription factor cooperativity in bacterial cells but failed to increase cooperativity in mammalian cells, possibly because the additional mammalian activation domain disrupted their oligomerization capability. Therefore, we chose a different type of oligomerized domain (CarHC), whose ability to oligomerize is not dependent on its C-terminal domains, to fuse with a transcription factor (RpaR) and activation domain (VTR3), forming a potential cooperative transcription activator RpaR-CarH-VTR3 for mammalian regulatory systems. Compared with RpaR-VTR3, the cooperativity of RpaR-CarH-VTR3 was significantly improved with higher Hill coefficient and a narrower input range in the inducible switch system in mammalian cells. Moreover, a mathematical model based on statistical mechanics model was developed and the simulation results supported the hypothesis that the tetramer of the CarH domain in mammalian cells was the reason for the cooperative capacity of RpaR-CarH-VTR3.

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

Statistical estimates of multiple transcription factors binding in the model plant genomes based on ChIP-seq data

The development of high-throughput genomic sequencing coupled with chromatin immunoprecipitation technologies allows studying the binding sites of the protein transcription factors (TF) in the genome scale. The growth of data volume on the experimentally determined binding sites raises qualitatively new problems for the analysis of gene expression regulation, prediction of transcription factors target genes, and regulatory gene networks reconstruction. Genome regulation remains an insufficiently studied though plants have complex molecular regulatory mechanisms of gene expression and response to environmental stresses. It is important to develop new software tools for the analysis of the TF binding sites location and their clustering in the plant genomes, visualization, and the following statistical estimates. This study presents application of the analysis of multiple TF binding profiles in three evolutionarily distant model plant organisms. The construction and analysis of non-random ChIP-seq binding clusters of the different TFs in mammalian embryonic stem cells were discussed earlier using similar bioinformatics approaches. Such clusters of TF binding sites may indicate the gene regulatory regions, enhancers and gene transcription regulatory hubs. It can be used for analysis of the gene promoters as well as a background for transcription networks reconstruction. We discuss the statistical estimates of the TF binding sites clusters in the model plant genomes. The distributions of the number of different TFs per binding cluster follow same power law distribution for all the genomes studied. The binding clusters in Arabidopsis thaliana genome were discussed here in detail.

Recurrent heterozygous PAX6 missense variants cause severe bilateral microphthalmia via predictable effects on DNA-protein interaction

PURPOSE

Most classical aniridia is caused by PAX6 haploinsufficiency. PAX6 missense variants can be hypomorphic or mimic haploinsufficiency. We hypothesized that missense variants also cause previously undescribed disease by altering the affinity and/or specificity of PAX6 genomic interactions.

METHODS

We screened PAX6 in 372 individuals with bilateral microphthalmia, anophthalmia, or coloboma (MAC) from the Medical Research Council Human Genetics Unit eye malformation cohort (HGUeye) and reviewed data from the Deciphering Developmental Disorders study. We performed cluster analysis on PAX6-associated ocular phenotypes by variant type and molecular modeling of the structural impact of 86 different PAX6 causative missense variants.

RESULTS

Eight different PAX6 missense variants were identified in 17 individuals (15 families) with MAC, accounting for 4% (15/372) of our cohort. Seven altered the paired domain (p.[Arg26Gln]x1, p.[Gly36Val]x1, p.[Arg38Trp]x2, p.[Arg38Gln]x1, p.[Gly51Arg]x2, p.[Ser54Arg]x2, p.[Asn124Lys]x5) and one the homeodomain (p.[Asn260Tyr]x1). p.Ser54Arg and p.Asn124Lys were exclusively associated with severe bilateral microphthalmia. MAC-associated variants were predicted to alter but not ablate DNA interaction, consistent with the electrophoretic mobility shifts observed using mutant paired domains with well-characterized PAX6-binding sites. We found no strong evidence for novel PAX6-associated extraocular disease.

CONCLUSION

Altering the affinity and specificity of PAX6-binding genome-wide provides a plausible mechanism for the worse-than-null effects of MAC-associated missense variants.

Mechanistic insights into transcription factor cooperativity and its impact on protein-phenotype interactions

Recent high-throughput transcription factor (TF) binding assays revealed that TF cooperativity is a widespread phenomenon. However, a global mechanistic and functional understanding of TF cooperativity is still lacking. To address this, here we introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data. We identify DNA shape as driver for cooperativity, with a particularly strong effect for Forkhead-Ets pairs. Follow-up experiments reveal a local shape preference at the Ets-DNA-Forkhead interface and decreased cooperativity upon loss of the interaction. Additionally, we discover many functional associations for cooperatively bound TFs. Examination of the link between FOXO1:ETV6 and lymphomas reveals that their joint expression levels improve patient clinical outcome stratification. Altogether, our results demonstrate that inter-family cooperative TF binding is driven by position-specific DNA readout mechanisms, which provides an additional regulatory layer for downstream biological functions.

Although transcription factor (TF) cooperativity is widespread, a global mechanistic understanding of the role of TF cooperativity is still lacking. Here the authors introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data and provide mechanistic insights into TF cooperativity and its impact on protein-phenotype interactions.

DNA-mediated dimerization on a compact sequence signature controls enhancer engagement and regulation by FOXA1

FOXA1 is a transcription factor capable to bind silenced chromatin to direct context-dependent cell fate conversion. Here, we demonstrate that a compact palindromic DNA element (termed 'DIV' for its diverging half-sites) induces the homodimerization of FOXA1 with strongly positive cooperativity. Alternative structural models are consistent with either an indirect DNA-mediated cooperativity or a direct protein-protein interaction. The cooperative homodimer formation is strictly constrained by precise half-site spacing. Re-analysis of chromatin immunoprecipitation sequencing data indicates that the DIV is effectively targeted by FOXA1 in the context of chromatin. Reporter assays show that FOXA1-dependent transcriptional activity declines when homodimeric binding is disrupted. In response to phosphatidylinositol-3 kinase inhibition DIV sites pre-bound by FOXA1 such as at the PVT1/MYC locus exhibit a strong increase in accessibility suggesting a role of the DIV configuration in the chromatin closed-open dynamics. Moreover, several disease-associated single nucleotide polymorphisms map to DIV elements and show allelic differences in FOXA1 homodimerization, reporter gene expression and are annotated as quantitative trait loci. This includes the rs541455835 variant at the MAPT locus encoding the Tau protein associated with Parkinson's disease. Collectively, the DIV guides chromatin engagement and regulation by FOXA1 and its perturbation could be linked to disease etiologies.

OCT4 and PAX6 determine the dual function of SOX2 in human ESCs as a key pluripotent or neural factor

Background

Sox2 is a well-established pluripotent transcription factor that plays an essential role in establishing and maintaining pluripotent stem cells (PSCs). It is also thought to be a linage specifier that governs PSC neural lineage specification upon their exiting the pluripotent state. However, the exact role of SOX2 in human PSCs was still not fully understood. In this study, we studied the role of SOX2 in human embryonic stem cells (hESCs) by gain- and loss-of-function approaches and explored the possible underlying mechanisms.

Results

We demonstrate that knockdown of SOX2 induced hESC differentiation to endoderm-like cells, whereas overexpression of SOX2 in hESCs enhanced their pluripotency under self-renewing culture conditions but promoted their neural differentiation upon replacing the culture to non-self-renewal conditions. We show that this culture-dependent dual function of SOX2 was probably attributed to its interaction with different transcription factors predisposed by the culture environments. Whilst SOX2 interacts with OCT4 under self-renewal conditions, we found that, upon neural differentiation, PAX6, a key neural transcription factor, is upregulated and shows interaction with SOX2. The SOX2-PAX6 complex has different gene regulation pattern from that of SOX2-OCT4 complex.

Conclusions

Our work provides direct evidence that SOX2 is necessarily required for hESC pluripotency; however, it can also function as a neural factor, depending on the environmental input. OCT4 and PAX6 might function as key SOX2-interacting partners that determine the function of SOX2 in hESCs.

Electronic supplementary material

The online version of this article (10.1186/s13287-019-1228-7) contains supplementary material, which is available to authorized users.

Allosterism and signal transfer in DNA

We analysed the basic mechanisms of signal transmission in DNA and the origins of the allostery exhibited by systems such as the ternary complex BAMHI-DNA-GRDBD. We found that perturbation information generated by a primary protein binding event travels as a wave to distant regions of DNA following a hopping mechanism. However, such a structural perturbation is transient and does not lead to permanent changes in the DNA geometry and interaction properties at the secondary binding site. The BAMHI-DNA-GRDBD allosteric mechanism does not occur through any traditional models: direct (protein-protein), indirect (reorganization of the secondary site) readout or solvent-release. On the contrary, it is generated by a subtle and less common entropy-mediated mechanism, which might have an important role to explain other DNA-mediated cooperative effects.

Detection of cooperatively bound transcription factor pairs using ChIP-seq peak intensities and expectation maximization

Transcription factors (TFs) often work cooperatively, where the binding of one TF to DNA enhances the binding affinity of a second TF to a nearby location. Such cooperative binding is important for activating gene expression from promoters and enhancers in both prokaryotic and eukaryotic cells. Existing methods to detect cooperative binding of a TF pair rely on analyzing the sequence that is bound. We propose a method that uses, instead, only ChIP-seq peak intensities and an expectation maximization (CPI-EM) algorithm. We validate our method using ChIP-seq data from cells where one of a pair of TFs under consideration has been genetically knocked out. Our algorithm relies on our observation that cooperative TF-TF binding is correlated with weak binding of one of the TFs, which we demonstrate in a variety of cell types, including E. coli, S. cerevisiae and M. musculus cells. We show that this method performs significantly better than a predictor based only on the ChIP-seq peak distance of the TFs under consideration. This suggests that peak intensities contain information that can help detect the cooperative binding of a TF pair. CPI-EM also outperforms an existing sequence-based algorithm in detecting cooperative binding. The CPI-EM algorithm is available at https://github.com/vishakad/cpi-em.

Identification of a Wells-Dawson polyoxometalate-based AP-2γ inhibitor with pro-apoptotic activity

AP-2 gamma (AP-2γ) is a transcription factor that plays pivotal roles in breast cancer biology. To search for small molecule inhibitors of AP-2γ, we performed a high-throughput fluorescence anisotropy screen and identified a polyoxometalate compound with Wells-Dawson structure K₆[P₂Mo₁₈O₆₂] (Dawson-POM) that blocks the DNA-binding activity of AP-2γ. We showed that this blocking activity is due to the direct binding of Dawson-POM to AP-2γ. We also provided evidence to show that Dawson-POM decreases AP-2γ-dependent transcription similar to silencing the gene. Finally, we demonstrated that Dawson-POM contains anti-proliferative and pro-apoptotic effects in breast cancer cells. In summary, we identified the first small molecule inhibitor of AP-2γ and showed Dawson-POM-mediated inhibition of AP-2γ as a potential avenue for cancer therapy.

Diversity among POU transcription factors in chromatin recognition and cell fate reprogramming

The POU (Pit-Oct-Unc) protein family is an evolutionary ancient group of transcription factors (TFs) that bind specific DNA sequences to direct gene expression programs. The fundamental importance of POU TFs to orchestrate embryonic development and to direct cellular fate decisions is well established, but the molecular basis for this activity is insufficiently understood. POU TFs possess a bipartite 'two-in-one' DNA binding domain consisting of two independently folding structural units connected by a poorly conserved and flexible linker. Therefore, they represent a paradigmatic example to study the molecular basis for the functional versatility of TFs. Their modular architecture endows POU TFs with the capacity to accommodate alternative composite DNA sequences by adopting different quaternary structures. Moreover, associations with partner proteins crucially influence the selection of their DNA binding sites. The plentitude of DNA binding modes confers the ability to POU TFs to regulate distinct genes in the context of different cellular environments. Likewise, different binding modes of POU proteins to DNA could trigger alternative regulatory responses in the context of different genomic locations of the same cell. Prominent POU TFs such as Oct4, Brn2, Oct6 and Brn4 are not only essential regulators of development but have also been successfully employed to reprogram somatic cells to pluripotency and neural lineages. Here we review biochemical, structural, genomic and cellular reprogramming studies to examine how the ability of POU TFs to select regulatory DNA, alone or with partner factors, is tied to their capacity to epigenetically remodel chromatin and drive specific regulatory programs that give cells their identities.

Cellular, ultrastructural and molecular analyses of epidermal cell development in the planarian Schmidtea mediterranea

The epidermis is essential for animal survival, providing both a protective barrier and cellular sensor to external environments. The generally conserved embryonic origin of the epidermis, but the broad morphological and functional diversity of this organ across animals is puzzling. We define the transcriptional regulators underlying epidermal lineage differentiation in the planarian Schmidtea mediterranea, an invertebrate organism that, unlike fruitflies and nematodes, continuously replaces its epidermal cells. We find that Smed-p53, Sox and Pax transcription factors are essential regulators of epidermal homeostasis, and act cooperatively to regulate genes associated with early epidermal precursor cell differentiation, including a tandemly arrayed novel gene family (prog) of secreted proteins. Additionally, we report on the discovery of distinct and previously undescribed secreted organelles whose production is dependent on the transcriptional activity of soxP-3, and which we term Hyman vesicles.

Coop-Seq Analysis Demonstrates that Sox2 Evokes Latent Specificities in the DNA Recognition by Pax6

Sox2 and Pax6 co-regulate genes in neural lineages and the lens by forming a ternary complex likely facilitated allosterically through DNA. We used the quantitative and scalable cooperativity-by-sequencing (Coop-seq) approach to interrogate Sox2/Pax6 dimerization on a DNA library where five positions of the Pax6 half-site were randomized yielding 1024 cooperativity factors. Consensus positions normally required for the high-affinity DNA binding by Pax6 need to be mutated for effective dimerization with Sox2. Out of the five randomized bases, a 5' thymidine is present in most of the top ranking elements. However, this thymidine maps to a region outside of the Pax half site and is not expected to directly interact with Pax6 in known binding modes suggesting structural reconfigurations. Re-analysis of ChIP-seq data identified several genomic regions where the cooperativity promoting sequence pattern is co-bound by Sox2 and Pax6. A highly conserved Sox2/Pax6 bound site near the Sprouty2 locus was verified to promote cooperative dimerization designating Sprouty2 as a potential target reliant on Sox2/Pax6 cooperativity in several neural cell types. Collectively, the functional interplay of Sox2 and Pax6 demands the relaxation of high-affinity binding sites and is enabled by alternative DNA sequences. We conclude that this binding mode evolved to warrant that a subset of target genes is only regulated in the presence of suitable partner factors.

Direct reprogramming with SOX factors: masters of cell fate

Over the last decade significant advances have been made toward reprogramming the fate of somatic cells, typically by overexpression of cell lineage-determinant transcription factors. As key regulators of cell fate, the SOX family of transcription factors has emerged as potent drivers of direct somatic cell reprogramming into multiple lineages, in some cases as the sole overexpressed factor. The vast capacity of SOX factors, especially those of the SOXB1, E and F subclasses, to reprogram cell fate is enlightening our understanding of organismal development, cancer and disease, and offers tremendous potential for regenerative medicine and cell-based therapies. Understanding the molecular mechanisms through which SOX factors reprogram cell fate is essential to optimize the development of novel somatic cell transdifferentiation strategies.

Quantitative profiling of selective Sox/POU pairing on hundreds of sequences in parallel by Coop-seq

Cooperative binding of transcription factors is known to be important in the regulation of gene expression programs conferring cellular identities. However, current methods to measure cooperativity parameters have been laborious and therefore limited to studying only a few sequence variants at a time. We developed Coop-seq (cooperativity by sequencing) that is capable of efficiently and accurately determining the cooperativity parameters for hundreds of different DNA sequences in a single experiment. We apply Coop-seq to 12 dimer pairs from the Sox and POU families of transcription factors using 324 unique sequences with changed half-site orientation, altered spacing and discrete randomization within the binding elements. The study reveals specific dimerization profiles of different Sox factors with Oct4. By contrast, Oct4 and the three neural class III POU factors Brn2, Brn4 and Oct6 assemble with Sox2 in a surprisingly indistinguishable manner. Two novel half-site configurations can support functional Sox/Oct dimerization in addition to known composite motifs. Moreover, Coop-seq uncovers a nucleotide switch within the POU half-site when spacing is altered, which is mirrored in genomic loci bound by Sox2/Oct4 complexes.

Pax6 associates with H3K4-specific histone methyltransferases Mll1, Mll2, and Set1a and regulates H3K4 methylation at promoters and enhancers

BACKGROUND

Pax6 is a key regulator of the entire cascade of ocular lens formation through specific binding to promoters and enhancers of batteries of target genes. The promoters and enhancers communicate with each other through DNA looping mediated by multiple protein-DNA and protein-protein interactions and are marked by specific combinations of histone posttranslational modifications (PTMs). Enhancers are distinguished from bulk chromatin by specific modifications of core histone H3, including H3K4me1 and H3K27ac, while promoters show increased H3K4me3 PTM. Previous studies have shown the presence of Pax6 in as much as 1/8 of lens-specific enhancers but a much smaller fraction of tissue-specific promoters. Although Pax6 is known to interact with EP300/p300 histone acetyltransferase responsible for generation of H3K27ac, a potential link between Pax6 and histone H3K4 methylation remains to be established.

RESULTS

Here we show that Pax6 co-purifies with H3K4 methyltransferase activity in lens cell nuclear extracts. Proteomic studies show that Pax6 immunoprecipitates with Set1a, Mll1, and Mll2 enzymes, and their associated proteins, i.e., Wdr5, Rbbp5, Ash2l, and Dpy30. ChIP-seq studies using chromatin prepared from mouse lens and cultured lens cells demonstrate that Pax6-bound regions are mostly enriched with H3K4me2 and H3K4me1 in enhancers and promoters, though H3K4me3 marks only Pax6-containing promoters. The shRNA-mediated knockdown of Pax6 revealed down-regulation of a set of direct target genes, including Cap2, Farp1, Pax6, Plekha1, Prox1, Tshz2, and Zfp536. Pax6 knockdown was accompanied by reduced H3K4me1 at enhancers and H3K4me3 at promoters, with little or no changes of the H3K4me2 modifications. These changes were prominent in Plekha1, a gene regulated by Pax6 in both lens and retinal pigmented epithelium.

CONCLUSIONS

Our study supports a general model of Pax6-mediated recruitment of histone methyltransferases Mll1 and Mll2 to lens chromatin, especially at distal enhancers. Genome-wide data in lens show that Pax6 binding correlates with H3K4me2, consistent with the idea that H3K4me2 PTMs correlate with the binding of transcription factors. Importantly, partial reduction of Pax6 induces prominent changes in local H3K4me1 and H3K4me3 modification. Together, these data open the field to mechanistic studies of Pax6, Mll1, Mll2, and H3K4me1/2/3 dynamics at distal enhancers and promoters of developmentally controlled genes.

Examining cooperative binding of Sox2 on DC5 regulatory element upon complex formation with Pax6 through excess electron transfer assay

Functional cooperativity among transcription factors on regulatory genetic elements is pivotal for milestone decision-making in various cellular processes including mammalian development. However, their molecular interaction during the cooperative binding cannot be precisely understood due to lack of efficient tools for the analyses of protein-DNA interaction in the transcription complex. Here, we demonstrate that photoinduced excess electron transfer assay can be used for analysing cooperativity of proteins in transcription complex using cooperative binding of Pax6 to Sox2 on the regulatory DNA element (DC5 enhancer) as an example. In this assay, (Br)U-labelled DC5 was introduced for the efficient detection of transferred electrons from Sox2 and Pax6 to the DNA, and guanine base in the complementary strand was replaced with hypoxanthine (I) to block intra-strand electron transfer at the Sox2-binding site. By examining DNA cleavage occurred as a result of the electron transfer process, from tryptophan residues of Sox2 and Pax6 to DNA after irradiation at 280 nm, we not only confirmed their binding to DNA but also observed their increased occupancy on DC5 with respect to that of Sox2 and Pax6 alone as a result of their cooperative interaction.

Mechanisms of Specificity for Hox Factor Activity

Metazoans encode clusters of paralogous Hox genes that are critical for proper development of the body plan. However, there are a number of unresolved issues regarding how paralogous Hox factors achieve specificity to control distinct cell fates. First, how do Hox paralogs, which have very similar DNA binding preferences in vitro, drive different transcriptional programs in vivo? Second, the number of potential Hox binding sites within the genome is vast compared to the number of sites bound. Hence, what determines where in the genome Hox factors bind? Third, what determines whether a Hox factor will activate or repress a specific target gene? Here, we review the current evidence that is beginning to shed light onto these questions. In particular, we highlight how cooperative interactions with other transcription factors (especially PBC and HMP proteins) and the sequences of cis-regulatory modules provide a basis for the mechanisms of Hox specificity. We conclude by integrating a number of the concepts described throughout the review in a case study of a highly interrogated Drosophila cis-regulatory module named “The Distal-less Conserved Regulatory Element” (DCRE).

PAX6: 25th anniversary and more to learn

The DNA-binding transcription factor PAX6 was cloned 25 years ago by multiple teams pursuing identification of human and mouse eye disease causing genes, cloning vertebrate homologues of pattern-forming regulatory genes identified in Drosophila, or abundant eye-specific transcripts. Since its discovery in 1991, genetic, cellular, molecular and evolutionary studies on Pax6 mushroomed in the mid 1990s leading to the transformative thinking regarding the genetic program orchestrating both early and late stages of eye morphogenesis as well as the origin and evolution of diverse visual systems. Since Pax6 is also expressed outside of the eye, namely in the central nervous system and pancreas, a number of important insights into the development and function of these organs have been amassed. In most recent years, genome-wide technologies utilizing massively parallel DNA sequencing have begun to provide unbiased insights into the regulatory hierarchies of specification, determination and differentiation of ocular cells and neurogenesis in general. This review is focused on major advancements in studies on mammalian eye development driven by studies of Pax6 genes in model organisms and future challenges to harness the technology-driven opportunities to reconstruct, step-by-step, the transition from naïve ectoderm, neuroepithelium and periocular mesenchyme/neural crest cells into the three-dimensional architecture of the eye.

Glucose-dependent downregulation of glucagon gene expression mediated by selective interactions between ALX3 and PAX6 in mouse alpha cells

AIMS/HYPOTHESIS

The stimulation of glucagon secretion in response to decreased glucose levels has been studied extensively. In contrast, little is known about the regulation of glucagon gene expression in response to fluctuations in glucose concentration. Paired box 6 (PAX6) is a key transcription factor that regulates the glucagon promoter by binding to the G1 and G3 elements. Here, we investigated the role of the transcription factor aristaless-like homeobox 3 (ALX3) as a glucose-dependent modulator of PAX6 activity in alpha cells.

METHODS

Experiments were performed in wild-type or Alx3-deficient islets and alphaTC1 cells. We used chromatin immunoprecipitations and electrophoretic mobility shift assays for DNA binding, immunoprecipitations and pull-down assays for protein interactions, transfected cells for promoter activity, and small interfering RNA and quantitative RT-PCR for gene expression.

RESULTS

Elevated glucose concentration resulted in stimulated expression of Alx3 and decreased glucagon gene expression in wild-type islets. In ALX3-deficient islets, basal glucagon levels were non-responsive to changes in glucose concentration. In basal conditions ALX3 bound to the glucagon promoter at G3, but not at G1. ALX3 could form heterodimers with PAX6 that were permissive for binding to G3 but not to G1. Thus, increasing the levels of ALX3 in response to glucose resulted in the sequestration of PAX6 by ALX3 for binding to G1, thus reducing glucagon promoter activation and glucagon gene expression.

CONCLUSIONS/INTERPRETATION

Glucose-stimulated expression of ALX3 in alpha cells provides a regulatory mechanism for the downregulation of glucagon gene expression by interfering with PAX6-mediated transactivation on the glucagon G1 promoter element.

Structure and decoy-mediated inhibition of the SOX18/Prox1-DNA interaction

The transcription factor (TF) SOX18 drives lymphatic vessel development in both embryogenesis and tumour-induced neo-lymphangiogenesis. Genetic disruption of Sox18 in a mouse model protects from tumour metastasis and established the SOX18 protein as a molecular target. Here, we report the crystal structure of the SOX18 DNA binding high-mobility group (HMG) box bound to a DNA element regulating Prox1 transcription. The crystals diffracted to 1.75Å presenting the highest resolution structure of a SOX/DNA complex presently available revealing water structure, structural adjustments at the DNA contact interface and non-canonical conformations of the DNA backbone. To explore alternatives to challenging small molecule approaches for targeting the DNA-binding activity of SOX18, we designed a set of five decoys based on modified Prox1-DNA. Four decoys potently inhibited DNA binding of SOX18 in vitro and did not interact with non-SOX TFs. Serum stability, nuclease resistance and thermal denaturation assays demonstrated that a decoy circularized with a hexaethylene glycol linker and terminal phosphorothioate modifications is most stable. This SOX decoy also interfered with the expression of a luciferase reporter under control of a SOX18-dependent VCAM1 promoter in COS7 cells. Collectively, we propose SOX decoys as potential strategy for inhibiting SOX18 activity to disrupt tumour-induced neo-lymphangiogenesis.

Mapping gene regulatory circuitry of Pax6 during neurogenesis

Pax6 is a highly conserved transcription factor among vertebrates and is important in various aspects of the central nervous system development. However, the gene regulatory circuitry of Pax6 underlying these functions remains elusive. We find that Pax6 targets a large number of promoters in neural progenitors cells. Intriguingly, many of these sites are also bound by another progenitor factor, Sox2, which cooperates with Pax6 in gene regulation. A combinatorial analysis of Pax6-binding data set with transcriptome changes in Pax6-deficient neural progenitors reveals a dual role for Pax6, in which it activates the neuronal (ectodermal) genes while concurrently represses the mesodermal and endodermal genes, thereby ensuring the unidirectionality of lineage commitment towards neuronal differentiation. Furthermore, Pax6 is critical for inducing activity of transcription factors that elicit neurogenesis and repress others that promote non-neuronal lineages. In addition to many established downstream effectors, Pax6 directly binds and activates a number of genes that are specifically expressed in neural progenitors but have not been previously implicated in neurogenesis. The in utero knockdown of one such gene, Ift74, during brain development impairs polarity and migration of newborn neurons. These findings demonstrate new aspects of the gene regulatory circuitry of Pax6, revealing how it functions to control neuronal development at multiple levels to ensure unidirectionality and proper execution of the neurogenic program.

Structural and conformational insights into SOX2/OCT4-bound enhancer DNA: a computational perspective

The roles of SOX2 and OCT4 are critical in stem cell maintenance either in the context of iPSCs generation or cancer stem cell growth; therefore, it is imperative to study their cooperative binding and SOX2/OCT4-induced DNA conformational switching.

Pax factors in transcription and epigenetic remodelling

The nine Pax transcription factors that constitute the mammalian family of paired domain (PD) factors play key roles in many developmental processes. As DNA binding transcription factors, they exhibit tremendous variability and complexity in their DNA recognition patterns. This is ascribed to the presence of multiple DNA binding structural domains, namely helix-turn-helix (HTH) domains. The PD contains two HTH subdomains and four of the nine Pax factors have an additional HTH domain, the homeodomain (HD). We now review these diverse DNA binding modalities together with their properties as transcriptional activators and repressors. The action of Pax factors on gene expression is also exerted through recruitment of chromatin remodelling complexes that introduce either activating or repressive chromatin marks. Interestingly, the recent demonstration that Pax7 has pioneer activity, the unique property to "open" chromatin, further underlines the mechanistic versatility and the developmental importance of these factors.

Cooperative DNA Recognition Modulated by an Interplay between Protein-Protein Interactions and DNA-Mediated Allostery

Highly specific transcriptional regulation depends on the cooperative association of transcription factors into enhanceosomes. Usually, their DNA-binding cooperativity originates from either direct interactions or DNA-mediated allostery. Here, we performed unbiased molecular simulations followed by simulations of protein-DNA unbinding and free energy profiling to study the cooperative DNA recognition by OCT4 and SOX2, key components of enhanceosomes in pluripotent cells. We found that SOX2 influences the orientation and dynamics of the DNA-bound configuration of OCT4. In addition SOX2 modifies the unbinding free energy profiles of both DNA-binding domains of OCT4, the POU specific and POU homeodomain, despite interacting directly only with the first. Thus, we demonstrate that the OCT4-SOX2 cooperativity is modulated by an interplay between protein-protein interactions and DNA-mediated allostery. Further, we estimated the change in OCT4-DNA binding free energy due to the cooperativity with SOX2, observed a good agreement with experimental measurements, and found that SOX2 affects the relative DNA-binding strength of the two OCT4 domains. Based on these findings, we propose that available interaction partners in different biological contexts modulate the DNA exploration routes of multi-domain transcription factors such as OCT4. We consider the OCT4-SOX2 cooperativity as a paradigm of how specificity of transcriptional regulation is achieved through concerted modulation of protein-DNA recognition by different types of interactions.

Pluripotent stem cells can give rise to all somatic lineages. When taken out of the context of the embryo they can be maintained and for this a core transcriptional regulatory circuitry is crucial. OCT4 and SOX2, two factors of this network, are also critical for the induction of pluripotency in somatic cells. In pluripotent cells, OCT4 and SOX2 associate on DNA regulatory regions, enhancing or modifying each other's sequence specificity. In contrast, in the early stages during induction of pluripotency, it was proposed that OCT4 explores the genome independent of SOX2. Here we report the mechanism by which SOX2 influences the orientation, dynamics, and unbinding free energy profile of OCT4. This involves an interplay of protein-protein interactions and DNA-mediated allostery. We consider that this mechanism enables OCT4 to use its DNA binding domains and the interaction partners available in a certain biological context to access alternative genome exploration routes. This study enhances the understanding of the context specific function of OCT4 and provides a general perspective on how DNA-binding cooperativity is modulated by different types of interactions.

SOXE transcription factors form selective dimers on non-compact DNA motifs through multifaceted interactions between dimerization and high-mobility group domains

The SOXE transcription factors SOX8, SOX9 and SOX10 are master regulators of mammalian development directing sex determination, gliogenesis, pancreas specification and neural crest development. We identified a set of palindromic SOX binding sites specifically enriched in regulatory regions of melanoma cells. SOXE proteins homodimerize on these sequences with high cooperativity. In contrast to other transcription factor dimers, which are typically rigidly spaced, SOXE group proteins can bind cooperatively at a wide range of dimer spacings. Using truncated forms of SOXE proteins, we show that a single dimerization (DIM) domain, that precedes the DNA binding high mobility group (HMG) domain, is sufficient for dimer formation, suggesting that DIM : HMG rather than DIM:DIM interactions mediate the dimerization. All SOXE members can also heterodimerize in this fashion, whereas SOXE heterodimers with SOX2, SOX4, SOX6 and SOX18 are not supported. We propose a structural model where SOXE-specific intramolecular DIM:HMG interactions are allosterically communicated to the HMG of juxtaposed molecules. Collectively, SOXE factors evolved a unique mode to combinatorially regulate their target genes that relies on a multifaceted interplay between the HMG and DIM domains. This property potentially extends further the diversity of target genes and cell-specific functions that are regulated by SOXE proteins.