DNA looping induced by a transcriptional enhancer in vivo

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.

Localization of the Drosophila pioneer factor GAF to subnuclear foci is driven by DNA binding and required to silence satellite repeat expression

The eukaryotic genome is organized to enable the precise regulation of gene expression. This organization is established as the embryo transitions from a fertilized gamete to a totipotent zygote. To understand the factors and processes that drive genomic organization, we focused on the pioneer factor GAGA factor (GAF) that is required for early development in Drosophila. GAF transcriptionally activates the zygotic genome and is localized to subnuclear foci. This non-uniform distribution is driven by binding to highly abundant GA repeats. At GA repeats, GAF is necessary to form heterochromatin and silence transcription. Thus, GAF is required to establish both active and silent regions. We propose that foci formation enables GAF to have opposing transcriptional roles within a single nucleus. Our data support a model in which the subnuclear concentration of transcription factors acts to organize the nucleus into functionally distinct domains essential for the robust regulation of gene expression.

GAGA-associated factor fosters loop formation in the Drosophila genome

The impact of genome organization on the control of gene expression persists as a major challenge in regulatory biology. Most efforts have focused on the role of CTCF-enriched boundary elements and TADs, which enable long-range DNA-DNA associations via loop extrusion processes. However, there is increasing evidence for long-range chromatin loops between promoters and distal enhancers formed through specific DNA sequences, including tethering elements, which bind the GAGA-associated factor (GAF). Previous studies showed that GAF possesses amyloid properties in vitro, bridging separate DNA molecules. In this study, we investigated whether GAF functions as a looping factor in Drosophila development. We employed Micro-C assays to examine the impact of defined GAF mutants on genome topology. These studies suggest that the N-terminal POZ/BTB oligomerization domain is important for long-range associations of distant GAGA-rich tethering elements, particularly those responsible for promoter-promoter interactions that coordinate the activities of distant paralogous genes.

ANT1 overexpression models: Some similarities with facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal dominant disorder characterized by progressive muscle weakness. Adenine nucleotide translocator 1 (ANT1), the only 4q35 gene involved in mitochondrial function, is strongly expressed in FSHD skeletal muscle biopsies. However, its role in FSHD is unclear. In this study, we evaluated ANT1 overexpression effects in primary myoblasts from healthy controls and during Xenopus laevis organogenesis. We also compared ANT1 overexpression effects with the phenotype of FSHD muscle cells and biopsies. Here, we report that the ANT1 overexpression-induced phenotype presents some similarities with FSHD muscle cells and biopsies. ANT1-overexpressing muscle cells showed disorganized morphology, altered cytoskeletal arrangement, enhanced mitochondrial respiration/glycolysis, ROS production, oxidative stress, mitochondrial fragmentation and ultrastructure alteration, as observed in FSHD muscle cells. ANT1 overexpression in Xenopus laevis embryos affected skeletal muscle development, impaired skeletal muscle, altered mitochondrial ultrastructure and led to oxidative stress as observed in FSHD muscle biopsies. Moreover, ANT1 overexpression in X. laevis embryos affected heart structure and mitochondrial ultrastructure leading to cardiac arrhythmia, as described in some patients with FSHD. Overall our data suggest that ANT1 could contribute to mitochondria dysfunction and oxidative stress in FSHD muscle cells by modifying their bioenergetic profile associated with ROS production. Such interplay between energy metabolism and ROS production in FSHD will be of significant interest for future prospects.

DNA Bending Force Facilitates Z-DNA Formation under Physiological Salt Conditions

Z-DNA, a noncanonical helical structure of double-stranded DNA (dsDNA), plays pivotal roles in various biological processes, including transcription regulation. Mechanical stresses on dsDNA, such as twisting and stretching, help to form Z-DNA. However, the effect of DNA bending, one of the most common dsDNA deformations, on Z-DNA formation is utterly unknown. Here, we show that DNA bending induces the formation of Z-DNA, that is, more Z-DNA is formed as the bending force becomes stronger. We regulated the bending force on dsDNA by using D-shaped DNA nanostructures. The B-Z transition was observed by single-molecule fluorescence resonance energy transfer. We found that as the bending force became stronger, Z-DNA was formed at lower Mg²⁺ concentrations. When dsDNA contained cytosine methylations, the B-Z transition occurred at 78 mM Mg²⁺ (midpoint) in the absence of the bending force. However, the B-Z transition occurred at a 28-fold lower Mg²⁺ concentration (2.8 mM) in the presence of the bending force. Monte Carlo simulation suggested that the B-Z transition stabilizes the bent form via the formation of the B-Z junction with base extrusion, which effectively releases the bending stress on DNA. Our results clearly show that the bending force facilitates the B-Z transition under physiological salt conditions.

The regulatory pattern of target gene expression by aberrant enhancer methylation in glioblastoma

Background

Glioblastoma multiforme (GBM) is the most common and aggressive primary malignant brain tumor with grim prognosis. Aberrant DNA methylation is an epigenetic mechanism that promotes GBM carcinogenesis, while the function of DNA methylation at enhancer regions in GBM remains poorly described.

Results

