Enhancer pleiotropy, gene expression, and the architecture of human enhancer-gene interactions

Genome-Wide Silencer Screening Reveals Key Silencer Modulating Reprogramming Efficiency in Mouse Induced Pluripotent Stem Cells

The majority of the mouse genome is composed of non-coding regions, which harbor numerous regulatory sequences essential for gene regulation. While extensive research focuses on enhancers that activate gene expression, the role of silencers that repress gene expression remains less explored. In this study, the first genome-wide identification of silencers in the mouse genome is conducted. In mouse embryonic fibroblasts (MEFs) and embryonic stem cells (mESCs), 89 596 and 115 165 silencers are identified, respectively. These silencers are ubiquitously distributed across the genome and are predominantly associated with low-expression genes. Additionally, these silencers are mainly cell-specific and function by binding to repressive transcription factors (TFs). Further, these silencers are notably enriched with the histone modification H3K9me3. It is observed that the transformation between dual-function silencers and enhancers is correlated with intracellular transcription factor concentrations, accompanied by changes in epigenetic modifications. In terms of biological effects, we have identified silencers that can enhance the induction efficiency of MEFs and influence the pluripotency of mESCs. Collectively, this work offers the first comprehensive silencer landscape in the mouse genome and provides strong evidence for the role of silencers in the induction of induced pluripotent stem cells (iPSCs).

Pleiotropic expression quantitative trait loci are enriched in enhancers and transcription factor binding sites and impact more genes

Integrating expression quantitative trait loci (eQTL) data with genome-wide association studies (GWAS) enables the discovery of pleiotropic gene regulatory variants that influence a wide range of traits and disease susceptibilities. However, a comprehensive understanding of the distribution of pleiotropic QTLs across the genome and their phenotypic associations remain limited. In this study, we systematically annotated genetic variants associated with both trait variation and gene expression changes, focusing specifically on the unique characteristics of pleiotropic eQTLs. By integrating data from 127 eQTL studies and 417 traits from the IEU Open GWAS Project, we identified 476 pleiotropic eQTL variants affecting two or more distinct traits. Our analysis highlighted 5345 eQTL candidates potentially linked to gene expression changes across 293 GWAS traits. Notably, the 476 pleiotropic eQTLs associated with multiple trait categories were localized within a cumulative 2.5 Mbp genomic region. These pleiotropic eQTLs were enriched in enhancer regions and CTCF loops, influencing a larger number of genes in closer genomic proximity. Our findings reveal that pleiotropic eQTLs are concentrated within a small fraction of the genome and exhibit distinct molecular features. Colocalization results are accessible through an interactive web application and UCSC genome browser tracks at https://gwas2eqtl.tagc.univ-amu.fr/gwas2eqtl, facilitating the exploration of pleiotropic eQTLs and their roles in gene regulation and disease susceptibility.

MicroRNA-mediated network redundancy is constrained by purifying selection and contributes to expression robustness in Drosophila melanogaster

MicroRNAs (miRNAs) are post-transcriptional, non-coding regulatory RNAs that function coordinately with transcription factors (TFs) in gene regulatory networks. TFs and their targets are often co-regulated by miRNAs, forming composite feedforward circuits (cFFCs) with varying degrees of redundancy, primarily mediated by miRNAs. However, the maintenance of miRNA-mediated regulatory redundancy and its impact on gene expression evolution remain elusive. By integrating ChIP-seq data from ENCODE and miRNA targeting from TargetScanFly, we quantified miRNA-mediated cFFC redundancy in Drosophila melanogaster embryos and larvae, revealing more than three quarters of miRNA targets are involved in redundant cFFCs. Higher cFFC redundancy, where more miRNAs target the same gene within a cFFC, is correlated with stronger purifying selection, reduced expression divergence between species, and increased expression stability under heat shock stress. Redundant cFFCs primarily regulate older or broadly expressed young genes. These findings highlight the role of miRNA-mediated cFFC redundancy in enhancing gene expression robustness through natural selection.

Characterization of Single-Cell Cis-regulatory Elements Informs Implications for Cell Differentiation

