The histone reader PHF7 cooperates with the SWI/SNF complex at cardiac super enhancers to promote direct reprogramming.

Direct cardiac reprogramming of fibroblasts to cardiomyocytes presents an attractive therapeutic strategy to restore cardiac function following injury. Cardiac reprogramming was initially achieved through overexpression of the transcription factors Gata4, Mef2c and Tbx5; later, Hand2 and Akt1 were found to further enhance this process1-5. Yet, staunch epigenetic barriers severely limit the ability of these cocktails to reprogramme adult fibroblasts6,7. We undertook a screen of mammalian gene regulatory factors to discover novel regulators of cardiac reprogramming in adult fibroblasts and identified the histone reader PHF7 as the most potent activating factor8. Mechanistically, PHF7 localizes to cardiac super enhancers in fibroblasts, and through cooperation with the SWI/SNF complex, it increases chromatin accessibility and transcription factor binding at these sites. Furthermore, PHF7 recruits cardiac transcription factors to activate a positive transcriptional autoregulatory circuit in reprogramming. Importantly, PHF7 achieves efficient reprogramming in the absence of Gata4. Here, we highlight the underexplored necessity of cardiac epigenetic readers, such as PHF7, in harnessing chromatin remodelling and transcriptional complexes to overcome critical barriers to direct cardiac reprogramming.

Characterising open chromatin in chick embryos identifies cis-regulatory elements important for paraxial mesoderm formation and axis extension.

Somites arising from paraxial mesoderm are a hallmark of the segmented vertebrate body plan. They form sequentially during axis extension and generate musculoskeletal cell lineages. How paraxial mesoderm becomes regionalised along the axis and how this correlates with dynamic changes of chromatin accessibility and the transcriptome remains unknown. Here, we report a spatiotemporal series of ATAC-seq and RNA-seq along the chick embryonic axis. Footprint analysis shows differential coverage of binding sites for several key transcription factors, including CDX2, LEF1 and members of HOX clusters. Associating accessible chromatin with nearby expressed genes identifies cis-regulatory elements (CRE) for TCF15 and MEOX1. We determine their spatiotemporal activity and evolutionary conservation in Xenopus and human. Epigenome silencing of endogenous CREs disrupts TCF15 and MEOX1 gene expression and recapitulates phenotypic abnormalities of anterior-posterior axis extension. Our integrated approach allows dissection of paraxial mesoderm regulatory circuits in vivo and has implications for investigating gene regulatory networks.

Efficient chromatin accessibility mapping in situ by nucleosome-tethered tagmentation.

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. A popular example is ATAC-seq, whereby Tn5 transposase inserts sequencing adapters into accessible DNA ('tagmentation'). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here, we show that by simply modifying the tagmentation conditions for histone H3K4me2 or H3K4me3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce chromatin accessibility maps that are indistinguishable from the best ATAC-seq maps. Thus, chromatin accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

Cells keep their DNA tidy by wrapping it into structures called nucleosomes. Each of these structures contains a short section of DNA wound around a cluster of proteins called histones. Not only do nucleosomes keep the genetic code organized, they also control whether the proteins that can switch genes on or off have access to the DNA. When genes turn on, the nucleosomes unwrap, exposing sections of genetic code called 'gene regulatory elements'. These elements attract the proteins that help read and copy nearby genes so the cell can make new proteins. Determining which regulatory elements are exposed at any given time can provide useful information about what is happening inside a cell, but the procedure can be expensive. The most popular way to map which regulatory elements are exposed is using a technique called Assay for Transposase-Accessible Chromatin using sequencing, or ATAC-seq for short. The 'transposase' in the acronym is an enzyme that cuts areas of DNA that are not wound around histones and prepares them for detection by DNA sequencing. Unfortunately, the data from ATAC-seq are often noisy (there are random factors that produce a signal that is detected but is not a 'real' result), so more sequencing is required to differentiate between real signal and noise, increasing the expense of ATAC-seq experiments. Furthermore, although ATAC-seq can identify unspooled sections of DNA, it cannot provide a direct connection between active genes and unwrapped DNA. To find the link between unspooled DNA and active genes, Henikoff et al. adapted a technique called CUT&Tag. Like ATAC-seq, it also uses transposases to cut the genome, but it allows more control over where the cuts occur. When genes are switched on, the proteins reading them leave chemical marks on the histones they pass. CUT&Tag attaches a transposase to a molecule that recognizes and binds to those marks. This allowed Henikoff et al. to guide the transposases to unspooled regions of DNA bordering active genes. The maps of gene regulatory elements produced using this method were the same as the best ATAC-seq maps. And, because the transposases could only access gaps near active genes, the data provided evidence that genes switching on leads to regulatory elements in the genome unwrapping. This new technique is simple enough that Henikoff et al. were able to perform it from home on the countertop of a laundry room. By tethering the transposases to histone marks it was possible to detect unspooled DNA that was active more efficiently than with ATAC-seq. This lowers laboratory costs by reducing the cost of DNA sequencing, and may also improve the detection of gaps between nucleosomes in single cells.