We integrated multi-omics data to identify differential methylation enhancer region (DMER)-genes and revealed global enhancer hypomethylation in GBM. In addition, a DMER-mediated target genes regulatory network and functional enrichment analysis of target genes that might be regulated by hypomethylation enhancer regions showed that aberrant enhancer regions could contribute to tumorigenesis and progression in GBM. Further, we identified 22 modules in which lncRNAs and mRNAs synergistically competed with each other. Finally, through the construction of drug-target association networks, our study identified potential small-molecule drugs for GBM treatment.

Conclusions

Our study provides novel insights for understanding the regulation of aberrant enhancer region methylation and developing methylation-based biomarkers for the diagnosis and treatment of GBM.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04345-8.

GAF is essential for zygotic genome activation and chromatin accessibility in the early Drosophila embryo

Following fertilization, the genomes of the germ cells are reprogrammed to form the totipotent embryo. Pioneer transcription factors are essential for remodeling the chromatin and driving the initial wave of zygotic gene expression. In Drosophila melanogaster, the pioneer factor Zelda is essential for development through this dramatic period of reprogramming, known as the maternal-to-zygotic transition (MZT). However, it was unknown whether additional pioneer factors were required for this transition. We identified an additional maternally encoded factor required for development through the MZT, GAGA Factor (GAF). GAF is necessary to activate widespread zygotic transcription and to remodel the chromatin accessibility landscape. We demonstrated that Zelda preferentially controls expression of the earliest transcribed genes, while genes expressed during widespread activation are predominantly dependent on GAF. Thus, progression through the MZT requires coordination of multiple pioneer-like factors, and we propose that as development proceeds control is gradually transferred from Zelda to GAF.

Identifying Improved Sites for Heterologous Gene Integration Using ATAC-seq

Constructing efficient cellular factories often requires integration of heterologous pathways for synthesis of novel compounds and improved cellular productivity. Few genomic sites are routinely used, however, for efficient integration and expression of heterologous genes, especially in nonmodel hosts. Here, a data-guided framework for informing suitable integration sites for heterologous genes based on ATAC-seq was developed in the nonmodel yeast Komagataella phaffii. Single-copy GFP constructs were integrated using CRISPR/Cas9 into 38 intergenic regions (IGRs) to evaluate the effects of IGR size, intensity of ATAC-seq peaks, and orientation and expression of adjacent genes. Only the intensity of accessibility peaks was observed to have a significant effect, with higher expression observed from IGRs with low- to moderate-intensity peaks than from high-intensity peaks. This effect diminished for tandem, multicopy integrations, suggesting that the additional copies of exogenous sequence buffered the transcriptional unit of the transgene against effects from endogenous sequence context. The approach developed from these results should provide a basis for nominating suitable IGRs in other eukaryotic hosts from an annotated genome and ATAC-seq data.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.

Polymer models for the mechanisms of chromatin 3D folding: review and perspective

In this perspective paper, classical physical models for mammalian interphase chromatin folding are reviewed.

Greater genetic and regulatory plasticity of retained duplicates in Epichloë endophytic fungi

Gene duplicates can act as a source of genetic material from which new functions arise. Most duplicated genes revert to single copy genes and only a small proportion are retained. However, it remains unclear why some duplicate genes persist in the genome for an extended time. We investigate this question by analysing retained gene duplicates in the fungal genus Epichloë, ascomycete fungi that form close endophytic symbioses with their host grasses. Retained duplicates within this genus have two independent origins, but both long pre-date the origin and diversification of the genus Epichloë. We find that loss of retained duplicates within the genus is frequent and often associated with speciation. Retained duplicates have faster evolutionary rates (Ka) and show relaxed selection (Ka/Ks) compared to single copy genes. Both features are time-dependent. Through comparison of conspecific strains, we find greater evolutionary rates in coding regions and sequence divergence in regulatory regions of retained duplicates than single copy genes, with this pattern more pronounced for strains adapted to different grass host species. Consistent with this sequence divergence in regulatory regions, transcriptome analyses show greater expression variation of retained duplicates than single copy genes. This suggest that cis-regulatory changes make important contributions to the expression patterns of retained duplicates. Coupled with supporting observations from the model yeast Saccharomyces cerevisiae, these data suggest that genetic robustness and regulatory plasticity are common drivers behind the retention of duplicated genes in fungi.

Direct Single-Molecule Observation of Sequential DNA Bending Transitions by the Sox2 HMG Box

Sox2 is a pioneer transcription factor that initiates cell fate reprogramming through locus-specific differential regulation. Mechanistically, it was assumed that Sox2 achieves its regulatory diversity via heterodimerization with partner transcription factors. Here, utilizing single-molecule fluorescence spectroscopy, we show that Sox2 alone can modulate DNA structural landscape in a dosage-dependent manner. We propose that such stoichiometric tuning of regulatory DNAs is crucial to the diverse biological functions of Sox2, and represents a generic mechanism of conferring functional plasticity and multiplicity to transcription factors.

Correlation between HDAC4 enhancer DNA methylation and mRNA expression during palatal fusion induced by all-trans retinoic acid

