Transcriptional enhancers and their communication with gene promoters

Physical and functional interaction among Irf8 enhancers during dendritic cell differentiation

The production of type 1 conventional dendritic cells (cDC1s) requires high expression of the transcription factor IRF8. Three enhancers at the Irf8 3' region function in a differentiation stage-specific manner. However, whether and how these enhancers interact physically and functionally remains unclear. Here, we show that the Irf8 3' enhancers directly interact with each other and contact the Irf8 gene body during cDC1 differentiation. The +56 kb enhancer, which functions from multipotent progenitor stages, activates the other 3' enhancers through an IRF8-dependent transcription factor program, that is, in trans. Then, the +32 kb enhancer, which operates in cDC1-committed cells, reversely acts in cis on the other 3' enhancers to maintain the high expression of Irf8. Indeed, mice with compound heterozygous deletion of the +56 and +32 kb enhancers are unable to generate cDC1s. These results illustrate how multiple enhancers cooperate to induce a lineage-determining transcription factor gene during cell differentiation.

Dose rate dependent reduction in chromatin accessibility at transcriptional start sites long time after exposure to gamma radiation

Ionizing radiation (IR) impact cellular and molecular processes that require chromatin remodelling relevant for cellular integrity. However, the cellular implications of ionizing radiation (IR) delivered per time unit (dose rate) are still debated. This study investigates whether the dose rate is relevant for inflicting changes to the epigenome, represented by chromatin accessibility, or whether it is the total dose that is decisive. CBA/CaOlaHsd mice were whole-body exposed to either chronic low dose rate (2.5 mGy/h for 54 d) or the higher dose rates (10 mGy/h for 14 d and 100 mGy/h for 30 h) of gamma radiation (⁶⁰Co, total dose: 3 Gy). Chromatin accessibility was analysed in liver tissue samples using Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-Seq), both one day after and over three months post-radiation (>100 d). The results show that the dose rate contributes to radiation-induced epigenomic changes in the liver at both sampling timepoints. Interestingly, chronic low dose rate exposure to a high total dose (3 Gy) did not inflict long-term changes to the epigenome. In contrast to the acute high dose rate given to the same total dose, reduced accessibility at transcriptional start sites (TSS) was identified in genes relevant for the DNA damage response and transcriptional activity. Our findings link dose rate to essential biological mechanisms that could be relevant for understanding long-term changes after ionizing radiation exposure. However, future studies are needed to comprehend the biological consequence of these findings.

Clustered and diverse transcription factor binding underlies cell type specificity of enhancers for housekeeping genes

Housekeeping genes are considered to be regulated by common enhancers across different tissues. Here we report that most of the commonly expressed mouse or human genes across different cell types, including more than half of the previously identified housekeeping genes, are associated with cell type-specific enhancers. Furthermore, the binding of most transcription factors (TFs) is cell type-specific. We reason that these cell type specificities are causally related to the collective TF recruitment at regulatory sites, as TFs tend to bind to regions associated with many other TFs and each cell type has a unique repertoire of expressed TFs. Based on binding profiles of hundreds of TFs from HepG2, K562, and GM12878 cells, we show that 80% of all TF peaks overlapping H3K27ac signals are in the top 20,000-23,000 most TF-enriched H3K27ac peak regions, and approximately 12,000-15,000 of these peaks are enhancers (nonpromoters). Those enhancers are mainly cell type-specific and include those linked to the majority of commonly expressed genes. Moreover, we show that the top 15,000 most TF-enriched regulatory sites in HepG2 cells, associated with about 200 TFs, can be predicted largely from the binding profile of as few as 30 TFs. Through motif analysis, we show that major enhancers harbor diverse and clustered motifs from a combination of available TFs uniquely present in each cell type. We propose a mechanism that explains how the highly focused TF binding at regulatory sites results in cell type specificity of enhancers for housekeeping and commonly expressed genes.

Boundary bypass activity in the abdominal-B region of the Drosophila bithorax complex is position dependent and regulated

Expression of Abdominal-B (Abd-B) in abdominal segments A5-A8 is controlled by four regulatory domains, iab-5-iab-8. Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements (Mcp, Fab-7, Fab-7 and Fab-8). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can 'jump over' intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Characterization of an eye field-like state during optic vesicle organoid development

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Multiple Fra-1-bound enhancers showing different molecular and functional features can cooperate to repress gene transcription

BACKGROUND

