Global genome decompaction leads to stochastic activation of gene expression as a first step toward fate commitment in human hematopoietic cells.

When human cord blood-derived CD34+ cells are induced to differentiate, they undergo rapid and dynamic morphological and molecular transformations that are critical for fate commitment. In particular, the cells pass through a transitory phase known as "multilineage-primed" state. These cells are characterized by a mixed gene expression profile, different in each cell, with the coexpression of many genes characteristic for concurrent cell lineages. The aim of our study is to understand the mechanisms of the establishment and the exit from this transitory state. We investigated this issue using single-cell RNA sequencing and ATAC-seq. Two phases were detected. The first phase is a rapid and global chromatin decompaction that makes most of the gene promoters in the genome accessible for transcription. It results 24 h later in enhanced and pervasive transcription of the genome leading to the concomitant increase in the cell-to-cell variability of transcriptional profiles. The second phase is the exit from the multilineage-primed phase marked by a slow chromatin closure and a subsequent overall down-regulation of gene transcription. This process is selective and results in the emergence of coherent expression profiles corresponding to distinct cell subpopulations. The typical time scale of these events spans 48 to 72 h. These observations suggest that the nonspecificity of genome decompaction is the condition for the generation of a highly variable multilineage expression profile. The nonspecific phase is followed by specific regulatory actions that stabilize and maintain the activity of key genes, while the rest of the genome becomes repressed again by the chromatin recompaction. Thus, the initiation of differentiation is reminiscent of a constrained optimization process that associates the spontaneous generation of gene expression diversity to subsequent regulatory actions that maintain the activity of some genes, while the rest of the genome sinks back to the repressive closed chromatin state.

Transcription regulates the spatio-temporal dynamics of genes through micro-compartmentalization.

Although our understanding of the involvement of heterochromatin architectural factors in shaping nuclear organization is improving, there is still ongoing debate regarding the role of active genes in this process. In this study, we utilize publicly-available Micro-C data from mouse embryonic stem cells to investigate the relationship between gene transcription and 3D gene folding. Our analysis uncovers a nonmonotonic - globally positive - correlation between intragenic contact density and Pol II occupancy, independent of cohesin-based loop extrusion. Through the development of a biophysical model integrating the role of transcription dynamics within a polymer model of chromosome organization, we demonstrate that Pol II-mediated attractive interactions with limited valency between transcribed regions yield quantitative predictions consistent with chromosome-conformation-capture and live-imaging experiments. Our work provides compelling evidence that transcriptional activity shapes the 4D genome through Pol II-mediated micro-compartmentalization.

The human NUPR1/P8 gene is transcriptionally activated by transforming growth factor ß via the SMAD signalling pathway.

NUPR1 (nuclear protein 1), also called P8 (molecular mass 8 kDa) or COM1 (candidate of metastasis 1), is involved in the stress response and in cancer progression. In the present study, we investigated whether human NUPR1 expression was regulated by TGFß (transforming growth factor ß), a secreted polypeptide largely involved in tumorigenesis. We demonstrate that the expression of NUPR1 was activated by TGFß at the transcriptional level. We show that this activation is mediated by the SMAD proteins, which are transcription factors specifically involved in the signalling of TGFß superfamily members. NUPR1 promoter analysis reveals the presence of a functional TGFß-response element binding the SMAD proteins located in the genomic DNA region corresponding to the 5'-UTR (5'-untranslated region). Altogether, the molecular results of the present study, which demonstrate the existence of a TGFß/SMAD/NUPR1 activation cascade, open the way to consider and investigate further a new mechanism enabling TGFß to promote tumorigenesis by inducing stress resistance.

CLLD8/KMT1F is a lysine methyltransferase that is important for chromosome segregation.

Proteins bearing a SET domain have been shown to methylate lysine residues in histones and contribute to chromatin architecture. Methylation of histone H3 at lysine 9 (H3K9) has emerged as an important player in the formation of heterochromatin, chromatin condensation, and transcriptional repression. Here, we have characterized a previously undescribed member of the histone H3K9 methyltransferase family named CLLD8 (or SETDB2 or KMT1F). This protein contributes to the trimethylation of both interspersed repetitive elements and centromere-associated repeats and participates in the recruitment of heterochromatin protein 1 to centromeres. Consistently, depletion in CLLD8/KMT1F coincides with a loss of CENP proteins and delayed mitosis, suggesting that this protein participates in chromosome condensation and segregation. Altogether, our results provide evidence that CLLD8/KMT1F is recruited to heterochromatin regions and contributes in vivo to the deposition of trimethyl marks in concert with SUV39H1/KMT1A.