Epithelial-mesenchymal transformation of the medial edge epithelium is the most crucial process in embryonic palatal fusion. This study aimed to explore the relationship and potential mechanism between enhancer DNA methylation and mRNA expression of histone deacetylase 4 (HDAC4) during palatal fusion induced by maternal exposure to all-trans retinoic acid (ATRA). Pregnant mice were administered ATRA (70 mg/kg) by gavage at embryonic gestation day 10.5 (E10.5) to establish a cleft palate (CP) model in C57BL/6J mice. Control groups were given an equivalent volume of corn oil. Pregnant mice were dissected at E14.5 (n = 6) to obtain embryonic palates. HDAC4 enhancer DNA methylation data were obtained from a previous MethylRAD-seq. Methylation-specific polymerase chain reaction (MSP) and real-time quantitative PCR were used to quantify enhancer methylation and the mRNA expression level of HDAC4. Enhancer DNA methylation at a non-CpG site within the HDAC4 gene was hyper-methylated at E14.5 (P: 0.011, log₂ FC:1.67). The MSP results indicated a similar trend, in agreement with the MethylRAD-seq results. The change in the HDAC4 expression level was negatively correlated with its enhancer DNA methylation level, at the non-CpG site, during palatal fusion induced by ATRA. Enhancer DNA methylation of HDAC4 might play an important regulatory role during palatogenesis, especially in embryonic palatal fusion at E 14.5, and may facilitate the development of novel epigenetic biomarkers in the treatment of CP.

CRISPR-Mediated Reorganization of Chromatin Loop Structure

Recent studies have clearly shown that long-range, three-dimensional chromatin looping interactions play a significant role in the regulation of gene expression, but whether looping is responsible for or a result of alterations in gene expression is still unknown. Until recently, how chromatin looping affects the regulation of gene activity and cellular function has been relatively ambiguous, and limitations in existing methods to manipulate these structures prevented in-depth exploration of these interactions. To resolve this uncertainty, we engineered a method for selective and reversible chromatin loop re-organization using CRISPR-dCas9 (CLOuD9). The dynamism of the CLOuD9 system has been demonstrated by successful localization of CLOuD9 constructs to target genomic loci to modulate local chromatin conformation. Importantly, the ability to reverse the induced contact and restore the endogenous chromatin conformation has also been confirmed. Modulation of gene expression with this method establishes the capacity to regulate cellular gene expression and underscores the great potential for applications of this technology in creating stable de novo chromatin loops that markedly affect gene expression in the contexts of cancer and development.

Genome-Wide DNA Methylation Analysis During Palatal Fusion Reveals the Potential Mechanism of Enhancer Methylation Regulating Epithelial Mesenchyme Transformation

Epithelial mesenchyme transformation (EMT) of the medial edge epithelium (MEE) is the crucial process during palatal fusion. This work is aimed to elucidate the enhancer regulatory mechanism by genome-wide DNA methylation analysis of EMT during palatal fusion. Over 800 million clean reads, 325 million enzyme reads, and 234 million mapping reads were generated. The mapping rate was 68.85-74.32%, which included differentially methylated 17299 CCGG sites and 2363 CCWGG sites (p < 0.05, log₂FC >1). Methylated sites in intron and intergenic regions were more compared to other regions of all DNA elements. GO and KEGG analysis indicated that differential methylation sites related to embryonic palatal fusion genes (HDAC4, TCF7L2, and PDGFRB) at the enhancer were located on CCWGG region of non-CpG islands. In addition, the results showed that the enhancer for HDAC4 was hypermethylated, whereas the enhancers for TCF7L2 and PDGFRB were hypomethylated. The methylation status of enhancer regions of HDAC4, PDGFRB, and TCF7L2, involved in the regulation of the EMT during palatal fusion, may enlighten the development of novel epigenetic biomarkers in the treatment of cleft palate.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.

Manipulation of nuclear architecture through CRISPR-mediated chromosomal looping

Chromatin looping is key to gene regulation, yet no broadly applicable methods to selectively modify chromatin loops have been described. We have engineered a method for chromatin loop reorganization using CRISPR-dCas9 (CLOuD9) to selectively and reversibly establish chromatin loops. We demonstrate the power of this technology to selectively modulate gene expression at targeted loci.

Genetic-epigenetic interactions in cis: a major focus in the post-GWAS era

Studies on genetic-epigenetic interactions, including the mapping of methylation quantitative trait loci (mQTLs) and haplotype-dependent allele-specific DNA methylation (hap-ASM), have become a major focus in the post-genome-wide-association-study (GWAS) era. Such maps can nominate regulatory sequence variants that underlie GWAS signals for common diseases, ranging from neuropsychiatric disorders to cancers. Conversely, mQTLs need to be filtered out when searching for non-genetic effects in epigenome-wide association studies (EWAS). Sequence variants in CCCTC-binding factor (CTCF) and transcription factor binding sites have been mechanistically linked to mQTLs and hap-ASM. Identifying these sites can point to disease-associated transcriptional pathways, with implications for targeted treatment and prevention.

Transcriptional regulatory logic of the diurnal cycle in the mouse liver