How transcription factors (TFs) down-regulate gene expression remains ill-understood, especially when they bind to multiple enhancers contacting the same gene promoter. In particular, it is not known whether they exert similar or significantly different molecular effects at these enhancers.

RESULTS

To address this issue, we used a particularly well-suited study model consisting of the down-regulation of the TGFB2 gene by the TF Fra-1 in Fra-1-overexpressing cancer cells, as Fra-1 binds to multiple enhancers interacting with the TGFB2 promoter. We show that Fra-1 does not repress TGFB2 transcription via reducing RNA Pol II recruitment at the gene promoter but by decreasing the formation of its transcription-initiating form. This is associated with complex long-range chromatin interactions implicating multiple molecularly and functionally heterogeneous Fra-1-bound transcriptional enhancers distal to the TGFB2 transcriptional start site. In particular, the latter display differential requirements upon the presence and the activity of the lysine acetyltransferase p300/CBP. Furthermore, the final transcriptional output of the TGFB2 gene seems to depend on a balance between the positive and negative effects of Fra-1 at these enhancers.

CONCLUSION

Our work unveils complex molecular mechanisms underlying the repressive actions of Fra-1 on TGFB2 gene expression. This has consequences for our general understanding of the functioning of the ubiquitous transcriptional complex AP-1, of which Fra-1 is the most documented component for prooncogenic activities. In addition, it raises the general question of the heterogeneity of the molecular functions of TFs binding to different enhancers regulating the same gene.

Defining the fine structure of promoter activity on a genome-wide scale with CISSECTOR

Classic promoter mutagenesis strategies can be used to study how proximal promoter regions regulate the expression of particular genes of interest. This is a laborious process, in which the smallest sub-region of the promoter still capable of recapitulating expression in an ectopic setting is first identified, followed by targeted mutation of putative transcription factor binding sites. Massively parallel reporter assays such as survey of regulatory elements (SuRE) provide an alternative way to study millions of promoter fragments in parallel. Here we show how a generalized linear model (GLM) can be used to transform genome-scale SuRE data into a high-resolution genomic track that quantifies the contribution of local sequence to promoter activity. This coefficient track helps identify regulatory elements and can be used to predict promoter activity of any sub-region in the genome. It thus allows in silico dissection of any promoter in the human genome to be performed. We developed a web application, available at cissector.nki.nl, that lets researchers easily perform this analysis as a starting point for their research into any promoter of interest.

Boundary Bypass Activity in the Abdominal-B Region of the Drosophila Bithorax Complex is Position Dependent and Regulated

Expression of Abdominal-B ( Abd-B ) in abdominal segments A5 - A8 is controlled by four regulatory domains, iab-5 - iab-8 . Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements ( Mcp , Fab-7 , Fab-7 and Fab-8 ). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can "jump over" intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Summary Statement

Boundaries separating Abd-B regulatory domains block crosstalk between domains and mediate their interactions with Abd-B . The latter function is location but not orientation dependent.

Strategies for activity analysis of single nucleotide polymorphisms associated with human diseases

Genome-wide association studies (GWAS) have identified a large number of single nucleotide polymorphism (SNP) sites associated with human diseases. In the annotation of human diseases, especially cancers, SNPs, as an important component of genetic factors, have gained increasing attention. Given that most of the SNPs are located in non-coding regions, the functional verification of these SNPs is a great challenge. The key to functional annotation for risk SNPs is to screen SNPs with regulatory activity from thousands of disease associated-SNPs. In this review, we systematically recapitulate the characteristics and functional roles of SNP sites, discuss three parallel reporter screening strategies in detail based on barcode tag classification, and recommend the common in silico strategies to help supplement the annotation of SNP sites with epigenetic activity analysis, prediction of target genes and trans-acting factors. We hope that this review will contribute to this exuberant research field by providing robust activity analysis strategies that can facilitate the translation of GWAS results into personalized diagnosis and prevention measures for human diseases.

Epigenetically programmed resistance to chemo- and immuno-therapies

Resistance to cancer treatments remains a major barrier in developing cancer cures. While promising combination chemotherapy treatments and novel immunotherapies have improved patient outcomes, resistance to these treatments remains poorly understood. New insights into the dysregulation of the epigenome show how it promotes tumor growth and resistance to therapy. By altering control of gene expression, tumor cells can evade immune cell recognition, ignore apoptotic cues, and reverse DNA damage induced by chemotherapies. In this chapter, we summarize the data on epigenetic remodeling during cancer progression and treatment that enable cancer cell survival and describe how these epigenetic changes are being targeted clinically to overcome resistance.