A transcriptomic analysis of human centromeric and pericentric sequences in normal and tumor cells.

Although there is now evidence that the expression of centromeric (CT) and pericentric (PCT) sequences are key players in major genomic functions, their transcriptional status in human cells is still poorly known. The main reason for this lack of data is the complexity and high level of polymorphism of these repeated sequences, which hampers straightforward analyses by available transcriptomic approaches. Here a transcriptomic macro-array dedicated to the analysis of CT and PCT expression is developed and validated in heat-shocked (HS) HeLa cells. For the first time, the expression status of CT and PCT sequences is analyzed in a series of normal and cancer human cells and tissues demonstrating that they are repressed in all normal tissues except in the testis, where PCT transcripts are found. Moreover, PCT sequences are specifically expressed in HS cells in a Heat-Shock Factor 1 (HSF1)-dependent fashion, and we show here that another independent pathway, involving DNA hypo-methylation, can also trigger their expression. Interestingly, CT and PCT were found illegitimately expressed in somatic cancer samples, whereas PCT were repressed in testis cancer, suggesting that the expression of CT and PCT sequences may represent a good indicator of epigenetic deregulations occurring in response to environmental changes or in cell transformation.

Tubulin mutations in neurodevelopmental disorders as a tool to decipher microtubule function.

Malformations of cortical development (MCDs) are a group of severe brain malformations associated with intellectual disability and refractory childhood epilepsy. Human missense heterozygous mutations in the 9 α-tubulin and 10 ß-tubulin isoforms forming the heterodimers that assemble into microtubules (MTs) were found to cause MCDs. However, how a single mutated residue in a given tubulin isoform can perturb the entire microtubule population in a neuronal cell remains a crucial question. Here, we examined 85 MCD-associated tubulin mutations occurring in TUBA1A, TUBB2, and TUBB3 and their location in a three-dimensional (3D) microtubule cylinder. Mutations hitting residues exposed on the outer microtubule surface are likely to alter microtubule association with partners, while alteration of intradimer contacts may impair dimer stability and straightness. Other types of mutations are predicted to alter interdimer and lateral contacts, which are responsible for microtubule cohesion, rigidity, and dynamics. MCD-associated tubulin mutations surprisingly fall into all categories, thus providing unexpected insights into how a single mutation may impair microtubule function and elicit dominant effects in neurons.

Dynamics as a cause for the nanoscale organization of the genome.

Chromatin 'blobs' were recently identified by live super-resolution imaging of labeled nucleosomes as pervasive but fleeting structural entities. However, the mechanisms leading to the formation of these blobs and their functional implications are unknown. We explore here whether causal relationships exist between parameters that characterize the chromatin blob dynamics and structure, by adapting a framework for spatio-temporal Granger-causality inference. Our analysis reveals that chromatin dynamics is a key determinant for both blob area and local density. Such causality, however, could be demonstrated only in 10-20% of the nucleus, suggesting that chromatin dynamics and structure at the nanometer scale are dominated by stochasticity. We show that the theory of active semiflexible polymers can be invoked to provide potential mechanisms leading to the organization of chromatin into blobs. Our results represent a first step toward elucidating the mechanisms that govern the dynamic and stochastic organization of chromatin in the cell nucleus.

Differential requirement of DNA replication factors for subtelomeric ARS consensus sequence protosilencers in Saccharomyces cerevisiae.

The establishment of silent chromatin requires passage through S-phase, but not DNA replication per se. Nevertheless, many proteins that affect silencing are bona fide DNA replication factors. It is not clear if mutations in these replication factors affect silencing directly or indirectly via deregulation of S-phase or DNA replication. Consequently, the relationship between DNA replication and silencing remains an issue of debate. Here we analyze the effect of mutations in DNA replication factors (mcm5-461, mcm5-1, orc2-1, orc5-1, cdc45-1, cdc6-1, and cdc7-1) on the silencing of a group of reporter constructs, which contain different combinations of "natural" subtelomeric elements. We show that the mcm5-461, mcm5-1, and orc2-1 mutations affect silencing through subtelomeric ARS consensus sequences (ACS), while cdc6-1 affects silencing independently of ACS. orc5-1, cdc45-1, and cdc7-1 affect silencing through ACS, but also show ACS-independent effects. We also demonstrate that isolated nontelomeric ACS do not recapitulate the same effects when inserted in the telomere. We propose a model that defines the modes of action of MCM5 and CDC6 in silencing.

A methyltransferase targeting assay reveals silencer-telomere interactions in budding yeast.