Cis-regulatory elements govern the specific patterns and dynamics of gene expression in cells during development, which are the fundamental mechanisms behind cell differentiation. However, the genomic characteristics of single-cell cis-regulatory elements closely linked to cell differentiation during development remain unclear. To explore this, we systematically analyzed ∼250,000 putative single-cell cis-regulatory elements obtained from snATAC-seq analysis of the developing mouse cerebellum. We found that over 80% of these single-cell cis-regulatory elements show pleiotropic effects, being active in 2 or more cell types. The pleiotropic degrees of proximal and distal single-cell cis-regulatory elements are positively correlated with the density and diversity of transcription factor binding motifs and GC content. There is a negative correlation between the pleiotropic degrees of single-cell cis-regulatory elements and their distances to the nearest transcription start sites, and proximal single-cell cis-regulatory elements display higher relevance strengths than distal ones. Furthermore, both proximal and distal single-cell cis-regulatory elements related to cell differentiation exhibit enhanced sequence-level evolutionary conservation, increased density and diversity of transcription factor binding motifs, elevated GC content, and greater distances from their nearest genes. Together, our findings reveal the general genomic characteristics of putative single-cell cis-regulatory elements and provide insights into the genomic and evolutionary mechanisms by which single-cell cis-regulatory elements regulate cell differentiation during development.

Evidence for compensatory evolution within pleiotropic regulatory elements

Pleiotropy, measured as expression breadth across tissues, is one of the best predictors for protein sequence and expression conservation. In this study, we investigated its effect on the evolution of cis-regulatory elements (CREs). To this end, we carefully reanalyzed the Epigenomics Roadmap data for nine fetal tissues, assigning a measure of pleiotropic degree to nearly half a million CREs. To assess the functional conservation of CREs, we generated ATAC-seq and RNA-seq data from humans and macaques. We found that more pleiotropic CREs exhibit greater conservation in accessibility, and the mRNA expression levels of the associated genes are more conserved. This trend of higher conservation for higher degrees of pleiotropy persists when analyzing the transcription factor binding repertoire. In contrast, simple DNA sequence conservation of orthologous sites between species tends to be even lower for pleiotropic CREs than for species-specific CREs. Combining various lines of evidence, we propose that the lack of sequence conservation in functionally conserved pleiotropic CREs is owing to within-element compensatory evolution. In summary, our findings suggest that pleiotropy is also a good predictor for the functional conservation of CREs, even though this is not reflected in the sequence conservation of pleiotropic CREs.

Enhancer-mediated FOXO3 expression promotes MSC adipogenic differentiation by activating autophagy

BACKGROUND

Mesenchymal stem cells (MSCs) are pluripotent stem cells capable of differentiating into osteocytes, adipocytes and chondrocytes. However, in osteoporosis, the balance of differentiation is tipped toward adipogenesis and the key mechanism is controversial. Researches have shown that, as upstream regulatory elements of gene expression, enhancers ar involved in the expression of identity genes. In this study, we identified enhancers-mediated gene FOXO3 promoting MSC adipogenic differentiation by activating autophagy.

METHODS

We integrated data of RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq) and ATAC-sequencing (ATAC-seq) to find the identity gene FOXO3. The expression of FOXO3 protein, adipogenic transcription factors and the substrate of autophagy were measured by western blotting. The Oil Red O (ORO) staining was used to visualize the adipogenesis of MSCs. Immunohistochemistry was used to visualize the FOXO3 expression in adipocytes in bone marrow. Immunofluorescence was used to detect the expression of PPARγ and LC3B.

RESULTS

During adipogenesis, enhancers redistribute to genes associated with adipogenic differentiation, among which we identified the pivotal identity gene FOXO3. FOXO3 could promote the expression of the adipogenic transcription factors PPARγ, CEBPα, and CEBPβ during adipogenic differentiation, while PPARγ, CEBPα, and CEBPβ could in turn bind to FOXO3 and continue to promote FOXO3 expression to form a positive feedback loop. Consistently elevated FOXO3 expression promotes autophagy by activating the PI3K-AKT pathway which mediates adipogenic differentiation.

CONCLUSIONS

Pivotal identity gene FOXO3 promotes autophagy by activating PI3K-AKT pathway, which provokes adipogenic differentiation of MSCs. Enhancer-regulated adipogenic identity gene FOXO3 could be an attractive treatment for osteoporosis.

A body map of super-enhancers and their function in pig

Introduction

