Developmental Enhancers Research Articles

VISTA Enhancer browser: an updated database of tissue-specific developmental enhancers.

Regulatory elements (enhancers) are major drivers of gene expression in mammals and harbor many genetic variants associated with human diseases. Here, we present an updated VISTA Enhancer Browser (https://enhancer.lbl.gov), a database of transgenic enhancer assays conducted in developing mouse embryos in vivo. Since the original publication in 2007, the database grew nearly 20-fold from 250 to over 4500 experiments and currently harbors over 23500 images. The updated database provides structured information on experiments conducted at different stages of embryonic development, including enhancer activities of human pathogenic and synthetic variants and sequences derived from a variety of species. In addition to manually curated results of thousands of individual experiments, the new database also features hundreds of manually curated comparisons between alleles. The VISTA Enhancer Browser provides a crucial resource for study of human genetic variation, gene regulation and developmental biology.

Developmental and housekeeping transcriptional programs display distinct modes of enhancer-enhancer cooperativity in Drosophila

Genomic enhancers are key transcriptional regulators which, upon the binding of sequence-specific transcription factors, activate their cognate target promoters. Although enhancers have been extensively studied in isolation, a substantial number of genes have more than one simultaneously active enhancer, and it remains unclear how these cooperate to regulate transcription. Using Drosophila melanogaster S2 cells as a model, we assay the activities of more than a thousand individual enhancers and about a million enhancer pairs toward housekeeping and developmental core promoters with STARR-seq. We report that housekeeping and developmental enhancers show distinct modes of enhancer-enhancer cooperativity: while housekeeping enhancers are additive such that their combined activity mirrors the sum of their individual activities, developmental enhancers are super-additive and combine multiplicatively. Super-additivity between developmental enhancers is promiscuous and neither depends on the enhancers’ endogenous genomic contexts nor on specific transcription factor motif signatures. However, it can be further boosted by Twist and Trl motifs and saturates for the highest levels of enhancer activity. These results have important implications for our understanding of gene regulation in complex multi-enhancer developmental loci and genomically clustered housekeeping genes, providing a rationale to interpret the transcriptional impact of non-coding mutations at different loci.

A suboptimal OCT4-SOX2 binding site facilitates the naïve-state specific function of a Klf4 enhancer.

Enhancers have critical functions in the precise, spatiotemporal control of transcription during development. It is thought that enhancer grammar, or the characteristics and arrangements of transcription factor binding sites, underlie the specific functions of developmental enhancers. In this study, we sought to identify grammatical constraints that direct enhancer activity in the naïve state of pluripotency, focusing on the enhancers for the naïve-state specific gene, Klf4. Using a combination of biochemical tests, reporter assays, and endogenous mutations in mouse embryonic stem cells, we have studied the binding sites for the transcription factors OCT4 and SOX2. We have found that the three Klf4 enhancers contain suboptimal OCT4-SOX2 composite binding sites. Substitution with a high-affinity OCT4-SOX2 binding site in Klf4 enhancer E2 rescued enhancer function and Klf4 expression upon loss of the ESRRB and STAT3 binding sites. We also observed that the low-affinity of the OCT4-SOX2 binding site is crucial to drive the naïve-state specific activities of Klf4 enhancer E2. Altogether, our work suggests that the affinity of OCT4-SOX2 binding sites could facilitate enhancer functions in specific states of pluripotency.

Mutagenesis Sensitivity Mapping of Human Enhancers In Vivo.

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

FOXA2 rewires AP-1 for transcriptional reprogramming and lineage plasticity in prostate cancer

FOXA family proteins act as pioneer factors by remodeling compact chromatin structures. FOXA1 is crucial for the chromatin binding of the androgen receptor (AR) in both normal prostate epithelial cells and the luminal subtype of prostate cancer (PCa). Recent studies have highlighted the emergence of FOXA2 as an adaptive response to AR signaling inhibition treatments. However, the role of the FOXA1 to FOXA2 transition in regulating cancer lineage plasticity remains unclear. Our study demonstrates that FOXA2 binds to distinct classes of developmental enhancers in multiple AR-independent PCa subtypes, with its binding depending on LSD1. Moreover, we reveal that FOXA2 collaborates with JUN at chromatin and promotes transcriptional reprogramming of AP-1 in lineage-plastic cancer cells, thereby facilitating cell state transitions to multiple lineages. Overall, our findings underscore the pivotal role of FOXA2 as a pan-plasticity driver that rewires AP-1 to induce the differential transcriptional reprogramming necessary for cancer cell lineage plasticity.

Increased enhancer-promoter interactions during developmental enhancer activation in mammals.

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

The Binding Occupancy of Pro-Hypertrophic Transcription Factors among Cardiac Developmental Enhancers

Objective: To check which of the eight pro-hypertrophic transcription factors bind in abundance to a set of 140 functionally confirmed cardiac enhancers aiming to identify potential therapeutic targets for delaying the onset of cardiac hypertrophy. Study Design: Bioinformatics, in-Silico Study.Place and Duration of Study: The study was conducted at the Department of Biological Sciences, National University of Medical Sciences (NUMS), Rawalpindi, Pakistan, from January 2022 to December 2022.Methods: For experimentally verified human and mouse noncoding regions with gene enhancer activity measured in transgenic mice, VISTA Enhancer Browser was used to obtain 140 cardiac-specific functionally confirmed enhancers. Among these enhancers, transcription factor binding sites were computationally mapped to assess transcription factor occupancy. The binding sites for all eight transcription factors were screened in the collective set of 140 cardiac enhancers via an in-house script. Results: In this study, among a small set of eight pro-hypertrophic transcription factors, occupancy of transcription factor HAND (Heart- And Neural Crest Derivatives-Expressed Protein) was found to be maximum in the collected pool of functionally confirmed cardiac enhancers. Conclusion: The transcription factor HAND was found to have maximum binding occupancy among a set of experimentally confirmed cardiac enhancers.How to cite this: Raza RZ, Hussain I, Akhtar F, Mumtaz S. The Binding Occupancy of Pro-Hypertrophic Transcription Factors among Cardiac Developmental Enhancers. Life and Science. 2024; 5(1): 22-26. doi: http://doi.org/10.37185/LnS.1.1.464

Epigenetic reprogramming of a distal developmental enhancer cluster drives SOX2 overexpression in breast and lung adenocarcinoma.

Enhancer reprogramming has been proposed as a key source of transcriptional dysregulation during tumorigenesis, but the molecular mechanisms underlying this process remain unclear. Here, we identify an enhancer cluster required for normal development that is aberrantly activated in breast and lung adenocarcinoma. Deletion of the SRR124-134 cluster disrupts expression of the SOX2 oncogene, dysregulates genome-wide transcription and chromatin accessibility and reduces the ability of cancer cells to form colonies in vitro. Analysis of primary tumors reveals a correlation between chromatin accessibility at this cluster and SOX2 overexpression in breast and lung cancer patients. We demonstrate that FOXA1 is an activator and NFIB is a repressor of SRR124-134 activity and SOX2 transcription in cancer cells, revealing a co-opting of the regulatory mechanisms involved in early development. Notably, we show that the conserved SRR124 and SRR134 regions are essential during mouse development, where homozygous deletion results in the lethal failure of esophageal-tracheal separation. These findings provide insights into how developmental enhancers can be reprogrammed during tumorigenesis and underscore the importance of understanding enhancer dynamics during development and disease.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/− 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Breaking enhancers to gain insights into developmental defects.

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

Mutational scans reveal differential evolvability of Drosophila promoters and enhancers.

Rapid enhancer and slow promoter evolution have been demonstrated through comparative genomics. However, it is not clear how this information is encoded genetically and if this can be used to place evolution in a predictive context. Part of the challenge is that our understanding of the potential for regulatory evolution is biased primarily toward natural variation or limited experimental perturbations. Here, to explore the evolutionary capacity of promoter variation, we surveyed an unbiased mutation library for three promoters in Drosophila melanogaster. We found that mutations in promoters had limited to no effect on spatial patterns of gene expression. Compared to developmental enhancers, promoters are more robust to mutations and have more access to mutations that can increase gene expression, suggesting that their low activity might be a result of selection. Consistent with these observations, increasing the promoter activity at the endogenous locus of shavenbaby led to increased transcription yet limited phenotypic changes. Taken together, developmental promoters may encode robust transcriptional outputs allowing evolvability through the integration of diverse developmental enhancers. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Interdisciplinary approaches to predicting evolutionary biology.

Interdisciplinary approaches to predicting evolutionary biology.

Dynamic interplay between non-coding enhancer transcription and gene activity in development

Non-coding transcription at the intergenic regulatory regions is a prevalent feature of metazoan genomes, but its biological function remains uncertain. Here, we devise a live-imaging system that permits simultaneous visualization of gene activity along with intergenic non-coding transcription at single-cell resolution in Drosophila. Quantitative image analysis reveals that elongation of RNA polymerase II across the internal core region of enhancers leads to suppression of transcriptional bursting from linked genes. Super-resolution imaging and genome-editing analysis further demonstrate that enhancer transcription antagonizes molecular crowding of transcription factors, thereby interrupting the formation of a transcription hub at the gene locus. We also show that a certain class of developmental enhancers are structurally optimized to co-activate gene transcription together with non-coding transcription effectively. We suggest that enhancer function is flexibly tunable through the modulation of hub formation via surrounding non-coding transcription during development.

Enhancer architecture and chromatin accessibility constrain phenotypic space during Drosophila development.

Developmental enhancers bind transcription factors and dictate patterns of gene expression during development. Their molecular evolution can underlie phenotypical evolution, but the contributions of the evolutionary pathways involved remain little understood. Here, using mutation libraries in Drosophila melanogaster embryos, we observed that most point mutations in developmental enhancers led to changes in gene expression levels but rarely resulted in novel expression outside of the native pattern. In contrast, random sequences, often acting as developmental enhancers, drove expression across a range of cell types; random sequences including motifs for transcription factors with pioneer activity acted as enhancers even more frequently. Our findings suggest that the phenotypic landscapes of developmental enhancers are constrained by enhancer architecture and chromatin accessibility. We propose that the evolution of existing enhancers is limited in its capacity to generate novel phenotypes, whereas the activity of de novo elements is a primary source of phenotypic novelty.

Active DNA demethylation of developmental cis-regulatory regions predates vertebrate origins.

DNA methylation [5-methylcytosine (5mC)] is a repressive gene-regulatory mark required for vertebrate embryogenesis. Genomic 5mC is tightly regulated through the action of DNA methyltransferases, which deposit 5mC, and ten-eleven translocation (TET) enzymes, which participate in its active removal through the formation of 5-hydroxymethylcytosine (5hmC). TET enzymes are essential for mammalian gastrulation and activation of vertebrate developmental enhancers; however, to date, a clear picture of 5hmC function, abundance, and genomic distribution in nonvertebrate lineages is lacking. By using base-resolution 5mC and 5hmC quantification during sea urchin and lancelet embryogenesis, we shed light on the roles of nonvertebrate 5hmC and TET enzymes. We find that these invertebrate deuterostomes use TET enzymes for targeted demethylation of regulatory regions associated with developmental genes and show that the complement of identified 5hmC-regulated genes is conserved to vertebrates. This work demonstrates that active 5mC removal from regulatory regions is a common feature of deuterostome embryogenesis suggestive of an unexpected deep conservation of a major gene-regulatory module.

Integration of multimodal data in the developing tooth reveals candidate regulatory loci driving human odontogenic phenotypes.

Human odontogenic aberrations such as abnormal tooth number and delayed tooth eruption can occur as a symptom of rare syndromes or, more commonly, as nonsyndromic phenotypes. These phenotypes can require extensive and expensive dental treatment, posing a significant burden. While many dental phenotypes are heritable, most nonsyndromic cases have not been linked to causal genes. We demonstrate the novel finding that common sequence variants associated with human odontogenic phenotypes are enriched in developmental craniofacial enhancers conserved between human and mouse. However, the bulk nature of these samples obscures if this finding is due to the tooth itself or the surrounding tissues. We therefore sought to identify enhancers specifically active in the tooth anlagen and quantify their contribution to the observed genetic enrichments. We systematically identified 22,001 conserved enhancers active in E13.5 mouse incisors using ChIP-seq and machine learning pipelines and demonstrated biologically relevant enrichments in putative target genes, transcription factor binding motifs, and in vivo activity. Multi-tissue comparisons of human and mouse enhancers revealed that these putative tooth enhancers had the strongest enrichment of odontogenic phenotype-associated variants, suggesting a role for dysregulation of tooth developmental enhancers in human dental phenotypes. The large number of these regions genome-wide necessitated prioritization of enhancer loci for future investigations. As enhancers modulate gene expression, we prioritized regions based on enhancers' putative target genes. We predicted these target genes and prioritized loci by integrating chromatin state, bulk gene expression and coexpression, GWAS variants, and cell type resolved gene expression to generate a prioritized list of putative odontogenic phenotype-driving loci active in the developing tooth. These genomic regions are of particular interest for downstream experiments determining the role of specific dental enhancer:gene pairs in odontogenesis.

Sequence divergence and retrotransposon insertion underlie interspecific epigenetic differences in primates.

Changes in the epigenome can affect the phenotype without the presence of changes in the genomic sequence. Given the high identity of the human and chimpanzee genome sequences, a substantial portion of their phenotypic divergence likely arises from epigenomic differences between the two species. In this study, the transcriptome and epigenome were determined for induced pluripotent stem cells (iPSCs) generated from human and chimpanzee individuals. The transcriptome and epigenomes for trimethylated histone H3 at lysine-4 (H3K4me3) and at lysine-27 (H3K27me3) showed high levels of similarity between the two species. However, there were some differences in histone modifications. Although such regions, in general, did not show significant enrichment of interspecies nucleotide variations, gains in binding motifs for pluripotency-related transcription factors, especially POU5F1 and SOX2, were frequently found in species-specific H3K4me3 regions. We also revealed that species-specific insertions of retrotransposons, including the LTR5_Hs subfamily in human and a newly identified LTR5_Pt subfamily in chimpanzee, created species-specific H3K4me3 regions associated with increased expression of nearby genes. Human iPSCs have more species-specific H3K27me3 regions, resulting in more abundant bivalent domains. Only a limited number of these species-specific H3K4me3 and H3K27me3 regions overlap with species-biased enhancers in cranial neural crest cells, suggesting that differences in the epigenetic state of developmental enhancers appear late in development. Therefore, iPSCs serve as a suitable starting material for studying evolutionary changes in epigenome dynamics during development.

Regeneration and developmental enhancers are differentially compatible with minimal promoters

Enhancers and promoters are cis-regulatory elements that control gene expression. Enhancers are activated in a cell type-, tissue-, and condition-specific manner to stimulate promoter function and transcription. Zebrafish have emerged as a powerful animal model for examining the activities of enhancers derived from various species through transgenic enhancer assays, in which an enhancer is coupled with a minimal promoter. However, the efficiency of minimal promoters and their compatibility with multiple developmental and regeneration enhancers have not been systematically tested in zebrafish. Thus, we assessed the efficiency of six minimal promoters and comprehensively interrogated the compatibility of the promoters with developmental and regeneration enhancers. We found that the fos minimal promoter and Drosophila synthetic core promoter (DSCP) yielded high rates of leaky expression that may complicate the interpretation of enhancer assays. Notably, the adenovirus E1b promoter, the zebrafish lepb 0.8-kb (P0.8) and lepb 2-kb (P2) promoters, and a new zebrafish synthetic promoter (ZSP) that combines elements of the E1b and P0.8 promoters drove little or no ectopic expression, making them suitable for transgenic assays. We also found significant differences in compatibility among specific combinations of promoters and enhancers, indicating the importance of promoters as key regulatory elements determining the specificity of gene expression. Our study provides guidelines for transgenic enhancer assays in zebrafish to aid in the discovery of functional enhancers regulating development and regeneration.

DeepSTARR predicts enhancer activity from DNA sequence and enables the de novo design of synthetic enhancers.

Enhancer sequences control gene expression and comprise binding sites (motifs) for different transcription factors (TFs). Despite extensive genetic and computational studies, the relationship between DNA sequence and regulatory activity is poorly understood, and de novo enhancer design has been challenging. Here, we built a deep-learning model, DeepSTARR, to quantitatively predict the activities of thousands of developmental and housekeeping enhancers directly from DNA sequence in Drosophila melanogaster S2 cells. The model learned relevant TF motifs and higher-order syntax rules, including functionally nonequivalent instances of the same TF motif that are determined by motif-flanking sequence and intermotif distances. We validated these rules experimentally and demonstrated that they can be generalized to humans by testing more than 40,000 wildtype and mutant Drosophila and human enhancers. Finally, we designed and functionally validated synthetic enhancers with desired activities de novo.