HTSstation: a web application and open-access libraries for high-throughput sequencing data analysis

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

K29-linked free polyubiquitin chains affect ribosome biogenesis and direct ribosomal proteins to the intranuclear quality control compartment

Ribosome assembly requires precise coordination between the production and assembly of ribosomal components. Mutations in ribosomal proteins that inhibit the assembly process or ribosome function are often associated with ribosomopathies, some of which are linked to defects in proteostasis. In this study, we examine the interplay between several yeast proteostasis enzymes, including deubiquitylases (DUBs) Ubp2 and Ubp14, and E3 ligases Ufd4 and Hul5, and we explore their roles in the regulation of the cellular levels of K29-linked unanchored polyubiquitin (polyUb) chains. Accumulating K29-linked unanchored polyUb chains associate with maturing ribosomes to disrupt their assembly, activate the ribosome assembly stress response (RASTR), and lead to the sequestration of ribosomal proteins at the intranuclear quality control compartment (INQ). These findings reveal the physiological relevance of INQ and provide insights into mechanisms of cellular toxicity associated with ribosomopathies.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

K29-linked unanchored polyubiquitin chains disrupt ribosome biogenesis and direct ribosomal proteins to the Intranuclear Quality control compartment (INQ)

Ribosome assembly requires precise coordination between the production and assembly of ribosomal components. Mutations in ribosomal proteins that inhibit the assembly process or ribosome function are often associated with Ribosomopathies, some of which are linked to defects in proteostasis. In this study, we examine the interplay between several yeast proteostasis enzymes, including deubiquitylases (DUBs), Ubp2 and Ubp14, and E3 ligases, Ufd4 and Hul5, and we explore their roles in the regulation of the cellular levels of K29-linked unanchored polyubiquitin (polyUb) chains. Accumulating K29-linked unanchored polyUb chains associate with maturing ribosomes to disrupt their assembly, activate the Ribosome assembly stress response (RASTR), and lead to the sequestration of ribosomal proteins at the Intranuclear Quality control compartment (INQ). These findings reveal the physiological relevance of INQ and provide insights into mechanisms of cellular toxicity associated with Ribosomopathies.

The Drosophila Fab-7 boundary modulates Abd-B gene activity by guiding an inversion of collinear chromatin organization and alternate promoter use

Hox genes encode transcription factors that specify segmental identities along the anteroposterior body axis. These genes are organized in clusters, where their order corresponds to their activity along the body axis, a feature known as collinearity. In Drosophila, the BX-C cluster contains the three most posterior Hox genes, where their collinear activation incorporates progressive changes in histone modifications, chromatin architecture, and use of boundary elements and cis-regulatory regions. To dissect functional hierarchies, we compare chromatin organization in cell lines and larvae, with a focus on the Abd-B gene. Our work establishes the importance of the Fab-7 boundary for insulation between 3D domains carrying different histone modifications. Interestingly, we detect a non-canonical inversion of collinear chromatin dynamics at Abd-B, with the domain of active histone modifications progressively decreasing in size. This dynamic chromatin organization differentially activates the alternative promoters of the Abd-B gene, thereby expanding the possibilities for fine-tuning of transcriptional output.

Integration of 3D genome topology and local chromatin features uncovers enhancers underlying craniofacial-specific cartilage defects

Aberrations in tissue-specific enhancers underlie many developmental defects. Disrupting a noncoding region distal from the human SOX9 gene causes the Pierre Robin sequence (PRS) characterized by the undersized lower jaw. Such a craniofacial-specific defect has been previously linked to enhancers transiently active in cranial neural crest cells (CNCCs). We demonstrate that the PRS region also strongly regulates Sox9 in CNCC-derived Meckel's cartilage (MC), but not in limb cartilages, even after decommissioning of CNCC enhancers. Such an MC-specific regulatory effect correlates with the MC-specific chromatin contacts between the PRS region and Sox9, highlighting the importance of lineage-dependent chromatin topology in instructing enhancer usage. By integrating the enhancer signatures and chromatin topology, we uncovered >10,000 enhancers that function differentially between MC and limb cartilages and demonstrated their association with human diseases. Our findings provide critical insights for understanding the choreography of gene regulation during development and interpreting the genetic basis of craniofacial pathologies.

The Transcriptional Regulator Prdm1 Is Essential for the Early Development of the Sensory Whisker Follicle and Is Linked to the Beta-Catenin First Dermal Signal

Prdm1 mutant mice are one of the rare mutant strains that do not develop whisker hair follicles while still displaying a pelage. Here, we show that Prdm1 is expressed at the earliest stage of whisker development in clusters of mesenchymal cells before placode formation. Its conditional knockout in the murine soma leads to the loss of expression of Bmp2, Shh, Bmp4, Krt17, Edar, and Gli1, though leaving the β-catenin-driven first dermal signal intact. Furthermore, we show that Prdm1 expressing cells not only act as a signaling center but also as a multipotent progenitor population contributing to the several lineages of the adult whisker. We confirm by genetic ablation experiments that the absence of macro vibrissae reverberates on the organization of nerve wiring in the mystacial pads and leads to the reorganization of the barrel cortex. We demonstrate that Lef1 acts upstream of Prdm1 and identify a primate-specific deletion of a Lef1 enhancer named Leaf. This loss may have been significant in the evolutionary process, leading to the progressive defunctionalization and disappearance of vibrissae in primates.

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

Maintenance of mitochondrial integrity in midbrain dopaminergic neurons governed by a conserved developmental transcription factor

Progressive degeneration of dopaminergic (DA) neurons in the substantia nigra is a hallmark of Parkinson’s disease (PD). Dysregulation of developmental transcription factors is implicated in dopaminergic neurodegeneration, but the underlying molecular mechanisms remain largely unknown. Drosophila Fer2 is a prime example of a developmental transcription factor required for the birth and maintenance of midbrain DA neurons. Using an approach combining ChIP-seq, RNA-seq, and genetic epistasis experiments with PD-linked genes, here we demonstrate that Fer2 controls a transcriptional network to maintain mitochondrial structure and function, and thus confers dopaminergic neuroprotection against genetic and oxidative insults. We further show that conditional ablation of Nato3, a mouse homolog of Fer2, in differentiated DA neurons causes mitochondrial abnormalities and locomotor impairments in aged mice. Our results reveal the essential and conserved role of Fer2 homologs in the mitochondrial maintenance of midbrain DA neurons, opening new perspectives for modeling and treating PD.

Mitochondrial dysfunction in dopaminergic neurons is a pathological hallmark of Parkinson’s disease. Here, the authors find a conserved mechanism by which a single transcription factor controls mitochondrial health in dopaminergic neurons during the aging process.

