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

Optimization of the Irf8 +32-kb enhancer disrupts dendritic cell lineage segregation

Autoactivation of lineage-determining transcription factors mediates bistable expression, generating distinct cell phenotypes essential for complex body plans. Classical type 1 dendritic cell (cDC1) and type 2 dendritic cell (cDC2) subsets provide nonredundant functions for defense against distinct immune challenges. Interferon regulatory factor 8 (IRF8), the cDC1 lineage-determining transcription factor, undergoes autoactivation in cDC1 progenitors to establish cDC1 identity, yet its expression is downregulated during cDC2 differentiation by an unknown mechanism. This study reveals that the Irf8 +32-kb enhancer, responsible for IRF8 autoactivation, is naturally suboptimized with low-affinity IRF8 binding sites. Introducing multiple high-affinity IRF8 sites into the Irf8 +32-kb enhancer causes a gain-of-function effect, leading to erroneous IRF8 autoactivation in specified cDC2 progenitors, redirecting them toward cDC1 and a novel hybrid DC subset with mixed-lineage phenotypes. Further, this also causes a loss-of-function effect, reducing Irf8 expression in cDC1s. These developmental alterations critically impair both cDC1-dependent and cDC2-dependent arms of immunity. Collectively, our findings underscore the significance of enhancer suboptimization in the developmental segregation of cDCs required for normal immune function.

Anatomy of a superenhancer

Interferon regulatory factor-8 (IRF8) is the lineage determining transcription factor for the type one classical dendritic cell (cDC1) subset, a terminal selector for plasmacytoid dendritic cells and important for the function of monocytes. Studies of Irf8 gene regulation have identified several enhancers controlling its activity during development of progenitors in the bone marrow that precisely regulate expression at distinct developmental stages. Each enhancer responds to distinct transcription factors that are expressed at each stage. IRF8 is first expressed in early progenitors that form the monocyte dendritic cell progenitor (MDP) in response to induction of the transcription factor CCAAT/enhancer-binding protein alpha (C/EBPα) acting at the Irf8 +56 kb enhancer. IRF8 levels increase further as the MDP transits into the common dendritic cell progenitor (CDP) in response to E protein activity at the Irf8 +41 kb enhancer. Upon Nfil3-induction in CDPs leading to specification of the cDC1 progenitor, abrupt induction of BATF3 forms the JUN/BATF3/IRF8 heterotrimer that activates the Irf8 +32 kb enhancer that sustains Irf8 autoactivation throughout the cDC1 lifetime. Deletions of each of these enhancers has revealed their stage dependent activation. Surprisingly, studies of compound heterozygotes for each combination of enhancer deletions revealed that activation of each subsequent enhancer requires the successful activation of the previous enhancer in strictly cis-dependent mechanism. Successful progression of enhancer activation is finely tuned to alter the functional accessibility of subsequent enhancers to factors active in the next stage of development. The molecular basis for these phenomenon is still obscure but could have implications for genomic regulation in a broader developmental context.

Lens placode modulates extracellular matrix formation during early eye development

The role extracellular matrix (ECM) in multiple events of morphogenesis has been well described, little is known about its specific role in early eye development. One of the first morphogenic events in lens development is placodal thickening, which converts the presumptive lens ectoderm from cuboidal to pseudostratified epithelium. This process occurs in the anterior pre-placodal ectoderm when the optic vesicle approaches the cephalic ectoderm and is regulated by transcription factor Pax6 and secreted BMP4. Since cells and ECM have a dynamic relationship of interdependence and modulation, we hypothesized that the ECM evolves with cell shape changes during lens placode formation. This study investigates changes in optic ECM including both protein distribution deposition, extracellular gelatinase activity and gene expression patterns during early optic development using chicken and mouse models. In particular, the expression of Timp2, a metalloprotease inhibitor, corresponds with a decrease in gelatinase activity within the optic ECM. Furthermore, we demonstrate that optic ECM remodeling depends on BMP signaling in the placode. Together, our findings suggest that the lens placode plays an active role in remodeling the optic ECM during early eye development.

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Overlapping binding sites underlie TF genomic occupancy

Sequence-specific DNA binding by transcription factors (TFs) is a crucial step in gene regulation. However, current high-throughput in vitro approaches cannot reliably detect lower affinity TF-DNA interactions, which play key roles in gene regulation. Here, we developed PADIT-seq ( p rotein a ffinity to D NA by in vitro transcription and RNA seq uencing) to assay TF binding preferences to all 10-bp DNA sequences at far greater sensitivity than prior approaches. The expanded catalogs of low affinity DNA binding sites for the human TFs HOXD13 and EGR1 revealed that nucleotides flanking high affinity DNA binding sites create overlapping lower affinity sites that together modulate TF genomic occupancy in vivo . Formation of such extended recognition sequences stems from an inherent property of TF binding sites to interweave each other and expands the genomic sequence space for identifying noncoding variants that directly alter TF binding.

One-Sentence Summary

Overlapping DNA binding sites underlie TF genomic occupancy through their inherent propensity to interweave each other.

Lens Placode Modulates Extracellular Matrix Formation During Early Eye Development

The role extracellular matrix (ECM) in multiple events of morphogenesis has been well described, little is known about its specific role in early eye development. One of the first morphogenic events in lens development is placodal thickening, which converts the presumptive lens ectoderm from cuboidal to pseudostratified epithelium. This process occurs in the anterior pre-placodal ectoderm when the optic vesicle approaches the cephalic ectoderm. Since cells and ECM have a dynamic relationship of interdependence and modulation, we hypothesized that the ECM evolves with cell shape changes during lens placode formation. This study investigates changes in optic ECM including both protein distribution deposition, extracellular gelatinase activity and gene expression patterns during early optic development using chicken and mouse models. In particular, the expression of Timp2 , a metalloprotease inhibitor, corresponds with a decrease in gelatinase activity within the optic ECM. Furthermore, we demonstrate that optic ECM remodeling depends on BMP signaling in the placode. Together, our findings suggest that the lens placode plays an active role in remodeling the optic ECM during early eye development.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia

BACKGROUND

Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus.

METHODS

Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region.

RESULTS

Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements.

CONCLUSION

Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

Japanese medaka Olpax6.1 mutant as a potential model for spondylo-ocular syndrome

