Precise temporal control of the eye regulatory gene Pax6 via enhancer-binding site affinity

Pax6- and Six3-mediated induction of lens cell fate in mouse and human ES cells

Embryonic stem (ES) cells provide a potentially useful in vitro model for the study of in vivo tissue differentiation. We used mouse and human ES cells to investigate whether the lens regulatory genes Pax6 and Six3 could induce lens cell fate in vitro. To help assess the onset of lens differentiation, we derived a new mES cell line (Pax6-GFP mES) that expresses a GFP reporter under the control of the Pax6 P0 promoter and lens ectoderm enhancer. Pax6 or Six3 expression vectors were introduced into mES or hES cells by transfection or lentiviral infection and the differentiating ES cells analyzed for lens marker expression. Transfection of mES cells with Pax6 or Six3 but not with other genes induced the expression of lens cell markers and up-regulated GFP reporter expression in Pax6-GFP mES cells by 3 days post-transfection. By 7 days post-transfection, mES cell cultures exhibited a>10-fold increase over controls in the number of colonies expressing γA-crystallin, a lens fiber cell differentiation marker. RT-PCR and immunostaining revealed induction of additional lens epithelial or fiber cell differentiation markers including Foxe3, Prox1, α- and β-crystallins, and Tdrd7. Moreover, γA-crystallin- or Prox1-expressing lentoid bodies formed by 30 days in culture. In hES cells, Pax6 or Six3 lentiviral vectors also induced lens marker expression. mES cells that express lens markers reside close to but are distinct from the Pax6 or Six3 transduced cells, suggesting that the latter induce nearby undifferentiated ES cells to adopt a lens fate by non-cell autonomous mechanisms. In sum, we describe a novel mES cell GFP reporter line that is useful for monitoring induction of lens fate, and demonstrate that Pax6 or Six3 is sufficient to induce ES cells to adopt a lens fate, potentially via non-cell autonomous mechanisms. These findings should facilitate investigations of lens development.

The cellular and molecular mechanisms of vertebrate lens development

The ocular lens is a model system for understanding important aspects of embryonic development, such as cell specification and the spatiotemporally controlled formation of a three-dimensional structure. The lens, which is characterized by transparency, refraction and elasticity, is composed of a bulk mass of fiber cells attached to a sheet of lens epithelium. Although lens induction has been studied for over 100 years, recent findings have revealed a myriad of extracellular signaling pathways and gene regulatory networks, integrated and executed by the transcription factor Pax6, that are required for lens formation in vertebrates. This Review summarizes recent progress in the field, emphasizing the interplay between the diverse regulatory mechanisms employed to form lens progenitor and precursor cells and highlighting novel opportunities to fill gaps in our understanding of lens tissue morphogenesis.

Modular evolution of DNA-binding preference of a Tbrain transcription factor provides a mechanism for modifying gene regulatory networks

Gene regulatory networks (GRNs) describe the progression of transcriptional states that take a single-celled zygote to a multicellular organism. It is well documented that GRNs can evolve extensively through mutations to cis-regulatory modules (CRMs). Transcription factor proteins that bind these CRMs may also evolve to produce novelty. Coding changes are considered to be rarer, however, because transcription factors are multifunctional and hence are more constrained to evolve in ways that will not produce widespread detrimental effects. Recent technological advances have unearthed a surprising variation in DNA-binding abilities, such that individual transcription factors may recognize both a preferred primary motif and an additional secondary motif. This provides a source of modularity in function. Here, we demonstrate that orthologous transcription factors can also evolve a changed preference for a secondary binding motif, thereby offering an unexplored mechanism for GRN evolution. Using protein-binding microarray, surface plasmon resonance, and in vivo reporter assays, we demonstrate an important difference in DNA-binding preference between Tbrain protein orthologs in two species of echinoderms, the sea star, Patiria miniata, and the sea urchin, Strongylocentrotus purpuratus. Although both orthologs recognize the same primary motif, only the sea star Tbr also has a secondary binding motif. Our in vivo assays demonstrate that this difference may allow for greater evolutionary change in timing of regulatory control. This uncovers a layer of transcription factor binding divergence that could exist for many pairs of orthologs. We hypothesize that this divergence provides modularity that allows orthologous transcription factors to evolve novel roles in GRNs through modification of binding to secondary sites.

Absence of a simple code: how transcription factors read the genome

Transcription factors (TFs) influence cell fate by interpreting the regulatory DNA within a genome. TFs recognize DNA in a specific manner; the mechanisms underlying this specificity have been identified for many TFs based on 3D structures of protein-DNA complexes. More recently, structural views have been complemented with data from high-throughput in vitro and in vivo explorations of the DNA-binding preferences of many TFs. Together, these approaches have greatly expanded our understanding of TF-DNA interactions. However, the mechanisms by which TFs select in vivo binding sites and alter gene expression remain unclear. Recent work has highlighted the many variables that influence TF-DNA binding, while demonstrating that a biophysical understanding of these many factors will be central to understanding TF function.

On the identification of potential regulatory variants within genome wide association candidate SNP sets

BACKGROUND

Genome wide association studies (GWAS) are a population-scale approach to the identification of segments of the genome in which genetic variations may contribute to disease risk. Current methods focus on the discovery of single nucleotide polymorphisms (SNPs) associated with disease traits. As there are many SNPs within identified risk loci, and the majority of these are situated within non-coding regions, a key challenge is to identify and prioritize variants affecting regulatory sequences that are likely to contribute to the phenotype assessed.