Many organisms exhibit temporal rhythms in gene expression that propel diurnal cycles in physiology. In the liver of mammals, these rhythms are controlled by transcription-translation feedback loops of the core circadian clock and by feeding-fasting cycles. To better understand the regulatory interplay between the circadian clock and feeding rhythms, we mapped DNase I hypersensitive sites (DHSs) in the mouse liver during a diurnal cycle. The intensity of DNase I cleavages cycled at a substantial fraction of all DHSs, suggesting that DHSs harbor regulatory elements that control rhythmic transcription. Using chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq), we found that hypersensitivity cycled in phase with RNA polymerase II (Pol II) loading and H3K27ac histone marks. We then combined the DHSs with temporal Pol II profiles in wild-type (WT) and Bmal1-/- livers to computationally identify transcription factors through which the core clock and feeding-fasting cycles control diurnal rhythms in transcription. While a similar number of mRNAs accumulated rhythmically in Bmal1-/- compared to WT livers, the amplitudes in Bmal1-/- were generally lower. The residual rhythms in Bmal1-/- reflected transcriptional regulators mediating feeding-fasting responses as well as responses to rhythmic systemic signals. Finally, the analysis of DNase I cuts at nucleotide resolution showed dynamically changing footprints consistent with dynamic binding of CLOCK:BMAL1 complexes. Structural modeling suggested that these footprints are driven by a transient heterotetramer binding configuration at peak activity. Together, our temporal DNase I mappings allowed us to decipher the global regulation of diurnal transcription rhythms in the mouse liver.

The origins of developmental gene regulation

The leap from simple unicellularity to complex multicellularity remains one of life's major enigmas. The origins of metazoan developmental gene regulatory mechanisms are sought by analyzing gene regulation in extant eumetazoans, sponges, and unicellular organisms. The main hypothesis of this manuscript is that, developmental enhancers evolved from unicellular inducible promoters that diversified the expression of regulatory genes during metazoan evolution. Promoters and enhancers are functionally similar; both can regulate the transcription of distal promoters and both direct local transcription. Additionally, enhancers have experimentally characterized structural features that reveal their origin from inducible promoters. The distal co-operative regulation among promoters identified in unicellular opisthokonts possibly represents the precursor of distal regulation of promoters by enhancers. During metazoan evolution, constitutive-type promoters of regulatory genes would have acquired novel receptivity to distal regulatory inputs from promoters of inducible genes that eventually specialized as enhancers. The novel regulatory interactions would have caused constitutively expressed genes controlling differential gene expression in unicellular organisms to become themselves differentially expressed. The consequence of the novel regulatory interactions was that regulatory pathways of unicellular organisms became interlaced and ultimately evolved into the intricate developmental gene regulatory networks (GRNs) of extant metazoans.

iRegNet3D: three-dimensional integrated regulatory network for the genomic analysis of coding and non-coding disease mutations

The mechanistic details of most disease-causing mutations remain poorly explored within the context of regulatory networks. We present a high-resolution three-dimensional integrated regulatory network (iRegNet3D) in the form of a web tool, where we resolve the interfaces of all known transcription factor (TF)-TF, TF-DNA and chromatin-chromatin interactions for the analysis of both coding and non-coding disease-associated mutations to obtain mechanistic insights into their functional impact. Using iRegNet3D, we find that disease-associated mutations may perturb the regulatory network through diverse mechanisms including chromatin looping. iRegNet3D promises to be an indispensable tool in large-scale sequencing and disease association studies.

Plant Enhancers: A Call for Discovery

Higher eukaryotes typically contain many different cell types, displaying different cellular functions that are influenced by biotic and abiotic cues. The different functions are characterized by specific gene expression patterns mediated by regulatory sequences such as transcriptional enhancers. Recent genome-wide approaches have identified thousands of enhancers in animals, reviving interest in enhancers in gene regulation. Although the regulatory roles of plant enhancers are as crucial as those in animals, genome-wide approaches have only very recently been applied to plants. Here we review characteristics of enhancers at the DNA and chromatin level in plants and other species, their similarities and differences, and techniques widely used for genome-wide discovery of enhancers in animal systems that can be implemented in plants.

How far from the SNP may the causative genes be?

While GWAS identify many disease-associated SNPs, using them to decipher disease mechanisms is hindered by the difficulty in mapping SNPs to genes. Most SNPs are in non-coding regions and it is often hard to identify the genes they implicate. To explore how far the SNP may be from the affected genes we used a pathway-based approach. We found that affected genes are often up to 2 Mbps away from the associated SNP, and are not necessarily the closest genes to the SNP. Existing approaches for mapping SNPs to genes leave many SNPs unmapped to genes and reveal only 86 significant phenotype-pathway associations for all known GWAS hits combined. Using the pathway-based approach we propose here allows mapping of virtually all SNPs to genes and reveals 435 statistically significant phenotype-pathway associations. In search for mechanisms that may explain the relationships between SNPs and distant genes, we found that SNPs that are mapped to distant genes have significantly more large insertions/deletions around them than other SNPs, suggesting that these SNPs may sometimes be markers for large insertions/deletions that may affect large genomic regions.

Precision Modulation of Neurodegenerative Disease-Related Gene Expression in Human iPSC-Derived Neurons

The ability to reprogram adult somatic cells into induced pluripotent stem cells (iPSCs) and the subsequent development of protocols for their differentiation into disease-relevant cell types have enabled in-depth molecular analyses of multiple disease states as hitherto impossible. Neurons differentiated from patient-specific iPSCs provide a means to recapitulate molecular phenotypes of neurodegenerative diseases in vitro. However, it remains challenging to conduct precise manipulations of gene expression in iPSC-derived neurons towards modeling complex human neurological diseases. The application of CRISPR/Cas9 to mammalian systems is revolutionizing the utilization of genome editing technologies in the study of molecular contributors to the pathogenesis of numerous diseases. Here, we demonstrate that CRISPRa and CRISPRi can be used to exert precise modulations of endogenous gene expression in fate-committed iPSC-derived neurons. This highlights CRISPRa/i as a major technical advancement in accessible tools for evaluating the specific contributions of critical neurodegenerative disease-related genes to neuropathogenesis.