Robust inference of positive selection on regulatory sequences in the human brain.

A longstanding hypothesis is that divergence between humans and chimpanzees might have been driven more by regulatory level adaptations than by protein sequence adaptations. This has especially been suggested for regulatory adaptations in the evolution of the human brain. We present a new method to detect positive selection on transcription factor binding sites on the basis of measuring predicted affinity change with a machine learning model of binding. Unlike other methods, this approach requires neither defining a priori neutral sites nor detecting accelerated evolution, thus removing major sources of bias. We scanned the signals of positive selection for CTCF binding sites in 29 human and 11 mouse tissues or cell types. We found that human brain-related cell types have the highest proportion of positive selection. This result is consistent with the view that adaptive evolution to gene regulation has played an important role in evolution of the human brain.

Dynamics of sequestering the limiting p300/CBP, viral cis-regulatory elements, and disease.

Many studies showed that the p300/CBP coactivator is limiting. Here we review three studies that showed how transcription complexes formed on viral cis-regulator elements compete with cellular transcription complexes by sequestrating the p300/CBP coactivator. According to the microcompetition model, this sequestering can cause disease. We use the microcompetition model to explain how a specific type of sequestering, caused by a latent virus that has an active cis-regulatory element in its promoter/enhancer that binds the transcription complex p300/CBP ·GABP can cause diseases such as cancer, atherosclerosis, diabetes, and certain autoimmune diseases.

DNA methylation in the 5-HTT regulatory region is associated with CO2-induced fear in panic disorder patients.

A polymorphism in the gene encoding the serotonin (5-HT) transporter (5-HTT) has been shown to moderate the response to CO2 inhalation, an experimental model for panic attacks (PAs). Recurrent, unpredictable PAs represent, together with anticipatory anxiety of recurring attacks, the core feature of panic disorder (PD) and significantly interfere with patients' daily life. In addition to genetic components, accumulating evidence suggests that epigenetic mechanisms, which regulate gene expression by modifying chromatin structure, also play a fundamental role in the etiology of mental disorders. However, in PD, epigenetic mechanisms have barely been examined to date. In the present study, we investigated the relationship between methylation at the regulatory region of the gene encoding the 5-HTT and the reactivity to a 35% CO2 inhalation in PD patients. We focused on four specific CpG sites and found a significant association between the methylation level of one of these CpG sites and the fear response. This suggests that the emotional response to CO2 inhalation might be moderated by an epigenetic mechanism, and underlines the implication of the 5-HT system in PAs. Future studies are needed to further investigate epigenetic alterations in PD and their functional consequences. These insights can increase our understanding of the underlying pathophysiology and support the development of new treatment strategies.

Dynamics of Global Gene Expression and Regulatory Elements in Growing Brachypodium Root System.