METHODS

We focused investigation on SNPs within lung and breast cancer GWAS loci that reached genome-wide significance for potential roles in gene regulation with a specific focus on SNPs likely to disrupt transcription factor binding sites. Within risk loci, the regulatory potential of sub-regions was classified using relevant open chromatin and epigenetic high throughput sequencing data sets from the ENCODE project in available cancer and normal cell lines. Furthermore, transcription factor affinity altering variants were predicted by comparison of position weight matrix scores between disease and reference alleles. Lastly, ChIP-seq data of transcription associated factors and topological domains were included as binding evidence and potential gene target inference.

RESULTS

The sets of SNPs, including both the disease-associated markers and those in high linkage disequilibrium with them, were significantly over-represented in regulatory sequences of cancer and/or normal cells; however, over-representation was generally not restricted to disease-relevant tissue specific regions. The calculated regulatory potential, allelic binding affinity scores and ChIP-seq binding evidence were the three criteria used to prioritize candidates. Fitting all three criteria, we highlighted breast cancer susceptibility SNPs and a borderline lung cancer relevant SNP located in cancer-specific enhancers overlapping multiple distinct transcription associated factor ChIP-seq binding sites.

CONCLUSION

Incorporating high throughput sequencing epigenetic and transcription factor data sets from both cancer and normal cells into cancer genetic studies reveals potential functional SNPs and informs subsequent characterization efforts.

Biochemistry of the tale transcription factors PREP, MEIS, and PBX in vertebrates

TALE (three amino acids loop extension) homeodomain transcription factors are required in various steps of embryo development, in many adult physiological functions, and are involved in important pathologies. This review focuses on the PREP, MEIS, and PBX sub-families of TALE factors and aims at giving information on their biochemical properties, i.e., structure, interactors, and interaction surfaces. Members of the three sets of protein form dimers in which the common partner is PBX but they can also directly interact with other proteins forming higher-order complexes, in particular HOX. Finally, recent advances in determining the genome-wide DNA-binding sites of PREP1, MEIS1, and PBX1, and their partial correspondence with the binding sites of some HOX proteins, are reviewed. These studies have generated a few general rules that can be applied to all members of the three gene families. PREP and MEIS recognize slightly different consensus sequences: PREP prefers to bind to promoters and to have PBX as a DNA-binding partner; MEIS prefers HOX as partner, and both PREP and MEIS drive PBX to their own binding sites. This outlines the clear individuality of the PREP and MEIS proteins, the former mostly devoted to basic cellular functions, the latter more to developmental functions.

Protein-DNA binding: complexities and multi-protein codes

Binding of proteins to particular DNA sites across the genome is a primary determinant of specificity in genome maintenance and gene regulation. DNA-binding specificity is encoded at multiple levels, from the detailed biophysical interactions between proteins and DNA, to the assembly of multi-protein complexes. At each level, variation in the mechanisms used to achieve specificity has led to difficulties in constructing and applying simple models of DNA binding. We review the complexities in protein–DNA binding found at multiple levels and discuss how they confound the idea of simple recognition codes. We discuss the impact of new high-throughput technologies for the characterization of protein–DNA binding, and how these technologies are uncovering new complexities in protein–DNA recognition. Finally, we review the concept of multi-protein recognition codes in which new DNA-binding specificities are achieved by the assembly of multi-protein complexes.

Low-affinity transcription factor binding sites shape morphogen responses and enhancer evolution

In the era of functional genomics, the role of transcription factor (TF)-DNA binding affinity is of increasing interest: for example, it has recently been proposed that low-affinity genomic binding events, though frequent, are functionally irrelevant. Here, we investigate the role of binding site affinity in the transcriptional interpretation of Hedgehog (Hh) morphogen gradients. We noted that enhancers of several Hh-responsive Drosophila genes have low predicted affinity for Ci, the Gli family TF that transduces Hh signalling in the fly. Contrary to our initial hypothesis, improving the affinity of Ci/Gli sites in enhancers of dpp, wingless and stripe, by transplanting optimal sites from the patched gene, did not result in ectopic responses to Hh signalling. Instead, we found that these enhancers require low-affinity binding sites for normal activation in regions of relatively low signalling. When Ci/Gli sites in these enhancers were altered to improve their binding affinity, we observed patterning defects in the transcriptional response that are consistent with a switch from Ci-mediated activation to Ci-mediated repression. Synthetic transgenic reporters containing isolated Ci/Gli sites confirmed this finding in imaginal discs. We propose that the requirement for gene activation by Ci in the regions of low-to-moderate Hh signalling results in evolutionary pressure favouring weak binding sites in enhancers of certain Hh target genes.

Discriminative motif analysis of high-throughput dataset

MOTIVATION

High-throughput ChIP-seq studies typically identify thousands of peaks for a single transcription factor (TF). It is common for traditional motif discovery tools to predict motifs that are statistically significant against a naïve background distribution but are of questionable biological relevance.

RESULTS

We describe a simple yet effective algorithm for discovering differential motifs between two sequence datasets that is effective in eliminating systematic biases and scalable to large datasets. Tested on 207 ENCODE ChIP-seq datasets, our method identifies correct motifs in 78% of the datasets with known motifs, demonstrating improvement in both accuracy and efficiency compared with DREME, another state-of-art discriminative motif discovery tool. More interestingly, on the remaining more challenging datasets, we identify common technical or biological factors that compromise the motif search results and use advanced features of our tool to control for these factors. We also present case studies demonstrating the ability of our method to detect single base pair differences in DNA specificity of two similar TFs. Lastly, we demonstrate discovery of key TF motifs involved in tissue specification by examination of high-throughput DNase accessibility data.