Super-enhancers (SEs) are clusters of enhancers that act synergistically to drive the high-level expression of genes involved in cell identity and function. Although SEs have been extensively investigated in humans and mice, they have not been well characterized in pigs.

Methods

Here, we identified 42,380 SEs in 14 pig tissues using chromatin immunoprecipitation sequencing, and statistics of its overall situation, studied the composition and characteristics of SE, and explored the influence of SEs characteristics on gene expression.

Results

We observed that approximately 40% of normal enhancers (NEs) form SEs. Compared to NEs, we found that SEs were more likely to be enriched with an activated enhancer and show activated functions. Interestingly, SEs showed X chromosome depletion and short interspersed nuclear element enrichment, implying that SEs play an important role in sex traits and repeat evolution. Additionally, SE-associated genes exhibited higher expression levels and stronger conservation than NE-associated genes. However, genes with the largest SEs had higher expression levels than those with the smallest SEs, indicating that SE size may influence gene expression. Moreover, we observed a negative correlation between SE gene distance and gene expression, indicating that the proximity of SEs can affect gene activity. Gene ontology enrichment and motif analysis revealed that SEs have strong tissue-specific activity. For example, the CORO2B gene with a brain-specific SE shows strong brain-specific expression, and the phenylalanine hydroxylase gene with liver-specific SEs shows strong liver-specific expression.

Discussion

In this study, we illustrated a body map of SEs and explored their functions in pigs, providing information on the composition and tissue-specific patterns of SEs. This study can serve as a valuable resource of gene regulatory and comparative analyses to the scientific community and provides a theoretical reference for genetic control mechanisms of important traits in pigs.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data

MOTIVATION

Identifying regulatory regions in the genome is of great interest for understanding the epigenomic landscape in cells. One fundamental challenge in this context is to find the target genes whose expression is affected by the regulatory regions. A recent successful method is the Activity-By-Contact (ABC) model which scores enhancer-gene interactions based on enhancer activity and the contact frequency of an enhancer to its target gene. However, it describes regulatory interactions entirely from a gene's perspective, and does not account for all the candidate target genes of an enhancer. In addition, the ABC model requires two types of assays to measure enhancer activity, which limits the applicability. Moreover, there is neither implementation available that could allow for an integration with transcription factor (TF) binding information nor an efficient analysis of single-cell data.

RESULTS

We demonstrate that the ABC score can yield a higher accuracy by adapting the enhancer activity according to the number of contacts the enhancer has to its candidate target genes and also by considering all annotated transcription start sites of a gene. Further, we show that the model is comparably accurate with only one assay to measure enhancer activity. We combined our generalized ABC model with TF binding information and illustrated an analysis of a single-cell ATAC-seq dataset of the human heart, where we were able to characterize cell type-specific regulatory interactions and predict gene expression based on TF affinities. All executed processing steps are incorporated into our new computational pipeline STARE.

AVAILABILITY AND IMPLEMENTATION

The software is available at https://github.com/schulzlab/STARE.

CONTACT

marcel.schulz@em.uni-frankfurt.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Site pleiotropy of a stickleback Bmp6 enhancer

Development and regeneration are orchestrated by gene regulatory networks that operate in part through transcriptional enhancers. Although many enhancers are pleiotropic and are active in multiple tissues, little is known about whether enhancer pleiotropy is due to 1) site pleiotropy, in which individual transcription factor binding sites (TFBS) are required for activity in multiple tissues, or 2) multiple distinct sites that regulate expression in different tissues. Here, we investigated the pleiotropy of an intronic enhancer of the stickleback Bone morphogenetic protein 6 (Bmp6) gene. This enhancer was previously shown to regulate evolved changes in tooth number and tooth regeneration, and is highly pleiotropic, with robust activity in both fins and teeth throughout embryonic, larval, and adult life, and in the heart and kidney in adult fish. We tested the hypothesis that the pleiotropy of this enhancer is due to site pleiotropy of an evolutionarily conserved predicted Foxc1 TFBS. Transgenic analysis and site-directed mutagenesis experiments both deleting and scrambling this predicted Foxc1 TFBS revealed that the binding site is required for enhancer activity in both teeth and fins throughout embryonic, larval, and adult development, and in the heart and kidney in adult fish. Collectively these data support a model where the pleiotropy of this Bmp6 enhancer is due to site pleiotropy and this putative binding site is required for enhancer activity in multiple anatomical sites from the embryo to the adult.

