Enhancers predominantly regulate gene expression during differentiation via transcription initiation.

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and ß-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

Absolute Quantification of Transcription Factors Reveals Principles of Gene Regulation in Erythropoiesis.

Dynamic cellular processes such as differentiation are driven by changes in the abundances of transcription factors (TFs). However, despite years of studies, our knowledge about the protein copy number of TFs in the nucleus is limited. Here, by determining the absolute abundances of 103 TFs and co-factors during the course of human erythropoiesis, we provide a dynamic and quantitative scale for TFs in the nucleus. Furthermore, we establish the first gene regulatory network of cell fate commitment that integrates temporal protein stoichiometry data with mRNA measurements. The model revealed quantitative imbalances in TFs' cross-antagonistic relationships that underlie lineage determination. Finally, we made the surprising discovery that, in the nucleus, co-repressors are dramatically more abundant than co-activators at the protein level, but not at the RNA level, with profound implications for understanding transcriptional regulation. These analyses provide a unique quantitative framework to understand transcriptional regulation of cell differentiation in a dynamic context.

Defining genome architecture at base-pair resolution.

In higher eukaryotes, many genes are regulated by enhancers that are 104-106 base pairs (bp) away from the promoter. Enhancers contain transcription-factor-binding sites (which are typically around 7-22 bp), and physical contact between the promoters and enhancers is thought to be required to modulate gene expression. Although chromatin architecture has been mapped extensively at resolutions of 1 kilobase and above; it has not been possible to define physical contacts at the scale of the proteins that determine gene expression. Here we define these interactions in detail using a chromosome conformation capture method (Micro-Capture-C) that enables the physical contacts between different classes of regulatory elements to be determined at base-pair resolution. We find that highly punctate contacts occur between enhancers, promoters and CCCTC-binding factor (CTCF) sites and we show that transcription factors have an important role in the maintenance of the contacts between enhancers and promoters. Our data show that interactions between CTCF sites are increased when active promoters and enhancers are located within the intervening chromatin. This supports a model in which chromatin loop extrusion1 is dependent on cohesin loading at active promoters and enhancers, which explains the formation of tissue-specific chromatin domains without changes in CTCF binding.

Impaired DNA demethylation of C/EBP sites causes premature aging.

Changes in DNA methylation are among the best-documented epigenetic alterations accompanying organismal aging. However, whether and how altered DNA methylation is causally involved in aging have remained elusive. GADD45α (growth arrest and DNA damage protein 45A) and ING1 (inhibitor of growth family member 1) are adapter proteins for site-specific demethylation by TET (ten-eleven translocation) methylcytosine dioxygenases. Here we show that Gadd45a/Ing1 double-knockout mice display segmental progeria and phenocopy impaired energy homeostasis and lipodystrophy characteristic of Cebp (CCAAT/enhancer-binding protein) mutants. Correspondingly, GADD45α occupies C/EBPß/δ-dependent superenhancers and, cooperatively with ING1, promotes local DNA demethylation via long-range chromatin loops to permit C/EBPß recruitment. The results indicate that enhancer methylation can affect aging and imply that C/EBP proteins play an unexpected role in this process. Our study suggests a causal nexus between DNA demethylation, metabolism, and organismal aging.

Natural and Experimental Rewiring of Gene Regulatory Regions.

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

LanceOtron: a deep learning peak caller for genome sequencing experiments.

MOTIVATION: Genome sequencing experiments have revolutionized molecular biology by allowing researchers to identify important DNA-encoded elements genome wide. Regions where these elements are found appear as peaks in the analog signal of an assay's coverage track, and despite the ease with which humans can visually categorize these patterns, the size of many genomes necessitates algorithmic implementations. Commonly used methods focus on statistical tests to classify peaks, discounting that the background signal does not completely follow any known probability distribution and reducing the information-dense peak shapes to simply maximum height. Deep learning has been shown to be highly accurate for many pattern recognition tasks, on par or even exceeding human capabilities, providing an opportunity to reimagine and improve peak calling. RESULTS: We present the peak calling framework LanceOtron, which combines deep learning for recognizing peak shape with multifaceted enrichment calculations for assessing significance. In benchmarking ATAC-seq, ChIP-seq and DNase-seq, LanceOtron outperforms long-standing, gold-standard peak callers through its improved selectivity and near-perfect sensitivity. AVAILABILITY AND IMPLEMENTATION: A fully featured web application is freely available from LanceOtron.molbiol.ox.ac.uk, command line interface via python is pip installable from PyPI at https://pypi.org/project/lanceotron/, and source code and benchmarking tests are available at https://github.com/LHentges/LanceOtron. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic and functional insights into CDA-I prevalence and pathogenesis.