Root systems are dynamic and adaptable organs that play critical roles in plant development. However, how roots grow and accumulate biomass during plant life cycle and in relation to shoot growth phenology remains understudied. A comprehensive time-dependent root morphological analysis integrated with molecular signatures is then required to advance our understanding of root growth and development. Here we studied Brachypodium distachyon rooting process by monitoring root morphology, biomass production, and C/N ratios during developmental stages. To provide insight into gene regulation that accompanies root growth, we generated comprehensive transcript profiles of Brachypodium whole-root system at four developmental stages. Our data analysis revealed that multiple biological processes including trehalose metabolism and various families of transcription factors (TFs) were differentially expressed in root system during plant development. In particular, the AUX/IAA, ERFs, WRKY, NAC, and MADS TF family members were upregulated as plant entered the booting/heading stage, while ARFs and GRFs were downregulated suggesting these TF families as important factors involved in specific phases of rooting, and possibly in regulation of transition to plant reproductive stages. We identified several Brachypodium candidate root biomass-promoting genes and cis-regulatory elements for further functional validations and root growth improvements in grasses.

RUNX2 co-operates with EGR1 to regulate osteogenic differentiation through Htra1 enhancers.

Runt-related transcription factor 2 (Runx2) has been shown to regulate osteoblast differentiation by directly or indirectly regulating numerous osteoblast-related genes. However, our understanding of the transcriptional mechanisms of RUNX2 is mainly restricted to its transactivation, while the mechanism underlying its inhibitory effect during osteoblast differentiation remains largely unknown. Here, we incorporated the anti-RUNX2 chromatin immunoprecipitation (ChIP) sequencing in MC3T3-E1 cells and RNA-sequencing of parietal bone from Runx2 heterozygous mutant mice, to identify the putative genes negatively regulated by RUNX2. We identified HtrA serine peptidase 1 (Htra1) as a target gene and found ten candidate Htra1 enhancers potentially regulated by RUNX2, among which seven were verified by dual-luciferase assays. Furthermore, we investigated the motifs in the vicinity of RUNX2-binding sites and identified early growth response 1 (EGR1) as a potential partner transcription factor (TF) potentially regulating Htra1 expression, which was subsequently confirmed by Re-ChIP assays. RUNX2 and EGR1 co-repressed Htra1 and increased the expression levels of other osteoblast marker genes, such as osterix, osteocalcin, and osteoprotegerin at the messenger RNA and protein level. Moreover, Alizarin red staining combined with alkaline phosphatase (ALP) staining showed decreased calcified nodules and ALP activity in the siRUNX2+siEGR1 group compared with siRUNX2 group. Our findings revealed the detailed mechanism of the inhibitory function of RUNX2 towards its downstream genes, along with its partner TFs, to promote osteoblast differentiation.

Splicing of enhancer-associated lincRNAs contributes to enhancer activity.

Transcription is common at active mammalian enhancers sometimes giving rise to stable enhancer-associated long intergenic noncoding RNAs (elincRNAs). Expression of elincRNA is associated with changes in neighboring gene product abundance and local chromosomal topology, suggesting that transcription at these loci contributes to gene expression regulation in cis Despite the lack of evidence supporting sequence-dependent functions for most elincRNAs, splicing of these transcripts is unexpectedly common. Whether elincRNA splicing is a mere consequence of cognate enhancer activity or if it directly impacts enhancer function remains unresolved. Here, we investigate the association between elincRNA splicing and enhancer activity in mouse embryonic stem cells. We show that multi-exonic elincRNAs are enriched at conserved enhancers, and the efficient processing of elincRNAs is strongly associated with their cognate enhancer activity. This association is supported by their enrichment in enhancer-specific chromatin signatures; elevated binding of co-transcriptional regulators; increased local intra-chromosomal DNA contacts; and strengthened cis-regulation on target gene expression. Our results support the role of efficient RNA processing of enhancer-associated transcripts to cognate enhancer activity.

Hominin-specific regulatory elements selectively emerged in oligodendrocytes and are disrupted in autism patients.