Enhancer/gene relationships: Need for more reliable genome-wide reference sets

Differences in cells' functions arise from differential activity of regulatory elements, including enhancers. Enhancers are cis-regulatory elements that cooperate with promoters through transcription factors to activate the expression of one or several genes by getting physically close to them in the 3D space of the nucleus. There is increasing evidence that genetic variants associated with common diseases are enriched in enhancers active in cell types relevant to these diseases. Identifying the enhancers associated with genes and conversely, the sets of genes activated by each enhancer (the so-called enhancer/gene or E/G relationships) across cell types, can help understanding the genetic mechanisms underlying human diseases. There are three broad approaches for the genome-wide identification of E/G relationships in a cell type: 1) genetic link methods or eQTL, 2) functional link methods based on 1D functional data such as open chromatin, histone mark or gene expression and 3) spatial link methods based on 3D data such as HiC. Since 1) and 3) are costly, the current strategy is to develop functional link methods and to use data from 1) and 3) as reference to evaluate them. However, there is still no consensus on the best functional link method to date, and method comparison remain seldom. Here, we compared the relative performances of three recent methods for the identification of enhancer-gene links, TargetFinder, Average-Rank, and the ABC model, using the three latest benchmarks from the field: a reference that combines 3D and eQTL data, called BENGI, and two genetic screening references, called CRiFF and CRiSPRi. Overall, none of the three methods performed best on the three references. CRiFF and CRISPRi reference sets are likely more reliable, but CRiFF is not genome-wide and CRiFF and CRISPRi are mostly available on the K562 cancer cell line. The BENGI reference set is genome-wide but likely contains many false positives. This study therefore calls for new reliable and genome-wide E/G reference data rather than new functional link E/G identification methods.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Atlas of regulated target genes of transcription factors (ART-TF) in human ES cells

Background

Transcription factors (TFs) play central roles in maintaining “stemness” of embryonic stem (ES) cells and their differentiation into several hundreds of adult cell types. The regulatory competence of TFs is routinely assessed by detecting target genes to which they bind. However, these data do not indicate which target genes are activated, repressed, or not affected by the change of TF abundance. There is a lack of large-scale studies that compare the genome binding of TFs with the expression change of target genes after manipulation of each TF.

Results

In this paper we associated human TFs with their target genes by two criteria: binding to genes, evaluated from published ChIP-seq data (n = 1868); and change of target gene expression shortly after induction of each TF in human ES cells. Lists of direction- and strength-specific regulated target genes are generated for 311 TFs (out of 351 TFs tested) with expected proportion of false positives less than or equal to 0.30, including 63 new TFs not present in four existing databases of target genes. Our lists of direction-specific targets for 152 TFs (80.0%) are larger that in the TRRUST database. In average, 30.9% of genes that respond greater than or equal to twofold to the induction of TFs are regulated targets. Regulated target genes indicate that the majority of TFs are either strong activators or strong repressors, whereas sets of genes that responded greater than or equal to twofold to the induction of TFs did not show strong asymmetry in the direction of expression change. The majority of human TFs (82.1%) regulated their target genes primarily via binding to enhancers. Repression of target genes is more often mediated by promoter-binding than activation of target genes. Enhancer-promoter loops are more abundant among strong activator and repressor TFs.

Conclusions

We developed an atlas of regulated targets of TFs (ART-TF) in human ES cells by combining data on TF binding with data on gene expression change after manipulation of individual TFs. Sets of regulated gene targets were identified with a controlled rate of false positives. This approach contributes to the understanding of biological functions of TFs and organization of gene regulatory networks. This atlas should be a valuable resource for ES cell-based regenerative medicine studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04924-3.

Erythroid Cell Research: 3D Chromatin, Transcription Factors and Beyond

Studies of the regulatory networks and signals controlling erythropoiesis have brought important insights in several research fields of biology and have been a rich source of discoveries with far-reaching implications beyond erythroid cells biology. The aim of this review is to highlight key recent discoveries and show how studies of erythroid cells bring forward novel concepts and refine current models related to genome and 3D chromatin organization, signaling and disease, with broad interest in life sciences.

Epigenetic Reprogramming in Host-Parasite Coevolution: The Toxoplasma Paradigm