A fast Myosin super enhancer dictates muscle fiber phenotype through competitive interactions with Myosin genes

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here, we identify by snATAC-seq a 42 kb super-enhancer at the locus regrouping the fast Myosin genes. By 4C-seq we show that active fast Myosin promoters interact with this super-enhancer by DNA looping, leading to the activation of a single promoter per nucleus. A rainbow mouse transgenic model of the locus including the super-enhancer recapitulates the endogenous spatio-temporal expression of adult fast Myosin genes. In situ deletion of the super-enhancer by CRISPR/Cas9 editing demonstrates its major role in the control of associated fast Myosin genes, and deletion of two fast Myosin genes at the locus reveals an active competition of the promoters for the shared super-enhancer. Last, by disrupting the organization of fast Myosin, we uncover positional heterogeneity within limb skeletal muscles that may underlie selective muscle susceptibility to damage in certain myopathies.

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here the authors show that a super enhancer controls the spatiotemporal expression of the genes at the fast myosin heavy chain locus by DNA looping and that this expression profile is recapitulated in a rainbow transgenic mouse model of the locus.

Detection of Allele-Specific 3D Chromatin Interactions Using High-Resolution In-Nucleus 4C-seq

Chromosome conformation capture techniques are a set of methods used to determine 3D genome organization through the capture and identification of physical contacts between pairs of genomic loci. Among them, 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing) allows for the identification and quantification of the sequences interacting with a preselected locus of interest. 4C-seq has been widely used in the literature, mainly to study chromatin loops between enhancers and promoters or between CTCF binding sites and to identify chromatin domain boundaries. As 3D-contacts may be established in an allele-specific manner, we describe an up-to-date allele-specific 4C-seq protocol, starting from the selection of allele-specific viewpoints to Illumina sequencing. This protocol has mainly been optimized for cultured mammalian cells, but can be adapted for other cell types with relatively minor changes in initial steps.

Neurofibromin 1 in mushroom body neurons mediates circadian wake drive through activating cAMP-PKA signaling

Various behavioral and cognitive states exhibit circadian variations in animals across phyla including Drosophila melanogaster, in which only ~0.1% of the brain's neurons contain circadian clocks. Clock neurons transmit the timing information to a plethora of non-clock neurons via poorly understood mechanisms. Here, we address the molecular underpinning of this phenomenon by profiling circadian gene expression in non-clock neurons that constitute the mushroom body, the center of associative learning and sleep regulation. We show that circadian clocks drive rhythmic expression of hundreds of genes in mushroom body neurons, including the Neurofibromin 1 (Nf1) tumor suppressor gene and Pka-C1. Circadian clocks also drive calcium rhythms in mushroom body neurons via NF1-cAMP/PKA-C1 signaling, eliciting higher mushroom body activity during the day than at night, thereby promoting daytime wakefulness. These findings reveal the pervasive, non-cell-autonomous circadian regulation of gene expression in the brain and its role in sleep.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

Loss of EZH2-like or SU(VAR)3-9-like proteins causes simultaneous perturbations in H3K27 and H3K9 tri-methylation and associated developmental defects in the fungus Podospora anserina

BACKGROUND

Selective gene silencing is key to development. It is generally accepted that H3K27me3-enriched heterochromatin maintains transcriptional repression established during early development and regulates cell fate. Conversely, H3K9me3-enriched heterochromatin prevents differentiation but constitutes protection against transposable elements. We exploited the fungus Podospora anserina, a valuable alternative to higher eukaryote models, to question the biological relevance and functional interplay of these two distinct heterochromatin conformations.

RESULTS

We established genome-wide patterns of H3K27me3 and H3K9me3 modifications, and found these marks mutually exclusive within gene-rich regions but not within repeats. We generated the corresponding histone methyltransferase null mutants and showed an interdependence of H3K9me3 and H3K27me3 marks. Indeed, removal of the PaKmt6 EZH2-like enzyme resulted not only in loss of H3K27me3 but also in significant H3K9me3 reduction. Similarly, removal of PaKmt1 SU(VAR)3-9-like enzyme caused loss of H3K9me3 and substantial decrease of H3K27me3. Removal of the H3K9me binding protein PaHP1 provided further support to the notion that each type of heterochromatin requires the presence of the other. We also established that P. anserina developmental programs require H3K27me3-mediated silencing, since loss of the PaKmt6 EZH2-like enzyme caused severe defects in most aspects of the life cycle including growth, differentiation processes and sexual reproduction, whereas loss of the PaKmt1 SU(VAR)3-9-like enzyme resulted only in marginal defects, similar to loss of PaHP1.

CONCLUSIONS

Our findings support a conserved function of the PRC2 complex in fungal development. However, we uncovered an intriguing evolutionary fluidity in the repressive histone deposition machinery, which challenges canonical definitions of constitutive and facultative heterochromatin.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

During tetrapod limb development, the HOXA13 and HOXD13 transcription factors are critical for the emergence and organization of the autopod, the most distal aspect where digits will develop. Since previous work had suggested that the Dbx2 gene is a target of these factors, we set up to analyze in detail this potential regulatory interaction.

RESULTS

We show that HOX13 proteins bind to mammalian-specific sequences at the vicinity of the Dbx2 locus that have enhancer activity in developing digits. However, the functional inactivation of the DBX2 protein did not elicit any particular phenotype related to Hox genes inactivation in digits, suggesting either redundant or compensatory mechanisms. We report that the neighboring Nell2 and Ano6 genes are also expressed in distal limb buds and are in part controlled by the same Dbx2 enhancers despite being localized into two different topologically associating domains (TADs) flanking the Dbx2 locus.

CONCLUSIONS

We conclude that Hoxa13 and Hoxd genes cooperatively activate Dbx2 expression in developing digits through binding to mammalian specific regulatory sequences in the Dbx2 neighborhood. Furthermore, these enhancers can overcome TAD boundaries in either direction to co-regulate a set of genes located in distinct chromatin domains.

Loop extrusion as a mechanism for formation of DNA damage repair foci

The repair of DNA double-strand breaks (DSBs) is essential for safeguarding genome integrity. When a DSB forms, the PI3K-related ATM kinase rapidly triggers the establishment of megabase-sized, chromatin domains decorated with phosphorylated histone H2AX (γH2AX), which act as seeds for the formation of DNA-damage response foci1. It is unclear how these foci are rapidly assembled to establish a 'repair-prone' environment within the nucleus. Topologically associating domains are a key feature of 3D genome organization that compartmentalize transcription and replication, but little is known about their contribution to DNA repair processes2,3. Here we show that topologically associating domains are functional units of the DNA damage response, and are instrumental for the correct establishment of γH2AX-53BP1 chromatin domains in a manner that involves one-sided cohesin-mediated loop extrusion on both sides of the DSB. We propose a model in which H2AX-containing nucleosomes are rapidly phosphorylated as they actively pass by DSB-anchored cohesin. Our work highlights the importance of chromosome conformation in the maintenance of genome integrity and demonstrates the establishment of a chromatin modification by loop extrusion.

Oscillating and stable genome topologies underlie hepatic physiological rhythms during the circadian cycle