BACKGROUND: Congenital dyserythropoietic anaemia type I (CDA-I) is a hereditary anaemia caused by biallelic mutations in the widely expressed genes CDAN1 and C15orf41. Little is understood about either protein and it is unclear in which cellular pathways they participate. METHODS: Genetic analysis of a cohort of patients with CDA-I identifies novel pathogenic variants in both known causative genes. We analyse the mutation distribution and the predicted structural positioning of amino acids affected in Codanin-1, the protein encoded by CDAN1. Using western blotting, immunoprecipitation and immunofluorescence, we determine the effect of particular mutations on both proteins and interrogate protein interaction, stability and subcellular localisation. RESULTS: We identify six novel CDAN1 mutations and one novel mutation in C15orf41 and uncover evidence of further genetic heterogeneity in CDA-I. Additionally, population genetics suggests that CDA-I is more common than currently predicted. Mutations are enriched in six clusters in Codanin-1 and tend to affect buried residues. Many missense and in-frame mutations do not destabilise the entire protein. Rather C15orf41 relies on Codanin-1 for stability and both proteins, which are enriched in the nucleolus, interact to form an obligate complex in cells. CONCLUSION: Stability and interaction data suggest that C15orf41 may be the key determinant of CDA-I and offer insight into the mechanism underlying this disease. Both proteins share a common pathway likely to be present in a wide variety of cell types; however, nucleolar enrichment may provide a clue as to the erythroid specific nature of CDA-I. The surprisingly high predicted incidence of CDA-I suggests that better ascertainment would lead to improved patient care.

Recapitulation of erythropoiesis in congenital dyserythropoietic anaemia type I (CDA-I) identifies defects in differentiation and nucleolar abnormalities.

The investigation of inherited disorders of erythropoiesis has elucidated many of the principles underlying the production of normal red blood cells and how this is perturbed in human disease. Congenital Dyserythropoietic Anaemia type 1 (CDA-I) is a rare form of anaemia caused by mutations in two genes of unknown function: CDAN1 and CDIN1 (previously called C15orf41), whilst in some cases, the underlying genetic abnormality is completely unknown. Consequently, the pathways affected in CDA-I remain to be discovered. To enable detailed analysis of this rare disorder we have validated a culture system which recapitulates all of the cardinal haematological features of CDA-I, including the formation of the pathognomonic 'spongy' heterochromatin seen by electron microscopy. Using a variety of cell and molecular biological approaches we discovered that erythroid cells in this condition show a delay during terminal erythroid differentiation, associated with increased proliferation and widespread changes in chromatin accessibility. We also show that the proteins encoded by CDAN1 and CDIN1 are enriched in nucleoli which are structurally and functionally abnormal in CDA-I. Together these findings provide important pointers to the pathways affected in CDA-I which for the first time can now be pursued in the tractable culture system utilised here.

Functional characterisation of cis-regulatory elements governing dynamic Eomes expression in the early mouse embryo.

The T-box transcription factor (TF) Eomes is a key regulator of cell fate decisions during early mouse development. The cis-acting regulatory elements that direct expression in the anterior visceral endoderm (AVE), primitive streak (PS) and definitive endoderm (DE) have yet to be defined. Here, we identified three gene-proximal enhancer-like sequences (PSE_a, PSE_b and VPE) that faithfully activate tissue-specific expression in transgenic embryos. However, targeted deletion experiments demonstrate that PSE_a and PSE_b are dispensable, and only VPE is required for optimal Eomes expression in vivo Embryos lacking this enhancer display variably penetrant defects in anterior-posterior axis orientation and DE formation. Chromosome conformation capture experiments reveal VPE-promoter interactions in embryonic stem cells (ESCs), prior to gene activation. The locus resides in a large (500 kb) pre-formed compartment in ESCs and activation during DE differentiation occurs in the absence of 3D structural changes. ATAC-seq analysis reveals that VPE, PSE_a and four additional putative enhancers display increased chromatin accessibility in DE that is associated with Smad2/3 binding coincident with transcriptional activation. By contrast, activation of the Eomes target genes Foxa2 and Lhx1 is associated with higher order chromatin reorganisation. Thus, diverse regulatory mechanisms govern activation of lineage specifying TFs during early development.