pax6 is a canonic master gene for eye formation. Knockout of pax6 affects the development of craniofacial skeleton and eye in mice. Whether pax6 affects the development of spinal bone has not been reported yet. In the present study, we used CRISPR/Cas9 system to generate Olpax6.1 mutant in Japanese medaka. Phenotype analysis showed that ocular mutation caused by the Olpax6.1 mutation occurred in the homozygous mutant. The phenotype of heterozygotes is not significantly different from that of wild-type. In addition, knockout Olpax6.1 resulted in severe curvature of the spine in the homozygous F2 generation. Comparative transcriptome analysis and qRT-PCR revealed that the defective Olpax6.1 protein caused a decrease in the expression level of sp7, col10a1a, and bglap, while the expression level of xylt2 did not change significantly. The functional enrichment of differentially expressed genes (DEGs) using the Kyoto Encyclopedia of Genes and Genomes database showed that the DEGs between Olpax6.1 mutation and wild-type were enriched in p53 signaling pathway, extracellular matrix (ECM) -receptor interaction, et al. Our results indicated that the defective Olpax6.1 protein results in the reduction of sp7 expression level and the activation of p53 signaling pathway, which leads to a decrease in the expression of genes encoding ECM protein, such as collagen protein family and bone gamma-carboxyglutamate protein, which further inhibits bone development. Based on the phenotype and molecular mechanism of ocular mutation and spinal curvature induced by Olpax6.1 knockout, we believe that the Olpax6.1-/- mutant could be a potential model for the study of spondylo-ocular syndrome.

High-Throughput Transcriptomics of Celf1 Conditional Knockout Lens Identifies Downstream Networks Linked to Cataract Pathology

Defects in the development of the ocular lens can cause congenital cataracts. To understand the various etiologies of congenital cataracts, it is important to characterize the genes linked to this developmental defect and to define their downstream pathways that are relevant to lens biology and pathology. Deficiency or alteration of several RNA-binding proteins, including the conserved RBP Celf1 (CUGBP Elav-like family member 1), has been described to cause lens defects and early onset cataracts in animal models and/or humans. Celf1 is involved in various aspects of post-transcriptional gene expression control, including regulation of mRNA stability/decay, alternative splicing and translation. Celf1 germline knockout mice and lens conditional knockout (Celf1cKO) mice develop fully penetrant cataracts in early postnatal stages. To define the genome-level changes in RNA transcripts that result from Celf1 deficiency, we performed high-throughput RNA-sequencing of Celf1cKO mouse lenses at postnatal day (P) 0. Celf1cKO lenses exhibit 987 differentially expressed genes (DEGs) at cut-offs of >1.0 log2 counts per million (CPM), ≥±0.58 log2 fold-change and <0.05 false discovery rate (FDR). Of these, 327 RNAs were reduced while 660 were elevated in Celf1cKO lenses. The DEGs were subjected to various downstream analyses including iSyTE lens enriched-expression, presence in Cat-map, and gene ontology (GO) and representation of regulatory pathways. Further, a comparative analysis was done with previously generated microarray datasets on Celf1cKO lenses P0 and P6. Together, these analyses validated and prioritized several key genes mis-expressed in Celf1cKO lenses that are relevant to lens biology, including known cataract-linked genes (e.g., Cryab, Cryba2, Cryba4, Crybb1, Crybb2, Cryga, Crygb, Crygc, Crygd, Cryge, Crygf, Dnase2b, Bfsp1, Gja3, Pxdn, Sparc, Tdrd7, etc.) as well as novel candidates (e.g., Ell2 and Prdm16). Together, these data have defined the alterations in lens transcriptome caused by Celf1 deficiency, in turn uncovering downstream genes and pathways (e.g., structural constituents of eye lenses, lens fiber cell differentiation, etc.) associated with lens development and early-onset cataracts.

Knockout of signal peptide peptidase in the eye reduces HSV-1 replication and eye disease in ocularly infected mice

We previously reported that knocking out signal peptide peptidase (SPP), a glycoprotein K (gK) binding partner, in mouse peripheral sensory neurons reduced latency-reactivation in infected mice without affecting primary virus replication or eye disease. Since virus replication in the eye plays an essential role in eye disease, we generated a conditional knockout mouse lacking SPP expression in the eye by crossing Pax6 (paired box 6)-Cre mice that have intact Pax6 expression with SPPflox/flox mice. Significantly less SPP protein expression was detected in the eyes of Pax6-SPP-/- mice than in WT control mice. HSV-1 replication in the eyes of Pax6-SPP-/- mice was significantly lower than in WT control mice. Levels of gB, gK, and ICP0 transcripts in corneas, but not trigeminal ganglia (TG), of Pax6-SPP-/- infected mice were also significantly lower than in WT mice. Corneal scarring and angiogenesis were significantly lower in Pax6-SPP-/- mice than in WT control mice, while corneal sensitivity was significantly higher in Pax6-SPP-/- mice compared with WT control mice. During acute viral infection, absence of SPP in the eye did not affect CD4 expression but did affect CD8α and IFNγ expression in the eye. However, in the absence of SPP, latency-reactivation was similar in Pax6-SPP-/- and WT control groups. Overall, our results showed that deleting SPP expression in the eyes reduced primary virus replication in the eyes, reduced CD8α and IFNγ mRNA expression, reduced eye disease and reduced angiogenesis but did not alter corneal sensitivity or latency reactivation to HSV-1 infection. Thus, blocking gK binding to SPP in the eye may have therapeutic potential by reducing both virus replication in the eye and eye disease associated with virus replication.

"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes

Regulatory elements activate promoters by recruiting transcription factors (TFs) to specific motifs. Notably, TF-DNA interactions often depend on cooperativity with colocalized partners, suggesting an underlying cis-regulatory syntax. To explore TF cooperativity in mammals, we analyze ∼500 mouse and human primary cells by combining an atlas of TF motifs, footprints, ChIP-seq, transcriptomes, and accessibility. We uncover two TF groups that colocalize with most expressed factors, forming stripes in hierarchical clustering maps. The first group includes lineage-determining factors that occupy DNA elements broadly, consistent with their key role in tissue-specific transcription. The second one, dubbed universal stripe factors (USFs), comprises ∼30 SP, KLF, EGR, and ZBTB family members that recognize overlapping GC-rich sequences in all tissues analyzed. Knockouts and single-molecule tracking reveal that USFs impart accessibility to colocalized partners and increase their residence time. Mammalian cells have thus evolved a TF superfamily with overlapping DNA binding that facilitate chromatin accessibility.

Concentration-dependent transcriptional switching through a collective action of cis-elements

Gene expression specificity of homeobox transcription factors has remained paradoxical. WUSCHEL activates and represses CLAVATA3 transcription at lower and higher concentrations, respectively. We use computational modeling and experimental analysis to investigate the properties of the cis-regulatory module. We find that intrinsically each cis-element can only activate CLAVATA3 at a higher WUSCHEL concentration. However, together, they repress CLAVATA3 at higher WUSCHEL and activate only at lower WUSCHEL, showing that the concentration-dependent interactions among cis-elements regulate both activation and repression. Biochemical experiments show that two adjacent functional cis-elements bind WUSCHEL with higher affinity and dimerize at relatively lower levels. Moreover, increasing the distance between cis-elements prolongs WUSCHEL monomer binding window, resulting in higher CLAVATA3 activation. Our work showing a constellation of optimally spaced cis-elements of defined affinities determining activation and repression thresholds in regulating CLAVATA3 transcription provides a previously unknown mechanism of cofactor-independent regulation of transcription factor binding in mediating gene expression specificity.