The circadian clock drives extensive temporal gene expression programs controlling daily changes in behavior and physiology. In mouse liver, transcription factors dynamics, chromatin modifications, and RNA Polymerase II (PolII) activity oscillate throughout the 24-hour (24h) day, regulating the rhythmic synthesis of thousands of transcripts. Also, 24h rhythms in gene promoter-enhancer chromatin looping accompany rhythmic mRNA synthesis. However, how chromatin organization impinges on temporal transcription and liver physiology remains unclear. Here, we applied time-resolved chromosome conformation capture (4C-seq) in livers of WT and arrhythmic Bmal1 knockout mice. In WT, we observed 24h oscillations in promoter-enhancer loops at multiple loci including the core-clock genes Period1, Period2 and Bmal1. In addition, we detected rhythmic PolII activity, chromatin modifications and transcription involving stable chromatin loops at clock-output gene promoters representing key liver function such as glucose metabolism and detoxification. Intriguingly, these contacts persisted in clock-impaired mice in which both PolII activity and chromatin marks no longer oscillated. Finally, we observed chromatin interaction hubs connecting neighbouring genes showing coherent transcription regulation across genotypes. Thus, both clock-controlled and clock-independent chromatin topology underlie rhythmic regulation of liver physiology.

Mechanisms coordinating ribosomal protein gene transcription in response to stress

While expression of ribosomal protein genes (RPGs) in the budding yeast has been extensively studied, a longstanding enigma persists regarding their co-regulation under fluctuating growth conditions. Most RPG promoters display one of two distinct arrangements of a core set of transcription factors (TFs) and are further differentiated by the presence or absence of the HMGB protein Hmo1. However, a third group of promoters appears not to be bound by any of these proteins, raising the question of how the whole suite of genes is co-regulated. We demonstrate here that all RPGs are regulated by two distinct, but complementary mechanisms driven by the TFs Ifh1 and Sfp1, both of which are required for maximal expression in optimal conditions and coordinated downregulation upon stress. At the majority of RPG promoters, Ifh1-dependent regulation predominates, whereas Sfp1 plays the major role at all other genes. We also uncovered an unexpected protein homeostasis-dependent binding property of Hmo1 at RPG promoters. Finally, we show that the Ifh1 paralog Crf1, previously described as a transcriptional repressor, can act as a constitutive RPG activator. Our study provides a more complete picture of RPG regulation and may serve as a paradigm for unravelling RPG regulation in multicellular eukaryotes.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

Three-dimensional chromatin interactions remain stable upon CAG/CTG repeat expansion

Expanded CAG/CTG repeats underlie 13 neurological disorders, including myotonic dystrophy type 1 (DM1) and Huntington's disease (HD). Upon expansion, disease loci acquire heterochromatic characteristics, which may provoke changes to chromatin conformation and thereby affect both gene expression and repeat instability. Here, we tested this hypothesis by performing 4C sequencing at the DMPK and HTT loci from DM1 and HD-derived cells. We find that allele sizes ranging from 15 to 1700 repeats displayed similar chromatin interaction profiles. This was true for both loci and for alleles with different DNA methylation levels and CTCF binding. Moreover, the ectopic insertion of an expanded CAG repeat tract did not change the conformation of the surrounding chromatin. We conclude that CAG/CTG repeat expansions are not enough to alter chromatin conformation in cis. Therefore, it is unlikely that changes in chromatin interactions drive repeat instability or changes in gene expression in these disorders.

A complex regulatory landscape involved in the development of mammalian external genitals

Developmental genes are often controlled by large regulatory landscapes matching topologically associating domains (TADs). In various contexts, the associated chromatin backbone is modified by specific enhancer-enhancer and enhancer-promoter interactions. We used a TAD flanking the mouse HoxD cluster to study how these regulatory architectures are formed and deconstructed once their function achieved. We describe this TAD as a functional unit, with several regulatory sequences acting together to elicit a transcriptional response. With one exception, deletion of these sequences didn't modify the transcriptional outcome, a result at odds with a conventional view of enhancer function. The deletion and inversion of a CTCF site located near these regulatory sequences did not affect transcription of the target gene. Slight modifications were nevertheless observed, in agreement with the loop extrusion model. We discuss these unexpected results considering both conventional and alternative explanations relying on the accumulation of poorly specific factors within the TAD backbone.

Lymphatic endothelial cells prime naïve CD8+ T cells into memory cells under steady-state conditions

Lymphatic endothelial cells (LECs) chemoattract naïve T cells and promote their survival in the lymph nodes, and can cross-present antigens to naïve CD8+ T cells to drive their proliferation despite lacking key costimulatory molecules. However, the functional consequence of LEC priming of CD8+ T cells is unknown. Here, we show that while many proliferating LEC-educated T cells enter early apoptosis, the remainders comprise a long-lived memory subset, with transcriptional, metabolic, and phenotypic features of central memory and stem cell-like memory T cells. In vivo, these memory cells preferentially home to lymph nodes and display rapid proliferation and effector differentiation following memory recall, and can protect mice against a subsequent bacterial infection. These findings introduce a new immunomodulatory role for LECs in directly generating a memory-like subset of quiescent yet antigen-experienced CD8+ T cells that are long-lived and can rapidly differentiate into effector cells upon inflammatory antigenic challenge.

Neutrophils suppress tumor-infiltrating T cells in colon cancer via matrix metalloproteinase-mediated activation of TGFβ

High T-cell infiltration in colorectal cancer (CRC) correlates with a favorable disease outcome and immunotherapy response. This, however, is only observed in a small subset of CRC patients. A better understanding of the factors influencing tumor T-cell responses in CRC could inspire novel therapeutic approaches to achieve broader immunotherapy responsiveness. Here, we investigated T cell-suppressive properties of different myeloid cell types in an inducible colon tumor mouse model. The most potent inhibitors of T-cell activity were tumor-infiltrating neutrophils. Gene expression analysis and combined in vitro and in vivo tests indicated that T-cell suppression is mediated by neutrophil-secreted metalloproteinase activation of latent TGFβ. CRC patient neutrophils similarly suppressed T cells via TGFβ in vitro, and public gene expression datasets suggested that T-cell activity is lowest in CRCs with combined neutrophil infiltration and TGFβ activation. Thus, the interaction of neutrophils with a TGFβ-rich tumor microenvironment may represent a conserved immunosuppressive mechanism in CRC.

CTCF modulates allele-specific sub-TAD organization and imprinted gene activity at the mouse Dlk1-Dio3 and Igf2-H19 domains

BACKGROUND

Genomic imprinting is essential for mammalian development and provides a unique paradigm to explore intra-cellular differences in chromatin configuration. So far, the detailed allele-specific chromatin organization of imprinted gene domains has mostly been lacking. Here, we explored the chromatin structure of the two conserved imprinted domains controlled by paternal DNA methylation imprints-the Igf2-H19 and Dlk1-Dio3 domains-and assessed the involvement of the insulator protein CTCF in mouse cells.

