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

RLS-associated MEIS transcription factors control distinct processes in human neural stem cells

MEIS1 and MEIS2 encode highly conserved homeodomain transcription factors crucial for developmental processes in a wide range of tissues, including the brain. They can execute redundant functions when co-expressed in the same cell types, but their roles during early stages of neural differentiation have not been systematically compared. By separate knockout and overexpression of MEIS1 and MEIS2 in human neural stem cells, we find they control specific sets of target genes, associated with distinct biological processes. Integration of DNA binding sites with differential transcriptomics implicates MEIS1 to co-regulate gene expression by interaction with transcription factors of the SOX and FOX families. MEIS1 harbors the strongest risk factor for restless legs syndrome (RLS). Our data suggest that MEIS1 can directly regulate the RLS-associated genes NTNG1, MDGA1 and DACH1, constituting new approaches to study the elusive pathomechanism or RLS.

Decoding the epigenetics and chromatin loop dynamics of androgen receptor-mediated transcription

Androgen receptor (AR)-mediated transcription plays a critical role in development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigate the AR-driven epigenetic and chromosomal chromatin looping changes by generating a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we find that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we show that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then propose and experimentally validate an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide insights into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Genetic analysis of medaka fish illuminates conserved and divergent roles of Pax6 in vertebrate eye development

Landmark discovery of eye defects caused by Pax6 gene mutations in humans, rodents, and even fruit flies combined with Pax6 gene expression studies in various phyla, led to the master control gene hypothesis postulating that the gene is required almost universally for animal visual system development. However, this assumption has not been broadly tested in genetically trackable organisms such as vertebrates. Here, to determine the functional role of the fish orthologue of mammalian Pax6 in eye development we analyzed mutants in medaka Pax6.1 gene generated by genome editing. We found that transcription factors implicated in vertebrate lens development (Prox1a, MafB, c-Maf, FoxE3) failed to initiate expression in the presumptive lens tissue of Pax6.1 mutant fish resulting in aphakia, a phenotype observed previously in Pax6 mutant mice. Surprisingly, the overall differentiation potential of Pax6.1-deficient retinal progenitor cells (RPCs) is not severely compromised, and the only cell types affected by the absence of Pax6.1 transcription factor are retinal ganglion cells. This is in stark contrast to the situation in mice where the Pax6 gene is required cell-autonomously for the expansion of RPCs, and the differentiation of all retina cell types. Our results provide novel insight into the conserved and divergent roles of Pax6 gene orthologues in vertebrate eye development indicating that the lens-specific role is more evolutionarily conserved than the role in retina differentiation.

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

Multiple cis-regulatory elements control prox1a expression in distinct lymphatic vascular beds

During embryonic development, lymphatic endothelial cell (LEC) precursors are distinguished from blood endothelial cells by the expression of Prospero-related homeobox 1 (Prox1), which is essential for lymphatic vasculature formation in mouse and zebrafish. Prox1 expression initiation precedes LEC sprouting and migration, serving as the marker of specified LECs. Despite its crucial role in lymphatic development, Prox1 upstream regulation in LECs remains to be uncovered. SOX18 and COUP-TFII are thought to regulate Prox1 in mice by binding its promoter region. However, the specific regulation of Prox1 expression in LECs remains to be studied in detail. Here, we used evolutionary conservation and chromatin accessibility to identify enhancers located in the proximity of zebrafish prox1a active in developing LECs. We confirmed the functional role of the identified sequences through CRISPR/Cas9 mutagenesis of a lymphatic valve enhancer. The deletion of this region results in impaired valve morphology and function. Overall, our results reveal an intricate control of prox1a expression through a collection of enhancers. Ray-finned fish-specific distal enhancers drive pan-lymphatic expression, whereas vertebrate-conserved proximal enhancers refine expression in functionally distinct subsets of lymphatic endothelium.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Decoding the Epigenetics and Chromatin Loop Dynamics of Androgen Receptor-Mediated Transcription

Androgen receptor (AR)-mediated transcription plays a critical role in normal prostate development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigated the AR-driven epigenetic and chromosomal chromatin looping changes. We collected a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we found that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we showed that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then proposed and experimentally validated an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide new insight into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Unique activities of two overlapping PAX6 retinal enhancers

Enhancers play a critical role in development by precisely modulating spatial, temporal, and cell type-specific gene expression. Sequence variants in enhancers have been implicated in diseases; however, establishing the functional consequences of these variants is challenging because of a lack of understanding of precise cell types and developmental stages where the enhancers are normally active. PAX6 is the master regulator of eye development, with a regulatory landscape containing multiple enhancers driving the expression in the eye. Whether these enhancers perform additive, redundant or distinct functions is unknown. Here, we describe the precise cell types and regulatory activity of two PAX6 retinal enhancers, HS5 and NRE. Using a unique combination of live imaging and single-cell RNA sequencing in dual enhancer-reporter zebrafish embryos, we uncover differences in the spatiotemporal activity of these enhancers. Our results show that although overlapping, these enhancers have distinct activities in different cell types and therefore likely nonredundant functions. This work demonstrates that unique cell type-specific activities can be uncovered for apparently similar enhancers when investigated at high resolution in vivo.

Conserved enhancers control notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

The interaction between enhancer variants and environmental factors as an overlooked aetiological paradigm in human complex disease