Evolution of transcription factor binding through sequence variations and turnover of binding sites

Variations in noncoding regulatory sequences play a central role in evolution. Interpreting such variations, however, remains difficult even in the context of defined attributes such as transcription factor (TF) binding sites. Here, we systematically link variations in cis-regulatory sequences to TF binding by profiling the allele-specific binding of 27 TFs expressed in a yeast hybrid, in which two related genomes are present within the same nucleus. TFs localize preferentially to sites containing their known consensus motifs but occupy only a small fraction of the motif-containing sites available within the genomes. Differential binding of TFs to the orthologous alleles was well explained by variations that alter motif sequence, whereas differences in chromatin accessibility between alleles were of little apparent effect. Motif variations that abolished binding when present in only one allele were still bound when present in both alleles, suggesting evolutionary compensation, with a potential role for sequence conservation at the motif's vicinity. At the level of the full promoter, we identify cases of binding-site turnover, in which binding sites are reciprocally gained and lost, yet most interspecific differences remained uncompensated. Our results show the flexibility of TFs to bind imprecise motifs and the fast evolution of TF binding sites between related species.

Androgen receptor positively regulates gonadotropin-releasing hormone receptor in pituitary gonadotropes

Within pituitary gonadotropes, the gonadotropin-releasing hormone receptor (GnRHR) receives hypothalamic input from GnRH neurons that is critical for reproduction. Previous studies have suggested that androgens may regulate GnRHR, although the mechanisms remain unknown. In this study, we demonstrated that androgens positively regulate Gnrhr mRNA in mice. We then investigated the effects of androgens and androgen receptor (AR) on Gnrhr promoter activity in immortalized mouse LβT2 cells, which represent mature gonadotropes. We found that AR positively regulates the Gnrhr proximal promoter, and that this effect requires a hormone response element (HRE) half site at -159/-153 relative to the transcription start site. We also identified nonconsensus, full-length HREs at -499/-484 and -159/-144, which are both positively regulated by androgens on a heterologous promoter. Furthermore, AR associates with the Gnrhr promoter in ChIP. Altogether, we report that GnRHR is positively regulated by androgens through recruitment of AR to the Gnrhr proximal promoter.

Temporal Control of the TGF-β Signaling Network by Mouse ESC MicroRNA Targets of Different Affinities

Although microRNAs (miRNAs) function in the control of embryonic stem cell (ESC) pluripotency, a systems-level understanding is still being developed. Through the analysis of progressive Argonaute (Ago)-miRNA depletion and rescue, including stable Ago knockout mouse ESCs, we uncover transforming growth factor beta (TGF-β) pathway activation as a direct and early response to ESC miRNA reduction. Mechanistically, we link the derepression of weaker miRNA targets, including TGF-β receptor 1 (Tgfbr1), to the sensitive TGF-β pathway activation. In contrast, stronger miRNA targets impart a more robust repression, which dampens concurrent transcriptional activation. We verify such dampened induction for TGF-β antagonist Lefty. We find that TGF-β pathway activation contributes to the G1 cell-cycle accumulation of miRNA-deficient ESCs. We propose that miRNA target affinity is a determinant of the temporal response to miRNA changes, which enables the coordination of gene network responses.

Kelly et al. report the transcriptional and post-transcriptional dynamics that occur with loss of Argonaute proteins in embryonic stem cells. They find that Argonaute proteins are not required for ESC viability, function to control the transforming growth factor beta (TGF-β) pathway, and mediate temporal responses during changes in miRNA levels.

Robust and efficient gene regulation through localized nuclear microenvironments

Developmental enhancers drive gene expression in specific cell types during animal development. They integrate signals from many different sources mediated through the binding of transcription factors, producing specific responses in gene expression. Transcription factors often bind low-affinity sequences for only short durations. How brief, low-affinity interactions drive efficient transcription and robust gene expression is a central question in developmental biology. Localized high concentrations of transcription factors have been suggested as a possible mechanism by which to use these enhancer sites effectively. Here, we discuss the evidence for such transcriptional microenvironments, mechanisms for their formation and the biological consequences of such sub-nuclear compartmentalization for developmental decisions and evolution.

Concentration-dependent splicing is enabled by Rbfox motifs of intermediate affinity

The Rbfox family of splicing factors regulate alternative splicing during animal development and in disease, impacting thousands of exons in the maturing brain, heart, and muscle. Rbfox proteins have long been known to bind to the RNA sequence GCAUG with high affinity, but just half of Rbfox binding sites contain a GCAUG motif in vivo. We incubated recombinant RBFOX2 with over 60,000 mouse and human transcriptomic sequences to reveal substantial binding to several moderate-affinity, non-GCAYG sites at a physiologically relevant range of RBFOX concentrations. We find that many of these “secondary motifs” bind Rbfox robustly in cells and that several together can exert regulation comparable to GCAUG in a trichromatic splicing reporter assay. Furthermore, secondary motifs regulate RNA splicing in neuronal development and in neuronal subtypes where cellular Rbfox concentrations are highest, enabling a second wave of splicing changes as Rbfox levels increase.

Functional genomic approaches to elucidate the role of enhancers during development

Successful development depends on the precise tissue-specific regulation of genes by enhancers, genetic elements that act as switches to control when and where genes are expressed. Because enhancers are critical for development, and the majority of disease-associated mutations reside within enhancers, it is essential to understand which sequences within enhancers are important for function. Advances in sequencing technology have enabled the rapid generation of genomic data that predict putative active enhancers, but functionally validating these sequences at scale remains a fundamental challenge. Herein, we discuss the power of genome-wide strategies used to identify candidate enhancers, and also highlight limitations and misconceptions that have arisen from these data. We discuss the use of massively parallel reporter assays to test enhancers for function at scale. We also review recent advances in our ability to study gene regulation during development, including CRISPR-based tools to manipulate genomes and single-cell transcriptomics to finely map gene expression. Finally, we look ahead to a synthesis of complementary genomic approaches that will advance our understanding of enhancer function during development. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Developmental Biology > Developmental Processes in Health and Disease Laboratory Methods and Technologies > Genetic/Genomic Methods.

The Spectrum of PAX6 Mutations and Genotype-Phenotype Correlations in the Eye