Speciation is associated with substantial rewiring of the regulatory circuitry underlying the expression of genes. Determining which changes are relevant and underlie the emergence of the human brain or its unique susceptibility to neural disease has been challenging. Here we annotate changes to gene regulatory elements (GREs) at cell type resolution in the brains of multiple primate species spanning most of primate evolution. We identify a unique set of regulatory elements that emerged in hominins prior to the separation of humans and chimpanzees. We demonstrate that these hominin gains perferentially affect oligodendrocyte function postnatally and are preferentially affected in the brains of autism patients. This preference is also observed for human-specific GREs suggesting this system is under continued selective pressure. Our data provide a roadmap of regulatory rewiring across primate evolution providing insight into the genomic changes that underlie the emergence of the brain and its susceptibility to neural disease.

Regulation of OmpA Translation and Shigella dysenteriae Virulence by an RNA Thermometer.

RNA thermometers are cis-acting riboregulators that mediate the posttranscriptional regulation of gene expression in response to environmental temperature. Such regulation is conferred by temperature-responsive structural changes within the RNA thermometer that directly result in differential ribosomal binding to the regulated transcript. The significance of RNA thermometers in controlling bacterial physiology and pathogenesis is becoming increasingly clear. This study combines in silico, molecular genetics, and biochemical analyses to characterize both the structure and function of a newly identified RNA thermometer within the ompA transcript of Shigella dysenteriae First identified by in silico structural predictions, genetic analyses have demonstrated that the ompA RNA thermometer is a functional riboregulator sufficient to confer posttranscriptional temperature-dependent regulation, with optimal expression observed at the host-associated temperature of 37°C. Structural studies and ribosomal binding analyses have revealed both increased exposure of the ribosomal binding site and increased ribosomal binding to the ompA transcript at permissive temperatures. The introduction of site-specific mutations predicted to alter the temperature responsiveness of the ompA RNA thermometer has predictable consequences for both the structure and function of the regulatory element. Finally, in vitro tissue culture-based analyses implicate the ompA RNA thermometer as a bona fide S. dysenteriae virulence factor in this bacterial pathogen. Given that ompA is highly conserved among Gram-negative pathogens, these studies not only provide insight into the significance of riboregulation in controlling Shigella virulence, but they also have the potential to facilitate further understanding of the physiology and/or pathogenesis of a wide range of bacterial species.

Deep Conservation of cis-Element Variants Regulating Plant Hormonal Responses.

Phytohormones regulate many aspects of plant life by activating transcription factors (TFs) that bind sequence-specific response elements (REs) in regulatory regions of target genes. Despite their short length, REs are degenerate, with a core of just 3 to 4 bp. This degeneracy is paradoxical, as it reduces specificity and REs are extremely common in the genome. To study whether RE degeneracy might serve a biological function, we developed an algorithm for the detection of regulatory sequence conservation and applied it to phytohormone REs in 45 angiosperms. Surprisingly, we found that specific RE variants are highly conserved in core hormone response genes. Experimental evidence showed that specific variants act to regulate the magnitude and spatial profile of hormonal response in Arabidopsis (Arabidopsis thaliana) and tomato (Solanum lycopersicum). Our results suggest that hormone-regulated TFs bind a spectrum of REs, each coding for a distinct transcriptional response profile. Our approach has implications for precise genome editing and for rational promoter design.

The lineage-specific, intrinsically disordered N-terminal extension of monothiol glutaredoxin 1 from trypanosomes contains a regulatory region.

Glutaredoxins (Grx) are small proteins conserved throughout all the kingdoms of life that are engaged in a wide variety of biological processes and share a common thioredoxin-fold. Among them, class II Grx are redox-inactive proteins involved in iron-sulfur (FeS) metabolism. They contain a single thiol group in their active site and use low molecular mass thiols such as glutathione as ligand for binding FeS-clusters. In this study, we investigated molecular aspects of 1CGrx1 from the pathogenic parasite Trypanosoma brucei brucei, a mitochondrial class II Grx that fulfills an indispensable role in vivo. Mitochondrial 1CGrx1 from trypanosomes differs from orthologues in several features including the presence of a parasite-specific N-terminal extension (NTE) whose role has yet to be elucidated. Previously we have solved the structure of a truncated form of 1CGrx1 containing only the conserved glutaredoxin domain but lacking the NTE. Our aim here is to investigate the effect of the NTE on the conformation of the protein. We therefore solved the NMR structure of the full-length protein, which reveals subtle but significant differences with the structure of the NTE-less form. By means of different experimental approaches, the NTE proved to be intrinsically disordered and not involved in the non-redox dependent protein dimerization, as previously suggested. Interestingly, the portion comprising residues 65-76 of the NTE modulates the conformational dynamics of the glutathione-binding pocket, which may play a role in iron-sulfur cluster assembly and delivery. Furthermore, we disclosed that the class II-strictly conserved loop that precedes the active site is critical for stabilizing the protein structure. So far, this represents the first communication of a Grx containing an intrinsically disordered region that defines a new protein subgroup within class II Grx.