We have designed a modified version of the Dam identification technique and used it to probe higher-order chromatin structure in Saccharomyces cerevisiae. We fused the bacterial DNA methyltransferase Dam to the DNA-binding domain of TetR and targeted the resulting chimera to Tet operators inserted in the yeast genome at the repressed locus HML. We then monitored the methylation status of HML and other sequences by a quantitative technique combining methylation-sensitive restriction and real-time PCR. As expected, we found that TetR-Dam efficiently methylated HML in cis. More strikingly, when TetR-Dam was present at HML, we observed increased methylation in the III-L subtelomeric region but not in intervening sequences. This effect was lost when the HML silencers were inactivated by mutations. When the HM silencers and the Tet operators were transferred to a plasmid, strong methylation was clearly observed not only in the III-L subtelomeric region but also at other telomeres. These data indicate that HM silencers can specifically associate with telomeres, even those located on different chromosomes.

TIF1γ Suppresses Tumor Progression by Regulating Mitotic Checkpoints and Chromosomal Stability.

The transcription accessory factor TIF1γ/TRIM33/RFG7/PTC7/Ectodermin functions as a tumor suppressor that promotes development and cellular differentiation. However, its precise function in cancer has been elusive. In the present study, we report that TIF1γ inactivation causes cells to accumulate chromosomal defects, a hallmark of cancer, due to attenuations in the spindle assembly checkpoint and the post-mitotic checkpoint. TIF1γ deficiency also caused a loss of contact growth inhibition and increased anchorage-independent growth in vitro and in vivo. Clinically, reduced TIF1γ expression in human tumors correlated with an increased rate of genomic rearrangements. Overall, our work indicates that TIF1γ exerts its tumor-suppressive functions in part by promoting chromosomal stability.

Subtelomeric ACS-containing proto-silencers act as antisilencers in replication factors mutants in Saccharomyces cerevisiae.

Subtelomeric genes are either fully active or completely repressed and can switch their state about once per 20 generations. This meta-stable telomeric position effect is mediated by strong repression signals emitted by the telomere and relayed/enhanced by weaker repressor elements called proto-silencers. In addition, subtelomeric regions contain sequences with chromatin partitioning and antisilencing activities referred to as subtelomeric antisilencing regions. Using extensive mutational analysis of subtelomeric elements, we show that ARS consensus sequence (ACS)-containing proto-silencers convert to antisilencers in several replication factor mutants. We point out the significance of the B1 auxiliary sequence next to ACS in mediating these effects. In contrast, an origin-derived ACS does not convert to antisilencer in mutants and its B1 element has little bearing on silencing. These results are specific for the analyzed ACS and in addition to the effects of each mutation (relative to wild type) on global silencing. Another line of experiments shows that Mcm5p possesses antisilencing activity and is recruited to telomeres in an ACS-dependent manner. Mcm5p persists at this location at the late stages of S phase. We propose that telomeric ACS are not static proto-silencers but conduct finely tuned silencing and antisilencing activities mediated by ACS-bound factors.

Global analysis of DNA methylation and transcription of human repetitive sequences.

Half of the human genome consists of repetitive DNA sequences. Recent studies in various organisms highlight the role of chromatin regulation of repetitive DNA in gene regulation as well as in maintainance of chromosomes and genome integrity. Hence, repetitive DNA sequences might be potential "sensors" for chromatin changes associated with pathogenesis. Therefore, we developed a new genomic tool called RepArray. RepArray is a repeat-specific microarray composed of a representative set of human repeated sequences including transposon-derived repeats, simple sequences repeats, tandemly repeated sequences such as centromeres and telomeres. We showed that combined to anti-methylcytosine immunoprecipitation assay, the RepArray can be used to generate repeat-specific methylation maps. Using cell lines impaired chemically or genetically for DNA methyltransferases activities, we were able to distinguish different epigenomes demonstrating that repeats can be used as markers of genome-wide methylation changes. Besides, using a well-documented system model, the thermal stress, we demonstrated that RepArray is also a fast and reliable tool to obtain an overview of overall transcriptional activity on whole repetitive compartment in a given cell type. Thus, the RepArray represents the first valuable tool for systematic and genome-wide analyses of the methylation and transcriptional status of the repetitive counterpart of the human genome.

Protosilencers as building blocks for heterochromatin.