The transcription factor PAX6 is essential in ocular development in vertebrates, being considered the master regulator of the eye. During eye development, it is essential for the correct patterning and formation of the multi-layered optic cup and it is involved in the developing lens and corneal epithelium. In adulthood, it is mostly expressed in cornea, iris, and lens. PAX6 is a dosage-sensitive gene and it is highly regulated by several elements located upstream, downstream, and within the gene. There are more than 500 different mutations described to affect PAX6 and its regulatory regions, the majority of which lead to PAX6 haploinsufficiency, causing several ocular and systemic abnormalities. Aniridia is an autosomal dominant disorder that is marked by the complete or partial absence of the iris, foveal hypoplasia, and nystagmus, and is caused by heterozygous PAX6 mutations. Other ocular abnormalities have also been associated with PAX6 changes, and genotype-phenotype correlations are emerging. This review will cover recent advancements in PAX6 regulation, particularly the role of several enhancers that are known to regulate PAX6 during eye development and disease. We will also present an updated overview of the mutation spectrum, where an increasing number of mutations in the non-coding regions have been reported. Novel genotype-phenotype correlations will also be discussed.

Selective deployment of transcription factor paralogs with submaximal strength facilitates gene regulation in the immune system

In multicellular organisms, duplicated genes can diverge through tissue-specific gene expression patterns, as exemplified by highly regulated expression of Runx transcription factor paralogs with apparent functional redundancy. Here we asked what cell type-specific biologies might be supported by the selective expression of Runx paralogs during Langerhans cell and inducible regulatory T cell differentiation. We uncovered functional non-equivalence between Runx paralogs. Selective expression of native paralogs allowed integration of transcription factor activity with extrinsic signals, while non-native paralogs enforced differentiation even in the absence of exogenous inducers. DNA-binding affinity was controlled by divergent amino acids within the otherwise highly conserved RUNT domain, and evolutionary reconstruction suggested convergence of RUNT domain residues towards sub-maximal strength. Hence, the selective expression of gene duplicates in specialized cell types can synergize with the acquisition of functional differences to enable appropriate gene expression, lineage choice and differentiation in the mammalian immune system.

Altering the Temporal Regulation of One Transcription Factor Drives Evolutionary Trade-Offs between Head Sensory Organs

Size trade-offs of visual versus olfactory organs is a pervasive feature of animal evolution. This could result from genetic or functional constraints. We demonstrate that head sensory organ size trade-offs in Drosophila are genetically encoded and arise through differential subdivision of the head primordium into visual versus non-visual fields. We discover that changes in the temporal regulation of the highly conserved eyeless/Pax6 gene expression during development is a conserved mechanism for sensory trade-offs within and between Drosophila species. We identify a natural single nucleotide polymorphism in the cis-regulatory region of eyeless in a binding site of its repressor Cut that is sufficient to alter its temporal regulation and eye size. Because eyeless/Pax6 is a conserved regulator of head sensory placode subdivision, we propose that its temporal regulation is key to define the relative size of head sensory organs.

Multi-enhancer transcriptional hubs confer phenotypic robustness

We previously showed in Drosophila melanogaster embryos that low-affinity Ultrabithorax (Ubx)-responsive shavenbaby (svb) enhancers drive expression using localized transcriptional environments and that active svb enhancers on different chromosomes tended to colocalize (Tsai et al., 2017). Here, we test the hypothesis that these multi-enhancer 'hubs' improve phenotypic resilience to stress by buffering against decreases in transcription factor concentrations and transcriptional output. Deleting a redundant enhancer from the svb locus led to reduced trichome numbers in embryos raised at elevated temperatures. Using high-resolution fluorescence microscopy, we observed lower Ubx concentration and transcriptional output in this deletion allele. Transcription sites of the full svb cis-regulatory region inserted into a different chromosome colocalized with the svb locus, increasing Ubx concentration, the transcriptional output of svb, and partially rescuing the phenotype. Thus, multiple enhancers could reinforce a local transcriptional hub to buffer against environmental stresses and genetic perturbations, providing a mechanism for phenotypical robustness.

Profiling of chromatin accessibility and identification of general cis-regulatory mechanisms that control two ocular lens differentiation pathways

BACKGROUND

Promoters and enhancers are cis-regulatory DNA sequences that control specificity and quantity of transcription. Both are rich on clusters of cis-acting sites that interact with sequence-specific DNA-binding transcription factors (TFs). At the level of chromatin, these regions display increased nuclease sensitivity, reduced nucleosome density, including nucleosome-free regions, and specific combinations of posttranslational modifications of core histone proteins. Together, "open" and "closed" chromatins represent transcriptionally active and repressed states of individual genes, respectively. Cellular differentiation is marked by changes in local chromatin structure. Lens morphogenesis, regulated by TF Pax6, includes differentiation of epithelial precursor cells into lens fibers in parallel with differentiation of epithelial precursors into the mature lens epithelium.

RESULTS

Using ATAC-seq, we investigated dynamics of chromatin changes during mouse lens fibers and epithelium differentiation. Tissue-specific features of these processes are demonstrated via comparative studies of embryonic stem cells, forebrain, and liver chromatins. Unbiased analysis reveals cis-regulatory logic of lens differentiation through known (e.g., AP-1, Ets, Hsf4, Maf, and Pax6 sites) and novel (e.g., CTCF, Tead, and NF1) motifs. Twenty-six DNA-binding TFs, recognizing these cis-motifs, are markedly up-regulated in differentiating lens fibers. As specific examples, our ATAC-seq data uncovered both the regulatory regions and TF binding motifs in Foxe3, Prox1, and Mip loci that are consistent with previous, though incomplete, experimental data. A cross-examination of Pax6 binding with ATAC-seq data demonstrated that Pax6 bound to both open (H3K27ac and P300-enriched) and closed chromatin domains in lens and forebrain.

CONCLUSIONS

Our study has generated the first lens chromatin accessibility maps that support a general model of stage-specific chromatin changes associated with transcriptional activities of batteries of genes required for lens fiber cell formation. Analysis of active (or open) promoters and enhancers reveals important cis-DNA motifs that establish the molecular foundation for temporally and spatially regulated gene expression in lens. Together, our data and models open new avenues for the field to conduct mechanistic studies of transcriptional control regions, reconstruction of gene regulatory networks that govern lens morphogenesis, and identification of cataract-causing mutations in noncoding sequences.

Diversification of transcription factor-DNA interactions and the evolution of gene regulatory networks

Sequence-specific transcription factors (TFs) bind short DNA sequences in the genome to regulate the expression of target genes. In the last decade, numerous technical advances have enabled the determination of the DNA-binding specificities of many of these factors. Large-scale screens of many TFs enabled the creation of databases of TF DNA-binding specificities, typically represented as position weight matrices (PWMs). Although great progress has been made in determining and predicting binding specificities systematically, there are still many surprises to be found when studying a particular TF's interactions with DNA in detail. Paralogous TFs' binding specificities can differ in subtle ways, in a manner that is not immediately apparent from looking at their PWMs. These differences affect gene regulatory outputs and enable TFs to rewire transcriptional networks over evolutionary time. This review discusses recent observations made in the study of TF-DNA interactions that highlight the importance of continued in-depth analysis of TF-DNA interactions and their inherent complexity. This article is categorized under: Biological Mechanisms > Regulatory Biology.

Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