RESULTS

Both imprinted domains are located within overarching topologically associating domains (TADs) that are similar on both parental chromosomes. At each domain, a single differentially methylated region is bound by CTCF on the maternal chromosome only, in addition to multiple instances of bi-allelic CTCF binding. Combinations of allelic 4C-seq and DNA-FISH revealed that bi-allelic CTCF binding alone, on the paternal chromosome, correlates with a first level of sub-TAD structure. On the maternal chromosome, additional CTCF binding at the differentially methylated region adds a further layer of sub-TAD organization, which essentially hijacks the existing paternal-specific sub-TAD organization. Perturbation of maternal-specific CTCF binding site at the Dlk1-Dio3 locus, using genome editing, results in perturbed sub-TAD organization and bi-allelic Dlk1 activation during differentiation.

CONCLUSIONS

Maternal allele-specific CTCF binding at the imprinted Igf2-H19 and the Dlk1-Dio3 domains adds an additional layer of sub-TAD organization, on top of an existing three-dimensional configuration and prior to imprinted activation of protein-coding genes. We speculate that this allele-specific sub-TAD organization provides an instructive or permissive context for imprinted gene activation during development.

Conditional effects of the epigenetic regulator JUMONJI 14 in Arabidopsis root growth

Methylation of lysine 4 in histone 3 (H3K4) is a post-translational modification that promotes gene expression. H3K4 methylation can be reversed by specific demethylases with an enzymatic Jumonji C domain. In Arabidopsis thaliana, H3K4-specific JUMONJI (JMJ) proteins distinguish themselves by the association with an F/Y-rich (FYR) domain. Here, we report that jmj14 mutations partially suppress reduced root meristem size and growth vigor of brevis radix (brx) mutants. Similar to its close homologs, JMJ15, JMJ16 and JMJ18, the JMJ14 promoter confers expression in mature root vasculature. Yet, unlike jmj14, neither jmj16 nor jmj18 mutation markedly suppresses brx phenotypes. Domain-swapping experiments suggest that the specificity of JMJ14 function resides in the FYR domain. Despite JMJ14 promoter activity in the mature vasculature, jmj14 mutation affects root meristem size. However, JMJ14 protein is observed throughout the meristem, suggesting that the JMJ14 transcript region contributes substantially to the spatial aspect of JMJ14 expression. In summary, our data reveal a role for JMJ14 in root growth in sensitized genetic backgrounds that depends on its FYR domain and regulatory input from the JMJ14 cistron.

Fork pausing complex engages topoisomerases at the replisome

In this study, Shyian et al. set out to address mechanistically how the evolutionarily conserved fork pausing complex acts at proteinaceous replication fork barriers (RFBs) to promote fork passage and genome stability. Using several molecular and cell-based assays, the authors propose that forks pause at proteinaceous RFBs through a “sTOP” mechanism (“slowing down with topoisomerases I–II”), which also contributes to protecting cells from topoisomerase-blocking agents.

Replication forks temporarily or terminally pause at hundreds of hard-to-replicate regions around the genome. A conserved pair of budding yeast replisome components Tof1–Csm3 (fission yeast Swi1–Swi3 and human TIMELESS–TIPIN) act as a “molecular brake” and promote fork slowdown at proteinaceous replication fork barriers (RFBs), while the accessory helicase Rrm3 assists the replisome in removing protein obstacles. Here we show that the Tof1–Csm3 complex promotes fork pausing independently of Rrm3 helicase by recruiting topoisomerase I (Top1) to the replisome. Topoisomerase II (Top2) partially compensates for the pausing decrease in cells when Top1 is lost from the replisome. The C terminus of Tof1 is specifically required for Top1 recruitment to the replisome and fork pausing but not for DNA replication checkpoint (DRC) activation. We propose that forks pause at proteinaceous RFBs through a “sTOP” mechanism (“slowing down with topoisomerases I–II”), which we show also contributes to protecting cells from topoisomerase-blocking agents.

Opposing chromatin remodelers control transcription initiation frequency and start site selection

Precise nucleosome organization at eukaryotic promoters is thought to be generated by multiple chromatin remodeler (CR) enzymes and to affect transcription initiation. Using an integrated analysis of chromatin remodeler binding and nucleosome occupancy following rapid remodeler depletion, we investigated the interplay between these enzymes and their impact on transcription in yeast. We show that many promoters are affected by multiple CRs that operate in concert or in opposition to position the key transcription start site (TSS)-associated +1 nucleosome. We also show that nucleosome movement after CR inactivation usually results from the activity of another CR and that in the absence of any remodeling activity, +1 nucleosomes largely maintain their positions. Finally, we present functional assays suggesting that +1 nucleosome positioning often reflects a trade-off between maximizing RNA polymerase recruitment and minimizing transcription initiation at incorrect sites. Our results provide a detailed picture of fundamental mechanisms linking promoter nucleosome architecture to transcription initiation.

Impact of genome architecture on the functional activation and repression of Hox regulatory landscapes

BACKGROUND

The spatial organization of the mammalian genome relies upon the formation of chromatin domains of various scales. At the level of gene regulation in cis, collections of enhancer sequences define large regulatory landscapes that usually match with the presence of topologically associating domains (TADs). These domains often contain ranges of enhancers displaying similar or related tissue specificity, suggesting that in some cases, such domains may act as coherent regulatory units, with a global on or off state. By using the HoxD gene cluster, which specifies the topology of the developing limbs via highly orchestrated regulation of gene expression, as a paradigm, we investigated how the arrangement of regulatory domains determines their activity and function.

RESULTS

Proximal and distal cells in the developing limb express different levels of Hoxd genes, regulated by flanking 3' and 5' TADs, respectively. We characterized the effect of large genomic rearrangements affecting these two TADs, including their fusion into a single chromatin domain. We show that, within a single hybrid TAD, the activation of both proximal and distal limb enhancers globally occurred as when both TADs are intact. However, the activity of the 3' TAD in distal cells is generally increased in the fused TAD, when compared to wild type where it is silenced. Also, target gene activity in distal cells depends on whether or not these genes had previously responded to proximal enhancers, which determines the presence or absence of H3K27me3 marks. We also show that the polycomb repressive complex 2 is mainly recruited at the Hox gene cluster and can extend its coverage to far-cis regulatory sequences as long as confined to the neighboring TAD structure.

CONCLUSIONS

We conclude that antagonistic limb proximal and distal enhancers can exert their specific effects when positioned into the same TAD and in the absence of their genuine target genes. We also conclude that removing these target genes reduced the coverage of a regulatory landscape by chromatin marks associated with silencing, which correlates with its prolonged activity in time.

A ribosome assembly stress response regulates transcription to maintain proteome homeostasis