Evolutionary Changes in the Chromatin Landscape Contribute to Reorganization of a Developmental Gene Network During Rapid Life History Evolution in Sea Urchins

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

Temporal analysis of enhancers during mouse cerebellar development reveals dynamic and novel regulatory functions

We have identified active enhancers in the mouse cerebellum at embryonic and postnatal stages which provides a view of novel enhancers active during cerebellar development. The majority of cerebellar enhancers have dynamic activity between embryonic and postnatal development. Cerebellar enhancers were enriched for neural transcription factor binding sites with temporally specific expression. Putative gene targets displayed spatially restricted expression patterns, indicating cell-type specific expression regulation. Functional analysis of target genes indicated that enhancers regulate processes spanning several developmental epochs such as specification, differentiation and maturation. We use these analyses to discover one novel regulator and one novel marker of cerebellar development: Bhlhe22 and Pax3, respectively. We identified an enrichment of de novo mutations and variants associated with autism spectrum disorder in cerebellar enhancers. Furthermore, by comparing our data with relevant brain development ENCODE histone profiles and cerebellar single-cell datasets we have been able to generalize and expand on the presented analyses, respectively. We have made the results of our analyses available online in the Developing Mouse Cerebellum Enhancer Atlas (https://goldowitzlab.shinyapps.io/developing_mouse_cerebellum_enhancer_atlas/), where our dataset can be efficiently queried, curated and exported by the scientific community to facilitate future research efforts. Our study provides a valuable resource for studying the dynamics of gene expression regulation by enhancers in the developing cerebellum and delivers a rich dataset of novel gene-enhancer associations providing a basis for future in-depth studies in the cerebellum.

Pleiotropic Enhancers are Ubiquitous Regulatory Elements in the Human Genome

Enhancers are regulatory elements of genomes that determine spatio-temporal patterns of gene expression. The human genome contains a vast number of enhancers, which largely outnumber protein-coding genes. Historically, enhancers have been regarded as highly tissue-specific. However, recent evidence has demonstrated that many enhancers are pleiotropic, with activity in multiple developmental contexts. Yet, the extent and impact of pleiotropy remain largely unexplored. In this study we analyzed active enhancers across human organs based on the analysis of both eRNA transcription (FANTOM5 consortium data sets) and chromatin architecture (ENCODE consortium data sets). We show that pleiotropic enhancers are pervasive in the human genome and that most enhancers active in a particular organ are also active in other organs. In addition, our analysis suggests that the proportion of context-specific enhancers of a given organ is explained, at least in part, by the proportion of context-specific genes in that same organ. The notion that such a high proportion of human enhancers can be pleiotropic suggests that small regions of regulatory DNA contain abundant regulatory information and that these regions evolve under important evolutionary constraints.

Adaptive eQTLs reveal the evolutionary impacts of pleiotropy and tissue-specificity while contributing to health and disease

Large numbers of expression quantitative trait loci (eQTLs) have recently been identified in humans, and many of these regulatory variants have large allele frequency differences between populations. Here, we conducted genome-wide scans of selection to identify adaptive eQTLs (i.e., eQTLs with large population branch statistics). We then tested if tissue pleiotropy affects whether eQTLs are more or less likely to be adaptive and identified tissues that have been key targets of positive selection during the last 100,000 years. Top adaptive eQTL outliers include rs1043809, rs66899053, and rs2814778 (a SNP that is associated with malaria resistance). We found that effect sizes of eQTLs were negatively correlated with population branch statistics and that adaptive eQTLs affect two-thirds as many tissues as do non-adaptive eQTLs. Because the tissue breadth of an eQTL can be viewed as a measure of pleiotropy, these results imply that pleiotropy inhibits adaptation. The proportion of eQTLs that are adaptive varies by tissue, and we found that eQTLs that regulate expression in testis, thyroid, blood, or sun-exposed skin are enriched for signatures of positive selection. By contrast, eQTLs that regulate expression in the cerebrum or female-specific tissues have a relative lack of adaptive outliers. Scans of selections also reveal that many adaptive eQTLs are closely linked to disease-associated loci. Taken together, our results indicate that eQTLs have played an important role in recent human evolution.