Allosteric communication between DNA-binding and light-responsive domains of diatom class I aureochromes

The modular architecture of aureochrome blue light receptors, found in several algal groups including diatoms, is unique by having the LOV-type photoreceptor domain fused to the C-terminus of its putative effector, an N-terminal DNA-binding bZIP module. The structural and functional understanding of aureochromes’ light-dependent signaling mechanism is limited, despite their promise as an optogenetic tool. We show that class I aureochromes 1a and 1c from the diatom Phaeodactylum tricornutum are regulated in a light-independent circadian rhythm. These aureochromes are capable to form functional homo- and heterodimers, which recognize the ACGT core sequence within the canonical ‘aureo box’, TGACGT, in a light-independent manner. The bZIP domain holds a more folded and less flexible but extended conformation in the duplex DNA-bound state. FT-IR spectroscopy in the absence and the presence of DNA shows light-dependent helix unfolding in the LOV domain, which leads to conformational changes in the bZIP region. The solution structure of DNA bound to aureochrome points to a tilted orientation that was further validated by molecular dynamics simulations. We propose that aureochrome signaling relies on an allosteric pathway from LOV to bZIP that results in conformational changes near the bZIP-DNA interface without major effects on the binding affinity.

Architectural proteins Pita, Zw5,and ZIPIC contain homodimerization domain and support specific long-range interactions in Drosophila

According to recent models, as yet poorly studied architectural proteins appear to be required for local regulation of enhancer-promoter interactions, as well as for global chromosome organization. Transcription factors ZIPIC, Pita and Zw5 belong to the class of chromatin insulator proteins and preferentially bind to promoters near the TSS and extensively colocalize with cohesin and condensin complexes. ZIPIC, Pita and Zw5 are structurally similar in containing the N-terminal zinc finger-associated domain (ZAD) and different numbers of C2H2-type zinc fingers at the C-terminus. Here we have shown that the ZAD domains of ZIPIC, Pita and Zw5 form homodimers. In Drosophila transgenic lines, these proteins are able to support long-distance interaction between GAL4 activator and the reporter gene promoter. However, no functional interaction between binding sites for different proteins has been revealed, suggesting that such interactions are highly specific. ZIPIC facilitates long-distance stimulation of the reporter gene by GAL4 activator in yeast model system. Many of the genomic binding sites of ZIPIC, Pita and Zw5 are located at the boundaries of topologically associated domains (TADs). Thus, ZAD-containing zinc-finger proteins can be attributed to the class of architectural proteins.

Put your 3D glasses on: plant chromatin is on show

The three-dimensional organization of the eukaryotic nucleus and its chromosomal conformation have emerged as important features in the complex network of mechanisms behind gene activity and genome connectivity dynamics, which can be evidenced in the regionalized chromosomal spatial distribution and the clustering of diverse genomic regions with similar expression patterns. The development of chromatin conformation capture (3C) techniques has permitted the elucidation of commonalities between the eukaryotic phyla, as well as important differences among them. The growing number of studies in the field performed in plants has shed light on the structural and regulatory features of these organisms. For instance, it has been proposed that plant chromatin can be arranged into different conformations such as Rabl, Rosette-like, and Bouquet, and that both short- and long-range chromatin interactions occur in Arabidopsis. In this review, we compile the current knowledge about chromosome architecture characteristics in plants, as well as the molecular events and elements (including long non-coding RNAs, histone and DNA modifications, chromatin remodeling complexes, and transcription factors) shaping the genome three-dimensional conformation. Furthermore, we discuss the developmental outputs of genome topology-mediated gene expression regulation. It is becoming increasingly clear that new tools and techniques with higher resolution need to be developed and implemented in Arabidopsis and other model plants in order to better understand chromosome architecture dynamics, from an integrative perspective with other fields of plant biology such as development, stress biology, and finally agriculture.

Spatially coordinated replication and minimization of expression noise constrain three-dimensional organization of yeast genome

Despite recent advances, the underlying functional constraints that shape the three-dimensional organization of eukaryotic genome are not entirely clear. Through comprehensive multivariate analyses of genome-wide datasets, we show that cis and trans interactions in yeast genome have significantly distinct functional associations. In particular, (i) the trans interactions are constrained by coordinated replication and co-varying mutation rates of early replicating domains through interactions among early origins, while cis interactions are constrained by coordination of late replication through interactions among late origins; (ii)cis and trans interactions exhibit differential preference for nucleosome occupancy; (iii)cis interactions are also constrained by the essentiality and co-fitness of interacting genes. Essential gene clusters associate with high average interaction frequency, relatively short-range interactions of low variance, and exhibit less fluctuations in chromatin conformation, marking a physically restrained state of engaged loci that, we suggest, is important to mitigate the epigenetic errors by restricting the spatial mobility of loci. Indeed, the genes with lower expression noise associate with relatively short-range interactions of lower variance and exhibit relatively higher average interaction frequency, a property that is conserved across Escherichia coli,yeast, and mESCs. Altogether, our observations highlight the coordination of replication and the minimization of expression noise, not necessarily co-expression of genes, as potent evolutionary constraints shaping the spatial organization of yeast genome.