Ribosome biogenesis is a complex and energy-demanding process requiring tight coordination of ribosomal RNA (rRNA) and ribosomal protein (RP) production. Given the extremely high level of RP synthesis in rapidly growing cells, alteration of any step in the ribosome assembly process may impact growth by leading to proteotoxic stress. Although the transcription factor Hsf1 has emerged as a central regulator of proteostasis, how its activity is coordinated with ribosome biogenesis is unknown. Here, we show that arrest of ribosome biogenesis in the budding yeast Saccharomyces cerevisiae triggers rapid activation of a highly specific stress pathway that coordinately upregulates Hsf1 target genes and downregulates RP genes. Activation of Hsf1 target genes requires neo-synthesis of RPs, which accumulate in an insoluble fraction and presumably titrate a negative regulator of Hsf1, the Hsp70 chaperone. RP aggregation is also coincident with that of the RP gene activator Ifh1, a transcription factor that is rapidly released from RP gene promoters. Our data support a model in which the levels of newly synthetized RPs, imported into the nucleus but not yet assembled into ribosomes, work to continuously balance Hsf1 and Ifh1 activity, thus guarding against proteotoxic stress during ribosome assembly.

When yeast cells are growing at top speed, they can make 2,000 new ribosomes every minute. These enormous molecular assemblies are the protein-making machines of the cell. Building new ribosomes is one of the most energy-demanding parts of cell growth and, if the process goes wrong, the results can be catastrophic. The proteins that make up the ribosomes themselves are sticky. Left unattended, they start to form toxic clumps inside the compartment that houses most of the cell’s DNA, the nucleus.

A protein called Heat shock factor 1, or Hsf1 for short, plays an important role in the cell's quality control systems. It helps to manage sticky proteins by switching on genes that break down protein clumps and prevent new clumps from forming. Hsf1 levels start to rise whenever cells are struggling to keep up with protein production. If it is half-finished ribosomes that are causing the problem, cells can stop making ribosome proteins. The protein in charge of this in yeast is Ifh1. It is a transcription factor that sits at the front of the genes for ribosome proteins, switching them on. When yeast cells get stressed, Ifh1 drops away from the genes within minutes, switching them off again. Yet how this happens, and how it links to Hsf1, is a mystery.

To start to provide some answers, Albert et al. disrupted the production of ribosomes in yeast cells and examined the consequences. This revealed a new rescue response, that they named the “ribosome assembly stress response”. Both Hsf1 and Ifh1 are sensitive to the build-up of unfinished ribosomes in the nucleus. As expected, Hsf1 activated when ribosome proteins started to build up, and switched on the genes needed to manage the protein clumps. The effect on Isfh1, however, was unexpected. When the unassembled ribosome proteins started to build up, it was the clumps themselves that pulled the Ifh1 proteins off the genes. The unassembled ribosomes proteins seemed to be stopping their own production. Low levels of clumped ribosome proteins in the nuclei of unstressed cells also helped to keep Hsf1 active and pull Ifh1 off the ribosome genes. It is possible that this provides continual protection against a toxic protein build-up.

These findings are not only important for understanding yeast cells; cancer cells also need to produce ribosomes at a very high rate to sustain their rapid growth. They too might be prone to stresses that interrupt their ribosome assembly. As such, understanding more about this process could one day lead to new therapies to target cancer cells.

Sequence-Directed Action of RSC Remodeler and General Regulatory Factors Modulates +1 Nucleosome Position to Facilitate Transcription

Accessible chromatin is important for RNA polymerase II recruitment and transcription initiation at eukaryotic promoters. We investigated the mechanistic links between promoter DNA sequence, nucleosome positioning, and transcription. Our results indicate that positioning of the transcription start site-associated +1 nucleosome in yeast is critical for efficient TBP binding and is driven by two key factors, the essential chromatin remodeler RSC and a small set of ubiquitous general regulatory factors (GRFs). Our findings indicate that the strength and directionality of RSC action on promoter nucleosomes depends on the arrangement and proximity of two specific DNA motifs. This, together with the effect on nucleosome position observed in double depletion experiments, suggests that, despite their widespread co-localization, RSC and GRFs predominantly act through independent signals to generate accessible chromatin. Our results provide mechanistic insight into how the promoter DNA sequence instructs trans-acting factors to control nucleosome architecture and stimulate transcription initiation.

An extended regulatory landscape drives Tbx18 activity in a variety of prostate-associated cell lineages

The evolutionarily conserved transcription factor, Tbx18, is expressed in a dynamic pattern throughout embryonic and early postnatal life and plays crucial roles in the development of multiple organ systems. Previous studies have indicated that this dynamic function is controlled by an expansive regulatory structure, extending far upstream and downstream of the gene. With the goal of identifying elements that interact with the Tbx18 promoter in developing prostate, we coupled chromatin conformation capture (4C) and ATAC-seq from embryonic day 18.5 (E18.5) mouse urogenital sinus (UGS), where Tbx18 is highly expressed. The data revealed dozens of active chromatin elements distributed throughout a 1.5 million base pair topologically associating domain (TAD). To identify cell types contributing to this chromatin signal, we used lineage tracing methods with a Tbx18 Cre "knock-in" allele; these data show clearly that Tbx18-expressing precursors differentiate into wide array of cell types in multiple tissue compartments, most of which have not been previously reported. We also used a 209 kb Cre-expressing Tbx18 transgene, to partition enhancers for specific precursor types into two rough spatial domains. Within this central 209 kb compartment, we identified ECR1, previously described to regulate Tbx18 expression in ureter, as an active regulator of UGS expression. Together these data define the diverse fates of Tbx18+ precursors in prostate-associated tissues for the first time, and identify a highly active TAD controlling the gene's essential function in this tissue.

Noncoding transcription influences the replication initiation program through chromatin regulation

In eukaryotic organisms, replication initiation follows a temporal program. Among the parameters that regulate this program in Saccharomyces cerevisiae, chromatin structure has been at the center of attention without considering the contribution of transcription. Here, we revisit the replication initiation program in the light of widespread genomic noncoding transcription. We find that noncoding RNA transcription termination in the vicinity of autonomously replicating sequences (ARSs) shields replication initiation from transcriptional readthrough. Consistently, high natural nascent transcription correlates with low ARS efficiency and late replication timing. High readthrough transcription is also linked to increased nucleosome occupancy and high levels of H3K36me3. Moreover, forcing ARS readthrough transcription promotes these chromatin features. Finally, replication initiation defects induced by increased transcriptional readthrough are partially rescued in the absence of H3K36 methylation. Altogether, these observations indicate that natural noncoding transcription into ARSs influences replication initiation through chromatin regulation.

Similarities and differences in the regulation of HoxD genes during chick and mouse limb development