Functional analysis of a tomato (Solanum lycopersicum L.) rhamnogalacturonan lyase promoter.

The enzyme rhamnogalacturonan lyase (RGL) cleaves α-1,4 glycosidic bonds located between rhamnose and galacturonic acid residues in the main chain of rhamnogalacturonan-I (RG-I), a component of the plant cell wall polymer pectin. Although the mode of action of RGL is well known, its physiological functions associated with fruit biology are less understood. Here, we generated transgenic tomato plants expressing the ß-glucuronidase (GUS) reporter gene under the control of a -504 bp or a -776 bp fragment of the promoter of a tomato RGL gene, Solyc11g011300. GUS enzymatic activity and the expression levels of GUS and Solyc11g011300 were measured in a range of organs and fruit developmental stages. GUS staining was undetectable in leaves and roots, but high GUS enzymatic activity was detected in flowers and red ripe (RR) fruits. Maximal expression levels of Solyc11g011300 were detected at the RR developmental stage. GUS activity was 5-fold higher in flowers expressing GUS driven by the -504 bp RGL promoter fragment (RGFL3::GUS) than in the isogenic line, and 1.7-fold higher when GUS gene was driven by the -776 bp RGL promoter fragment (RGLF2::GUS) or the constitutive CaMV35S promoter. Quantitative real-time polymerase chain reaction analysis showed that the highest expression of GUS was in fruits at 40 days after anthesis, for both promoter fragments. The promoter of Solyc11g011300 is predicted to contain cis-acting elements, and to be active in pollen grains, pollen tubes, flowers and during tomato fruit ripening, suggesting that the Solyc11g011300 promoter is transcriptionally active and organ-specific.

Variation in the regulatory region of FZP causes increases in secondary inflorescence branching and grain yield in rice domestication.

Inflorescence branching is a key agronomic trait determining rice yield. The primary branch of the ancestral wild rice (Oryza rufipogon Griff.) bears few grains, due to minimal secondary branching. By contrast, Oryza sativa cultivars have been selected to produce large panicles with more secondary branches. Here we showed that the CONTROL OF SECONDARY BRANCH 1 (COS1) gene, which is identical to FRIZZY PANICLE (FZP), plays an important role in the key transition from few secondary branches in wild rice to more secondary branches in domesticated rice cultivars. A 4-bp tandem repeat deletion approximately 2.7 kb upstream of FZP may affect the binding activities of auxin response factors to the FZP promoter, decrease the expression level of FZP and significantly enhance the number of secondary branches and grain yield in cultivated rice. Functional analyses showed that NARROW LEAF 1 (NAL1), a trypsin-like serine and cysteine protease, interacted with FZP and promoted its degradation. Consistently, downregulating FZP expression or upregulating NAL1 expression in the commercial cultivar Zhonghua 17 increased the number of secondary branches per panicle, grain number per panicle and grain yield per plant. Our findings not only provide insights into the molecular mechanism of increasing grain number and yield during rice domestication, but also offer favorable genes for improving the grain yield of rice.

coTRaCTE predicts co-occurring transcription factors within cell-type specific enhancers.