Sea stars and sea urchins are model systems for interrogating the types of deep evolutionary changes that have restructured developmental gene regulatory networks (GRNs). Although cis-regulatory DNA evolution is likely the predominant mechanism of change, it was recently shown that Tbrain, a Tbox transcription factor protein, has evolved a changed preference for a low-affinity, secondary binding motif. The primary, high-affinity motif is conserved. To date, however, no genome-wide comparisons have been performed to provide an unbiased assessment of the evolution of GRNs between these taxa, and no study has attempted to determine the interplay between transcription factor binding motif evolution and GRN topology. The study here measures genome-wide binding of Tbrain orthologs by using ChIP-sequencing and associates these orthologs with putative target genes to assess global function. Targets of both factors are enriched for other regulatory genes, although nonoverlapping sets of functional enrichments in the two datasets suggest a much diverged function. The number of low-affinity binding motifs is significantly depressed in sea urchins compared with sea star, but both motif types are associated with genes from a range of functional categories. Only a small fraction (∼10%) of genes are predicted to be orthologous targets. Collectively, these data indicate that Tbr has evolved significantly different developmental roles in these echinoderms and that the targets and the binding motifs in associated cis-regulatory sequences are dispersed throughout the hierarchy of the GRN, rather than being biased toward terminal process or discrete functional blocks, which suggests extensive evolutionary tinkering.

A feed-forward relay integrates the regulatory activities of Bicoid and Orthodenticle via sequential binding to suboptimal sites

Datta et al. define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bicoid (Bcd) and Orthodenticle (Otd). The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins.

The K50 (lysine at amino acid position 50) homeodomain (HD) protein Orthodenticle (Otd) is critical for anterior patterning and brain and eye development in most metazoans. In Drosophila melanogaster, another K50HD protein, Bicoid (Bcd), has evolved to replace Otd's ancestral function in embryo patterning. Bcd is distributed as a long-range maternal gradient and activates transcription of a large number of target genes, including otd. Otd and Bcd bind similar DNA sequences in vitro, but how their transcriptional activities are integrated to pattern anterior regions of the embryo is unknown. Here we define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bcd and Otd. Class 1 enhancers are initially activated by Bcd, and activation is transferred to Otd via a feed-forward relay (FFR) that involves sequential binding of the two proteins to the same DNA motif. Class 2 enhancers are activated by Bcd and maintained by an Otd-independent mechanism. Class 3 enhancers are never bound by Bcd, but Otd binds and activates them in a second wave of zygotic transcription. The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins. Our results define specific patterning roles for Bcd and Otd and provide mechanisms for coordinating the precise timing of gene expression patterns during embryonic development.

Degenerate Pax2 and Senseless binding motifs improve detection of low-affinity sites required for enhancer specificity

Cells use thousands of regulatory sequences to recruit transcription factors (TFs) and produce specific transcriptional outcomes. Since TFs bind degenerate DNA sequences, discriminating functional TF binding sites (TFBSs) from background sequences represents a significant challenge. Here, we show that a Drosophila regulatory element that activates Epidermal Growth Factor signaling requires overlapping, low-affinity TFBSs for competing TFs (Pax2 and Senseless) to ensure cell- and segment-specific activity. Testing available TF binding models for Pax2 and Senseless, however, revealed variable accuracy in predicting such low-affinity TFBSs. To better define parameters that increase accuracy, we developed a method that systematically selects subsets of TFBSs based on predicted affinity to generate hundreds of position-weight matrices (PWMs). Counterintuitively, we found that degenerate PWMs produced from datasets depleted of high-affinity sequences were more accurate in identifying both low- and high-affinity TFBSs for the Pax2 and Senseless TFs. Taken together, these findings reveal how TFBS arrangement can be constrained by competition rather than cooperativity and that degenerate models of TF binding preferences can improve identification of biologically relevant low affinity TFBSs.

While all cells in an organism share a common genome, each cell type must express the appropriate combination of genes needed for its specific function. Cells activate and repress different parts of the genome using transcription factor proteins that bind regulatory regions known as enhancers. We currently have an incomplete view of how enhancers recruit transcription factors to yield accurate gene activation and repression. This problem is complicated by the fact that most animals contain over a thousand different transcription factors, and each can generally bind multiple DNA sequences. Thus, it is difficult to predict which transcription factors interact with which enhancers. To gain insights into this process, we focused on determining how an enhancer that activates a gene needed to make liver-like cells is regulated in a precise manner in the fruit-fly embryo. We demonstrate that the specific activity of this enhancer depends on weak and overlapping transcription factor binding sites. Furthermore, we demonstrate that computational models that include weak transcription factor interactions yield better predictive accuracy. These results shed light on how DNA sequences determine enhancer activity and the types of strategies that are most useful for predicting transcription factor binding sites in the genome.

The RNA-binding protein Celf1 post-transcriptionally regulates p27Kip1 and Dnase2b to control fiber cell nuclear degradation in lens development

Opacification of the ocular lens, termed cataract, is a common cause of blindness. To become transparent, lens fiber cells undergo degradation of their organelles, including their nuclei, presenting a fundamental question: does signaling/transcription sufficiently explain differentiation of cells progressing toward compromised transcriptional potential? We report that a conserved RNA-binding protein Celf1 post-transcriptionally controls key genes to regulate lens fiber cell differentiation. Celf1-targeted knockout mice and celf1-knockdown zebrafish and Xenopus morphants have severe eye defects/cataract. Celf1 spatiotemporally down-regulates the cyclin-dependent kinase (Cdk) inhibitor p27Kip1 by interacting with its 5' UTR and mediating translation inhibition. Celf1 deficiency causes ectopic up-regulation of p21Cip1. Further, Celf1 directly binds to the mRNA of the nuclease Dnase2b to maintain its high levels. Together these events are necessary for Cdk1-mediated lamin A/C phosphorylation to initiate nuclear envelope breakdown and DNA degradation in fiber cells. Moreover, Celf1 controls alternative splicing of the membrane-organization factor beta-spectrin and regulates F-actin-crosslinking factor Actn2 mRNA levels, thereby controlling fiber cell morphology. Thus, we illustrate new Celf1-regulated molecular mechanisms in lens development, suggesting that post-transcriptional regulatory RNA-binding proteins have evolved conserved functions to control vertebrate oculogenesis.