In all tetrapods examined thus far, the development and patterning of limbs require the activation of gene members of the HoxD cluster. In mammals, they are regulated by a complex bimodal process that controls first the proximal patterning and then the distal structure. During the shift from the former to the latter regulation, this bimodal regulatory mechanism allows the production of a domain with low Hoxd gene expression, at which both telomeric (T-DOM) and centromeric regulatory domains (C-DOM) are silent. These cells generate the future wrist and ankle articulations. We analyzed the implementation of this regulatory mechanism in chicken, i.e., in an animal for which large morphological differences exist between fore- and hindlimbs. We report that although this bimodal regulation is globally conserved between the mouse and the chick, some important modifications evolved at least between these two model systems, in particular regarding the activity of specific enhancers, the width of the TAD boundary separating the two regulations, and the comparison between the forelimb versus hindlimb regulatory controls. At least one aspect of these regulations seems to be more conserved between chick and bats than with mouse, which may relate to the extent to which forelimbs and hindlimbs of these various animals differ in their morphologies.

A comparison of Hox gene regulation during the development of limbs in birds and mammals reveals that whereas the characteristic bimodal regulatory system, based on large chromatin domains, is largely conserved between these morphologically distinct structures, some differences are revealed in the way this is implemented in various vertebrates.

The shapes of limbs vary greatly among tetrapod species, even between the forelimbs and hindlimbs of the same animal. Hox genes regulate the proper growth and patterning of tetrapod limbs. In order to evaluate whether variations in the complex regulation of a cluster of Hox genes—the Hoxd genes—during limb development contribute to the differences in limb shape, we compared their transcriptional control during limb bud development in the forelimbs and hindlimbs of mouse and chicken embryos. We found that the regulatory mechanism underlying Hoxd gene expression is highly conserved, but some clear differences exist. For instance, we observed a variation in the topologically associating domain (TAD; a self-interacting genomic region) boundary interval between the mouse and the chick, as well as differences in the activity of a conserved enhancer element situated within the telomeric regulatory domain. In contrast to the mouse, the chicken enhancer has a stronger activity in the forelimb buds than in the hindlimb buds, which is correlated with the striking differences in the mRNA levels of the genes. We conclude that differences in both the timing and duration of TAD activities and in the width of their boundary may parallel the important decrease in Hoxd gene transcription in chick hindlimb buds versus forelimb buds. These differences may also account for the slightly distinct regulatory strategies implemented by mammals and birds at this locus.

Oestrogen receptor α AF-1 and AF-2 domains have cell population-specific functions in the mammary epithelium

Oestrogen receptor α (ERα) is a transcription factor with ligand-independent and ligand-dependent activation functions (AF)-1 and -2. Oestrogens control postnatal mammary gland development acting on a subset of mammary epithelial cells (MECs), termed sensor cells, which are ERα-positive by immunohistochemistry (IHC) and secrete paracrine factors, which stimulate ERα-negative responder cells. Here we show that deletion of AF-1 or AF-2 blocks pubertal ductal growth and subsequent development because both are required for expression of essential paracrine mediators. Thirty percent of the luminal cells are ERα-negative by IHC but express Esr1 transcripts. This low level ERα expression through AF-2 is essential for cell expansion during puberty and growth-inhibitory during pregnancy. Cell-intrinsic ERα is not required for cell proliferation nor for secretory differentiation but controls transcript levels of cell motility and cell adhesion genes and a stem cell and epithelial mesenchymal transition (EMT) signature identifying ERα as a key regulator of mammary epithelial cell plasticity.

Oestrogen receptors α (ERα) are expressed in a subset of mammary epithelial cells. Here, the authors identify cells with low-ERα protein levels and show that distinct cell populations have distinct requirements for the AF1 and AF2 domains of the ERα, and ERα acts in a biphasic manner dependent on developmental stage.

TFAP2 transcription factors are regulators of lipid droplet biogenesis

How trafficking pathways and organelle abundance adapt in response to metabolic and physiological changes is still mysterious, although a few transcriptional regulators of organellar biogenesis have been identified in recent years. We previously found that the Wnt signaling directly controls lipid droplet formation, linking the cell storage capacity to the established functions of Wnt in development and differentiation. In the present paper, we report that Wnt-induced lipid droplet biogenesis does not depend on the canonical TCF/LEF transcription factors. Instead, we find that TFAP2 family members mediate the pro-lipid droplet signal induced by Wnt3a, leading to the notion that the TFAP2 transcription factor may function as a 'master' regulator of lipid droplet biogenesis.

Heterogeneous combinatorial expression of Hoxd genes in single cells during limb development

BACKGROUND

Global analyses of gene expression during development reveal specific transcription patterns associated with the emergence of various cell types, tissues, and organs. These heterogeneous patterns are instrumental to ensure the proper formation of the different parts of our body, as shown by the phenotypic effects generated by functional genetic approaches. However, variations at the cellular level can be observed within each structure or organ. In the developing mammalian limbs, expression of Hox genes from the HoxD cluster is differentially controlled in space and time, in cells that will pattern the digits and the forearms. While the Hoxd genes broadly share a common regulatory landscape and large-scale analyses have suggested a homogenous Hox gene transcriptional program, it has not previously been clear whether Hoxd genes are expressed together at the same levels in the same cells.

RESULTS

We report a high degree of heterogeneity in the expression of the Hoxd11 and Hoxd13 genes. We analyzed single-limb bud cell transcriptomes and show that Hox genes are expressed in specific combinations that appear to match particular cell types. In cells giving rise to digits, we find that the expression of the five relevant Hoxd genes (Hoxd9 to Hoxd13) is unbalanced, despite their control by known global enhancers. We also report that specific combinatorial expression follows a pseudo-time sequence, which is established based on the transcriptional diversity of limb progenitors.

CONCLUSIONS

Our observations reveal the existence of distinct combinations of Hoxd genes at the single-cell level during limb development. In addition, we document that the increasing combinatorial expression of Hoxd genes in this developing structure is associated with specific transcriptional signatures and that these signatures illustrate a temporal progression in the differentiation of these cells.

Distinct patterns of histone acetyltransferase and Mediator deployment at yeast protein-coding genes

Here, Bruzzone et al. explore the relative contributions of the transcriptional coactivators Mediator and two histone acetyltransferase (HAT) complexes, NuA4 and SAGA, to RNA polymerase II association at specific genes and gene classes by rapid nuclear depletion of key complex subunits. They show that the NuA4 HAT Esa1 differentially affects certain groups of genes, whereas the SAGA HAT Gcn5 has a weaker but more uniform effect, and their findings suggest that at least three distinct combinations of coactivator deployment are used to generate moderate or high transcription levels.

The transcriptional coactivators Mediator and two histone acetyltransferase (HAT) complexes, NuA4 and SAGA, play global roles in transcriptional activation. Here we explore the relative contributions of these factors to RNA polymerase II association at specific genes and gene classes by rapid nuclear depletion of key complex subunits. We show that the NuA4 HAT Esa1 differentially affects certain groups of genes, whereas the SAGA HAT Gcn5 has a weaker but more uniform effect. Relative dependence on Esa1 and Tra1, a shared component of NuA4 and SAGA, distinguishes two large groups of coregulated growth-promoting genes. In contrast, we show that the activity of Mediator is particularly important at a separate, small set of highly transcribed TATA-box-containing genes. Our analysis indicates that at least three distinct combinations of coactivator deployment are used to generate moderate or high transcription levels and suggests that each may be associated with distinct forms of regulation.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