The interactions between genetic and environmental risk factors contribute to the aetiology of complex human diseases. Genome-wide association studies (GWAS) have revealed that most of the genetic variants associated with complex diseases are located in the non-coding part of the genome, preferentially within enhancers. Enhancers are distal cis-regulatory elements composed of clusters of transcription factors binding sites that positively regulate the expression of their target genes. The generation of genome-wide maps for histone marks (e.g., H3K27ac), chromatin accessibility and transcription factor and coactivator (e.g., p300) binding profiles have enabled the identification of enhancers across many human cell types and tissues. Nonetheless, the functional and pathological consequences of the majority of disease-associated genetic variants located within enhancers seem to be rather minor under normal conditions, thus questioning their medical relevance. Here we propose that, due to the prevalence of enhancer redundancy, the pathological effects of many disease-associated non-coding genetic variants might be preferentially (or even only) manifested under environmental stress.

Rho enhancers play unexpectedly minor roles in Rhodopsin transcription and rod cell integrity

Enhancers function with a basal promoter to control the transcription of target genes. Enhancer regulatory activity is often studied using reporter-based transgene assays. However, unmatched results have been reported when selected enhancers are silenced in situ. In this study, using genomic deletion analysis in mice, we investigated the roles of two previously identified enhancers and the promoter of the Rho gene that codes for the visual pigment rhodopsin. The Rho gene is robustly expressed by rod photoreceptors of the retina, and essential for the subcellular structure and visual function of rod photoreceptors. Mutations in RHO cause severe vision loss in humans. We found that each Rho regulatory region can independently mediate local epigenomic changes, but only the promoter is absolutely required for establishing active Rho chromatin configuration and transcription and maintaining the cell integrity and function of rod photoreceptors. To our surprise, two Rho enhancers that enable strong promoter activation in reporter assays are largely dispensable for Rho expression in vivo. Only small and age-dependent impact is detectable when both enhancers are deleted. Our results demonstrate context-dependent roles of enhancers and highlight the importance of studying functions of cis-regulatory regions in the native genomic context.

Regulatory changes associated with the head to trunk developmental transition

BACKGROUND

Development of vertebrate embryos is characterized by early formation of the anterior tissues followed by the sequential extension of the axis at their posterior end to build the trunk and tail structures, first by the activity of the primitive streak and then of the tail bud. Embryological, molecular and genetic data indicate that head and trunk development are significantly different, suggesting that the transition into the trunk formation stage involves major changes in regulatory gene networks.

RESULTS

We explored those regulatory changes by generating differential interaction networks and chromatin accessibility profiles from the posterior epiblast region of mouse embryos at embryonic day (E)7.5 and E8.5. We observed changes in various cell processes, including several signaling pathways, ubiquitination machinery, ion dynamics and metabolic processes involving lipids that could contribute to the functional switch in the progenitor region of the embryo. We further explored the functional impact of changes observed in Wnt signaling associated processes, revealing a switch in the functional relevance of Wnt molecule palmitoleoylation, essential during gastrulation but becoming differentially required for the control of axial extension and progenitor differentiation processes during trunk formation. We also found substantial changes in chromatin accessibility at the two developmental stages, mostly mapping to intergenic regions and presenting differential footprinting profiles to several key transcription factors, indicating a significant switch in the regulatory elements controlling head or trunk development. Those chromatin changes are largely independent of retinoic acid, despite the key role of this factor in the transition to trunk development. We also tested the functional relevance of potential enhancers identified in the accessibility assays that reproduced the expression profiles of genes involved in the transition. Deletion of these regions by genome editing had limited effect on the expression of those genes, suggesting the existence of redundant enhancers that guarantee robust expression patterns.

CONCLUSIONS

This work provides a global view of the regulatory changes controlling the switch into the axial extension phase of vertebrate embryonic development. It also revealed mechanisms by which the cellular context influences the activity of regulatory factors, channeling them to implement one of several possible biological outputs.

Congenital aniridia beyond black eyes: From phenotype and novel genetic mechanisms to innovative therapeutic approaches

Congenital PAX6-aniridia, initially characterized by the absence of the iris, has progressively been shown to be associated with other developmental ocular abnormalities and systemic features making congenital aniridia a complex syndromic disorder rather than a simple isolated disease of the iris. Moreover, foveal hypoplasia is now recognized as a more frequent feature than complete iris hypoplasia and a major visual prognosis determinant, reversing the classical clinical picture of this disease. Conversely, iris malformation is also a feature of various anterior segment dysgenesis disorders caused by PAX6-related developmental genes, adding a level of genetic complexity for accurate molecular diagnosis of aniridia. Therefore, the clinical recognition and differential genetic diagnosis of PAX6-related aniridia has been revealed to be much more challenging than initially thought, and still remains under-investigated. Here, we update specific clinical features of aniridia, with emphasis on their genotype correlations, as well as provide new knowledge regarding the PAX6 gene and its mutational spectrum, and highlight the beneficial utility of clinically implementing targeted Next-Generation Sequencing combined with Whole-Genome Sequencing to increase the genetic diagnostic yield of aniridia. We also present new molecular mechanisms underlying aniridia and aniridia-like phenotypes. Finally, we discuss the appropriate medical and surgical management of aniridic eyes, as well as innovative therapeutic options. Altogether, these combined clinical-genetic approaches will help to accelerate time to diagnosis, provide better determination of the disease prognosis and management, and confirm eligibility for future clinical trials or genetic-specific therapies.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

Dynamic changes in whole genome DNA methylation, chromatin and gene expression during mouse lens differentiation