AVAILABILITY

The motifRG package is publically available via the bioconductor repository.

CONTACT

yzizhen@fhcrc.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pax6 downregulation mediates abnormal lineage commitment of the ocular surface epithelium in aqueous-deficient dry eye disease

Keratinizing squamous metaplasia (SQM) of the ocular surface is a blinding consequence of systemic autoimmune disease and there is no cure. Ocular SQM is traditionally viewed as an adaptive tissue response during chronic keratoconjunctivitis sicca (KCS) that provokes pathological keratinization of the corneal epithelium and fibrosis of the corneal stroma. Recently, we established the autoimmune regulator-knockout (Aire KO) mouse as a model of autoimmune KCS and identified an essential role for autoreactive CD4+ T cells in SQM pathogenesis. In subsequent studies, we noted the down-regulation of paired box gene 6 (Pax6) in both human patients with chronic KCS associated with Sjögren’s syndrome and Aire KO mice. Pax6 encodes a pleiotropic transcription factor guiding eye morphogenesis during development. While the postnatal function of Pax6 is largely unknown, we hypothesized that its role in maintaining ocular surface homeostasis was disrupted in the inflamed eye and that loss of Pax6 played a functional role in the initiation and progression of SQM. Adoptive transfer of autoreactive T cells from Aire KO mice to immunodeficient recipients confirmed CD4+ T cells as the principal downstream effectors promoting Pax6 downregulation in Aire KO mice. CD4+ T cells required local signaling via Interleukin-1 receptor (IL-1R1) to provoke Pax6 loss, which prompted a switch from corneal-specific cytokeratin, CK12, to epidermal-specific CK10. The functional role of Pax6 loss in SQM pathogenesis was indicated by the reversal of SQM and restoration of ocular surface homeostasis following forced expression of Pax6 in corneal epithelial cells using adenovirus. Thus, tissue-restricted restoration of Pax6 prevented aberrant epidermal-lineage commitment suggesting adjuvant Pax6 gene therapy may represent a novel therapeutic approach to prevent SQM in patients with chronic inflammatory diseases of the ocular surface.

Genome-wide analyses of Shavenbaby target genes reveals distinct features of enhancer organization

Background

Developmental programs are implemented by regulatory interactions between Transcription Factors (TFs) and their target genes, which remain poorly understood. While recent studies have focused on regulatory cascades of TFs that govern early development, little is known about how the ultimate effectors of cell differentiation are selected and controlled. We addressed this question during late Drosophila embryogenesis, when the finely tuned expression of the TF Ovo/Shavenbaby (Svb) triggers the morphological differentiation of epidermal trichomes.

Results

We defined a sizeable set of genes downstream of Svb and used in vivo assays to delineate 14 enhancers driving their specific expression in trichome cells. Coupling computational modeling to functional dissection, we investigated the regulatory logic of these enhancers. Extending the repertoire of epidermal effectors using genome-wide approaches showed that the regulatory models learned from this first sample are representative of the whole set of trichome enhancers. These enhancers harbor remarkable features with respect to their functional architectures, including a weak or non-existent clustering of Svb binding sites. The in vivo function of each site relies on its intimate context, notably the flanking nucleotides. Two additional cis-regulatory motifs, present in a broad diversity of composition and positioning among trichome enhancers, critically contribute to enhancer activity.

Conclusions

Our results show that Svb directly regulates a large set of terminal effectors of the remodeling of epidermal cells. Further, these data reveal that trichome formation is underpinned by unexpectedly diverse modes of regulation, providing fresh insights into the functional architecture of enhancers governing a terminal differentiation program.

Why detailed model gene studies in higher eukaryotes are still necessary

Alterations in gene expression programmes are controlled by sequence-specific DNA-binding proteins that interact with the epigenetic regulatory machinery. The sum of such processes comprises a gene regulatory network and differentiation processes involve transitions between such networks. However, while great progress has been made to identify network components, this list is not complete, and we still do not fully understand how they work together. In this article, I argue that one reason for this lack of knowledge is the fact that we still do not understand what controls the cell stage and cell state-specific regulation of individual genes and review examples highlighting this notion.

Evaluating thermodynamic models of enhancer activity on cellular resolution gene expression data

With the advent of high throughput sequencing and high resolution transcriptomic technologies, there exists today an unprecedented opportunity to understand gene regulation at a quantitative level. State of the art models of the relationship between regulatory sequence and gene expression have shown great promise, but also suffer from some major shortcomings. In this paper, we identify and address methodological challenges pertaining to quantitative modeling of gene expression from sequence, and test our models on the anterior-posterior patterning system in the Drosophila embryo. We first develop a framework to process cellular resolution three-dimensional gene expression data from the Drosophila embryo and create data sets on which quantitative models can be trained. Next we propose a new score, called 'weighted pattern generating potential' (w-PGP), to evaluate model predictions, and show its advantages over the two most common scoring schemes in use today. The model building exercise uses w-PGP as the evaluation score and adopts a systematic strategy to increase a model's complexity while guarding against over-fitting. Our model identifies three transcription factors--ZELDA, SLOPPY-PAIRED, and NUBBIN--that have not been previously incorporated in quantitative models of this system, as having significant regulatory influence. Finally, we show how fitting quantitative models on data sets comprising a handful of enhancers, as reported in earlier work, may lead to unreliable models.