Cross-Resistance of UV- or Chlorine Dioxide-Resistant Echovirus 11 to Other Disinfectants

The emergence of waterborne viruses with resistance to disinfection has been demonstrated in the laboratory and in the environment. Yet, the implications of such resistance for virus control remain obscure. In this study we investigate if viruses with resistance to a given disinfection method exhibit cross-resistance to other disinfectants. Chlorine dioxide (ClO₂)- or UV-resistant populations of echovirus 11 were exposed to five inactivating treatments (free chlorine, ClO₂, UV radiation, sunlight, and heat), and the extent of cross-resistance was determined. The ClO₂-resistant population exhibited cross-resistance to free chlorine, but to none of the other inactivating treatments tested. We furthermore demonstrated that ClO₂ and free chlorine act by a similar mechanism, in that they mainly inhibit the binding of echovirus 11 to its host cell. As such, viruses with host binding mechanisms that can withstand ClO₂ treatment were also better able to withstand oxidation by free chlorine. Conversely, the UV-resistant population was not significantly cross-resistant to any other disinfection treatment. Overall, our results indicate that viruses with resistance to multiple disinfectants exist, but that they can be controlled by inactivating methods that operate by a distinctly different mechanism. We therefore suggest to utilize two disinfection barriers that act by different mechanisms in order to control disinfection-resistant viruses.

Large scale genomic reorganization of topological domains at the HoxD locus

BACKGROUND

The transcriptional activation of HoxD genes during mammalian limb development involves dynamic interactions with two topologically associating domains (TADs) flanking the HoxD cluster. In particular, the activation of the most posterior HoxD genes in developing digits is controlled by regulatory elements located in the centromeric TAD (C-DOM) through long-range contacts.

RESULTS

To assess the structure-function relationships underlying such interactions, we measured compaction levels and TAD discreteness using a combination of chromosome conformation capture (4C-seq) and DNA FISH. We assessed the robustness of the TAD architecture by using a series of genomic deletions and inversions that impact the integrity of this chromatin domain and that remodel long-range contacts. We report multi-partite associations between HoxD genes and up to three enhancers. We find that the loss of native chromatin topology leads to the remodeling of TAD structure following distinct parameters.

CONCLUSIONS

Our results reveal that the recomposition of TAD architectures after large genomic re-arrangements is dependent on a boundary-selection mechanism in which CTCF mediates the gating of long-range contacts in combination with genomic distance and sequence specificity. Accordingly, the building of a recomposed TAD at this locus depends on distinct functional and constitutive parameters.

The SAGA complex, together with transcription factors and the endocytic protein Rvs167p, coordinates the reprofiling of gene expression in response to changes in sterol composition in Saccharomyces cerevisiae

The SAGA complex, together with transcription factors and Rvs167p, coordinates sterol-dependent transcription changes. In ergosterol mutants the SAGA complex increases its occupancy on ergosterol biosynthesis and anaerobic gene promoters, recruits the SWI/SNF complex, and binds to transcription factors and Rvs167p. Genes encoding stress proteins and basic amino acid synthesis are also affected even though promoter occupancy is not changed.

Changes in cellular sterol species and concentrations can have profound effects on the transcriptional profile. In yeast, mutants defective in sterol biosynthesis show a wide range of changes in transcription, including a coinduction of anaerobic genes and ergosterol biosynthesis genes, biosynthesis of basic amino acids, and several stress genes. However the mechanisms underlying these changes are unknown. We identified mutations in the SAGA complex, a coactivator of transcription, which abrogate the ability to carry out most of these sterol-dependent transcriptional changes. In the erg3 mutant, the SAGA complex increases its occupancy time on many of the induced ergosterol and anaerobic gene promoters, increases its association with several relevant transcription factors and the SWI/SNF chromatin remodeling complex, and surprisingly, associates with an endocytic protein, Rvs167p, suggesting a moonlighting function for this protein in the sterol-regulated induction of the heat shock protein, HSP42 and HSP102, mRNAs.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

The formation of multicellular organisms from the egg to the adult stage is largely under genetic control. The activation of specific genes is governed by regulatory DNA sequences present nearby on the chromosome. Most of these sequences promote activation and are called enhancers. In this paper, we study two enhancers governing the expression of a gene involved in the formation of the posterior brain in vertebrates. One of these enhancers is involved in a positive feedback loop: it is itself activated by the protein product of the gene that it regulates. The other enhancer was thought to be only involved in the initial accumulation of the protein, necessary for the subsequent activation of the feedback loop. Here we show that the second enhancer directly cooperates with the autoregulatory enhancer to increase its accessibility and its activity. Our work uncovers a novel, long-range mode of cooperation between enhancers that restricts the domain of action of autoregulatory enhancers within embryos and might be essential to avoid their inappropriate activation.

Chromosomal contacts connect loci associated with autism, BMI and head circumference phenotypes

Copy number variants (CNVs) are major contributors to genomic imbalance disorders. Phenotyping of 137 unrelated deletion and reciprocal duplication carriers of the distal 16p11.2 220 kb BP2-BP3 interval showed that these rearrangements are associated with autism spectrum disorders and mirror phenotypes of obesity/underweight and macrocephaly/microcephaly. Such phenotypes were previously associated with rearrangements of the non-overlapping proximal 16p11.2 600 kb BP4-BP5 interval. These two CNV-prone regions at 16p11.2 are reciprocally engaged in complex chromatin looping, as successfully confirmed by 4C-seq, fluorescence in situ hybridization and Hi-C, as well as coordinated expression and regulation of encompassed genes. We observed that genes differentially expressed in 16p11.2 BP4-BP5 CNV carriers are concomitantly modified in their chromatin interactions, suggesting that disruption of chromatin interplays could participate in the observed phenotypes. We also identified cis- and trans-acting chromatin contacts to other genomic regions previously associated with analogous phenotypes. For example, we uncovered that individuals with reciprocal rearrangements of the trans-contacted 2p15 locus similarly display mirror phenotypes on head circumference and weight. Our results indicate that chromosomal contacts' maps could uncover functionally and clinically related genes.

Transcription facilitated genome-wide recruitment of topoisomerase I and DNA gyrase