BACKGROUND

Cellular differentiation is marked by temporally and spatially coordinated gene expression regulated at multiple levels. DNA methylation represents a universal mechanism to control chromatin organization and its accessibility. Cytosine methylation of CpG dinucleotides regulates binding of methylation-sensitive DNA-binding transcription factors within regulatory regions of transcription, including promoters and distal enhancers. Ocular lens differentiation represents an advantageous model system to examine these processes as lens comprises only two cell types, the proliferating lens epithelium and postmitotic lens fiber cells all originating from the epithelium.

RESULTS

Using whole genome bisulfite sequencing (WGBS) and microdissected lenses, we investigated dynamics of DNA methylation and chromatin changes during mouse lens fiber and epithelium differentiation between embryos (E14.5) and newborns (P0.5). Histone H3.3 variant chromatin landscapes were also generated for both P0.5 lens epithelium and fibers by chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Tissue-specific features of DNA methylation patterns are demonstrated via comparative studies with embryonic stem (ES) cells and neural progenitor cells (NPCs) at Nanog, Pou5f1, Sox2, Pax6 and Six3 loci. Comparisons with ATAC-seq and RNA-seq data demonstrate that reduced methylation is associated with increased expression of fiber cell abundant genes, including crystallins, intermediate filament (Bfsp1 and Bfsp2) and gap junction proteins (Gja3 and Gja8), marked by high levels of histone H3.3 within their transcribed regions. Interestingly, Pax6-binding sites exhibited predominantly DNA hypomethylation in lens chromatin. In vitro binding of Pax6 proteins showed Pax6's ability to interact with sites containing one or two methylated CpG dinucleotides.

CONCLUSIONS

Our study has generated the first data on methylation changes between two different stages of mammalian lens development and linked these data with chromatin accessibility maps, presence of histone H3.3 and gene expression. Reduced DNA methylation correlates with expression of important genes involved in lens morphogenesis and lens fiber cell differentiation.

Monoallelic variants resulting in substitutions of MAB21L1 Arg51 Cause Aniridia and microphthalmia

Classical aniridia is a congenital and progressive panocular disorder almost exclusively caused by heterozygous loss-of-function variants at the PAX6 locus. We report nine individuals from five families with severe aniridia and/or microphthalmia (with no detectable PAX6 mutation) with ultrarare monoallelic missense variants altering the Arg51 codon of MAB21L1. These mutations occurred de novo in 3/5 families, with the remaining families being compatible with autosomal dominant inheritance. Mice engineered to carry the p.Arg51Leu change showed a highly-penetrant optic disc anomaly in heterozygous animals with severe microphthalmia in homozygotes. Substitutions of the same codon (Arg51) in MAB21L2, a close homolog of MAB21L1, cause severe ocular and skeletal malformations in humans and mice. The predicted nucleotidyltransferase function of MAB21L1 could not be demonstrated using purified protein with a variety of nucleotide substrates and oligonucleotide activators. Induced expression of GFP-tagged wildtype and mutant MAB21L1 in human cells caused only modest transcriptional changes. Mass spectrometry of immunoprecipitated protein revealed that both mutant and wildtype MAB21L1 associate with transcription factors that are known regulators of PAX6 (MEIS1, MEIS2 and PBX1) and with poly(A) RNA binding proteins. Arg51 substitutions reduce the association of wild-type MAB21L1 with TBL1XR1, a component of the NCoR complex. We found limited evidence for mutation-specific interactions with MSI2/Musashi-2, an RNA-binding proteins with effects on many different developmental pathways. Given that biallelic loss-of-function variants in MAB21L1 result in a milder eye phenotype we suggest that Arg51-altering monoallelic variants most plausibly perturb eye development via a gain-of-function mechanism.

Paralogous Genes Involved in Embryonic Development: Lessons from the Eye and Other Tissues

During evolution, gene duplications lead to a naturally increased gene dosage. Duplicated genes can be further retained or eliminated over time by purifying selection pressure. The retention probability is increased by functional diversification and by the acquisition of novel functions. Interestingly, functionally diverged paralogous genes can maintain a certain level of functional redundancy and at least a partial ability to replace each other. In such cases, diversification probably occurred at the level of transcriptional regulation. Nevertheless, some duplicated genes can maintain functional redundancy after duplication and the ability to functionally compensate for the loss of each other. Many of them are involved in proper embryonic development. The development of particular tissues/organs and developmental processes can be more or less sensitive to the overall gene dosage. Alterations in the gene dosage or a decrease below a threshold level may have dramatic phenotypic consequences or even lead to embryonic lethality. The number of functional alleles of particular paralogous genes and their mutual cooperation and interactions influence the gene dosage, and therefore, these factors play a crucial role in development. This review will discuss individual interactions between paralogous genes and gene dosage sensitivity during development. The eye was used as a model system, but other tissues are also included.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

Multiple roles of Pax6 in postnatal cornea development