Perturbing Chromatin Structure to Understand Mechanisms of Gene Expression

The study of nuclear structure falls between the fields of cell biology and molecular biology and draws on techniques from both fields. In recent years, many exciting advances have been made in these areas, including single-molecule and superresolution imaging and the development of chromosome conformation capture (3C)-based technologies, which have brought the advent of genome-wide analysis of chromatin structure and contacts. However, many questions remain as to the function of nuclear structures, in particular their influence on transcription. Here we describe studies that have directly manipulated nuclear architecture at various levels and thus have clarified the causal influence of structure on transcription. We will also highlight open questions in the field, most notably regarding our understanding of the dynamics and variability in nuclear structure and its influence on gene expression.

ZNF555 protein binds to transcriptional activator site of 4qA allele and ANT1: potential implication in Facioscapulohumeral dystrophy

Facioscapulohumeral dystrophy (FSHD) is an epi/genetic satellite disease associated with at least two satellite sequences in 4q35: (i) D4Z4 macrosatellite and (ii) β-satellite repeats (BSR), a prevalent part of the 4qA allele. Most of the recent FSHD studies have been focused on a DUX4 transcript inside D4Z4 and its tandem contraction in FSHD patients. However, the D4Z4-contraction alone is not pathological, which would also require the 4qA allele. Since little is known about BSR, we investigated the 4qA BSR functional role in the transcriptional control of the FSHD region 4q35. We have shown that an individual BSR possesses enhancer activity leading to activation of the Adenine Nucleotide Translocator 1 gene (ANT1), a major FSHD candidate gene. We have identified ZNF555, a previously uncharacterized protein, as a putative transcriptional factor highly expressed in human primary myoblasts that interacts with the BSR enhancer site and impacts the ANT1 promoter activity in FSHD myoblasts. The discovery of the functional role of the 4qA allele and ZNF555 in the transcriptional control of ANT1 advances our understanding of FSHD pathogenesis and provides potential therapeutic targets.

Quantitation of interactions between two DNA loops demonstrates loop domain insulation in E. coli cells

Eukaryotic gene regulation involves complex patterns of long-range DNA-looping interactions between enhancers and promoters, but how these specific interactions are achieved is poorly understood. Models that posit other DNA loops--that aid or inhibit enhancer-promoter contact--are difficult to test or quantitate rigorously in eukaryotic cells. Here, we use the well-characterized DNA-looping proteins Lac repressor and phage λ CI to measure interactions between pairs of long DNA loops in E. coli cells in the three possible topological arrangements. We find that side-by-side loops do not affect each other. Nested loops assist each other's formation consistent with their distance-shortening effect. In contrast, alternating loops, where one looping element is placed within the other DNA loop, inhibit each other's formation, thus providing clear support for the loop domain model for insulation. Modeling shows that combining loop assistance and loop interference can provide strong specificity in long-range interactions.

Mechanisms and proteins involved in long-distance interactions

Due to advances in genome-wide technologies, consistent distant interactions within chromosomes of higher eukaryotes have been revealed. In particular, it has been shown that enhancers can specifically and directly interact with promoters by looping out intervening sequences, which can be up to several hundred kilobases long. This review is focused on transcription factors that are supposed to be involved in long-range interactions. Available data are in agreement with the model that several known transcription factors and insulator proteins belong to an abundant but poorly studied class of proteins that are responsible for chromosomal architecture.

Manipulating nuclear architecture

The eukaryotic genome is highly organized in the nucleus. Genes can be localized to specific nuclear compartments in a manner reflecting their activity. A plethora of recent reports has described multiple levels of chromosomal folding that can be related to gene-specific expression states. Here we discuss studies designed to probe the causal impact of genome organization on gene expression. The picture that emerges is that of a reciprocal relationship in which nuclear organization is not only shaped by gene expression states but also directly influences them.

Chromatin insulators and long-distance interactions in Drosophila

Data on long-distance enhancer-mediated activation of gene promoters and complex regulation of gene expression by multiple enhancers have prompted the hypothesis that the action of enhancers is restricted by insulators. Studies with transgenic lines have shown that insulators are responsible for establishing proper local interactions between regulatory elements, but not for defining independent transcriptional domains that restrict the activity of enhancers. It has also become apparent that enhancer blocking is only one of several functional activities of known insulator proteins, which also contribute to the organization of chromosome architecture and the integrity of regulatory elements.

Developmental regulation of chromatin conformation by Hox proteins in Drosophila

We present a strategy to examine the chromatin conformation of individual loci in specific cell types during Drosophila embryogenesis. Regulatory DNA is tagged with binding sites (lacO) for LacI, which is used to immunoprecipitate the tagged chromatin from specific cell types. We applied this approach to Distalless (Dll), a gene required for limb development in Drosophila. We show that the local chromatin conformation at Dll depends on the cell type: in cells that express Dll, the 5' regulatory region is in close proximity to the Dll promoter. In Dll-nonexpressing cells this DNA is in a more extended configuration. In addition, transcriptional activators and repressors are bound to Dll regulatory DNA in a cell type-specific manner. The pattern of binding by GAGA factor and the variant histone H2Av suggest that they play a role in the regulation of Dll chromatin conformation in expressing and nonexpressing cell types, respectively.

