Transcription factors mediate long-range enhancer-promoter interactions

Genetic factors mediating long-range enhancer-promoter communication in mammalian development

Enhancers are remotely located noncoding DNA sequences that regulate gene expression in response to developmental, homeostatic, and environmental cues. Canonical short-range enhancers located <50 kb from their cognate promoters function by binding transcription factors, coactivators, and chromatin modifiers. In this review, we discuss recent evidence that medium-range (50-400 kb) and long-range (>400 kb) enhancers rely on additional mechanisms, including cohesin, CCCTC-binding factor, and high-affinity protein-protein interactions. These mechanisms are crucial for establishing the physical proximity and interaction between enhancers and their target promoters over extended genomic distances and ensuring robust gene activation during mammalian development. Future studies will be critical to unravel their prevalence and evolutionary significance across various genomic loci, cell types, and species.

CRISPR-mediated iron and folate biofortification in crops: advances and perspectives

Micronutrient deficiency conditions, such as anemia, are the most prevalent global health problem due to inadequate iron and folate in dietary sources. Biofortification advancements can propel the rapid amelioration of nutritionally beneficial components in crops that are required to combat the adverse effects of micronutrient deficiencies on human health. To date, several strategies have been proposed to increase micronutrients in plants to improve food quality, but very few approaches have intrigued `clustered regularly interspaced short palindromic repeats' (CRISPR) modules for the enhancement of iron and folate concentration in the edible parts of plants. In this review, we discuss two important approaches to simultaneously enhance the bioavailability of iron and folate concentrations in rice endosperms by utilizing advanced CRISPR-Cas9-based technology. This includes the 'tuning of cis-elements' and 'enhancer re-shuffling' in the regulatory components of genes that play a vital role in iron and folate biosynthesis/transportation pathways. In particular, base-editing and enhancer re-installation in native promoters of selected genes can lead to enhanced accumulation of iron and folate levels in the rice endosperm. The re-distribution of micronutrients in specific plant organs can be made possible using the above-mentioned contemporary approaches. Overall, the present review discusses the possible approaches for synchronized iron and folate biofortification through modification in regulatory gene circuits employing CRISPR-Cas9 technology.

Fluctuating Chromatin Facilitates Enhancer-Promoter Communication by Regulating Transcriptional Clustering Dynamics

Enhancers regulate gene expression by forming contacts with distant promoters. Phase-separated condensates or clusters formed by transcription factors (TFs) and cofactors are thought to facilitate these enhancer-promoter (E-P) interactions. Using polymer physics, we developed distinct coarse-grained chromatin models that produce similar ensemble-averaged Hi-C maps but with "stable" and "dynamic" characteristics. Our findings, consistent with recent experiments, reveal a multistep E-P communication process. The dynamic model facilitates E-P proximity by enhancing TF clustering and subsequently promotes direct E-P interactions by destabilizing the TF clusters through chain flexibility. Our study promotes physical understanding of the molecular mechanisms governing E-P communication in transcriptional regulation.

The EGR1-Artemin Axis in Keratinocytes Enhances the Innervation of Epidermal Sensory Neurons during Skin Inflammation Induced by House Dust Mite Extract from Dermatophagoidesfarinae

Epidermal hyperinnervation is a critical feature of pruritus during skin inflammation. However, the mechanisms underlying epidermal hyperinnervation are unclear. This study investigates the role of the transcription factor EGR1 in epidermal innervation by utilizing wild-type (Egr1+/+) and Egr1-null (Egr1‒/‒) mice topically applied Dermatophagoides farinae extract from dust mite. Our findings revealed that Egr1‒/‒ mice exhibited reduced scratching behaviors and decreased density of epidermal innervation compared with Egr1+/+ mice. Furthermore, we identified artemin, a neurotrophic factor, as an EGR1 target responsible for Dermatophagoides farinae extract-induced hyperinnervation. It has been demonstrated that Dermatophagoides farinae extract stimulates toll-like receptors in keratinocytes. To elucidate the cellular mechanism, we stimulated keratinocytes with Pam3CSK4, a toll-like receptor 1/2 ligand. Pam3CSK4 triggered a toll-like receptor 1/2-mediated signaling cascade involving IRAK4, IκB kinase, MAPKs, ELK1, EGR1, and artemin, leading to increased neurite outgrowth and neuronal migration. In addition, increased expression of EGR1 and artemin was observed in the skin tissues of patients with atopic dermatitis. These findings highlight the significance of the EGR1-artemin axis in keratinocytes, promoting the process of epidermal innervation and suggesting it as a potential therapeutic target for alleviating itch and pain associated with house dust mite-induced skin inflammation.

Functional Identification of Malus halliana MhbZIP23 Gene Demonstrates That It Enhances Saline-Alkali Stress Tolerance in Arabidopsis thaliana

Saline-alkali stress is a significant abiotic stress that restricts plant growth globally. Basic region leucine zipper (bZIP) transcription factor proteins are widely involved in plants in response to abiotic stress such as saline-alkali stress. Based on transcriptome and quantitative real-time PCR (qRT-PCR), we found that the MhbZIP23 gene could respond to saline-alkali stress. Despite this discovery, the underlying mechanism by which the MhbZIP23 transcription factor responds to saline-alkaline stress remains unexplored. To address this gap in knowledge, we successfully cloned the MhbZIP23 (MD05G1121500) gene from Malus halliana for heterologous expression in Arabidopsis thaliana, facilitating the investigation of its functional role in stress response. Compared to the wild type (WT), Arabidopsis plants demonstrated enhanced growth and a lower degree of wilting when subjected to saline-alkali stress. Furthermore, several physiological indices of the plants altered under such stress conditions. The transgenic Arabidopsis plants (OE-5, 6, and 8), which grew normally, exhibited a higher chlorophyll content and had greater root length in comparison to the control check (CK). MhbZIP23 effectively regulated the levels of the osmoregulatory substance proline (Pro), enhanced the activities of antioxidant enzymes such as peroxidase (POD) and superoxide dismutase (SOD), and reduced the levels of malondialdehyde (MDA) and relative conductivity (REC). These actions improved the ability of plant cells in transgenic Arabidopsis to counteract ROS, as evidenced by the decreased accumulation of O2- and hydrogen peroxide (H2O2). In summary, the MhbZIP23 gene demonstrated effectiveness in alleviating saline-alkali stress in M. halliana, presenting itself as an outstanding resistance gene for apples to combat saline-alkali stress.

Non-canonical functions of enhancers: regulation of RNA polymerase III transcription, DNA replication, and V(D)J recombination

Enhancers are the key regulators of other DNA-based processes by virtue of their unique ability to generate nucleosome-depleted regions in a highly regulated manner. Enhancers regulate cell-type-specific transcription of tRNA genes by RNA polymerase III (Pol III). They are also responsible for the binding of the origin replication complex (ORC) to DNA replication origins, thereby regulating origin utilization, replication timing, and replication-dependent chromosome breaks. Additionally, enhancers regulate V(D)J recombination by increasing access of the recombination-activating gene (RAG) recombinase to target sites and by generating non-coding enhancer RNAs and localized regions of trimethylated histone H3-K4 recognized by the RAG2 PHD domain. Thus, enhancers represent the first step in decoding the genome, and hence they regulate biological processes that, unlike RNA polymerase II (Pol II) transcription, do not have dedicated regulatory proteins.

The Blueprint of Logical Decisions in a NF-κB Signaling System

Nearly identical cells can exhibit substantially different responses to the same stimulus that causes phenotype diversity. Such interplay between phenotype diversity and the architecture of regulatory circuits is crucial since it determines the state of a biological cell. Here, we theoretically analyze how the circuit blueprints of NF-κB in cellular environments are formed and their role in determining the cells' metabolic state. The NF-κB is a collective name for a developmental conserved family of five different transcription factors that can form homodimers or heterodimers and often promote DNA looping to reprogram the inflammatory gene response. The NF-κB controls many biological functions, including cellular differentiation, proliferation, migration, and survival. Our model shows that nuclear localization of NF-κB differentially promotes logic operations such as AND, NAND, NOR, and OR in its regulatory network. Through the quantitative thermodynamic model of transcriptional regulation and systematic variation of promoter-enhancer interaction modes, we can account for the origin of various logic gates as formed in the NF-κB system. We further show that the interconversion or switching of logic gates yielded under systematic variations of the stimuli activity and DNA looping parameters. Such computation occurs in regulatory and signaling pathways in individual cells at a molecular scale, which one can exploit to design a biomolecular computer.

Polymer Collapse & Liquid-Liquid Phase-Separation are Coupled in a Generalized Prewetting Transition

The three-dimensional organization of chromatin is thought to play an important role in controlling gene expression. Specificity in expression is achieved through the interaction of transcription factors and other nuclear proteins with particular sequences of DNA. At unphysiological concentrations many of these nuclear proteins can phase-separate in the absence of DNA, and it has been hypothesized that, in vivo, the thermodynamic forces driving these phases help determine chromosomal organization. However it is unclear how DNA, itself a long polymer subject to configurational transitions, interacts with three-dimensional protein phases. Here we show that a long compressible polymer can be coupled to interacting protein mixtures, leading to a generalized prewetting transition where polymer collapse is coincident with a locally stabilized liquid droplet. We use lattice Monte-Carlo simulations and a mean-field theory to show that these phases can be stable even in regimes where both polymer collapse and coexisting liquid phases are unstable in isolation, and that these new transitions can be either abrupt or continuous. For polymers with internal linear structure we further show that changes in the concentration of bulk components can lead to changes in three-dimensional polymer structure. In the nucleus there are many distinct proteins that interact with many different regions of chromatin, potentially giving rise to many different Prewet phases. The simple systems we consider here highlight chromatin's role as a lower-dimensional surface whose interactions with proteins are required for these novel phases.

Polymer Collapse & Liquid-Liquid Phase-Separation are Coupled in a Generalized Prewetting Transition

The three-dimensional organization of chromatin is thought to play an important role in controlling gene expression. Specificity in expression is achieved through the interaction of transcription factors and other nuclear proteins with particular sequences of DNA. At unphysiological concentrations many of these nuclear proteins can phase-separate in the absence of DNA, and it has been hypothesized that, in vivo, the thermodynamic forces driving these phases help determine chromosomal organization. However it is unclear how DNA, itself a long polymer subject to configurational transitions, interacts with three-dimensional protein phases. Here we show that a long compressible polymer can be coupled to interacting protein mixtures, leading to a generalized prewetting transition where polymer collapse is coincident with a locally stabilized liquid droplet. We use lattice Monte-Carlo simulations and a mean-field theory to show that these phases can be stable even in regimes where both polymer collapse and coexisting liquid phases are unstable in isolation, and that these new transitions can be either abrupt or continuous. For polymers with internal linear structure we further show that changes in the concentration of bulk components can lead to changes in three-dimensional polymer structure. In the nucleus there are many distinct proteins that interact with many different regions of chromatin, potentially giving rise to many different Prewet phases. The simple systems we consider here highlight chromatin's role as a lower-dimensional surface whose interactions with proteins are required for these novel phases.

E2F1 promotes cell migration in hepatocellular carcinoma via FNDC3B

FNDC3B (fibronectin type III domain containing 3B) is highly expressed in hepatocellular carcinoma (HCC) and other cancer types, and fusion genes involving FNDC3B have been identified in HCC and leukemia. Growing evidence suggests the significance of FNDC3B in tumorigenesis, particularly in cell migration and tumor metastasis. However, its regulatory mechanisms remain elusive. In this study, we employed bioinformatic, gene regulation, and protein-DNA interaction screening to investigate the transcription factors (TFs) involved in regulating FNDC3B. Initially, 338 candidate TFs were selected based on previous chromatin immunoprecipitation (ChIP)-seq experiments available in public domain databases. Through TF knockdown screening and ChIP coupled with Droplet Digital PCR assays, we identified that E2F1 (E2F transcription factor 1) is crucial for the activation of FNDC3B. Overexpression or knockdown of E2F1 significantly impacts the expression of FNDC3B. In conclusion, our study elucidated the mechanistic link between FNDC3B and E2F1. These findings contribute to a better understanding of FNDC3B in tumorigenesis and provide insights into potential therapeutic targets for cancer treatment.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

Enhancer-MDLF: a novel deep learning framework for identifying cell-specific enhancers

Enhancers, noncoding DNA fragments, play a pivotal role in gene regulation, facilitating gene transcription. Identifying enhancers is crucial for understanding genomic regulatory mechanisms, pinpointing key elements and investigating networks governing gene expression and disease-related mechanisms. Existing enhancer identification methods exhibit limitations, prompting the development of our novel multi-input deep learning framework, termed Enhancer-MDLF. Experimental results illustrate that Enhancer-MDLF outperforms the previous method, Enhancer-IF, across eight distinct human cell lines and exhibits superior performance on generic enhancer datasets and enhancer-promoter datasets, affirming the robustness of Enhancer-MDLF. Additionally, we introduce transfer learning to provide an effective and potential solution to address the prediction challenges posed by enhancer specificity. Furthermore, we utilize model interpretation to identify transcription factor binding site motifs that may be associated with enhancer regions, with important implications for facilitating the study of enhancer regulatory mechanisms. The source code is openly accessible at https://github.com/HaoWuLab-Bioinformatics/Enhancer-MDLF.

Remote regulation of rs80245547 and rs72673891 mediated by transcription factors C-Jun and CREB1 affect GSTCD expression

Chronic obstructive pulmonary disease (COPD), the third leading cause of death worldwide, is influenced by genetic factors. The genetic signal rs10516526 in the glutathione S-transferase C-terminal domain containing (GSTCD) gene is a highly significant and reproducible signal associated with lung function and COPD on chromosome 4q24. In this study, comprehensive bioinformatics analyses and experimental verifications were detailly implemented to explore the regulation mechanism of rs10516526 and GSTCD in COPD. The results suggested that low expression of GSTCD was associated with COPD (p = 0.010). And C-Jun and CREB1 transcription factors were found to be essential for the regulation of GSTCD by rs80245547 and rs72673891. Moreover, rs80245547T and rs72673891G had a stronger binding ability to these transcription factors, which may promote the allele-specific long-range enhancer-promoter interactions on GSTCD, thus making COPD less susceptible. Our study provides a new insight into the relationship between rs10516526, GSTCD, and COPD.

Rescuing biologically relevant consensus regions across replicated samples

BACKGROUND

Protein-DNA binding sites of ChIP-seq experiments are identified where the binding affinity is significant based on a given threshold. The choice of the threshold is a trade-off between conservative region identification and discarding weak, but true binding sites.

RESULTS

We rescue weak binding sites using MSPC, which efficiently exploits replicates to lower the threshold required to identify a site while keeping a low false-positive rate, and we compare it to IDR, a widely used post-processing method for identifying highly reproducible peaks across replicates. We observe several master transcription regulators (e.g., SP1 and GATA3) and HDAC2-GATA1 regulatory networks on rescued regions in K562 cell line.

CONCLUSIONS

We argue the biological relevance of weak binding sites and the information they add when rescued by MSPC. An implementation of the proposed extended MSPC methodology and the scripts to reproduce the performed analysis are freely available at https://genometric.github.io/MSPC/ ; MSPC is distributed as a command-line application and an R package available from Bioconductor ( https://doi.org/doi:10.18129/B9.bioc.rmspc ).

A survey on algorithms to characterize transcription factor binding sites

Transcription factors (TFs) are key regulatory proteins that control the transcriptional rate of cells by binding short DNA sequences called transcription factor binding sites (TFBS) or motifs. Identifying and characterizing TFBS is fundamental to understanding the regulatory mechanisms governing the transcriptional state of cells. During the last decades, several experimental methods have been developed to recover DNA sequences containing TFBS. In parallel, computational methods have been proposed to discover and identify TFBS motifs based on these DNA sequences. This is one of the most widely investigated problems in bioinformatics and is referred to as the motif discovery problem. In this manuscript, we review classical and novel experimental and computational methods developed to discover and characterize TFBS motifs in DNA sequences, highlighting their advantages and drawbacks. We also discuss open challenges and future perspectives that could fill the remaining gaps in the field.

Theoretical investigation of functional responses of bio-molecular assembly networks

Cooperative protein-protein and protein-DNA interactions form programmable complex assemblies, often performing non-linear gene regulatory operations involved in signal transductions and cell fate determination. The apparent structure of those complex assemblies is very similar, but their functional response strongly depends on the topology of the protein-DNA interaction networks. Here, we demonstrate how the coordinated self-assembly creates gene regulatory network motifs that corroborate the existence of a precise functional response at the molecular level using thermodynamic and dynamic analyses. Our theoretical and Monte Carlo simulations show that a complex network of interactions can form a decision-making loop, such as feedback and feed-forward circuits, only by a few molecular mechanisms. We characterize each possible network of interactions by systematic variations of free energy parameters associated with the binding among biomolecules and DNA looping. We also find that the higher-order networks exhibit alternative steady states from the stochastic dynamics of each network. We capture this signature by calculating stochastic potentials and attributing their multi-stability features. We validate our findings against the Gal promoter system in yeast cells. Overall, we show that the network topology is vital in phenotype diversity in regulatory circuits.

Single-cell multiomics reveals the complexity of TGFβ signalling to chromatin in iPSC-derived kidney organoids

TGFβ1 plays a regulatory role in the determination of renal cell fate and the progression of renal fibrosis. Here we show an association between SMAD3 and the histone methyltransferase, EZH2, during cell differentiation; ChIP-seq revealed that SMAD3 and EZH2 co-occupy the genome in iPSCs and in iPSC-derived nephron progenitors. Through integration of single cell gene expression and epigenome profiling, we identified de novo ACTA2^+ve/POSTN^+ve myofibroblasts in kidney organoids treated with TGFβ1, characterised by increased SMAD3-dependent cis chromatin accessibility and gene expression associated with fibroblast activation. We have identified fibrosis-associated regulons characterised by enrichment of SMAD3, AP1, the ETS family of transcription factors, and NUAK1, CREB3L1, and RARG, corresponding to enriched motifs at accessible loci identified by scATACseq. Treatment with the EZH2 specific inhibitor GSK343, blocked SMAD3-dependent cis co-accessibility and inhibited myofibroblast activation. This mechanism, through which TGFβ signals directly to chromatin, represents a critical determinant of fibrotic, differentiated states.

iPro-WAEL: a comprehensive and robust framework for identifying promoters in multiple species

Promoters are consensus DNA sequences located near the transcription start sites and they play an important role in transcription initiation. Due to their importance in biological processes, the identification of promoters is significantly important for characterizing the expression of the genes. Numerous computational methods have been proposed to predict promoters. However, it is difficult for these methods to achieve satisfactory performance in multiple species. In this study, we propose a novel weighted average ensemble learning model, termed iPro-WAEL, for identifying promoters in multiple species, including Human, Mouse, E.coli, Arabidopsis, B.amyloliquefaciens, B.subtilis and R.capsulatus. Extensive benchmarking experiments illustrate that iPro-WAEL has optimal performance and is superior to the current methods in promoter prediction. The experimental results also demonstrate a satisfactory prediction ability of iPro-WAEL on cross-cell lines, promoters annotated by other methods and distinguishing between promoters and enhancers. Moreover, we identify the most important transcription factor binding site (TFBS) motif in promoter regions to facilitate the study of identifying important motifs in the promoter regions. The source code of iPro-WAEL is freely available at https://github.com/HaoWuLab-Bioinformatics/iPro-WAEL.

Complex chromosomal rearrangements induced by transposons in maize

Eukaryotic genomes are large and complex, and gene expression can be affected by multiple regulatory elements and their positions within the dynamic chromatin architecture. Transposable elements are known to play important roles in genome evolution, yet questions remain as to how transposable elements alter genome structure and affect gene expression. Previous studies have shown that genome rearrangements can be induced by Reversed Ends Transposition involving termini of Activator and related transposable elements in maize and other plants. Here, we show that complex alleles can be formed by the rapid and progressive accumulation of Activator-induced duplications and rearrangements. The p1 gene enhancer in maize can induce ectopic expression of the nearby p2 gene in pericarp tissue when placed near it via different structural rearrangements. By screening for p2 expression, we identified and studied 5 cases in which multiple sequential transposition events occurred and increased the p1 enhancer copy number. We see active p2 expression due to multiple copies of the p1 enhancer present near p2 in all 5 cases. The p1 enhancer effects are confirmed by the observation that loss of p2 expression is correlated with transposition-induced excision of the p1 enhancers. We also performed a targeted Chromosome Conformation Capture experiment to test the physical interaction between the p1 enhancer and p2 promoter region. Together, our results show that transposon-induced rearrangements can accumulate rapidly and progressively increase genetic variation important for genomic evolution.

Identification of mechanism of cancer-cell-specific reactivation of hTERT offers therapeutic opportunities for blocking telomerase specifically in human colorectal cancer

Transcriptional reactivation of hTERT is the limiting step in tumorigenesis. While mutations in hTERT promoter present in 19% of cancers are recognized as key drivers of hTERT reactivation, mechanisms by which wildtype hTERT (WT-hTERT) promoter is reactivated, in majority of human cancers, remain unknown. Using primary colorectal cancers (CRC) we identified Tert INTeracting region 2 (T-INT2), the critical chromatin region essential for reactivating WT-hTERT promoter in CRCs. Elevated β-catenin and JunD level in CRC facilitates chromatin interaction between hTERT promoter and T-INT2 that is necessary to turn on hTERTexpression. Pharmacological screens uncovered salinomycin, which inhibits JunD mediated hTERT-T-INT2 interaction that is required for the formation of a stable transcription complex on the hTERT promoter. Our results showed for the first time how known CRC alterations, such as APC, lead to WT-hTERT promoter reactivation during stepwise-tumorigenesis and provide a new perspective for developing cancer-specific drugs.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Genome-wide screens identify specific drivers of mutant hTERT promoters

Mutations in hTERT promoter are seen in over 19% of human cancers, irrespective of the cancer type. Understanding the molecular players that regulate Mut-hTERT promoters may help the design of effective targeting strategies to inhibit telomerase reactivation specifically in cancer cells. Our work uses genome-wide functional screens to identify 30 specific regulators of Mut-hTERT promoters. These candidates identified from the screening serve as an excellent resource to understand how telomerase is reactivated and as targets for making inhibitors to telomerase, a key driver of cancer.

Cancer-specific hTERT promoter mutations reported in 19% of cancers result in enhanced telomerase activity. Understanding the distinctions between transcriptional regulation of wild-type (WT) and mutant (Mut) hTERT promoters may open up avenues for development of inhibitors which specially block hTERT expression in cancer cells. To comprehensively identify physiological regulators of WT- or Mut-hTERT promoters, we generated several isogenic reporter cells driven by endogenous hTERT loci. Genome-wide CRISPR-Cas9 and small interfering RNA screens using these isogenic reporter lines identified specific regulators of Mut-hTERT promoters. We validate and characterize one of these hits, namely, MED12, a kinase subunit of mediator complex. We demonstrate that MED12 specifically drives expression of hTERT from the Mut-hTERT promoter by mediating long-range chromatin interaction between the proximal Mut-hTERT promoter and T-INT1 distal regulatory region 260 kb upstream. Several hits identified in our screens could serve as potential therapeutic targets, inhibition of which may specifically block Mut-hTERT promoter driven telomerase reactivation in cancers.

ZNF382 controls mouse neuropathic pain via silencer-based epigenetic inhibition of Cxcl13 in DRG neurons

Nerve injury-induced changes of gene expression in dorsal root ganglion (DRG) are critical for neuropathic pain genesis. However, how these changes occur remains elusive. Here we report the down-regulation of zinc finger protein 382 (ZNF382) in injured DRG neurons after nerve injury. Rescuing this down-regulation attenuates nociceptive hypersensitivity. Conversely, mimicking this down-regulation produces neuropathic pain symptoms, which are alleviated by C-X-C motif chemokine 13 (CXCL13) knockdown or its receptor CXCR5 knockout. Mechanistically, an identified cis-acting silencer at distal upstream of the Cxcl13 promoter suppresses Cxcl13 transcription via binding to ZNF382. Blocking this binding or genetically deleting this silencer abolishes the ZNF382 suppression on Cxcl13 transcription and impairs ZNF382-induced antinociception. Moreover, ZNF382 down-regulation disrupts the repressive epigenetic complex containing histone deacetylase 1 and SET domain bifurcated 1 at the silencer-promoter loop, resulting in Cxcl13 transcriptional activation. Thus, ZNF382 down-regulation is required for neuropathic pain likely through silencer-based epigenetic disinhibition of CXCL13, a key neuropathic pain player, in DRG neurons.

Predict long-range enhancer regulation based on protein-protein interactions between transcription factors

Long-range regulation by distal enhancers plays critical roles in cell-type specific transcriptional programs. Computational predictions of genome-wide enhancer-promoter interactions are still challenging due to limited accuracy and the lack of knowledge on the molecular mechanisms. Based on recent biological investigations, the protein-protein interactions (PPIs) between transcription factors (TFs) have been found to participate in the regulation of chromatin loops. Therefore, we developed a novel predictive model for cell-type specific enhancer-promoter interactions by leveraging the information of TF PPI signatures. Evaluated by a series of rigorous performance comparisons, the new model achieves superior performance over other methods. The model also identifies specific TF PPIs that may mediate long-range regulatory interactions, revealing new mechanistic understandings of enhancer regulation. The prioritized TF PPIs are associated with genes in distinct biological pathways, and the predicted enhancer-promoter interactions are strongly enriched with cis-eQTLs. Most interestingly, the model discovers enhancer-mediated trans-regulatory links between TFs and genes, which are significantly enriched with trans-eQTLs. The new predictive model, along with the genome-wide analyses, provides a platform to systematically delineate the complex interplay among TFs, enhancers and genes in long-range regulation. The novel predictions also lead to mechanistic interpretations of eQTLs to decode the genetic associations with gene expression.

Chromatin Alterations in Neurological Disorders and Strategies of (Epi)Genome Rescue

Neurological disorders (NDs) comprise a heterogeneous group of conditions that affect the function of the nervous system. Often incurable, NDs have profound and detrimental consequences on the affected individuals' lives. NDs have complex etiologies but commonly feature altered gene expression and dysfunctions of the essential chromatin-modifying factors. Hence, compounds that target DNA and histone modification pathways, the so-called epidrugs, constitute promising tools to treat NDs. Yet, targeting the entire epigenome might reveal insufficient to modify a chosen gene expression or even unnecessary and detrimental to the patients' health. New technologies hold a promise to expand the clinical toolkit in the fight against NDs. (Epi)genome engineering using designer nucleases, including CRISPR-Cas9 and TALENs, can potentially help restore the correct gene expression patterns by targeting a defined gene or pathway, both genetically and epigenetically, with minimal off-target activity. Here, we review the implication of epigenetic machinery in NDs. We outline syndromes caused by mutations in chromatin-modifying enzymes and discuss the functional consequences of mutations in regulatory DNA in NDs. We review the approaches that allow modifying the (epi)genome, including tools based on TALENs and CRISPR-Cas9 technologies, and we highlight how these new strategies could potentially change clinical practices in the treatment of NDs.

Models behind the mystery of establishing enhancer-promoter interactions

Enhancers and promoters are transcriptional regulatory elements whose facilitated interactions increase gene expression. Enhancer DNA sequences can be located far away from the promoter sequences that they regulate. Currently, the mechanism facilitating the establishment of enhancer-promoter interactions remains unclear. However, mutations causing errors in these interactions have been linked to cancer and disease, further conveying the need to understand the full mechanism. This review discusses multiple models that have been proposed to describe how enhancers go the distance to interact with promoters. Evidence supporting loop formation models is reviewed in addition to more complex hypotheses involving aspects of 3D chromatin organization and phase separation.

The Causes and Consequences of Spatial Organization of the Genome in Regulation of Gene Expression

Regulation of gene expression in time, space and quantity is orchestrated by the functional interplay of cis-acting elements and trans-acting factors. Our current view postulates that transcription factors recognize enhancer DNA and read the transcriptional regulatory code by cooperative DNA binding to specific DNA motifs, thus instructing the recruitment of transcriptional regulatory complexes forming a plethora of higher-ordered multi-protein-DNA and protein-protein complexes. Here, we reviewed the formation of multi-dimensional chromatin assemblies implicated in gene expression with emphasis on the regulatory role of enhancer hubs as coordinators of stochastic gene expression. Enhancer hubs contain many interacting regulatory elements and represent a remarkably dynamic and heterogeneous network of multivalent interactions. A functional consequence of such complex interaction networks could be that individual enhancers function synergistically to ensure coordination, tight control and robustness in regulation of expression of spatially connected genes. In this review, we discuss fundamental paradigms of such inter- and intra- chromosomal associations both in the context of immune-related genes and beyond.

The loopometer: a quantitative in vivo assay for DNA-looping proteins

Proteins that can bring together separate DNA sites, either on the same or on different DNA molecules, are critical for a variety of DNA-based processes. However, there are no general and technically simple assays to detect proteins capable of DNA looping in vivo nor to quantitate their in vivo looping efficiency. Here, we develop a quantitative in vivo assay for DNA-looping proteins in Escherichia coli that requires only basic DNA cloning techniques and a LacZ assay. The assay is based on loop assistance, where two binding sites for the candidate looping protein are inserted internally to a pair of operators for the E. coli LacI repressor. DNA looping between the sites shortens the effective distance between the lac operators, increasing LacI looping and strengthening its repression of a lacZ reporter gene. Analysis based on a general model for loop assistance enables quantitation of the strength of looping conferred by the protein and its binding sites. We use this ‘loopometer’ assay to measure DNA looping for a variety of bacterial and phage proteins.

Engineering three-dimensional genome folding

Animal genomes are partitioned and folded at various scales that contribute distinctly to nuclear processes. While structural features have been disrupted either globally or at select loci in loss-of-function studies, gain-of-function studies that probe the role of genome architecture have lagged behind. Here we examine recent advances in experimentally creating chromatin loops, contact domains, boundaries and compartments. Furthermore, we explore parallels between this emerging theme and natural evolution of mammalian genomes with increasing architectural complexity. Finally, we provide a perspective on how insights arising from recent gain-of-function studies may inform future endeavors toward engineering the three-dimensional genome.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Single cell transcriptional and chromatin accessibility profiling redefine cellular heterogeneity in the adult human kidney

The integration of single cell transcriptome and chromatin accessibility datasets enables a deeper understanding of cell heterogeneity. We performed single nucleus ATAC (snATAC-seq) and RNA (snRNA-seq) sequencing to generate paired, cell-type-specific chromatin accessibility and transcriptional profiles of the adult human kidney. We demonstrate that snATAC-seq is comparable to snRNA-seq in the assignment of cell identity and can further refine our understanding of functional heterogeneity in the nephron. The majority of differentially accessible chromatin regions are localized to promoters and a significant proportion are closely associated with differentially expressed genes. Cell-type-specific enrichment of transcription factor binding motifs implicates the activation of NF-κB that promotes VCAM1 expression and drives transition between a subpopulation of proximal tubule epithelial cells. Our multi-omics approach improves the ability to detect unique cell states within the kidney and redefines cellular heterogeneity in the proximal tubule and thick ascending limb.

Using functional genomics to advance the understanding of psoriatic arthritis

Psoriatic arthritis (PsA) is a complex disease where susceptibility is determined by genetic and environmental risk factors. Clinically, PsA involves inflammation of the joints and the skin, and, if left untreated, results in irreversible joint damage. There is currently no cure and the few treatments available to alleviate symptoms do not work in all patients. Over the past decade, genome-wide association studies (GWAS) have uncovered a large number of disease-associated loci but translating these findings into functional mechanisms and novel targets for therapeutic use is not straightforward. Most variants have been predicted to affect primarily long-range regulatory regions such as enhancers. There is now compelling evidence to support the use of chromatin conformation analysis methods to discover novel genes that can be affected by disease-associated variants. Here, we will review the studies published in the field that have given us a novel understanding of gene regulation in the context of functional genomics and how this relates to the study of PsA and its underlying disease mechanism.

BET inhibition disrupts transcription but retains enhancer-promoter contact

Enhancers are DNA sequences that enable complex temporal and tissue-specific regulation of genes in higher eukaryotes. Although it is not entirely clear how enhancer-promoter interactions can increase gene expression, this proximity has been observed in multiple systems at multiple loci and is thought to be essential for the maintenance of gene expression. Bromodomain and Extra-Terminal domain (BET) and Mediator proteins have been shown capable of forming phase condensates and are thought to be essential for super-enhancer function. Here, we show that targeting of cells with inhibitors of BET proteins or pharmacological degradation of BET protein Bromodomain-containing protein 4 (BRD4) has a strong impact on transcription but very little impact on enhancer-promoter interactions. Dissolving phase condensates reduces BRD4 and Mediator binding at enhancers and can also strongly affect gene transcription, without disrupting enhancer-promoter interactions. These results suggest that activation of transcription and maintenance of enhancer-promoter interactions are separable events. Our findings further indicate that enhancer-promoter interactions are not dependent on high levels of BRD4 and Mediator, and are likely maintained by a complex set of factors including additional activator complexes and, at some sites, CTCF and cohesin.

Trans-Ancestral Fine-Mapping and Epigenetic Annotation as Tools to Delineate Functionally Relevant Risk Alleles at IKZF1 and IKZF3 in Systemic Lupus Erythematosus

Background: Prioritizing tag-SNPs carried on extended risk haplotypes at susceptibility loci for common disease is a challenge. Methods: We utilized trans-ancestral exclusion mapping to reduce risk haplotypes at IKZF1 and IKZF3 identified in multiple ancestries from SLE GWAS and ImmunoChip datasets. We characterized functional annotation data across each risk haplotype from publicly available datasets including ENCODE, RoadMap Consortium, PC Hi-C data from 3D genome browser, NESDR NTR conditional eQTL database, GeneCards Genehancers and TF (transcription factor) binding sites from Haploregv4. Results: We refined the 60 kb associated haplotype upstream of IKZF1 to just 12 tag-SNPs tagging a 47.7 kb core risk haplotype. There was preferential enrichment of DNAse I hypersensitivity and H3K27ac modification across the 3′ end of the risk haplotype, with four tag-SNPs sharing allele-specific TF binding sites with promoter variants, which are eQTLs for IKZF1 in whole blood. At IKZF3, we refined a core risk haplotype of 101 kb (27 tag-SNPs) from an initial extended haplotype of 194 kb (282 tag-SNPs), which had widespread DNAse I hypersensitivity, H3K27ac modification and multiple allele-specific TF binding sites. Dimerization of Fox family TFs bound at the 3′ and promoter of IKZF3 may stabilize chromatin looping across the locus. Conclusions: We combined trans-ancestral exclusion mapping and epigenetic annotation to identify variants at both IKZF1 and IKZF3 with the highest likelihood of biological relevance. The approach will be of strong interest to other complex trait geneticists seeking to attribute biological relevance to risk alleles on extended risk haplotypes in their disease of interest.

Transcriptional regulators of the Golli/myelin basic protein locus integrate additive and stealth activities

Myelin is composed of plasma membrane spirally wrapped around axons and compacted into dense sheaths by myelin-associated proteins. Myelin is elaborated by neuroepithelial derived oligodendrocytes in the central nervous system (CNS) and by neural crest derived Schwann cells in the peripheral nervous system (PNS). While some myelin proteins accumulate in only one lineage, myelin basic protein (Mbp) is expressed in both. Overlapping the Mbp gene is Golli, a transcriptional unit that is expressed widely both within and beyond the nervous system. A super-enhancer domain within the Golli/Mbp locus contains multiple enhancers shown previously to drive reporter construct expression specifically in oligodendrocytes or Schwann cells. In order to determine the contribution of each enhancer to the Golli/Mbp expression program, and to reveal if functional interactions occur among them, we derived mouse lines in which they were deleted, either singly or in different combinations, and relative mRNA accumulation was measured at key stages of early development and at maturity. Although super-enhancers have been shown previously to facilitate interaction among their component enhancers, the enhancers investigated here demonstrated largely additive relationships. However, enhancers demonstrating autonomous activity strictly in one lineage, when missing, were found to significantly reduce output in the other, thus revealing cryptic "stealth" activity. Further, in the absence of a key oligodendrocyte enhancer, Golli accumulation was markedly and uniformly attenuated in all cell types investigated. Our observations suggest a model in which enhancer-mediated DNA-looping and potential super-enhancer properties underlie Golli/Mbp regulatory organization.

Identification and Analysis of p53-Regulated Enhancers in Hepatic Carcinoma

Enhancers can act as cis-regulatory elements to control transcriptional regulation by recruiting transcription factors (TFs) in a distance and orientation-independent manner. However, it is still unclear how p53 participates in the enhancer network as TF in hepatic carcinoma under the condition of DNA damage. A total of 14,286 active enhancers were identified through the integration of stable and unstable enhancer RNAs (eRNAs) captured by CAGE and GRO-seq, respectively. Furthermore, 218 p53-bound enhancers (Enh_p53) were identified by analyzing p53 ChIP-seq in HepG2 cells after DNA damage. The results showed that the enhancer expression and histone markers of enhancers (H3K4me1, H3K4me2, H3K4me3, H3K9ac, and H3K27ac) revealed significantly higher level on Enh_p53 than Enh_no−p53 which suggested that p53 participated in regulating enhancer activity and chromatin structure. By analyzing 124 TFs ChIP-seq from ENCODE, 93 TFs were found significantly enriched on Enh_p53 such as GATA4, YY1, and CTCF, indicating p53 may co-regulate enhancers with TFs participation. Moreover, significantly differentially expressed 438 miRNAs and 1,264 mRNAs were identified by analyzing small RNA-seq and RNA-seq, and 26 Enh_p53-miRNAs and 145 Enh_p53-mRNA interactions were identified by the integration of 3D genome data and genomic distance. The functional enrichment analysis showed that these miRNA targets and mRNAs were significantly involved in tumor biological processes and signaling pathways such as DNA replication, p53 signaling pathway, hepatitis B, focal adhesion, etc. The above results indicated that p53 participated in regulating enhancer network in hepatic carcinoma and Enh_p53 exhibited significantly different characteristics with Enh_no−p53.

HiChIP-Peaks: a HiChIP peak calling algorithm

MOTIVATION

HiChIP is a powerful tool to interrogate 3D chromatin organization. Current tools to analyse chromatin looping mechanisms using HiChIP data require the identification of loop anchors to work properly. However, current approaches to discover these anchors from HiChIP data are not satisfactory, having either a very high false discovery rate or strong dependence on sequencing depth. Moreover, these tools do not allow quantitative comparison of peaks across different samples, failing to fully exploit the information available from HiChIP datasets.

RESULTS

We develop a new tool based on a representation of HiChIP data centred on the re-ligation sites to identify peaks from HiChIP datasets, which can subsequently be used in other tools for loop discovery. This increases the reliability of these tools and improves recall rate as sequencing depth is reduced. We also provide a method to count reads mapping to peaks across samples, which can be used for differential peak analysis using HiChIP data.

AVAILABILITY AND IMPLEMENTATION

HiChIP-Peaks is freely available at https://github.com/ChenfuShi/HiChIP_peaks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CDK-Mediator and FBXL19 prime developmental genes for activation by promoting atypical regulatory interactions

Appropriate developmental gene regulation relies on the capacity of gene promoters to integrate inputs from distal regulatory elements, yet how this is achieved remains poorly understood. In embryonic stem cells (ESCs), a subset of silent developmental gene promoters are primed for activation by FBXL19, a CpG island binding protein, through its capacity to recruit CDK-Mediator. How mechanistically these proteins function together to prime genes for activation during differentiation is unknown. Here we discover that in mouse ESCs FBXL19 and CDK-Mediator support long-range interactions between silent gene promoters that rely on FBXL19 for their induction during differentiation and gene regulatory elements. During gene induction, these distal regulatory elements behave in an atypical manner, in that the majority do not acquire histone H3 lysine 27 acetylation and no longer interact with their target gene promoter following gene activation. Despite these atypical features, we demonstrate by targeted deletions that these distal elements are required for appropriate gene induction during differentiation. Together these discoveries demonstrate that CpG-island associated gene promoters can prime genes for activation by communicating with atypical distal gene regulatory elements to achieve appropriate gene expression.

Chromatin topology reorganization and transcription repression by PML-RARα in acute promyeloid leukemia

BACKGROUND

Acute promyeloid leukemia (APL) is characterized by the oncogenic fusion protein PML-RARα, a major etiological agent in APL. However, the molecular mechanisms underlying the role of PML-RARα in leukemogenesis remain largely unknown.

RESULTS

Using an inducible system, we comprehensively analyze the 3D genome organization in myeloid cells and its reorganization after PML-RARα induction and perform additional analyses in patient-derived APL cells with native PML-RARα. We discover that PML-RARα mediates extensive chromatin interactions genome-wide. Globally, it redefines the chromatin topology of the myeloid genome toward a more condensed configuration in APL cells; locally, it intrudes RNAPII-associated interaction domains, interrupts myeloid-specific transcription factors binding at enhancers and super-enhancers, and leads to transcriptional repression of genes critical for myeloid differentiation and maturation.

CONCLUSIONS

Our results not only provide novel topological insights for the roles of PML-RARα in transforming myeloid cells into leukemia cells, but further uncover a topological framework of a molecular mechanism for oncogenic fusion proteins in cancers.

Promoter-Proximal Chromatin Domain Insulator Protein BEAF Mediates Local and Long-Range Communication with a Transcription Factor and Directly Activates a Housekeeping Promoter in Drosophila

BEAF (Boundary Element-Associated Factor) was originally identified as a Drosophila melanogaster chromatin domain insulator-binding protein, suggesting a role in gene regulation through chromatin organization and dynamics. Genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, suggesting a role in promoter function. This would be a nontraditional role for an insulator-binding protein. To gain insight into molecular mechanisms of BEAF function, we identified interacting proteins using yeast two-hybrid assays. Here, we focus on the transcription factor Serendipity δ (Sry-δ). Interactions were confirmed in pull-down experiments using bacterially expressed proteins, by bimolecular fluorescence complementation, and in a genetic assay in transgenic flies. Sry-δ interacted with promoter-proximal BEAF both when bound to DNA adjacent to BEAF or > 2-kb upstream to activate a reporter gene in transient transfection experiments. The interaction between BEAF and Sry-δ was detected using both a minimal developmental promoter (y) and a housekeeping promoter (RpS12), while BEAF alone strongly activated the housekeeping promoter. These two functions for BEAF implicate it in playing a direct role in gene regulation at hundreds of BEAF-associated promoters.

PAN RNA: transcriptional exhaust from a viral engine

Kaposi’s sarcoma-associated herpesvirus (KSHV), also designated human herpesvirus 8 (HHV-8), has been linked to Kaposi’s sarcoma, as well as to primary effusion lymphoma (PEL), and a subset of multicentric Castleman’s disease. KSHV genomes are maintained as episomes within infected cells and the virus exhibits a biphasic life cycle consisting of a life-long latent phase during which only a few viral genes are expressed and no viral progeny are produced and a transient lytic reactivation phase, in which a full repertoire of ~ 80 lytic genes are activated in a temporally regulated manner culminating in the release of new virions. Lytic replication is initiated by a single viral protein, K-Rta (ORF50), which activates more than 80 viral genes from multiple resident viral episomes (i.e., viral chromosomes). One of the major targets of K-Rta is a long non-coding nuclear RNA, PAN RNA (polyadenylated nuclear RNA), a lncRNA that accumulates to exceedingly high levels in the nucleus during viral reactivation. K-Rta directly binds to the PAN RNA promoter and robustly activates PAN RNA expression. Although PAN RNA has been known for over 20 years, its role in viral replication is still incompletely understood. In this perspective, we will briefly review the current understanding of PAN RNA and then describe our current working model of this RNA. The model is based on our observations concerning events that occur during KSHV lytic reactivation including (i) a marked accumulation of RNA Pol II at the PAN promoter, (ii) genomic looping emanating from the PAN locus, (iii) interaction of a second viral lytic protein (ORF57) with K-Rta, PAN RNA and RNA Pol II, (iv) the essential requirement for PAN RNA expression in cis for optimal transcriptional execution needed for the entire lytic program, and (v) ORF57 recruitment of RNA Pol II to the PAN genomic locus. Together our results generate a model in which the PAN locus serves as a hub for sequestration/trapping of the cellular transcriptional machinery proximal to viral episomes. Sequestration at the PAN locus facilitates high levels of viral transcription throughout the viral genome during lytic replication. ORF57 acts as a transcription-dependent transactivator at the PAN locus by binding to both Rta and PAN to locally trap RNA Pol II. The resulting accumulation of high levels of nuclear PAN RNA created by this process is an inducible enhancer-derived (eRNA) by-product that litters the infected cell nucleus.

Histone ChIP-Seq identifies differential enhancer usage during chondrogenesis as critical for defining cell-type specificity

Epigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterized the epigenome during the in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next‐generation sequencing (ChIP‐seq) was used to assess a range of N‐terminal posttranscriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3, and H3K36me3) in both hMSCs and differentiated chondrocytes. Chromatin states were characterized using histone ChIP‐seq and cis‐regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis‐related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de‐methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this study with epigenomes of cell‐types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell‐types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers. Luciferase reporter assays confirmed that chondrocyte enhancers characterized in this study exhibited enhancer activity which may be modulated by DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross‐talk between different epigenetic mechanisms during chondrocyte differentiation.

Whole-Genome Sequencing of Finnish Type 1 Diabetic Siblings Discordant for Kidney Disease Reveals DNA Variants associated with Diabetic Nephropathy

BACKGROUND

Several genetic susceptibility loci associated with diabetic nephropathy have been documented, but no causative variants implying novel pathogenetic mechanisms have been elucidated.

METHODS

We carried out whole-genome sequencing of a discovery cohort of Finnish siblings with type 1 diabetes who were discordant for the presence (case) or absence (control) of diabetic nephropathy. Controls had diabetes without complications for 15-37 years. We analyzed and annotated variants at genome, gene, and single-nucleotide variant levels. We then replicated the associated variants, genes, and regions in a replication cohort from the Finnish Diabetic Nephropathy study that included 3531 unrelated Finns with type 1 diabetes.

RESULTS

We observed protein-altering variants and an enrichment of variants in regions associated with the presence or absence of diabetic nephropathy. The replication cohort confirmed variants in both regulatory and protein-coding regions. We also observed that diabetic nephropathy-associated variants, when clustered at the gene level, are enriched in a core protein-interaction network representing proteins essential for podocyte function. These genes include protein kinases (protein kinase C isoforms ε and ι) and protein tyrosine kinase 2.

CONCLUSIONS

Our comprehensive analysis of a diabetic nephropathy cohort of siblings with type 1 diabetes who were discordant for kidney disease points to variants and genes that are potentially causative or protective for diabetic nephropathy. This includes variants in two isoforms of the protein kinase C family not previously linked to diabetic nephropathy, adding support to previous hypotheses that the protein kinase C family members play a role in diabetic nephropathy and might be attractive therapeutic targets.

The liquid nucleome - phase transitions in the nucleus at a glance

Cells organize membrane-less internal compartments through a process called liquid-liquid phase separation (LLPS) to create chemically distinct compartments, referred to as condensates, which emerge from interactions among biological macromolecules. These condensates include various cytoplasmic structures such as P-granules and stress granules. However, an even wider array of condensates subcompartmentalize the cell nucleus, forming liquid-like structures that range from nucleoli and Cajal bodies to nuclear speckles and gems. Phase separation provides a biophysical assembly mechanism underlying this non-covalent form of fluid compartmentalization and functionalization. In this Cell Science at a Glance article and the accompanying poster, we term these phase-separated liquids that organize the nucleus the liquid nucleome; we discuss examples of biological phase transitions in the nucleus, how the cell utilizes biophysical aspects of phase separation to form and regulate condensates, and suggest interpretations for the role of phase separation in nuclear organization and function.

Neurobiological functions of transcriptional enhancers

Transcriptional enhancers are regulatory DNA elements that underlie the specificity and dynamic patterns of gene expression. Over the past decade, large-scale functional genomics projects have driven transformative progress in our understanding of enhancers. These data have relevance for identifying mechanisms of gene regulation in the CNS, elucidating the function of non-coding regulatory sequences in neurobiology and linking sequence variation within enhancers to genetic risk for neurological and psychiatric disorders. However, the sheer volume and complexity of genomic data presents a challenge to interpreting enhancer function in normal and pathogenic neurobiological processes. Here, to advance the application of genome-scale enhancer data, we offer a primer on current models of enhancer function in the CNS, we review how enhancers regulate gene expression across the neuronal lifespan, and we suggest how emerging findings regarding the role of non-coding sequence variation offer opportunities for understanding brain disorders and developing new technologies for neuroscience.

Bioinformatics Discovery of Putative Enhancers within Mouse Odorant Receptor Gene Clusters

Olfactory neuronal function depends on the expression and proper regulation of odorant receptor (OR) genes. Previous studies have identified 54 putative intergenic enhancers within or flanking 40 mouse OR clusters. At least 2 of these putative enhancers have been shown to regulate the expression of a small subset of proximal OR genes. In recognition of the large size of the mouse OR gene family (~1400 OR genes distributed across multiple chromosomal loci), it is likely that there remain many additional not-as-yet discovered OR enhancers. We utilized 23 of the previously identified enhancers as a training set (TS) and designed an algorithm that combines a broad range of epigenetic criteria (histone-3-lysine-4 monomethylation, histone-3-lysine-79 trimethylation, histone-3-lysine-27 acetylation, and DNase hypersensitivity) and genetic criteria (cross-species sequence conservation and transcription-factor binding site enrichment) to more broadly search OR gene clusters for additional candidates. We identified 181 new candidate enhancers located at 58 (of 68) mouse OR loci, including 25 new candidates identified by stringent search criteria whose signal strengths are not significantly different from the 23 previously characterized OR enhancers used as the TS. Additionally, we compared OR enhancer versus generic enhancer features in order to evaluate likelihoods that new enhancer candidates specifically function in OR regulation. We found that features distinguishing OR-specific function are significantly more evident for enhancer candidates located within OR clusters as compared with those in flanking regions.

Pax3 cooperates with Ldb1 to direct local chromosome architecture during myogenic lineage specification

Chromatin looping allows enhancer-bound regulatory factors to influence transcription. Large domains, referred to as topologically associated domains, participate in genome organization. However, the mechanisms underlining interactions within these domains, which control gene expression, are not fully understood. Here we report that activation of embryonic myogenesis is associated with establishment of long-range chromatin interactions centered on Pax3-bound loci. Using mass spectrometry and genomic studies, we identify the ubiquitously expressed LIM-domain binding protein 1 (Ldb1) as the mediator of looping interactions at a subset of Pax3 binding sites. Ldb1 is recruited to Pax3-bound elements independently of CTCF-Cohesin, and is necessary for efficient deposition of H3K4me1 at these sites and chromatin looping. When Ldb1 is deleted in Pax3-expressing cells in vivo, specification of migratory myogenic progenitors is severely impaired. These results highlight Ldb1 requirement for Pax3 myogenic activity and demonstrate how transcription factors can promote formation of sub-topologically associated domain interactions involved in lineage specification.

A combination of transcription factors mediates inducible interchromosomal contacts

The genome forms specific three-dimensional contacts in response to cellular or environmental conditions. However, it remains largely unknown which proteins specify and mediate such contacts. Here we describe an assay, MAP-C (Mutation Analysis in Pools by Chromosome conformation capture), that simultaneously characterizes the effects of hundreds of cis or trans-acting mutations on a chromosomal contact. Using MAP-C, we show that inducible interchromosomal pairing between HAS1pr-TDA1pr alleles in saturated cultures of Saccharomyces yeast is mediated by three transcription factors, Leu3, Sdd4 (Ypr022c), and Rgt1. The coincident, combined binding of all three factors is strongest at the HAS1pr-TDA1pr locus and is also specific to saturated conditions. We applied MAP-C to further explore the biochemical mechanism of these contacts, and find they require the structured regulatory domain of Rgt1, but no known interaction partners of Rgt1. Altogether, our results demonstrate MAP-C as a powerful method for dissecting the mechanistic basis of chromosome conformation.

Hierarchical cooperation of transcription factors from integration analysis of DNA sequences, ChIP-Seq and ChIA-PET data

Background

Chromosomal architecture, which is constituted by chromatin loops, plays an important role in cellular functions. Gene expression and cell identity can be regulated by the chromatin loop, which is formed by proximal or distal enhancers and promoters in linear DNA (1D). Enhancers and promoters are fundamental non-coding elements enriched with transcription factors (TFs) to form chromatin loops. However, the specific cooperation of TFs involved in forming chromatin loops is not fully understood.

Results

Here, we proposed a method for investigating the cooperation of TFs in four cell lines by the integrative analysis of DNA sequences, ChIP-Seq and ChIA-PET data. Results demonstrate that the interaction of enhancers and promoters is a hierarchical and dynamic complex process with cooperative interactions of different TFs synergistically regulating gene expression and chromatin structure. The TF cooperation involved in maintaining and regulating the chromatin loop of cells can be regulated by epigenetic factors, such as other TFs and DNA methylation.

Conclusions

Such cooperation among TFs provides the potential features that can affect chromatin’s 3D architecture in cells. The regulation of chromatin 3D organization and gene expression is a complex process associated with the hierarchical and dynamic prosperities of TFs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5535-2) contains supplementary material, which is available to authorized users.

BRCA1 mutations attenuate super-enhancer function and chromatin looping in haploinsufficient human breast epithelial cells

Background

BRCA1-associated breast cancer originates from luminal progenitor cells. BRCA1 functions in multiple biological processes, including double-strand break repair, replication stress suppression, transcriptional regulation, and chromatin reorganization. While non-malignant cells carrying cancer-predisposing BRCA1 mutations exhibit increased genomic instability, it remains unclear whether BRCA1 haploinsufficiency affects transcription and chromatin dynamics in breast epithelial cells.

Methods

H3K27ac-associated super-enhancers were compared in primary breast epithelial cells from BRCA1 mutation carriers (BRCA1^mut/+) and non-carriers (BRCA1^+/+). Non-tumorigenic MCF10A breast epithelial cells with engineered BRCA1 haploinsufficiency were used to confirm the H3K27ac changes. The impact of BRCA1 mutations on enhancer function and enhancer-promoter looping was assessed in MCF10A cells.

Results

Here, we show that primary mammary epithelial cells from women with BRCA1 mutations display significant loss of H3K27ac-associated super-enhancers. These BRCA1-dependent super-enhancers are enriched with binding motifs for the GATA family. Non-tumorigenic BRCA1^mut/+ MCF10A cells recapitulate the H3K27ac loss. Attenuated histone mark and enhancer activity in these BRCA1^mut/+ MCF10A cells can be partially restored with wild-type BRCA1. Furthermore, chromatin conformation analysis demonstrates impaired enhancer-promoter looping in BRCA1^mut/+ MCF10A cells.

Conclusions

H3K27ac-associated super-enhancer loss is a previously unappreciated functional deficiency in ostensibly normal BRCA1 mutation-carrying breast epithelium. Our findings offer new mechanistic insights into BRCA1 mutation-associated transcriptional and epigenetic abnormality in breast epithelial cells and tissue/cell lineage-specific tumorigenesis.

Electronic supplementary material

The online version of this article (10.1186/s13058-019-1132-1) contains supplementary material, which is available to authorized users.