Natural variation in stochastic photoreceptor specification and color preference in Drosophila

Each individual perceives the world in a unique way, but little is known about the genetic basis of variation in sensory perception. In the fly eye, the random mosaic of color-detecting R7 photoreceptor subtypes is determined by stochastic on/off expression of the transcription factor Spineless (Ss). In a genome-wide association study, we identified a naturally occurring insertion in a regulatory DNA element in ss that lowers the ratio of Ss^ON to Ss^OFF cells. This change in photoreceptor fates shifts the innate color preference of flies from green to blue. The genetic variant increases the binding affinity for Klumpfuss (Klu), a zinc finger transcriptional repressor that regulates ss expression. Klu is expressed at intermediate levels to determine the normal ratio of Ss^ON to Ss^OFF cells. Thus, binding site affinity and transcription factor levels are finely tuned to regulate stochastic expression, setting the ratio of alternative fates and ultimately determining color preference.

Nuclear microenvironments modulate transcription from low-affinity enhancers

Transcription factors bind low-affinity DNA sequences for only short durations. It is not clear how brief, low-affinity interactions can drive efficient transcription. Here, we report that the transcription factor Ultrabithorax (Ubx) utilizes low-affinity binding sites in the Drosophila melanogaster shavenbaby (svb) locus and related enhancers in nuclear microenvironments of high Ubx concentrations. Related enhancers colocalize to the same microenvironments independently of their chromosomal location, suggesting that microenvironments are highly differentiated transcription domains. Manipulating the affinity of svb enhancers revealed an inverse relationship between enhancer affinity and Ubx concentration required for transcriptional activation. The Ubx cofactor, Homothorax (Hth), was co-enriched with Ubx near enhancers that require Hth, even though Ubx and Hth did not co-localize throughout the nucleus. Thus, microenvironments of high local transcription factor and cofactor concentrations could help low-affinity sites overcome their kinetic inefficiency. Mechanisms that generate these microenvironments could be a general feature of eukaryotic transcriptional regulation.

Signaling and Gene Regulatory Networks in Mammalian Lens Development

Ocular lens development represents an advantageous system in which to study regulatory mechanisms governing cell fate decisions, extracellular signaling, cell and tissue organization, and the underlying gene regulatory networks. Spatiotemporally regulated domains of BMP, FGF, and other signaling molecules in late gastrula-early neurula stage embryos generate the border region between the neural plate and non-neural ectoderm from which multiple cell types, including lens progenitor cells, emerge and undergo initial tissue formation. Extracellular signaling and DNA-binding transcription factors govern lens and optic cup morphogenesis. Pax6, c-Maf, Hsf4, Prox1, Sox1, and a few additional factors regulate the expression of the lens structural proteins, the crystallins. Extensive crosstalk between a diverse array of signaling pathways controls the complexity and order of lens morphogenetic processes and lens transparency.

Quantifying the Impact of Non-coding Variants on Transcription Factor-DNA Binding

Many recent studies have emphasized the importance of genetic variants and mutations in cancer and other complex human diseases. The overwhelming majority of these variants occur in non-coding portions of the genome, where they can have a functional impact by disrupting regulatory interactions between transcription factors (TFs) and DNA. Here, we present a method for assessing the impact of non-coding mutations on TF-DNA interactions, based on regression models of DNA-binding specificity trained on high-throughput in vitro data. We use ordinary least squares (OLS) to estimate the parameters of the binding model for each TF, and we show that our predictions of TF-binding changes due to DNA mutations correlate well with measured changes in gene expression. In addition, by leveraging distributional results associated with OLS estimation, for each predicted change in TF binding we also compute a normalized score (z-score) and a significance value (p-value) reflecting our confidence that the mutation affects TF binding. We use this approach to analyze a large set of pathogenic non-coding variants, and we show that these variants lead to significant differences in TF binding between alleles, compared to a control set of common variants. Thus, our results indicate that there is a strong regulatory component to the pathogenic non-coding variants identified thus far.

Systematic dissection of genomic features determining transcription factor binding and enhancer function

Enhancers regulate gene expression through the binding of sequence-specific transcription factors (TFs) to cognate motifs. Various features influence TF binding and enhancer function-including the chromatin state of the genomic locus, the affinities of the binding site, the activity of the bound TFs, and interactions among TFs. However, the precise nature and relative contributions of these features remain unclear. Here, we used massively parallel reporter assays (MPRAs) involving 32,115 natural and synthetic enhancers, together with high-throughput in vivo binding assays, to systematically dissect the contribution of each of these features to the binding and activity of genomic regulatory elements that contain motifs for PPARγ, a TF that serves as a key regulator of adipogenesis. We show that distinct sets of features govern PPARγ binding vs. enhancer activity. PPARγ binding is largely governed by the affinity of the specific motif site and higher-order features of the larger genomic locus, such as chromatin accessibility. In contrast, the enhancer activity of PPARγ binding sites depends on varying contributions from dozens of TFs in the immediate vicinity, including interactions between combinations of these TFs. Different pairs of motifs follow different interaction rules, including subadditive, additive, and superadditive interactions among specific classes of TFs, with both spatially constrained and flexible grammars. Our results provide a paradigm for the systematic characterization of the genomic features underlying regulatory elements, applicable to the design of synthetic regulatory elements or the interpretation of human genetic variation.

The Gene Regulatory Network of Lens Induction Is Wired through Meis-Dependent Shadow Enhancers of Pax6

Lens induction is a classical developmental model allowing investigation of cell specification, spatiotemporal control of gene expression, as well as how transcription factors are integrated into highly complex gene regulatory networks (GRNs). Pax6 represents a key node in the gene regulatory network governing mammalian lens induction. Meis1 and Meis2 homeoproteins are considered as essential upstream regulators of Pax6 during lens morphogenesis based on their interaction with the ectoderm enhancer (EE) located upstream of Pax6 transcription start site. Despite this generally accepted regulatory pathway, Meis1-, Meis2- and EE-deficient mice have surprisingly mild eye phenotypes at placodal stage of lens development. Here, we show that simultaneous deletion of Meis1 and Meis2 in presumptive lens ectoderm results in arrested lens development in the pre-placodal stage, and neither lens placode nor lens is formed. We found that in the presumptive lens ectoderm of Meis1/Meis2 deficient embryos Pax6 expression is absent. We demonstrate using chromatin immunoprecipitation (ChIP) that in addition to EE, Meis homeoproteins bind to a remote, ultraconserved SIMO enhancer of Pax6. We further show, using in vivo gene reporter analyses, that the lens-specific activity of SIMO enhancer is dependent on the presence of three Meis binding sites, phylogenetically conserved from man to zebrafish. Genetic ablation of EE and SIMO enhancers demostrates their requirement for lens induction and uncovers an apparent redundancy at early stages of lens development. These findings identify a genetic requirement for Meis1 and Meis2 during the early steps of mammalian eye development. Moreover, they reveal an apparent robustness in the gene regulatory mechanism whereby two independent "shadow enhancers" maintain critical levels of a dosage-sensitive gene, Pax6, during lens induction.