Mammalian corneal development is a multistep process, including formation of the corneal epithelium (CE), endothelium and stroma during embryogenesis, followed by postnatal stratification of the epithelial layers and continuous renewal of the epithelium to replace the outermost corneal cells. Here, we employed the Cre-loxP system to conditionally deplete Pax6 proteins in two domains of ocular cells, i.e., the ocular surface epithelium (cornea, limbus and conjunctiva) (OSE) or postnatal CE via K14-cre or Aldh3-cre, respectively. Earlier and broader inactivation of Pax6 in the OSE resulted in thickened OSE with CE and limbal cells adopting the conjunctival keratin expression pattern. More restricted depletion of Pax6 in postnatal CE resulted in an abnormal cornea marked by reduced epithelial thickness despite increased epithelial cell proliferation. Immunofluorescence studies revealed loss of intermediate filament Cytokeratin 12 and diffused expression of adherens junction components, together with reduced tight junction protein, Zonula occludens-1. Furthermore, the expression of Cytokeratin 14, a basal cell marker in apical layers, indicates impaired differentiation of CE cells. Collectively, our data demonstrate that Pax6 is essential for maintaining proper differentiation and strong intercellular adhesion in postnatal CE cells, whereas limbal Pax6 is required to prevent the outgrowth of conjunctival cells to the cornea.

Changes in DNA methylation hallmark alterations in chromatin accessibility and gene expression for eye lens differentiation

Background

Methylation at cytosines (mCG) is a well-known regulator of gene expression, but its requirements for cellular differentiation have yet to be fully elucidated. A well-studied cellular differentiation model system is the eye lens, consisting of a single anterior layer of epithelial cells that migrate laterally and differentiate into a core of fiber cells. Here, we explore the genome-wide relationships between mCG methylation, chromatin accessibility and gene expression during differentiation of eye lens epithelial cells into fiber cells.

Results

Whole genome bisulfite sequencing identified 7621 genomic loci exhibiting significant differences in mCG levels between lens epithelial and fiber cells. Changes in mCG levels were inversely correlated with the differentiation state-specific expression of 1285 genes preferentially expressed in either lens fiber or lens epithelial cells (Pearson correlation r = − 0.37, p < 1 × 10^–42). mCG levels were inversely correlated with chromatin accessibility determined by assay for transposase-accessible sequencing (ATAC-seq) (Pearson correlation r = − 0.86, p < 1 × 10^–300). Many of the genes exhibiting altered regions of DNA methylation, chromatin accessibility and gene expression levels in fiber cells relative to epithelial cells are associated with lens fiber cell structure, homeostasis and transparency. These include lens crystallins (CRYBA4, CRYBB1, CRYGN, CRYBB2), lens beaded filament proteins (BFSP1, BFSP2), transcription factors (HSF4, SOX2, HIF1A), and Notch signaling pathway members (NOTCH1, NOTCH2, HEY1, HES5). Analysis of regions exhibiting cell-type specific alterations in DNA methylation revealed an overrepresentation of consensus sequences of multiple transcription factors known to play key roles in lens cell differentiation including HIF1A, SOX2, and the MAF family of transcription factors.

Conclusions

Collectively, these results link DNA methylation with control of chromatin accessibility and gene expression changes required for eye lens differentiation. The results also point to a role for DNA methylation in the regulation of transcription factors previously identified to be important for lens cell differentiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00440-z.

Repurposing of the enhancer-promoter communication underlies the compensation of Mesp2 by Mesp1

Organisms are inherently equipped with buffering systems against genetic perturbations. Genetic compensation, the compensatory response by upregulating another gene or genes, is one such buffering mechanism. Recently, a well-conserved compensatory mechanism was proposed: transcriptional adaptation of homologs under the nonsense-mediated mRNA decay pathways. However, this model cannot explain the onset of all compensatory events. We report a novel genetic compensation mechanism operating over the Mesp gene locus. Mesp1 and Mesp2 are paralogs located adjacently in the genome. Mesp2 loss is partially rescued by Mesp1 upregulation in the presomitic mesoderm (PSM). Using a cultured PSM induction system, we reproduced the compensatory response in vitro and found that the Mesp2-enhancer is required to promote Mesp1. We revealed that the Mesp2-enhancer directly interacts with the Mesp1 promoter, thereby upregulating Mesp1 expression upon the loss of Mesp2. Of note, this interaction is established by genomic arrangement upon PSM development independently of Mesp2 disruption. We propose that the repurposing of this established enhancer-promoter communication is the mechanism underlying this compensatory response for the upregulation of the adjacent gene.

Genetic compensation, the compensatory response by upregulating another gene or genes, is one of the inherent mechanisms against gene disruption to confer cellular fitness. However, the regulatory mechanisms are largely unknown. Nonsense-mediated mutant mRNA degradation was recently proposed as a conserved mechanism across species to upregulate homologous genes to compensate for a disrupted gene, but this cannot explain compensation events with no mutant mRNA. This study investigated the compensation mechanism operating over adjacent paralogs, Mesp1 and Mesp2, in the genome. Mesp genes encode essential transcription factors in the presomitic mesoderm for development. In general, an enhancer is considered to activate a target gene when it physically interacts with the target. The communication of the Mesp2-enhancer with the Mesp1 promoter is established upon differentiation of the presomitic mesoderm, but this communication activates Mesp1 only when Mesp2 is disrupted, leading to compensation. We revealed a novel compensation mechanism depending on the repurposing of this enhancer-promoter communication by gene disruption. Our study also provides new insight into transcriptional regulation by providing the concept that an enhancer changes its target even among its physically interacting genes in a context-dependent manner.

No Need to Stick Together to Be Connected: Multiple Types of Enhancers' Networking

Simple Summary