An allelic series at the paired box gene 6 (Pax6) locus reveals the functional specificity of Pax genes

The advent of the ocular and nervous system in metazoan evolution coincides with the diversification of a single ancestral paired box (Pax) gene into Pax6, Pax6(5a), and Pax2. To investigate the role of these Pax genes in neural development, we have generated an allelic series of knock-in models at the Pax6 locus. We showed that although Pax6(5a) and Pax2 could not replace Pax6 for its autoregulation in lens induction or for neural differentiation in retina, Pax6(5a) was sufficient for corneal-lenticular detachment. In brain development, cell proliferation in the cerebral cortex and dorsoventral patterning of the telencephalon and neural tube were partially rescued in either knock-in mutant. Contrary to the previous belief, our genetic studies showed that the Pax6 isoform Pax6(5a) could potentially play a role in neuronal differentiation in brain development. Importantly, Pax2 showed greater rescue efficiency than Pax6(5a) in the telencephalon even though the latter was identical to Pax6 outside the paired domain. In studying Ngn2, a Pax6 direct target gene in telencephalon, we showed that the level of Ngn2 expression correlated with the in vitro binding of Pax2, Pax6, and Pax6(5a) paired domain on its enhancer. Our results show that Pax6 is uniquely required for eye development, but in brain development, Pax6 can be functionally substituted by related Pax family genes that share a similar paired domain binding specificity.

Neural-specific Sox2 input and differential Gli-binding affinity provide context and positional information in Shh-directed neural patterning

In the vertebrate neural tube, regional Sonic hedgehog (Shh) signaling invokes a time- and concentration-dependent induction of six different cell populations mediated through Gli transcriptional regulators. Elsewhere in the embryo, Shh/Gli responses invoke different tissue-appropriate regulatory programs. A genome-scale analysis of DNA binding by Gli1 and Sox2, a pan-neural determinant, identified a set of shared regulatory regions associated with key factors central to cell fate determination and neural tube patterning. Functional analysis in transgenic mice validates core enhancers for each of these factors and demonstrates the dual requirement for Gli1 and Sox2 inputs for neural enhancer activity. Furthermore, through an unbiased determination of Gli-binding site preferences and analysis of binding site variants in the developing mammalian CNS, we demonstrate that differential Gli-binding affinity underlies threshold-level activator responses to Shh input. In summary, our results highlight Sox2 input as a context-specific determinant of the neural-specific Shh response and differential Gli-binding site affinity as an important cis-regulatory property critical for interpreting Shh morphogen action in the mammalian neural tube.

Unraveling the KNOTTED1 regulatory network in maize meristems

KNOTTED1 (KN1)-like homeobox (KNOX) transcription factors function in plant meristems, self-renewing structures consisting of stem cells and their immediate daughters. We defined the KN1 cistrome in maize inflorescences and found that KN1 binds to several thousand loci, including 643 genes that are modulated in one or multiple tissues. These KN1 direct targets are strongly enriched for transcription factors (including other homeobox genes) and genes participating in hormonal pathways, most significantly auxin, demonstrating that KN1 plays a key role in orchestrating the upper levels of a hierarchical gene regulatory network that impacts plant meristem identity and function.

Reduction of Prep1 levels affects differentiation of normal and malignant B cells and accelerates Myc driven lymphomagenesis

The Prep1 homeodomain transcription factor has recently been recognized as a tumor suppressor. Among other features, haploinsufficiency of Prep1 is able to strongly accelerate the B-lymphomagenesis in EμMyc mice. Now we report that this occurs concomitantly with a change in the type of B-cell lymphomas generated by the Myc oncogene. Indeed, the tumors generated in the EμMyc-Prep1^+/− mice are much more immature, being mostly made up of Pro-B or Pre-B cells, while those in the EμMyc-Prep1^+/+ mice are more differentiated being invariably IgM⁺. Moreover, we show that Prep1 is in fact required for the differentiation of Pro-B and Pre-B cells into IgM⁺ lymphocytes and/or their proliferation, thus showing also how a normal function of Prep1 affects EμMyc lymphomagenesis. Finally, we show that the haploinsufficiency of Prep1 is accompanied with a major decrease of Myc-induced apoptosis and that the haploinsufficieny is sufficient for all these effects because the second allele of Prep1 is not lost even at late stages. Therefore, the tumor-suppressive activity of Prep1 is intertwined with both the interference with Myc-induced apoptosis as well as with natural developmental functions of the protein.

Mutations and binding sites of human transcription factors

Mutations in any genome may lead to phenotype characteristics that determine ability of an individual to cope with adaptation to environmental challenges. In studies of human biology, among the most interesting ones are phenotype characteristics that determine responses to drug treatments, response to infections, or predisposition to specific inherited diseases. Most of the research in this field has been focused on the studies of mutation effects on the final gene products, peptides, and their alterations. Considerably less attention was given to the mutations that may affect regulatory mechanism(s) of gene expression, although these may also affect the phenotype characteristics. In this study we make a pilot analysis of mutations observed in the regulatory regions of 24,667 human RefSeq genes. Our study reveals that out of eight studied mutation types, "insertions" are the only one that in a statistically significant manner alters predicted transcription factor binding sites (TFBSs). We also find that 25 families of TFBSs have been altered by mutations in a statistically significant manner in the promoter regions we considered. Moreover, we find that the related transcription factors are, for example, prominent in processes related to intracellular signaling; cell fate; morphogenesis of organs and epithelium; development of urogenital system, epithelium, and tube; neuron fate commitment. Our study highlights the significance of studying mutations within the genes regulatory regions and opens way for further detailed investigations on this topic, particularly on the downstream affected pathways.