Cell-type specific gene expression is regulated by the combinatorial action of transcription factors (TFs). In this study, we predict transcription factor (TF) combinations that cooperatively bind in a cell-type specific manner. We first divide DNase hypersensitive sites into cell-type specifically open vs. ubiquitously open sites in 64 cell types to describe possible cell-type specific enhancers. Based on the pattern contrast between these two groups of sequences we develop "co-occurring TF predictor on Cell-Type specific Enhancers" (coTRaCTE) - a novel statistical method to determine regulatory TF co-occurrences. Contrasting the co-binding of TF pairs between cell-type specific and ubiquitously open chromatin guarantees the high cell-type specificity of the predictions. coTRaCTE predicts more than 2000 co-occurring TF pairs in 64 cell types. The large majority (70%) of these TF pairs is highly cell-type specific and overlaps in TF pair co-occurrence are highly consistent among related cell types. Furthermore, independently validated co-occurring and directly interacting TFs are significantly enriched in our predictions. Focusing on the regulatory network derived from the predicted co-occurring TF pairs in embryonic stem cells (ESCs) we find that it consists of three subnetworks with distinct functions: maintenance of pluripotency governed by OCT4, SOX2 and NANOG, regulation of early development governed by KLF4, STAT3, ZIC3 and ZNF148 and general functions governed by MYC, TCF3 and YY1. In summary, coTRaCTE predicts highly cell-type specific co-occurring TFs which reveal new insights into transcriptional regulatory mechanisms.

The immunoglobulin heavy chain 3' regulatory region superenhancer controls mouse B1 B-cell fate and late VDJ repertoire diversity.

The immunoglobulin heavy chain (IgH) 3' regulatory region (3'RR) superenhancer controls B2 B-cell IgH transcription and cell fate at the mature stage but not early repertoire diversity. B1 B cells represent a small percentage of total B cells differing from B2 B cells by several points such as precursors, development, functions, and regulation. B1 B cells act at the steady state to maintain homeostasis in the organism and during the earliest phases of an immune response, setting them at the interface between innate and acquired immunity. We investigated the role of the 3'RR superenhancer on B1 B-cell fate. Similar to B2 B cells, the 3'RR controls µ transcription and cell fate in B1 B cells. In contrast to B2 B cells, 3'RR deletion affects B1 B-cell late repertoire diversity. Thus, differences exist for B1 and B2 B-cell 3'RR control during B-cell maturation. For the first time, these results highlight the contribution of the 3'RR superenhancer at this interface between innate and acquired immunity.

Structural and functional conservation of cis-acting RNA elements in coronavirus 5'-terminal genome regions.

Structure predictions suggest a partial conservation of RNA structure elements in coronavirus terminal genome regions. Here, we determined the structures of stem-loops (SL) 1 and 2 of two alphacoronaviruses, human coronavirus (HCoV) 229E and NL63, by RNA structure probing and studied the functional relevance of these putative cis-acting elements. HCoV-229E SL1 and SL2 mutants generated by reverse genetics were used to study the effects on viral replication of single-nucleotide substitutions predicted to destabilize the SL1 and SL2 structures. The data provide conclusive evidence for the critical role of SL1 and SL2 in HCoV-229E replication and, in some cases, revealed parallels with previously characterized betacoronavirus SL1 and SL2 elements. Also, we were able to rescue viable HCoV-229E mutants carrying replacements of SL2 with equivalent betacoronavirus structural elements. The data obtained in this study reveal a remarkable degree of structural and functional conservation of 5'-terminal RNA structural elements across coronavirus genus boundaries.

[The IgH 3'RR: Doctor Jekyll and Mister Hyde of B-cell maturation and lymphomagenesis]. / 3'RR - Docteur Jekyll et Mister Hyde de la lymphopoïèse/lymphomagenèse B.

The four transcriptional enhancers located in the 3' regulatory region (3'RR) of the IgH locus control the late phases of B-cell maturation, namely IgH locus transcription, somatic hypermutation and class switch recombination. Doctor Jekyll by nature, the 3'RR acts as Mister Hyde in case of oncogenic translocation at the IgH locus taking under its transcriptional control the translocated oncogene. The aim of this review is to show this duality on the basis of the latest scientific advances in the structure and function of the 3'RR and to hIghlight the targeting of the 3'RR as a potential therapeutic approach in mature B-cell lymphomas.