Transcription regulation programs require the functional interaction of distal and proximal regulatory regions, interacting by specific 3D chromatin configurations. Enhancers are cis-acting regulatory elements able to promote gene expression regardless their orientation and distance from the transcription starting site. Their systematic mapping by genome-wide chromatin profiling and chromosome conformation analysis, combined with the development of gene-editing approaches to modulate their function, revealed that many enhancers work together to fine-tune the expression of their target genes. This review aim to describe the functions of different types of enhancers and the modalities of enhancers’ interaction, focusing on their role in the regulation of complex biological processes like cancer development.

Abstract

The control of gene expression at a transcriptional level requires a widespread landscape of regulatory elements. Central to these regulatory circuits are enhancers (ENHs), which are defined as cis-acting DNA elements able to increase the transcription of a target gene in a distance- and orientation-independent manner. ENHs are not independent functional elements but work in a complex and dynamic cooperative network, constituting the building blocks of multimodular domains of gene expression regulation. The information from each of these elements converges on the target promoter, contributing to improving the precision and sharpness of gene modulation. ENHs’ interplay varies in its nature and extent, ranging from an additive to redundant effect depending on contexts. Moving from super-enhancers that drive the high expression levels of identity genes, to shadow-enhancers, whose redundant functions contribute to buffering the variation in gene expression, this review aims to describe the different modalities of ENHs’ interaction and their role in the regulation of complex biological processes like cancer development.

dbCNS: A New Database for Conserved Noncoding Sequences