DNA repetitions may provoke heterochromatinization. We explore here a model in which multiple cis-acting sequences that display no silencing activity on their own (protosilencers) may cooperate to establish and maintain a heterochromatin domain efficiently. Protosilencers, first defined in budding yeast, have now been found in a wide range of genomes where they appear to stabilize and to extend the propagation of heterochromatin domains. Strikingly, isolated or moderately repeated protosilencers can also be found in promoters where they participate in transcriptional activation and have insulation functions. This suggests that the proper juxtaposition of a threshold number of protosilencers converts them from neutral or transactivating elements into ones that nucleate heterochromatin. Interactions might be transient or permanent, and are likely to occur over distances by looping. This model provides a conceptual framework for as varied phenomena as telomere-driven silencing in Drosophila, X inactivation in mammals, and rDNA silencing in S. cerevisiae. It may also account for the silencing that occurs when multiple copies of a transgene are inserted in tandem.

Insulator dynamics and the setting of chromatin domains.

The early discovery of cis-regulatory elements able to promote transcription of genes over large distances led to the postulate that elements, termed insulators, should also exist that would limit the action of enhancers, LCRs and silencers to defined domains. Such insulators were indeed found during the past fifteen years in a wide range of organisms, from yeast to humans. Recent advances point to an important role of transcription factors in insulator activity and demonstrate that the operational observation of an insulator effect relies on a delicate balance between the "efficiency" of the insulator and that of the element to be counteracted. In addition, genuine insulator elements now appear less common than initially envisaged, and they are only found at loci displaying a high density of coding or regulatory information. Where this is not the case, chromatin domains of opposing properties are thought to confront each other at "fuzzy" boundaries. In this article, we propose models for both fixed and fuzzy boundaries that incorporate probabilistic and dynamic parameters.

General regulatory factors (GRFs) as genome partitioners.

Insulators are sequences that uncouple adjacent chromosome domains. Here we have shown that Saccharomyces cerevisiae Rap1p and Abf1p proteins are endowed with a potent insulating capacity. Insulating domains in Rap1p coincide with previously described transcription activation domains, whereas four adjacent subdomains spanning the whole of the Abf1p C terminus (440-731) were found to display autonomous insulating capacity. That both Rap1p and Abf1p silencing domains either contain or largely overlap with an insulating domain suggests that insulation conveys some undefined chromosome organization capacity that also contributes a function in silencing. Together with Reb1p and Tbf1p, previously involved in the activity of Saccharomyces cerevisiae subtelomeric insulators, insulating potential emerges as a supplementary common property of General Regulatory Factors (GRFs). Thus GRFs, which bind to sites scattered throughout the genome within promoters, would not only play a key role in regulating gene expression but also partition the genome in functionally independent domains.

The win locus involved in activation of the distal N-myc2 gene upon WHV integration in woodchuck liver tumors harbors S/MAR elements.

Woodchuck hepatitis virus (WHV) and the woodchuck (Marmota monax) are models for hepatocellular carcinoma (HCC) induced by hepatitis B virus (HBV). In woodchuck liver tumors, the N-myc2 proto-oncogene is frequently activated by WHV integration either close to the gene or in the b3n and win downstream loci, located 10 and 150 kb from N-myc2, respectively. A scaffold/matrix attachment region (S/MAR) regulative element was shown to be in b3n, possibly mediating activation of the upstream N-myc2 gene upon WHV integration. To investigate if S/MAR elements are in win too, a 17-kb DNA fragment corresponding to the major region of WHV insertion in this locus was cloned and sequenced. Overlapping subcloned fragments spanning candidate S/MARs predicted by sequence analysis were tested by standard in vitro binding assays. Results showed the presence of two S/MAR elements in win. The distribution of previously described WHV insertions relative to the S/MARs reinforces the hypothesis that S/MARs nearby distal WHV insertions might be involved in long-range activation of N-myc2.

Chromatin domain boundaries delimited by a histone-binding protein in yeast.

When located next to chromosomal elements such as telomeres, genes can be subjected to epigenetic silencing. In yeast, this is mediated by the propagation of the SIR proteins from telomeres toward more centromeric regions. Particular transcription factors can protect downstream genes from silencing when tethered between the gene and the telomere, and they may thus act as chromatin domain boundaries. Here we have studied one such transcription factor, CTF-1, that binds directly histone H3. A deletion mutagenesis localized the barrier activity to the CTF-1 histone-binding domain. A saturating point mutagenesis of this domain identified several amino acid substitutions that similarly inhibited the boundary and histone binding activities. Chromatin immunoprecipitation experiments indicated that the barrier protein efficiently prevents the spreading of SIR proteins, and that it separates domains of hypoacetylated and hyperacetylated histones. Together, these results suggest a mechanism by which proteins such as CTF-1 may interact directly with histone H3 to prevent the propagation of a silent chromatin structure, thereby defining boundaries of permissive and silent chromatin domains.