Evolved Repression Overcomes Enhancer Robustness

Biological systems display extraordinary robustness. Robustness of transcriptional enhancers results mainly from clusters of binding sites for the same transcription factor, and it is not clear how robust enhancers can evolve loss of expression through point mutations. Here, we report the high-resolution functional dissection of a robust enhancer of the shavenbaby gene that has contributed to morphological evolution. We found that robustness is encoded by many binding sites for the transcriptional activator Arrowhead and that, during evolution, some of these activator sites were lost, weakening enhancer activity. Complete silencing of enhancer function, however, required evolution of a binding site for the spatially restricted potent repressor Abrupt. These findings illustrate that recruitment of repressor binding sites can overcome enhancer robustness and may minimize pleiotropic consequences of enhancer evolution. Recruitment of repression may be a general mode of evolution to break robust regulatory linkages.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

Genetic Analysis of 'PAX6-Negative' Individuals with Aniridia or Gillespie Syndrome

We report molecular genetic analysis of 42 affected individuals referred with a diagnosis of aniridia who previously screened as negative for intragenic PAX6 mutations. Of these 42, the diagnoses were 31 individuals with aniridia and 11 individuals referred with a diagnosis of Gillespie syndrome (iris hypoplasia, ataxia and mild to moderate developmental delay). Array-based comparative genomic hybridization identified six whole gene deletions: four encompassing PAX6 and two encompassing FOXC1. Six deletions with plausible cis-regulatory effects were identified: five that were 3' (telomeric) to PAX6 and one within a gene desert 5' (telomeric) to PITX2. Sequence analysis of the FOXC1 and PITX2 coding regions identified two plausibly pathogenic de novo FOXC1 missense mutations (p.Pro79Thr and p.Leu101Pro). No intragenic mutations were detected in PITX2. FISH mapping in an individual with Gillespie-like syndrome with an apparently balanced X;11 reciprocal translocation revealed disruption of a gene at each breakpoint: ARHGAP6 on the X chromosome and PHF21A on chromosome 11. In the other individuals with Gillespie syndrome no mutations were identified in either of these genes, or in HCCS which lies close to the Xp breakpoint. Disruption of PHF21A has previously been implicated in the causation of intellectual disability (but not aniridia). Plausibly causative mutations were identified in 15 out of 42 individuals (12/32 aniridia; 3/11 Gillespie syndrome). Fourteen of these mutations presented in the known aniridia genes; PAX6, FOXC1 and PITX2. The large number of individuals in the cohort with no mutation identified suggests greater locus heterogeneity may exist in both isolated and syndromic aniridia than was previously appreciated.

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.

To Be Specific or Not: The Critical Relationship Between Hox And TALE Proteins

Hox proteins are key regulatory transcription factors that act in different tissues of the embryo to provide specific spatial and temporal coordinates to each cell. These patterning functions often depend on the presence of the TALE-homeodomain class cofactors, which form cooperative DNA-binding complexes with all Hox proteins. How this family of cofactors contributes to the highly diverse and specific functions of Hox proteins in vivo remains an important unsolved question. We review here the most recent advances in understanding the molecular mechanisms underlying Hox-TALE function. In particular, we discuss the role of DNA shape, DNA-binding affinity, and protein-protein interaction flexibility in dictating Hox-TALE specificity. We propose several models to explain how these mechanisms are integrated with each other in the context of the many distinct functions that Hox and TALE factors carry out in vivo.

Systems biology of lens development: A paradigm for disease gene discovery in the eye

Over the past several decades, the biology of the developing lens has been investigated using molecular genetics-based approaches in various vertebrate model systems. These efforts, involving target gene knockouts or knockdowns, have led to major advances in our understanding of lens morphogenesis and the pathological basis of cataracts, as well as of other lens related eye defects. In particular, we now have a functional understanding of regulators such as Pax6, Six3, Sox2, Oct1 (Pou2f1), Meis1, Pnox1, Zeb2 (Sip1), Mab21l1, Foxe3, Tfap2a (Ap2-alpha), Pitx3, Sox11, Prox1, Sox1, c-Maf, Mafg, Mafk, Hsf4, Fgfrs, Bmp7, and Tdrd7 in this tissue. However, whether these individual regulators interact or their targets overlap, and the significance of such interactions during lens morphogenesis, is not well defined. The arrival of high-throughput approaches for gene expression profiling (microarrays, RNA-sequencing (RNA-seq), etc.), which can be coupled with chromatin immunoprecipitation (ChIP) or RNA immunoprecipitation (RIP) assays, along with improved computational resources and publically available datasets (e.g. those containing comprehensive protein-protein, protein-DNA information), presents new opportunities to advance our understanding of the lens tissue on a global systems level. Such systems-level knowledge will lead to the derivation of the underlying lens gene regulatory network (GRN), defined as a circuit map of the regulator-target interactions functional in lens development, which can be applied to expedite cataract gene discovery. In this review, we cover the various systems-level approaches such as microarrays, RNA-seq, and ChIP that are already being applied to lens studies and discuss strategies for assembling and interpreting these vast amounts of high-throughput information for effective dispersion to the scientific community. In particular, we discuss strategies for effective interpretation of this new information in the context of the rich knowledge obtained through the application of traditional single-gene focused experiments on the lens. Finally, we discuss our vision for integrating these diverse high-throughput datasets in a single web-based user-friendly tool iSyTE (integrated Systems Tool for Eye gene discovery) - a resource that is already proving effective in the identification and characterization of genes linked to lens development and cataract. We anticipate that application of a similar approach to other ocular tissues such as the retina and the cornea, and even other organ systems, will significantly impact disease gene discovery.

The Soft Touch: Low-Affinity Transcription Factor Binding Sites in Development and Evolution

Transcription factor proteins regulate gene expression by binding to specific DNA regions. Most studies of transcription factor binding sites have focused on the highest affinity sites for each factor. There is abundant evidence, however, that binding sites with a range of affinities, including very low affinities, are critical to gene regulation. Here, we present the theoretical and experimental evidence for the importance of low-affinity sites in gene regulation and development. We also discuss the implications of the widespread use of low-affinity sites in eukaryotic genomes for robustness, precision, specificity, and evolution of gene regulation.

Loss of the Homeodomain Transcription Factor Prep1 Perturbs Adult Hematopoiesis in the Bone Marrow