Robust detection of chromosomal interactions from small numbers of cells using low-input Capture-C.

Chromosome conformation capture (3C) techniques are crucial to understanding tissue-specific regulation of gene expression, but current methods generally require large numbers of cells. This hampers the investigation of chromatin architecture in rare cell populations. We present a new low-input Capture-C approach that can generate high-quality 3C interaction profiles from 10 000-20 000 cells, depending on the resolution used for analysis. We also present a PCR-free, sequencing-free 3C technique based on NanoString technology called C-String. By comparing C-String and Capture-C interaction profiles we show that the latter are not skewed by PCR amplification. Furthermore, we demonstrate that chromatin interactions detected by Capture-C do not depend on the degree of cross-linking by performing experiments with varying formaldehyde concentrations.

Dual DNA binding and coactivator functions of Aspergillus nidulans†TamA, a Zn(II)2Cys6 transcription factor.

Transcription factors containing DNA binding domains generally regulate transcription by direct interaction with DNA. For most transcription factors, including the fungal Zn(II)2Cys6 zinc binuclear cluster transcription factors, the DNA binding motif is essential for function. However, Aspergillus nidulans TamA and the related Saccharomyces cerevisiae Dal81p protein contain Zn(II)2Cys6 motifs shown to be dispensable for function. TamA acts at several promoters as a coactivator of the global nitrogen GATA transcription factor AreA. We now show that TamA is the major transcriptional activator of gdhA, encoding the key nitrogen metabolism enzyme NADP-glutamate dehydrogenase. Moreover, activation of gdhA by TamA occurs primarily by a mechanism requiring the TamA DNA binding motif. We show that the TamA DNA binding motif is required for DNA binding of FLAG-epitope-tagged TamA to the gdhA promoter. We identify a conserved promoter element required for TamA activation, and show that TamA and AreA are reciprocally required for full binding at the gdhA promoter under conditions where AreA is inactive at most promoters but active at gdhA. Therefore TamA has dual functions as a DNA-binding transcription factor and a non-DNA-binding coactivator. Dual DNA-binding and coactivator functions provide an additional level of combinatorial control to mediate gene-specific expression.

Multiple nuclear localization signals mediate nuclear localization of the GATA transcription factor AreA.

The Aspergillus nidulans GATA transcription factor AreA activates transcription of nitrogen metabolic genes in response to nitrogen limitation and is known to accumulate in the nucleus during nitrogen starvation. Sequence analysis of AreA revealed multiple nuclear localization signals (NLSs), five putative classical NLSs conserved in fungal AreA orthologs but not in the Saccharomyces cerevisiae functional orthologs Gln3p and Gat1p, and one putative noncanonical RRX33RXR bipartite NLS within the DNA-binding domain. In order to identify the functional NLSs in AreA, we constructed areA mutants with mutations in individual putative NLSs or combinations of putative NLSs and strains expressing green fluorescent protein (GFP)-AreA NLS fusion genes. Deletion of all five classical NLSs individually or collectively did not affect utilization of nitrogen sources or AreA-dependent gene expression and did not prevent AreA nuclear localization. Mutation of the bipartite NLS conferred the inability to utilize alternative nitrogen sources and abolished AreA-dependent gene expression likely due to effects on DNA binding but did not prevent AreA nuclear localization. Mutation of all six NLSs simultaneously prevented AreA nuclear accumulation. The bipartite NLS alone strongly directed GFP to the nucleus, whereas the classical NLSs collaborated to direct GFP to the nucleus. Therefore, AreA contains multiple conserved NLSs, which show redundancy and together function to mediate nuclear import. The noncanonical bipartite NLS is conserved in GATA factors from Aspergillus, yeast, and mammals, indicating an ancient origin.

Regulation of the NADP-glutamate dehydrogenase gene gdhA in Aspergillus nidulans by the Zn(II)2Cys6 transcription factor LeuB.