Controlling long-range genomic interactions at a native locus by targeted tethering of a looping factor

Chromatin loops juxtapose distal enhancers with active promoters, but their molecular architecture and relationship with transcription remain unclear. In erythroid cells, the locus control region (LCR) and β-globin promoter form a chromatin loop that requires transcription factor GATA1 and the associated molecule Ldb1. We employed artificial zinc fingers (ZF) to tether Ldb1 to the β-globin promoter in GATA1 null erythroblasts, in which the β-globin locus is relaxed and inactive. Remarkably, targeting Ldb1 or only its self-association domain to the β-globin promoter substantially activated β-globin transcription in the absence of GATA1. Promoter-tethered Ldb1 interacted with endogenous Ldb1 complexes at the LCR to form a chromatin loop, causing recruitment and phosphorylation of RNA polymerase II. ZF-Ldb1 proteins were inactive at alleles lacking the LCR, demonstrating that their activities depend on long-range interactions. Our findings establish Ldb1 as a critical effector of GATA1-mediated loop formation and indicate that chromatin looping causally underlies gene regulation.

Small RNA and transcriptional upregulation

Small RNA molecules, such as microRNA and small interfering RNA, have emerged as master regulators of gene expression through their ability to suppress target genes in a phenomenon collectively called RNA interference (RNAi). There is growing evidence that small RNAs can also serve as activators of gene expression by targeting gene regulatory sequences. This novel mechanism, known as RNA activation (RNAa), appears to be conserved in at least mammalian cells and triggered by both endogenous and artificially designed small RNAs. RNAa depends on Argonaute proteins, but possesses kinetics distinct from that of RNAi. Epigenetic changes are associated with RNAa and may contribute to transcriptional activation of target genes, but the underlying mechanism remains elusive. Given the potential of RNAa as a molecular tool for studying gene function and as a therapeutic for disease, further research is needed to completely elucidate its molecular mechanism in order to refine the rules for target selection and improve strategies for exploiting it therapeutically. WIREs RNA 2011 2 748-760 DOI: 10.1002/wrna.90 For further resources related to this article, please visit the WIREs website.

Drosophila BTB/POZ domains of "ttk group" can form multimers and selectively interact with each other

The BTB (bric-a-brac, tramtrack and broad complex)/POZ (poxvirus and zinc finger) domain is a conserved protein-protein interaction motif contained in a variety of transcription factors involved in development, chromatin remodeling, insulator activity, and carcinogenesis. All well-studied mammalian BTB domains form obligate homodimers and, rarely, tetramers. Only the BTB domain of the Drosophila GAGA factor (GAF) has been shown to exist as higher-order multimers. The BTB domain of GAF belongs to the "ttk group" that contains several highly conserved sequences not found in other BTB domains. Here, we have shown by size-exclusion chromatography, chemical cross-linking, and nondenaturing PAGE that four additional BTB domains of the ttk group-Batman, Mod(mdg4), Pipsqueak, and Tramtrack-can form multimers, like GAF. Interestingly, the BTB domains of GAF and Batman have formed a wide range of complexes and interacted in the yeast two-hybrid assay with other BTB domains tested. In contrast, the BTB domains of Mod(mdg4), Pipsqueak, and Tramtrack have formed stable high-order multimer complexes and failed to interact with each other. The BTB domain of Drosophila CP190 protein does not belong to the ttk group. This BTB domain has formed stable dimers and has not interacted with domains of the ttk group. Previously, it was suggested that GAF oligomerization into higher-order complexes facilitates long-range activation by providing a protein bridge between an enhancer and a promoter. Unexpectedly, experiments in the Drosophila model system have not supported the role of GAF in organization of long-distance interaction between the yeast GAL4 activator and the white promoter.

Functional and mechanistic diversity of distal transcription enhancers

Biological differences among metazoans and between cell types in a given organism arise in large part due to differences in gene expression patterns. Gene-distal enhancers are key contributors to these expression patterns, exhibiting both sequence diversity and cell type specificity. Studies of long-range interactions indicate that enhancers are often important determinants of nuclear organization, contributing to a general model for enhancer function that involves direct enhancer-promoter contact. However, mechanisms for enhancer function are emerging that do not fit solely within such a model, suggesting that enhancers as a class of DNA regulatory element may be functionally and mechanistically diverse.

Alanine zipper-like coiled-coil domains are necessary for homotypic dimerization of plant GAGA-factors in the nucleus and nucleolus

GAGA-motif binding proteins control transcriptional activation or repression of homeotic genes. Interestingly, there are no sequence similarities between animal and plant proteins. Plant BBR/BPC-proteins can be classified into two distinct groups: Previous studies have elaborated on group I members only and so little is known about group II proteins. Here, we focused on the initial characterization of AtBPC6, a group II protein from Arabidopsis thaliana. Comparison of orthologous BBR/BPC sequences disclosed two conserved signatures besides the DNA binding domain. A first peptide signature is essential and sufficient to target AtBPC6-GFP to the nucleus and nucleolus. A second domain is predicted to form a zipper-like coiled-coil structure. This novel type of domain is similar to Leucine zippers, but contains invariant alanine residues with a heptad spacing of 7 amino acids. By yeast-2-hybrid and BiFC-assays we could show that this Alanine zipper domain is essential for homotypic dimerization of group II proteins in vivo. Interhelical salt bridges and charge-stabilized hydrogen bonds between acidic and basic residues of the two monomers are predicted to form an interaction domain, which does not follow the classical knobs-into-holes zipper model. FRET-FLIM analysis of GFP/RFP-hybrid fusion proteins validates the formation of parallel dimers in planta. Sequence comparison uncovered that this type of domain is not restricted to BBR/BPC proteins, but is found in all kingdoms.