Pax6: a multi-level regulator of ocular development

Eye development has been a paradigm for the study of organogenesis, from the demonstration of lens induction through epithelial tissue morphogenesis, to neuronal specification and differentiation. The transcription factor Pax6 has been shown to play a key role in each of these processes. Pax6 is required for initiation of developmental pathways, patterning of epithelial tissues, activation of tissue-specific genes and interaction with other regulatory pathways. Herein we examine the data accumulated over the last few decades from extensive analyses of biochemical modules and genetic manipulation of the Pax6 gene. Specifically, we describe the regulation of Pax6's expression pattern, the protein's DNA-binding properties, and its specific roles and mechanisms of action at all stages of lens and retinal development. Pax6 functions at multiple levels to integrate extracellular information and execute cell-intrinsic differentiation programs that culminate in the specification and differentiation of a distinct ocular lineage.

iSyTE: integrated Systems Tool for Eye gene discovery

PURPOSE

To facilitate the identification of genes associated with cataract and other ocular defects, the authors developed and validated a computational tool termed iSyTE (integrated Systems Tool for Eye gene discovery; http://bioinformatics.udel.edu/Research/iSyTE). iSyTE uses a mouse embryonic lens gene expression data set as a bioinformatics filter to select candidate genes from human or mouse genomic regions implicated in disease and to prioritize them for further mutational and functional analyses.

METHODS

Microarray gene expression profiles were obtained for microdissected embryonic mouse lens at three key developmental time points in the transition from the embryonic day (E)10.5 stage of lens placode invagination to E12.5 lens primary fiber cell differentiation. Differentially regulated genes were identified by in silico comparison of lens gene expression profiles with those of whole embryo body (WB) lacking ocular tissue.

RESULTS

Gene set analysis demonstrated that this strategy effectively removes highly expressed but nonspecific housekeeping genes from lens tissue expression profiles, allowing identification of less highly expressed lens disease-associated genes. Among 24 previously mapped human genomic intervals containing genes associated with isolated congenital cataract, the mutant gene is ranked within the top two iSyTE-selected candidates in approximately 88% of cases. Finally, in situ hybridization confirmed lens expression of several novel iSyTE-identified genes.

CONCLUSIONS

iSyTE is a publicly available Web resource that can be used to prioritize candidate genes within mapped genomic intervals associated with congenital cataract for further investigation. Extension of this approach to other ocular tissue components will facilitate eye disease gene discovery.

Genetic and physical interaction of Meis2, Pax3 and Pax7 during dorsal midbrain development

Background

During early stages of brain development, secreted molecules, components of intracellular signaling pathways and transcriptional regulators act in positive and negative feed-back or feed-forward loops at the mid-hindbrain boundary. These genetic interactions are of central importance for the specification and subsequent development of the adjacent mid- and hindbrain. Much less, however, is known about the regulatory relationship and functional interaction of molecules that are expressed in the tectal anlage after tectal fate specification has taken place and tectal development has commenced.

Results

Here, we provide experimental evidence for reciprocal regulation and subsequent cooperation of the paired-type transcription factors Pax3, Pax7 and the TALE-homeodomain protein Meis2 in the tectal anlage. Using in ovo electroporation of the mesencephalic vesicle of chick embryos we show that (i) Pax3 and Pax7 mutually regulate each other's expression in the mesencephalic vesicle, (ii) Meis2 acts downstream of Pax3/7 and requires balanced expression levels of both proteins, and (iii) Meis2 physically interacts with Pax3 and Pax7. These results extend our previous observation that Meis2 cooperates with Otx2 in tectal development to include Pax3 and Pax7 as Meis2 interacting proteins in the tectal anlage.

Conclusion

The results described here suggest a model in which interdependent regulatory loops involving Pax3 and Pax7 in the dorsal mesencephalic vesicle modulate Meis2 expression. Physical interaction with Meis2 may then confer tectal specificity to a wide range of otherwise broadly expressed transcriptional regulators, including Otx2, Pax3 and Pax7.

Genetic and physical interaction of Meis2, Pax3 and Pax7 during dorsal midbrain development

BACKGROUND

During early stages of brain development, secreted molecules, components of intracellular signaling pathways and transcriptional regulators act in positive and negative feed-back or feed-forward loops at the mid-hindbrain boundary. These genetic interactions are of central importance for the specification and subsequent development of the adjacent mid- and hindbrain. Much less, however, is known about the regulatory relationship and functional interaction of molecules that are expressed in the tectal anlage after tectal fate specification has taken place and tectal development has commenced.

RESULTS

Here, we provide experimental evidence for reciprocal regulation and subsequent cooperation of the paired-type transcription factors Pax3, Pax7 and the TALE-homeodomain protein Meis2 in the tectal anlage. Using in ovo electroporation of the mesencephalic vesicle of chick embryos we show that (i) Pax3 and Pax7 mutually regulate each other's expression in the mesencephalic vesicle, (ii) Meis2 acts downstream of Pax3/7 and requires balanced expression levels of both proteins, and (iii) Meis2 physically interacts with Pax3 and Pax7. These results extend our previous observation that Meis2 cooperates with Otx2 in tectal development to include Pax3 and Pax7 as Meis2 interacting proteins in the tectal anlage.

CONCLUSION

The results described here suggest a model in which interdependent regulatory loops involving Pax3 and Pax7 in the dorsal mesencephalic vesicle modulate Meis2 expression. Physical interaction with Meis2 may then confer tectal specificity to a wide range of otherwise broadly expressed transcriptional regulators, including Otx2, Pax3 and Pax7.

Differential regulation of mesodermal gene expression by Drosophila cell type-specific Forkhead transcription factors

A common theme in developmental biology is the repeated use of the same gene in diverse spatial and temporal domains, a process that generally involves transcriptional regulation mediated by multiple separate enhancers, each with its own arrangement of transcription factor (TF)-binding sites and associated activities. Here, by contrast, we show that the expression of the Drosophila Nidogen (Ndg) gene at different embryonic stages and in four mesodermal cell types is governed by the binding of multiple cell-specific Forkhead (Fkh) TFs – including Biniou (Bin), Checkpoint suppressor homologue (CHES-1-like) and Jumeau (Jumu) – to three functionally distinguishable Fkh-binding sites in the same enhancer. Whereas Bin activates the Ndg enhancer in the late visceral musculature, CHES-1-like cooperates with Jumu to repress this enhancer in the heart. CHES-1-like also represses the Ndg enhancer in a subset of somatic myoblasts prior to their fusion to form multinucleated myotubes. Moreover, different combinations of Fkh sites, corresponding to two different sequence specificities, mediate the particular functions of each TF. A genome-wide scan for the occurrence of both classes of Fkh domain recognition sites in association with binding sites for known cardiac TFs showed an enrichment of combinations containing the two Fkh motifs in putative enhancers found within the noncoding regions of genes having heart expression. Collectively, our results establish that different cell-specific members of a TF family regulate the activity of a single enhancer in distinct spatiotemporal domains, and demonstrate how individual binding motifs for a TF class can differentially influence gene expression.

Comparison between aniridia with and without PAX6 mutations: clinical and molecular analysis in 14 Korean patients with aniridia

PURPOSE

To describe clinical and molecular characteristics of Korean patients with aniridia and to compare the clinical phenotype between those having an identifiable PAX6 mutation and those not.

DESIGN

Comparative case series.

PARTICIPANTS

A total of 14 Korean patients from 10 families with aniridia.

METHODS

Complete ophthalmologic examinations were performed for all patients. PAX6 analysis included direct sequencing of all coding regions and multiplex ligation-dependent probe amplification (MLPA) to detect large deletions when the sequencing was negative. If the PAX6 analysis failed to reveal any identifiable mutations, genomic copy number variation analysis via array comparative genomic hybridization (CGH) and candidate gene PITX3 and FOXE3 sequencing were then performed.

MAIN OUTCOME MEASURES

Severity of ocular abnormalities and genetic findings.

RESULTS

Sequencing of PAX6 exhibited 5 different heterozygous mutations in 8 patients from 5 families; 2 (p.Ser43Phe, IVS8-9C>G) were novel, and 3 (p.Arg208Trp, p.Arg317X, and p.X423L) have been previously reported. Among the remaining 6 patients in whom the PAX6 sequencing was negative, MLPA identified large deletions in 2 sporadic patients. However, the array CGH and candidate gene sequencing found no genomic or genetic abnormalities. The mutation detection rate was therefore 70%. Patients harboring an identifiable mutation in PAX6 had either a severe or a mild variant phenotype depending on the type of mutations. Likewise, among patients without an identifiable PAX6 mutation, their phenotypes varied widely from severe to very mild.

CONCLUSIONS

This study adds 2 novel PAX6 mutations to those previously reported, providing further evidence for genetic and phenotypic heterogeneity in aniridic ocular malformation. There was no difference in the clinical phenotype between patients with and without detectable mutations in the PAX6 gene. The wide variability of ocular phenotype regardless of the presence or absence of PAX6 mutations calls for a further appreciation of the complexity in the molecular diagnosis of aniridia and suggests that this ocular malformation may be better regarded as a group of heterogeneous disorders, rather than a single disease entity, associated with mutations in PAX6 and/or other genes located elsewhere in the human genome.

FINANCIAL DISCLOSURE(S)

The authors have no proprietary or commercial interest in any of the materials discussed in this article.

Transcription factors involved in lens development from the preplacodal ectoderm

Lens development is a stepwise process accompanied by the sequential activation of transcription factors. Transcription factor genes can be classified into three groups according to their functions: the first group comprises preplacodal genes, which are implicated in the formation of the preplacodal ectoderm that serves as a common primordium for cranial sensory tissues, including the lens. The second group comprises lens-specification genes, which establish the lens-field within the preplacodal ectoderm. The third group comprises lens-differentiation genes, which promote lens morphogenesis after the optic vesicle makes contact with the presumptive lens ectoderm. Analyses of the regulatory interactions between these genes have provided an overview of lens development, highlighting crucial roles for positive cross-regulation in fate specification and for feed-forward regulation in the execution of terminal differentiation. This overview also sheds light upon the mechanisms of how preplacodal gene activities lead to the activation of genes involved in lens-specification.

The functional role of the Meis/Prep-binding elements in Pax6 locus during pancreas and eye development

Pax6 is an essential transcription factor for lens, lacrimal gland and pancreas development. Previous transgenic analyses have identified several Pax6 regulatory elements, but their functional significance and binding factors remain largely unknown. In this study, we generated two genomic truncations to delete three elements that were previously shown to bind to the Meis/Prep family homeoproteins. One 3.1 kb deletion (Pax6(∆DP/∆DP)) removed two putative pancreatic enhancers and a previously identified ectodermal enhancer, while a 450 bp sub-deletion (Pax6(∆PE/∆PE)) eliminated only the promoter-proximal pancreatic enhancer. Immunohistochemistry and quantitative RT-PCR showed that the Pax6(∆PE/∆PE) pancreata had a significant decrease in Pax6, glucagon, and insulin expression, while no further reductions were observed in the Pax6(∆DP/∆DP) mice, indicating that only the 450 bp region is required for pancreatic development. In contrast, Pax6(∆DP/∆DP), but not Pax6(∆PE/∆PE) mice, developed stunted lacrimal gland and lens hypoplasia which was significantly more severe than that reported when only the ectodermal enhancer was deleted. This result suggested that the ectodermal enhancer must cooperate with its neighboring sequences to regulate the Pax6 ectodermal expression. Finally, we generated conditional knockouts of Prep1 in the lens and pancreas, but surprisingly, did not observe any developmental defects. Together, these results provide functional evidence for the independent and synergistic roles of the Pax6 upstream enhancers, and they suggest the potential redundancy of Meis/Prep protein in Pax6 regulation.

A model of spatially restricted transcription in opposing gradients of activators and repressors

Morphogens control patterns of transcription in development, often by establishing concentration gradients of a single transcriptional activator. However, many morphogens, including Hedgehog, create opposing activator and repressor gradients (OARGs). In contrast to single activator gradients, it is not well understood how OARGs control transcriptional patterns. We present a general thermodynamic model that explains how spatial patterns of gene expression are established within OARGs. The model predicts that differences in enhancer binding site affinities for morphogen-responsive transcription factors (TFs) produce discrete transcriptional boundaries, but only when either activators or repressors bind cooperatively. This model quantitatively predicts the boundaries of gene expression within OARGs. When trained on experimental data, our model accounts for the counterintuitive observation that increasing the affinity of binding sites in enhancers of Hedgehog target genes produces more restricted transcription within Hedgehog gradients in Drosophila. Because our model is general, it may explain the role of low-affinity binding sites in many contexts, including mammalian Hedgehog gradients.

Non-DNA-binding cofactors enhance DNA-binding specificity of a transcriptional regulatory complex

Recruitment of cofactors to specific DNA sites is integral for specificity in gene regulation. As a model system, we examined how targeting and transcriptional control of the sulfur metabolism genes in Saccharomyces cerevisiae is governed by recruitment of the transcriptional co-activator Met4. We developed genome-scale approaches to measure transcription factor (TF) DNA-binding affinities and cofactor recruitment to >1300 genomic binding site sequences. We report that genes responding to the TF Cbf1 and cofactor Met28 contain a novel 'recruitment motif' (RYAAT), adjacent to Cbf1 binding sites, which enhances the binding of a Met4-Met28-Cbf1 regulatory complex, and that abrogation of this motif significantly reduces gene induction under low-sulfur conditions. Furthermore, we show that correct recognition of this composite motif requires both non-DNA-binding cofactors Met4 and Met28. Finally, we demonstrate that the presence of an RYAAT motif next to a Cbf1 site, rather than Cbf1 binding affinity, specifies Cbf1-dependent sulfur metabolism genes. Our results highlight the need to examine TF/cofactor complexes, as novel specificity can result from cofactors that lack intrinsic DNA-binding specificity.

The orchestration of mammalian tissue morphogenesis through a series of coherent feed-forward loops

Tissue morphogenesis requires intricate temporal and spatial control of gene expression that is executed through specific gene regulatory networks (GRNs). GRNs are comprised from individual subcircuits of different levels of complexity. An important question is to elucidate the mutual relationship between those genes encoding DNA-binding factors that trigger the subcircuit with those that play major "later" roles during terminal differentiation via expression of specific genes that constitute the phenotype of individual tissues. The ocular lens is a classical model system to study tissue morphogenesis. Pax6 is essential for both lens placode formation and subsequent stages of lens morphogenesis, whereas c-Maf controls terminal differentiation of lens fibers, including regulation of crystallins, key lens structural proteins required for its transparency and refraction. Here, we show that Pax6 directly regulates c-Maf expression during lens development. A 1.3-kb c-Maf promoter with a 1.6-kb upstream enhancer (CR1) recapitulated the endogenous c-Maf expression pattern in lens and retinal pigmented epithelium. ChIP assays revealed binding of Pax6 and c-Maf to multiple regions of the c-Maf locus in lens chromatin. To predict functional Pax6-binding sites, nine novel variants of Pax6 DNA-binding motifs were identified and characterized. Two of these motifs predicted a pair of Pax6-binding sites in the CR1. Mutagenesis of these Pax6-binding sites inactivated transgenic expression in the lens but not in retinal pigmented epithelium. These data establish a novel regulatory role for Pax6 during lens development, link together the Pax6/c-Maf/crystallin regulatory network, and suggest a novel type of GRN subcircuit that controls a major part of embryonic lens development.

An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming

Alternative splicing (AS) is a key process underlying the expansion of proteomic diversity and the regulation of gene expression. Here, we identify an evolutionarily conserved embryonic stem cell (ESC)-specific AS event that changes the DNA-binding preference of the forkhead family transcription factor FOXP1. We show that the ESC-specific isoform of FOXP1 stimulates the expression of transcription factor genes required for pluripotency, including OCT4, NANOG, NR5A2, and GDF3, while concomitantly repressing genes required for ESC differentiation. This isoform also promotes the maintenance of ESC pluripotency and contributes to efficient reprogramming of somatic cells into induced pluripotent stem cells. These results reveal a pivotal role for an AS event in the regulation of pluripotency through the control of critical ESC-specific transcriptional programs.

Rapid evolutionary rewiring of a structurally constrained eye enhancer

BACKGROUND

Enhancers are genomic cis-regulatory sequences that integrate spatiotemporal signals to control gene expression. Enhancer activity depends on the combination of bound transcription factors as well as-in some cases-the arrangement and spacing of binding sites for these factors. Here, we examine evolutionary changes to the sequence and structure of sparkling, a Notch/EGFR/Runx-regulated enhancer that activates the dPax2 gene in cone cells of the developing Drosophila eye.

RESULTS

Despite functional and structural constraints on its sequence, sparkling has undergone major reorganization in its recent evolutionary history. Our data suggest that the relative strengths of the various regulatory inputs into sparkling change rapidly over evolutionary time, such that reduced input from some factors is compensated by increased input from different regulators. These gains and losses are at least partly responsible for the changes in enhancer structure that we observe. Furthermore, stereotypical spatial relationships between certain binding sites ("grammar elements") can be identified in all sparkling orthologs-although the sites themselves are often recently derived. We also find that low binding affinity for the Notch-regulated transcription factor Su(H), a conserved property of sparkling, is required to prevent ectopic responses to Notch in noncone cells.

CONCLUSIONS

Rapid DNA sequence turnover does not imply either the absence of critical cis-regulatory information or the absence of structural rules. Our findings demonstrate that even a severely constrained cis-regulatory sequence can be significantly rewired over a short evolutionary timescale.

Beyond the balance of activator and repressor

Cell differentiation during development is controlled by extracellular morphogens, which induce responding cells to differentiate into distinct cell fates based on the dose of morphogen they receive. Genes that specify the distinct cell fates are differentially responsive to morphogens, and the extracellular morphogen gradient is converted in responding cells to graded activity of transcription factors. In the case of Hedgehog, the gradient is converted to opposing gradients of transcriptional activator and repressor forms of the transcription factor Cubitus interruptus (Ci). It has been generally assumed that the balance between activator and repressor determines target gene responses within this gradient. However, new evidence shows that enhancers can respond selectively to the activator and repressor forms of Ci, and that this selectivity is determined by the affinity of Ci sites within the enhancers.

The cis-regulatory logic of Hedgehog gradient responses: key roles for gli binding affinity, competition, and cooperativity

Gradients of diffusible signaling proteins control precise spatial patterns of gene expression in the developing embryo. Here, we use quantitative expression measurements and thermodynamic modeling to uncover the cis-regulatory logic underlying spatially restricted gene expression in a Hedgehog (Hh) gradient in Drosophila. When Hh signaling is low, the Hh effector Gli, known as Cubitus interruptus (Ci) in Drosophila, acts as a transcriptional repressor; when Hh signaling is high, Gli acts as a transcriptional activator. Counterintuitively and in contrast to previous models of Gli-regulated gene expression, we found that low-affinity binding sites for Ci were required for proper spatial expression of the Hh target gene decapentaplegic (dpp) in regions of low Hh signal. Three low-affinity Ci sites enabled expression of dpp in response to low signal; increasing the affinity of these sites restricted dpp expression to regions of maximal signaling. A model incorporating cooperative repression by Ci correctly predicted the in vivo expression of a reporter gene controlled by a single Ci site. Our work clarifies how transcriptional activators and repressors, competing for common binding sites, can transmit positional information to the genome. It also provides an explanation for the widespread presence of conserved, nonconsensus Gli binding sites in Hh target genes.

In vivo Atoh1 targetome reveals how a proneural transcription factor regulates cerebellar development

The proneural, basic helix-loop-helix transcription factor Atoh1 governs the development of numerous key neuronal subtypes, such as cerebellar granule and brainstem neurons, inner ear hair cells, and several neurons of the proprioceptive system, as well as diverse nonneuronal cell types, such as Merkel cells and intestinal secretory lineages. However, the mere handful of targets that have been identified barely begin to account for Atoh1's astonishing range of functions, which also encompasses seemingly paradoxical activities, such as promoting cell proliferation and medulloblastoma formation in the cerebellum and inducing cell cycle exit and suppressing tumorigenesis in the intestine. We used a multipronged approach to create a comprehensive, unbiased list of over 600 direct Atoh1 target genes in the postnatal cerebellum. We found that Atoh1 binds to a 10 nucleotide motif (AtEAM) to directly regulate genes involved in migration, cell adhesion, metabolism, and other previously unsuspected functions. This study expands current thinking about the transcriptional activities driving neuronal differentiation and provides a framework for further neurodevelopmental studies.