We developed dbCNS (http://yamasati.nig.ac.jp/dbcns), a new database for conserved noncoding sequences (CNSs). CNSs exist in many eukaryotes and are assumed to be involved in protein expression control. Version 1 of dbCNS, introduced here, includes a powerful and precise CNS identification pipeline for multiple vertebrate genomes. Mutations in CNSs may induce morphological changes and cause genetic diseases. For this reason, many vertebrate CNSs have been identified, with special reference to primate genomes. We integrated ∼6.9 million CNSs from many vertebrate genomes into dbCNS, which allows users to extract CNSs near genes of interest using keyword searches. In addition to CNSs, dbCNS contains published genome sequences of 161 species. With purposeful taxonomic sampling of genomes, users can employ CNSs as queries to reconstruct CNS alignments and phylogenetic trees, to evaluate CNS modifications, acquisitions, and losses, and to roughly identify species with CNSs having accelerated substitution rates. dbCNS also produces links to dbSNP for searching pathogenic single-nucleotide polymorphisms in human CNSs. Thus, dbCNS connects morphological changes with genetic diseases. A test analysis using 38 gnathostome genomes was accomplished within 30 s. dbCNS results can evaluate CNSs identified by other stand-alone programs using genome-scale data.

Enhancer redundancy in development and disease

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

Meis homeobox genes control progenitor competence in the retina

The vertebrate eye is derived from the neuroepithelium, surface ectoderm, and extracellular mesenchyme. The neuroepithelium forms an optic cup in which the spatial separation of three domains is established, namely, the region of multipotent retinal progenitor cells (RPCs), the ciliary margin zone (CMZ)-which possesses both a neurogenic and nonneurogenic potential-and the optic disk (OD), the interface between the optic stalk and the neuroretina. Here, we show by genetic ablation in the developing optic cup that Meis1 and Meis2 homeobox genes function redundantly to maintain the retinal progenitor pool while they simultaneously suppress the expression of genes characteristic of CMZ and OD fates. Furthermore, we demonstrate that Meis transcription factors bind regulatory regions of RPC-, CMZ-, and OD-specific genes, thus providing a mechanistic insight into the Meis-dependent gene regulatory network. Our work uncovers the essential role of Meis1 and Meis2 as regulators of cell fate competence, which organize spatial territories in the vertebrate eye.

Ciona Brachyury proximal and distal enhancers have different FGF dose-response relationships

Many genes are regulated by two or more enhancers that drive similar expression patterns. Evolutionary theory suggests that these seemingly redundant enhancers must have functionally important differences. In the simple ascidian chordate Ciona, the transcription factor Brachyury is induced exclusively in the presumptive notochord downstream of lineage specific regulators and FGF-responsive Ets family transcription factors. Here we exploit the ability to finely titrate FGF signaling activity via the MAPK pathway using the MEK inhibitor U0126 to quantify the dependence of transcription driven by different Brachyury reporter constructs on this direct upstream regulator. We find that the more powerful promoter-adjacent proximal enhancer and a weaker distal enhancer have fundamentally different dose-response relationships to MAPK inhibition. The Distal enhancer is more sensitive to MAPK inhibition but shows a less cooperative response, whereas the Proximal enhancer is less sensitive and more cooperative. A longer construct containing both enhancers has a complex dose-response curve that supports the idea that the proximal and distal enhancers are moderately super-additive. We show that the overall expression loss from intermediate doses of U0126 is not only a function of the fraction of cells expressing these reporters, but also involves graded decreases in expression at the single-cell level. Expression of the endogenous gene shows a comparable dose-response relationship to the full length reporter, and we find that different notochord founder cells are differentially sensitive to MAPK inhibition. Together, these results indicate that although the two Brachyury enhancers have qualitatively similar expression patterns, they respond to FGF in quantitatively different ways and act together to drive high levels of Brachyury expression with a characteristic input/output relationship. This indicates that they are fundamentally not equivalent genetic elements.

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.

Modulation of the Promoter Activation Rate Dictates the Transcriptional Response to Graded BMP Signaling Levels in the Drosophila Embryo

Morphogen gradients specify cell fates during development, with a classic example being the bone morphogenetic protein (BMP) gradient's conserved role in embryonic dorsal-ventral axis patterning. Here, we elucidate how the BMP gradient is interpreted in the Drosophila embryo by combining live imaging with computational modeling to infer transcriptional burst parameters at single-cell resolution. By comparing burst kinetics in cells receiving different levels of BMP signaling, we show that BMP signaling controls burst frequency by regulating the promoter activation rate. We provide evidence that the promoter activation rate is influenced by both enhancer and promoter sequences, whereas Pol II loading rate is primarily modulated by the enhancer. Consistent with BMP-dependent regulation of burst frequency, the numbers of BMP target gene transcripts per cell are graded across their expression domains. We suggest that graded mRNA output is a general feature of morphogen gradient interpretation and discuss how this can impact on cell-fate decisions.

The Atoh7 remote enhancer provides transcriptional robustness during retinal ganglion cell development

The retinal ganglion cell (RGC) competence factor ATOH7 is dynamically expressed during retinal histogenesis. ATOH7 transcription is controlled by a promoter-adjacent primary enhancer and a remote shadow enhancer (SE). Deletion of the ATOH7 human SE causes nonsyndromic congenital retinal nonattachment (NCRNA) disease, characterized by optic nerve aplasia and total blindness. We used genome editing to model NCRNA in mice. Deletion of the murine SE reduces Atoh7 messenger RNA (mRNA) fivefold but does not recapitulate optic nerve loss; however, SEdel/knockout (KO) trans heterozygotes have thin optic nerves. By analyzing Atoh7 mRNA and protein levels, RGC development and survival, and chromatin landscape effects, we show that the SE ensures robust Atoh7 transcriptional output. Combining SE deletion and KO and wild-type alleles in a genotypic series, we determined the amount of Atoh7 needed to produce a normal complement of adult RGCs, and the secondary consequences of graded reductions in Atoh7 dosage. Together, these data reveal the workings of an evolutionary fail-safe, a duplicate enhancer mechanism that is hard-wired in the machinery of vertebrate retinal ganglion cell genesis.

Neural crest cells require Meis2 for patterning the mandibular arch via the Sonic hedgehog pathway

Cranial neural crest cells (cNCCs) originate in the anterior neural tube and populate pharyngeal arches in which they contribute to formation of bone and cartilage. This cell population also provides molecular signals for the development of tissues of non-neural crest origin, such as the tongue muscles, teeth enamel or gland epithelium. Here we show that the transcription factor Meis2 is expressed in the oral region of the first pharyngeal arch (PA1) and later in the tongue primordium. Conditional inactivation of Meis2 in cNCCs resulted in loss of Sonic hedgehog signalling in the oropharyngeal epithelium and impaired patterning of PA1 along the lateral-medial and oral-aboral axis. Failure of molecular specification of PA1, illustrated by altered expression of Hand1/2, Dlx5, Barx1, Gsc and other markers, led to hypoplastic tongue and ectopic ossification of the mandible. Meis2-mutant mice thus display craniofacial defects that are reminiscent of several human syndromes and patients with mutations in the Meis2 gene.

Transcriptome analysis of neural progenitor cells derived from Lowe syndrome induced pluripotent stem cells: identification of candidate genes for the neurodevelopmental and eye manifestations

BACKGROUND

Lowe syndrome (LS) is caused by loss-of-function mutations in the X-linked gene OCRL, which codes for an inositol polyphosphate 5-phosphatase that plays a key role in endosome recycling, clathrin-coated pit formation, and actin polymerization. It is characterized by congenital cataracts, intellectual and developmental disability, and renal proximal tubular dysfunction. Patients are also at high risk for developing glaucoma and seizures. We recently developed induced pluripotent stem cell (iPSC) lines from three patients with LS who have hypomorphic variants affecting the 3' end of the gene, and their neurotypical brothers to serve as controls.

METHODS

In this study, we used RNA sequencing (RNA-seq) to obtain transcriptome profiles in LS and control neural progenitor cells (NPCs).

RESULTS

In a comparison of the patient and control NPCs (n = 3), we found 16 differentially expressed genes (DEGs) at the multiple test adjusted p value (padj) < 0.1, with nine at padj < 0.05. Using nominal p value < 0.05, 319 DEGs were detected. The relatively small number of DEGs could be due to the fact that OCRL is not a transcription factor per se, although it could have secondary effects on gene expression through several different mechanisms. Although the number of DEGs passing multiple test correction was small, those that were found are quite consistent with some of the known molecular effects of OCRL protein, and the clinical manifestations of LS. Furthermore, using gene set enrichment analysis (GSEA), we found that genes increased expression in the patient NPCs showed enrichments of several gene ontology (GO) terms (false discovery rate < 0.25): telencephalon development, pallium development, NPC proliferation, and cortex development, which are consistent with a condition characterized by intellectual disabilities and psychiatric manifestations. In addition, a significant enrichment among the nominal DEGs for genes implicated in autism spectrum disorder (ASD) was found (e.g., AFF2, DNER, DPP6, DPP10, RELN, CACNA1C), as well as several that are strong candidate genes for the development of eye problems found in LS, including glaucoma. The most notable example is EFEMP1, a well-known candidate gene for glaucoma and other eye pathologies.

CONCLUSION

Overall, the RNA-seq findings present several candidate genes that could help explain the underlying basis for the neurodevelopmental and eye problems seen in boys with LS.

Lampreys, the jawless vertebrates, contain three Pax6 genes with distinct expression in eye, brain and pancreas

The transcription factor Pax6 is crucial for the development of the central nervous system, eye, olfactory system and pancreas, and is implicated in human disease. While a single Pax6 gene exists in human and chicken, Pax6 occurs as a gene family in other vertebrates, with two members in elephant shark, Xenopus tropicalis and Anolis lizard and three members in teleost fish such as stickleback and medaka. However, the complement of Pax6 genes in jawless vertebrates (cyclostomes), the sister group of jawed vertebrates (gnathostomes), is unknown. Using a combination of BAC sequencing and genome analysis, we discovered three Pax6 genes in lampreys. Unlike the paired-less Pax6 present in some gnathostomes, all three lamprey Pax6 have a highly conserved full-length paired domain. All three Pax6 genes are expressed in the eye and brain, with variable expression in other tissues. Notably, lamprey Pax6α transcripts are found in the pancreas, a vertebrate-specific organ, indicating the involvement of Pax6 in development of the pancreas in the vertebrate ancestor. Multi-species sequence comparisons revealed only a single conserved non-coding element, in the lamprey Pax6β locus, with similarity to the PAX6 neuroretina enhancer. Using a transgenic zebrafish enhancer assay we demonstrate functional conservation of this element over 500 million years of vertebrate evolution.

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.

MEIS transcription factors in development and disease

MEIS transcription factors are key regulators of embryonic development and cancer. Research on MEIS genes in the embryo and in stem cell systems has revealed novel and surprising mechanisms by which these proteins control gene expression. This Primer summarizes recent findings about MEIS protein activity and regulation in development, and discusses new insights into the role of MEIS genes in disease, focusing on the pathogenesis of solid cancers.

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.

Dental Epithelial Stem Cells Express the Developmental Regulator Meis1

MEIS1 is a key developmental regulator of several organs and participates in stem cell maintenance in different niches. However, despite the murine continuously growing incisor being a well described model for the study of adult stem cells, Meis1 has not been investigated in a dental context. Here, we uncover that Meis1 expression in the tooth is confined to the epithelial compartment. Its expression arises during morphogenesis and becomes restricted to the mouse incisor epithelial stem cell niche, the labial cervical loop. Meis1 is specifically expressed by Sox2⁺ stem cells, which give rise to all dental epithelial cell lineages. Also, we have found that Meis1 in the incisor is coexpressed with potential binding partner Pbx1 during both embryonic and adult stages. Interestingly, Meis2 is present in different areas of the forming tooth and it is not expressed by dental epithelial stem cells, suggesting different roles for these two largely homologous genes. Additionally, we have established the expression patterns of Meis1 and Meis2 during tongue, hair, salivary gland and palate formation. Finally, analysis of Meis1-null allele mice indicated that, similarly, to SOX2, MEIS1 is not essential for tooth initiation, but might have a role during adult incisor renewal.

HOXA9 is a novel myopia risk gene

Purpose

A recent meta-analysis revealed PAX6 as a risk gene for myopia. There is a link between PAX6 and HOXA9. Furthermore, HOXA9 has been reported to activate TGF-β that is a risk factor for myopia. We speculate HOXA9 may participate in myopia development.

Methods

The Singapore GUSTO birth cohort provides data on children’s cycloplegic refraction measured at age of 3 years and their methylation profile based on the umbilical cord DNA. The HOXA9 expression levels were measured in the eyes of mono-ocular form deprivation myopia in mice. The plasmid with the mouse HOXA9 cDNA was constructed and then transfected to mouse primary retinal pigment epithelial (RPE) cells. The expression levels of myopia-related genes and cell proliferation were measured in the HOXA9-overexpressed RPE cells.

Results

A total of 519 children had data on methylation profile and cycloplegic refraction. The mean spherical equivalent refraction (SE) was 0.90D. Among 8 SE outliers (worse than -2D), 7 children had HOXA9 hypomethylation. The HOXA9 levels in the retina of myopic eyes was 2.65-fold (p = 0.029; paired t-test) higher than the uncovered fellow eyes. When HOXA9 was over-expressed in the RPE cells, TGF-β, MMP2, FGF2 and IGF1R expression levels were dose-dependently increased by HOXA9. However, over-expression of HOXA9 had no significant influence on IGF1 or HGF expression. In addition, HOXA9 also increased RPE proliferation.

Conclusion

Based on the human, animal and cellular data, the transcription factor HOXA9 may promote the expression of pro-myopia genes and RPE proliferation, which eventually contribute to myopia development.

Promoter-enhancer looping and shadow enhancers of the mouse αA-crystallin locus

Gene regulation by enhancers is important for precise temporal and spatial gene expression. Enhancers can drive gene expression regardless of their location, orientation or distance from the promoter. Changes in chromatin conformation and chromatin looping occur to bring the promoter and enhancers into close proximity. αA-crystallin ranks among one of the most abundantly expressed genes and proteins in the mammalian lens. The αA-crystallin locus is characterized by a 16 kb chromatin domain marked by two distal enhancers, 5′ DCR1 and 3′ DCR3. Here we used chromatin conformation capture (3C) analysis and transgenic approaches to analyze temporal control of the mouse αA-crystallin gene. We find that DCR1 is necessary, but not sufficient alone to drive expression at E10.5 in the mouse lens pit. Chromatin looping revealed interaction between the promoter and the region 3′ to DCR1, identifying a novel enhancer region in the αA-crystallin locus. We determined that this novel enhancer region, DCR1S, recapitulates the temporal control by DCR1. Acting as shadow enhancers, DCR1 and DCR1S are able to control expression in the lens vesicle at E11.5. It remains to be elucidated however, which region of the αA-crystallin locus is responsible for expression in the lens pit at E10.5.

Summary: The αA-crystallin ranks amongst the most highly expressed tissue-specific genes. It is an advantageous model system to probe both promoter-enhancer looping and to identify distal enhancers and their temporal/spatial activities.

MEIS-1 level in unresectable hepatocellular carcinoma can predict the post-treatment outcomes of radiofrequency ablation

Radiofrequency ablation (RFA) is a local-ablative therapy for unresectable hepatocellular carcinoma (HCC). At present, there is no predictive marker for RFA treatment outcomes. This work aimed to valuate myeloid ecotropic viral integration site 1 (MEIS-1) in predicting post-RFA treatment outcomes of unresectable HCC patients. The time to progression (TTP) and overall survival (OS) of 81 HCC patients who received RFA treatment were measured. The protein level of MEIS-1 in tumor specimens was measured by western blot. The role of MEIS-1 in RFA-treating HCC in vivo growth nude mouse model was examined via PET/CT imaging. Higher level of MEIS-1 in tumor tissue is associated with better RFA treatment outcomes. The median TTP was 9.0 (95% confidence interval [CI]: 6.8-11.3) months in patients with high MEIS-1 expression (n = 43) versus 6.0 (95% CI: 4.6-7.4) months in patients with low MEIS-1 expression (n = 38). Moreover, in rodent HCC model we found overexpression of MEIS-1 enhanced the anti-tumor effect of RFA treatment. We conclude that high level of MEIS-1 expression predicts better RFA treatment outcome in HCC.

Enhancer redundancy provides phenotypic robustness in mammalian development

Distant-acting tissue-specific enhancers vastly outnumber protein-coding genes in mammalian genomes, but the functional significance of this regulatory complexity remains insufficiently understood^1,2. Here we show that the pervasive presence of multiple enhancers with similar activities near the same gene confers phenotypic robustness to loss-of-function mutations in individual enhancers. We used genome editing to create 23 mouse deletion lines and inter-crosses, including both single and combinatorial enhancer deletions at seven distinct loci required for limb development. Surprisingly, none of ten deletions of individual enhancers caused noticeable changes in limb morphology. In contrast, removal of pairs of limb enhancers near the same gene resulted in discernible phenotypes, indicating that enhancers function redundantly in establishing normal morphology. In a genetic background sensitized by reduced baseline expression of the target gene, even single enhancer deletions caused limb abnormalities, suggesting that functional redundancy is conferred by additive effects of enhancers on gene expression levels. A genome-wide analysis integrating epigenomic and transcriptomic data from 29 developmental mouse tissues revealed that mammalian genes are very commonly associated with multiple enhancers that have similar spatiotemporal activity. Systematic exploration of three representative developmental structures (limb, brain, heart) uncovered more than a thousand cases in which five or more enhancers with redundant activity patterns were found near the same gene. Taken together, our data indicate that enhancer redundancy is a remarkably widespread feature of mammalian genomes and provides an effective regulatory buffer preventing deleterious phenotypic consequences upon loss of individual enhancers.

Enhancer adoption caused by genomic insertion elicits interdigital Shh expression and syndactyly in mouse

In this study, we reexamined an old mouse mutant named Hammer toe (Hm), which arose spontaneously almost a half century ago and exhibits a limb phenotype with webbing. We revealed that a 150-kb noncoding genomic fragment that was originally located in chromosome 14 has been inserted into a genomic region proximal to Sonic hedgehog (Shh), located in chromosome 5. This inserted fragment possesses enhancer activity to induce Shh expression in the interdigital regions in Hm, which in turn down-regulates bone morphogenetic protein signaling and eventually results in syndactyly and web formation. Since the donor fragment residing in chromosome 14 has enhancer activity to induce interdigital gene expression, the Hm mutation appears to be an archetypal case of enhancer adoption.

Acquisition of new cis-regulatory elements (CREs) can cause alteration of developmental gene regulation and may introduce morphological novelty in evolution. Although structural variation in the genome generated by chromosomal rearrangement is one possible source of new CREs, only a few examples are known, except for cases of retrotransposition. In this study, we show the acquisition of novel regulatory sequences as a result of large genomic insertion in the spontaneous mouse mutation Hammer toe (Hm). Hm mice exhibit syndactyly with webbing, due to suppression of interdigital cell death in limb development. We reveal that, in the Hm genome, a 150-kb noncoding DNA fragment from chromosome 14 is inserted into the region upstream of the Sonic hedgehog (Shh) promoter in chromosome 5. Phenotyping of mouse embryos with a series of CRISPR/Cas9-aided partial deletion of the 150-kb insert clearly indicated that two different regions are necessary for the syndactyly phenotype of Hm. We found that each of the two regions contains at least one enhancer for interdigital regulation. These results show that a set of enhancers brought by the large genomic insertion elicits the interdigital Shh expression and the Hm phenotype. Transcriptome analysis indicates that ectopic expression of Shh up-regulates Chordin (Chrd) that antagonizes bone morphogenetic protein signaling in the interdigital region. Indeed, Chrd-overexpressing transgenic mice recapitulated syndactyly with webbing. Thus, the Hm mutation provides an insight into enhancer acquisition as a source of creation of novel gene regulation.

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.