Construction of cis-regulatory input functions of yeast promoters

Promoters contain a large number of binding sites for transcriptional factors transmitting signals from a variety of cellular pathways. The promoter processes these input signals and sets the level of gene expression, the output of the gene. Here, we describe how to design genetic constructs and measure gene expression to deliver data suitable for quantitative analysis. Synthetic genetic constructs are well suited to precisely control and measure gene expression to construct cis-regulatory input functions. These functions can be used to predict gene expression based on signal intensities transmitted to activators and repressors in the gene regulatory region. Simple models of gene expression are presented for competitive and noncompetitive repressions. Complex phenomena, exemplified by synergistic silencing, are modeled by reaction-diffusion equations.

Understanding distal transcriptional regulation from sequence, expression and interactome perspectives

Gene regulation in eukaryotes involves a complex interplay between the proximal promoter and distal genomic elements (such as enhancers) which work in concert to drive precise spatio-temporal gene expression. The experimental localization and characterization of gene regulatory elements is a very complex and resource-intensive process. The computational identification of regulatory regions that confer spatiotemporally specific tissue-restricted expression of a gene is thus an important challenge for computational biology. One of the most popular strategies for enhancer localization from DNA sequence is the use of conservation-based prefiltering and more recently, the use of canonical (transcription factor motifs) or de novo tissue-specific sequence motifs. However, there is an ongoing effort in the computational biology community to further improve the fidelity of enhancer predictions from sequence data by integrating other, complementary genomic modalities. In this work, we propose a framework that complements existing methodologies for prospective enhancer identification. The methods in this work are derived from two key insights: (i) that chromatin modification signatures can discriminate proximal and distally located regulatory regions and (ii) the notion of promoter-enhancer cross-talk (as assayed in 3C/5C experiments) might have implications in the search for regulatory sequences that co-operate with the promoter to yield tissue-restricted, gene-specific expression.

EZH2-dependent chromatin looping controls INK4a and INK4b, but not ARF, during human progenitor cell differentiation and cellular senescence

BACKGROUND

The INK4b-ARF-INK4a tumour suppressor locus controls the balance between progenitor cell renewal and cancer. In this study, we investigated how higher-order chromatin structure modulates differential expression of the human INK4b-ARF-INK4a locus during progenitor cell differentiation, cellular ageing and senescence of cancer cells.

RESULTS

We found that INK4b and INK4a, but not ARF, are upregulated following the differentiation of haematopoietic progenitor cells, in ageing fibroblasts and in senescing malignant rhabdoid tumour cells. To investigate the underlying molecular mechanism we analysed binding of polycomb group (PcG) repressive complexes (PRCs) and the spatial organization of the INK4b-ARF-INK4a locus. In agreement with differential derepression, PcG protein binding across the locus is discontinuous. As we described earlier, PcG repressors bind the INK4a promoter, but not ARF. Here, we identified a second peak of PcG binding that is located approximately 3 kb upstream of the INK4b promoter. During progenitor cell differentiation and ageing, PcG silencer EZH2 attenuates, causing loss of PRC binding and transcriptional activation of INK4b and INK4a. The expression pattern of the locus is reflected by its organization in space. In the repressed state, the PRC-binding regions are in close proximity, while the intervening chromatin harbouring ARF loops out. Down regulation of EZH2 causes release of the approximately 35 kb repressive chromatin loop and induction of both INK4a and INK4b, whereas ARF expression remains unaltered.

CONCLUSION

PcG silencers bind and coordinately regulate INK4b and INK4a, but not ARF, during a variety of physiological processes. Developmentally regulated EZH2 levels are one of the factors that can determine the higher order chromatin structure and expression pattern of the INK4b-ARF-INK4a locus, coupling human progenitor cell differentiation to proliferation control. Our results revealed a chromatin looping mechanism of long-range control and argue against models involving homogeneous spreading of PcG silencers across the INK4b-ARF-INK4a locus.

High local concentration: a fundamental strategy of life

Local increase in concentration is a basic principle of transcriptional control. Closer inspection reveals that it is a major force governing all interactions within and between proteins and DNA. Local increase in concentration acts on all levels of organization of living matter. The structures and functions of two central molecules of life-the linear polymers DNA and protein-are particularly illuminating examples. Local increase in concentration leads to cooperative or competitive interactions between molecules. It is an important principle of life.

Paused Pol II captures enhancer activity and acts as a potent insulator

Enhancers act over many kilobase pairs to activate target promoters, but their activity is constrained by insulator elements that prevent indiscriminate activation of nearby genes. In the July 1, 2009, issue of Genes & Development, Chopra and colleagues (pp. 1505-1509) report that promoters containing a stalled Pol II are activated by enhancers, but these promoters also serve as insulators that block enhancers from reaching more distal genes. This new class of insulators provide critical clues to regulatory mechanisms.