Prep1, a TALE-family homeodomain transcription factor, has been demonstrated to play a critical role in embryonic hematopoiesis, as its insufficiency caused late embryonic lethality associated with defective hematopoiesis and angiogenesis. In the present study, we generated hematopoietic- and endothelial cell-specific Prep1-deficient mice and demonstrated that expression of Prep1 in the hematopoietic cell compartment is not essential for either embryonic or adult hematopoiesis, although its absence causes significant hematopoietic abnormalities in the adult bone marrow. Loss of Prep1 promotes cell cycling of hematopoietic stem/progenitor cells (HSPC), leading to the expansion of the HSPC pool. Prep1 deficiency also results in the accumulation of lineage-committed progenitors, increased monocyte/macrophage differentiation and arrested erythroid maturation. Maturation of T cells and B cells is also perturbed in Prep-deficient mice. These findings provide novel insight into the pleiotropic roles of Prep1 in adult hematopoiesis that were unrecognized in previous studies using germline Prep1 hypomorphic mice.

Deficiency of the RNA binding protein caprin2 causes lens defects and features of Peters anomaly

BACKGROUND

It was recently demonstrated that deficiency of a conserved RNA binding protein (RBP) and RNA granule (RG) component Tdrd7 causes ocular defects including cataracts in human, mouse and chicken, indicating the importance of posttranscriptional regulation in eye development. Here we investigated the function of a second conserved RBP/RG component Caprin2 that is identified by the eye gene discovery tool iSyTE.

RESULTS

In situ hybridization, Western blotting and immunostaining confirmed highly enriched expression of Caprin2 mRNA and protein in mouse embryonic and postnatal lens. To gain insight into its function, lens-specific Caprin2 conditional knockout (cKO) mouse mutants were generated using a lens-Cre deleter line Pax6GFPCre. Phenotypic analysis of Caprin2(cKO/cKO) mutants revealed distinct eye defects at variable penetrance. Wheat germ agglutinin staining and scanning electron microscopy demonstrated that Caprin2(cKO/cKO) mutants have an abnormally compact lens nucleus, which is the core of the lens comprised of centrally located terminally differentiated fiber cells. Additionally, Caprin2(cKO/cKO) mutants also exhibited at 8% penetrance a developmental defect that resembles a human condition called Peters anomaly, wherein the lens and the cornea remain attached by a persistent stalk.

CONCLUSIONS

These data suggest that a conserved RBP Caprin2 functions in distinct morphological events in mammalian eye development.

Temporal hierarchy of gene expression mediated by transcription factor binding affinity and activation dynamics

Understanding cellular responses to environmental stimuli requires not only the knowledge of specific regulatory components but also the quantitative characterization of the magnitude and timing of regulatory events. The two-component system is one of the major prokaryotic signaling schemes and is the focus of extensive interest in quantitative modeling and investigation of signaling dynamics. Here we report how the binding affinity of the PhoB two-component response regulator (RR) to target promoters impacts the level and timing of expression of PhoB-regulated genes. Information content has often been used to assess the degree of conservation for transcription factor (TF)-binding sites. We show that increasing the information content of PhoB-binding sites in designed phoA promoters increased the binding affinity and that the binding affinity and concentration of phosphorylated PhoB (PhoB~P) together dictate the level and timing of expression of phoA promoter variants. For various PhoB-regulated promoters with distinct promoter architectures, expression levels appear not to be correlated with TF-binding affinities, in contrast to the intuitive and oversimplified assumption that promoters with higher affinity for a TF tend to have higher expression levels. However, the expression timing of the core set of PhoB-regulated genes correlates well with the binding affinity of PhoB~P to individual promoters and the temporal hierarchy of gene expression appears to be related to the function of gene products during the phosphate starvation response. Modulation of the information content and binding affinity of TF-binding sites may be a common strategy for temporal programming of the expression profile of RR-regulated genes.

IMPORTANCE

A single TF often orchestrates the expression of multiple genes in response to environmental stimuli. It is not clear how different TF-binding sites within the regulon dictate the expression profile. Our studies of Escherichia coli PhoB, a response regulator that controls expression of a core set of phosphate assimilation genes in response to phosphate starvation, showed that expression levels of PhoB-regulated genes are under sophisticated control and do not follow a simple correlation with the binding affinity of PhoB~P to individual promoters. However, the expression timing correlates with the PhoB-binding affinity and gene functions. Genes involved in direct Pi uptake contain high-affinity sites and are transcribed earlier than genes involved in phosphorus scavenging. This illustrates an elaborate mechanism of temporally programmed gene expression, even for nondevelopmental pathways.

ChIP-Seq and RNA-Seq analyses identify components of the Wnt and Fgf signaling pathways as Prep1 target genes in mouse embryonic stem cells

The Prep1 (Pknox1) homeodomain transcription factor is essential at multiple stages of embryo development. In the E11.5 embryo trunk, we previously estimated that Prep1 binds about 3,300 genomic sites at a highly specific decameric consensus sequence, mainly governing basal cellular functions. We now show that in embryonic stem (ES) cells Prep1 binding pattern only partly overlaps that of the embryo trunk, with about 2,000 novel sites. Moreover, in ES cells Prep1 still binds mostly to promoters, as in total embryo trunk but, among the peaks bound exclusively in ES cells, the percentage of enhancers was three-fold higher. RNA-seq identifies about 1800 genes down-regulated in Prep1-/- ES cells which belong to gene ontology categories not enriched in the E11.5 Prep1i/i differentiated embryo, including in particular essential components of the Wnt and Fgf pathways. These data agree with aberrant Wnt and Fgf expression levels in the Prep1-/- ES cells with a deficient embryoid bodies (EBs) formation and differentiation. Re-establishment of the Prep1 level rescues the phenotype.

Evolution of transcription factor function as a mechanism for changing metazoan developmental gene regulatory networks

The form that an animal takes during development is directed by gene regulatory networks (GRNs). Developmental GRNs interpret maternally deposited molecules and externally supplied signals to direct cell-fate decisions, which ultimately leads to the arrangements of organs and tissues in the organism. Genetically encoded modifications to these networks have generated the wide range of metazoan diversity that exists today. Most studies of GRN evolution focus on changes to cis-regulatory DNA, and it was historically theorized that changes to the transcription factors that bind to these cis-regulatory modules (CRMs) contribute to this process only rarely. A growing body of evidence suggests that changes to the coding regions of transcription factors play a much larger role in the evolution of developmental gene regulatory networks than originally imagined. Just as cis-regulatory changes make use of modular binding site composition and tissue-specific modules to avoid pleiotropy, transcription factor coding regions also predominantly evolve in ways that limit the context of functional effects. Here, we review the recent works that have led to this unexpected change in the field of Evolution and Development (Evo-Devo) and consider the implications these studies have had on our understanding of the evolution of developmental processes.

Low affinity binding site clusters confer hox specificity and regulatory robustness

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.