Movement of the transcription machinery along a template alters DNA topology resulting in the accumulation of supercoils in DNA. The positive supercoils generated ahead of transcribing RNA polymerase (RNAP) and the negative supercoils accumulating behind impose severe topological constraints impeding transcription process. Previous studies have implied the role of topoisomerases in the removal of torsional stress and the maintenance of template topology but the in vivo interaction of functionally distinct topoisomerases with heterogeneous chromosomal territories is not deciphered. Moreover, how the transcription-induced supercoils influence the genome-wide recruitment of DNA topoisomerases remains to be explored in bacteria. Using ChIP-Seq, we show the genome-wide occupancy profile of both topoisomerase I and DNA gyrase in conjunction with RNAP in Mycobacterium tuberculosis taking advantage of minimal topoisomerase representation in the organism. The study unveils the first in vivo genome-wide interaction of both the topoisomerases with the genomic regions and establishes that transcription-induced supercoils govern their recruitment at genomic sites. Distribution profiles revealed co-localization of RNAP and the two topoisomerases on the active transcriptional units (TUs). At a given locus, topoisomerase I and DNA gyrase were localized behind and ahead of RNAP, respectively, correlating with the twin-supercoiled domains generated. The recruitment of topoisomerases was higher at the genomic loci with higher transcriptional activity and/or at regions under high torsional stress compared to silent genomic loci. Importantly, the occupancy of DNA gyrase, sole type II topoisomerase in Mtb, near the Ter domain of the Mtb chromosome validates its function as a decatenase.

The generation of DNA topological constraint is intrinsic to transcription. Although in vitro studies indicated DNA gyrase and topoisomerase I action in the removal of excess supercoils, ahead and behind the transcribing RNA polymerase, in vivo recruitment and interaction of both topoisomerases with the genome has not been investigated. Using advanced sequencing, we have mapped the genome-wide footprints of topoisomerase I and DNA gyrase along with RNAP in deadly pathogen, Mycobacterium tuberculosis. We show that in vivo distribution of topoisomerases is guided by active transcription and both the enzymes co-occupy active transcription units. We establish their interaction with the regions of genome having propensity to accumulate negative and positive supercoiled domains, validating their role in managing the twin supercoiled domains.

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.

Integration of Shh and Fgf signaling in controlling Hox gene expression in cultured limb cells

During embryonic development, fields of progenitor cells form complex structures through dynamic interactions with external signaling molecules. How complex signaling inputs are integrated to yield appropriate gene expression responses is poorly understood. In the early limb bud, for instance, Sonic hedgehog (Shh) is expressed in the distal posterior mesenchyme, where it acts as a mediator of anterior to posterior (AP) patterning, whereas fibroblast growth factor 8 (Fgf8) is produced by the apical ectodermal ridge (AER) at the distal tip of the limb bud to direct outgrowth along the proximal to distal (PD) axis. Here we use cultured limb mesenchyme cells to assess the response of the target Hoxd genes to these two factors. We find that they act synergistically and that both factors are required to activate Hoxd13 in limb mesenchymal cells. However, the analysis of the enhancer landscapes flanking the HoxD cluster reveals that the bimodal regulatory switch observed in vivo is only partially achieved under these in vitro conditions, suggesting an additional requirement for other factors.

The Inosine Monophosphate Dehydrogenase, GuaB2, Is a Vulnerable New Bactericidal Drug Target for Tuberculosis

VCC234718, a molecule with growth inhibitory activity against Mycobacterium tuberculosis (Mtb), was identified by phenotypic screening of a 15344-compound library. Sequencing of a VCC234718-resistant mutant identified a Y487C substitution in the inosine monophosphate dehydrogenase, GuaB2, which was subsequently validated to be the primary molecular target of VCC234718 in Mtb. VCC234718 inhibits Mtb GuaB2 with a Ki of 100 nM and is uncompetitive with respect to IMP and NAD+. This compound binds at the NAD+ site, after IMP has bound, and makes direct interactions with IMP; therefore, the inhibitor is by definition uncompetitive. VCC234718 forms strong pi interactions with the Y487 residue side chain from the adjacent protomer in the tetramer, explaining the resistance-conferring mutation. In addition to sensitizing Mtb to VCC234718, depletion of GuaB2 was bactericidal in Mtb in vitro and in macrophages. When supplied at a high concentration (≥125 μM), guanine alleviated the toxicity of VCC234718 treatment or GuaB2 depletion via purine salvage. However, transcriptional silencing of guaB2 prevented Mtb from establishing an infection in mice, confirming that Mtb has limited access to guanine in this animal model. Together, these data provide compelling validation of GuaB2 as a new tuberculosis drug target.

A role for mitotic bookmarking of SOX2 in pluripotency and differentiation

Mitotic bookmarking transcription factors remain bound to chromosomes during mitosis and were proposed to regulate phenotypic maintenance of stem and progenitor cells at the mitosis-to-G1 (M-G1) transition. However, mitotic bookmarking remains largely unexplored in most stem cell types, and its functional relevance for cell fate decisions remains unclear. Here we screened for mitotic chromosome binding within the pluripotency network of embryonic stem (ES) cells and show that SOX2 and OCT4 remain bound to mitotic chromatin through their respective DNA-binding domains. Dynamic characterization using photobleaching-based methods and single-molecule imaging revealed quantitatively similar specific DNA interactions, but different nonspecific DNA interactions, of SOX2 and OCT4 with mitotic chromatin. Using ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) to assess the genome-wide distribution of SOX2 on mitotic chromatin, we demonstrate the bookmarking activity of SOX2 on a small set of genes. Finally, we investigated the function of SOX2 mitotic bookmarking in cell fate decisions and show that its absence at the M-G1 transition impairs pluripotency maintenance and abrogates its ability to induce neuroectodermal differentiation but does not affect reprogramming efficiency toward induced pluripotent stem cells. Our study demonstrates the mitotic bookmarking property of SOX2 and reveals its functional importance in pluripotency maintenance and ES cell differentiation.

The ChIP-Seq tools and web server: a resource for analyzing ChIP-seq and other types of genomic data

BACKGROUND

ChIP-seq and related high-throughput chromatin profilig assays generate ever increasing volumes of highly valuable biological data. To make sense out of it, biologists need versatile, efficient and user-friendly tools for access, visualization and itegrative analysis of such data.

RESULTS

Here we present the ChIP-Seq command line tools and web server, implementing basic algorithms for ChIP-seq data analysis starting with a read alignment file. The tools are optimized for memory-efficiency and speed thus allowing for processing of large data volumes on inexpensive hardware. The web interface provides access to a large database of public data. The ChIP-Seq tools have a modular and interoperable design in that the output from one application can serve as input to another one. Complex and innovative tasks can thus be achieved by running several tools in a cascade.

CONCLUSIONS

The various ChIP-Seq command line tools and web services either complement or compare favorably to related bioinformatics resources in terms of computational efficiency, ease of access to public data and interoperability with other web-based tools. The ChIP-Seq server is accessible at http://ccg.vital-it.ch/chipseq/ .

Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics

Background

Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors.

Methods

We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes.

Results

Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1^–/– mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1^–/– mouse model.

Conclusions

These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0359-z) contains supplementary material, which is available to authorized users.