NADP-dependent glutamate dehydrogenase (NADP-GDH) is a key enzyme in the assimilation of alternative nitrogen nutrient sources through ammonium in fungi. In Aspergillus nidulans, NADP-GDH is encoded by gdhA. Several transcription factors are known to regulate gdhA expression, including AreA, the major transcription activator of nitrogen metabolic genes, and TamA, a co-activator of AreA. TamA also interacts with LeuB, the regulator of leucine biosynthesis. We have investigated the effects of leucine biosynthesis on gdhA regulation, and found that leucine regulates the levels of NADP-GDH activity and gdhA expression. We show, using mutants with perturbed levels of α-isopropylmalate (α-IPM), that this leucine biosynthesis intermediate affects gdhA regulation. Leucine regulation of gdhA requires a functional LeuB with an intact Zn(II)2Cys6 DNA-binding domain. By analysing the prevalence of putative LeuB DNA-binding sites in promoters of gdhA orthologues we predict broad conservation of leucine regulation of NADP-GDH expression within ascomycetes except in the fusaria and fission yeasts. Using promoter mutations in gdhA-lacZ reporter genes we identified two sites of action for LeuB within the A. nidulans gdhA promoter. These two sites lack sequence identity, with one site conforming to the predicted LeuB DNA-binding site consensus motif, whereas the second site is a novel regulatory sequence element conserved in Aspergillus gdhA promoters. These data suggest that LeuB regulates NADP-GDH expression in response to leucine levels, which may act as an important sensor of nitrogen availability.

Inferring causal genes at type 2 diabetes GWAS loci through chromosome interactions in islet cells.

Background: Resolving causal genes for type 2 diabetes at loci implicated by genome-wide association studies (GWAS) requires integrating functional genomic data from relevant cell types. Chromatin features in endocrine cells of the pancreatic islet are particularly informative and recent studies leveraging chromosome conformation capture (3C) with Hi-C based methods have elucidated regulatory mechanisms in human islets. However, these genome-wide approaches are less sensitive and afford lower resolution than methods that target specific loci. Methods: To gauge the extent to which targeted 3C further resolves chromatin-mediated regulatory mechanisms at GWAS loci, we generated interaction profiles at 23 loci using next-generation (NG) capture-C in a human beta cell model (EndoC-ßH1) and contrasted these maps with Hi-C maps in EndoC-ßH1 cells and human islets and a promoter capture Hi-C map in human islets. Results: We found improvements in assay sensitivity of up to 33-fold and resolved ~3.6X more chromatin interactions. At a subset of 18 loci with 25 co-localised GWAS and eQTL signals, NG Capture-C interactions implicated effector transcripts at five additional genetic signals relative to promoter capture Hi-C through physical contact with gene promoters. Conclusions: High resolution chromatin interaction profiles at selectively targeted loci can complement genome- and promoter-wide maps.

Assessment of Multiway Interactions with Tri-C.

Tri-C is a chromosome conformation capture (3C) approach that can efficiently identify multiway chromatin interactions with viewpoints of interest. As opposed to pair-wise interactions identified in methods such as Hi-C, 4C, and Capture-C, the detection of multiway interactions allows researchers to investigate how multiple cis-regulatory elements interact together in higher-order structures in single nuclei and address questions regarding structural cooperation between these elements. Here, we describe the procedure for designing and performing a Tri-C experiment.

The chromatin remodeller ATRX facilitates diverse nuclear processes, in a stochastic manner, in both heterochromatin and euchromatin.

The chromatin remodeller ATRX interacts with the histone chaperone DAXX to deposit the histone variant H3.3 at sites of nucleosome turnover. ATRX is known to bind repetitive, heterochromatic regions of the genome including telomeres, ribosomal DNA and pericentric repeats, many of which are putative G-quadruplex forming sequences (PQS). At these sites ATRX plays an ancillary role in a wide range of nuclear processes facilitating replication, chromatin modification and transcription. Here, using an improved protocol for chromatin immunoprecipitation, we show that ATRX also binds active regulatory elements in euchromatin. Mutations in ATRX lead to perturbation of gene expression associated with a reduction in chromatin accessibility, histone modification, transcription factor binding and deposition of H3.3 at the sequences to which it normally binds. In erythroid cells where downregulation of α-globin expression is a hallmark of ATR-X syndrome, perturbation of chromatin accessibility and gene expression occurs in only a subset of cells. The stochastic nature of this process suggests that ATRX acts as a general facilitator of cell specific transcriptional and epigenetic programmes, both in heterochromatin and euchromatin.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture.

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.