Like many intracellular pathogens, the protozoan parasite Toxoplasma gondii has evolved sophisticated mechanisms to promote its transmission and persistence in a variety of hosts by injecting effector proteins that manipulate many processes in the cells it invades. Specifically, the parasite diverts host epigenetic modulators and modifiers from their native functions to rewire host gene expression to counteract the innate immune response and to limit its strength. The arms race between the parasite and its hosts has led to accelerated adaptive evolution of effector proteins and the unconventional secretion routes they use. This review provides an up-to-date overview of how T. gondii effectors, through the evolution of intrinsically disordered domains, the formation of supramolecular complexes, and the use of molecular mimicry, target host transcription factors that act as coordinating nodes, as well as chromatin-modifying enzymes, to control the fate of infected cells and ultimately the outcome of infection. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

3D chromatin architecture and transcription regulation in cancer

Chromatin has distinct three-dimensional (3D) architectures important in key biological processes, such as cell cycle, replication, differentiation, and transcription regulation. In turn, aberrant 3D structures play a vital role in developing abnormalities and diseases such as cancer. This review discusses key 3D chromatin structures (topologically associating domain, lamina-associated domain, and enhancer-promoter interactions) and corresponding structural protein elements mediating 3D chromatin interactions [CCCTC-binding factor, polycomb group protein, cohesin, and Brother of the Regulator of Imprinted Sites (BORIS) protein] with a highlight of their associations with cancer. We also summarise the recent development of technologies and bioinformatics approaches to study the 3D chromatin interactions in gene expression regulation, including crosslinking and proximity ligation methods in the bulk cell population (ChIA-PET and HiChIP) or single-molecule resolution (ChIA-drop), and methods other than proximity ligation, such as GAM, SPRITE, and super-resolution microscopy techniques.

Structural Variation at a Disease Mutation Hotspot: Strategies to Investigate Gene Regulation and the 3D Genome

A rare form of X-linked Charcot-Marie-Tooth neuropathy, CMTX3, is caused by an interchromosomal insertion occurring at chromosome Xq27.1. Interestingly, eight other disease phenotypes have been associated with insertions (or insertion-deletions) occurring at the same genetic locus. To date, the pathogenic mechanism underlying most of these diseases remains unsolved, although local gene dysregulation has clearly been implicated in at least two phenotypes. The challenges of accessing disease-relevant tissue and modelling these complex genomic rearrangements has led to this research impasse. We argue that recent technological advancements can overcome many of these challenges, particularly induced pluripotent stem cells (iPSC) and their capacity to provide access to patient-derived disease-relevant tissue. However, to date these valuable tools have not been utilized to investigate the disease-associated insertions at chromosome Xq27.1. Therefore, using CMTX3 as a reference disease, we propose an experimental approach that can be used to explore these complex mutations, as well as similar structural variants located elsewhere in the genome. The mutational hotspot at Xq27.1 is a valuable disease paradigm with the potential to improve our understanding of the pathogenic consequences of complex structural variation, and more broadly, refine our knowledge of the multifaceted process of long-range gene regulation. Intergenic structural variation is a critically understudied class of mutation, although it is likely to contribute significantly to unsolved genetic disease.

Prioritisation of Candidate Genes Underpinning COVID-19 Host Genetic Traits Based on High-Resolution 3D Chromosomal Topology

Genetic variants showing associations with specific biological traits and diseases detected by genome-wide association studies (GWAS) commonly map to non-coding DNA regulatory regions. Many of these regions are located considerable distances away from the genes they regulate and come into their proximity through 3D chromosomal interactions. We previously developed COGS, a statistical pipeline for linking GWAS variants with their putative target genes based on 3D chromosomal interaction data arising from high-resolution assays such as Promoter Capture Hi-C (PCHi-C). Here, we applied COGS to COVID-19 Host Genetic Consortium (HGI) GWAS meta-analysis data on COVID-19 susceptibility and severity using our previously generated PCHi-C results in 17 human primary cell types and SARS-CoV-2-infected lung carcinoma cells. We prioritise 251 genes putatively associated with these traits, including 16 out of 47 genes highlighted by the GWAS meta-analysis authors. The prioritised genes are expressed in a broad array of tissues, including, but not limited to, blood and brain cells, and are enriched for genes involved in the inflammatory response to viral infection. Our prioritised genes and pathways, in conjunction with results from other prioritisation approaches and targeted validation experiments, will aid in the understanding of COVID-19 pathology, paving the way for novel treatments.