A deep learning method that identifies cellular heterogeneity using nanoscale nuclear features.

Cellular phenotypic heterogeneity is an important hallmark of many biological processes and understanding its origins remains a substantial challenge. This heterogeneity often reflects variations in the chromatin structure, influenced by factors such as viral infections and cancer, which dramatically reshape the cellular landscape. To address the challenge of identifying distinct cell states, we developed artificial intelligence of the nucleus (AINU), a deep learning method that can identify specific nuclear signatures at the nanoscale resolution. AINU can distinguish different cell states based on the spatial arrangement of core histone H3, RNA polymerase II or DNA from super-resolution microscopy images. With only a small number of images as the training data, AINU correctly identifies human somatic cells, human-induced pluripotent stem cells, very early stage infected cells transduced with DNA herpes simplex virus type 1 and even cancer cells after appropriate retraining. Finally, using AI interpretability methods, we find that the RNA polymerase II localizations in the nucleoli aid in distinguishing human-induced pluripotent stem cells from their somatic cells. Overall, AINU coupled with super-resolution microscopy of nuclear structures provides a robust tool for the precise detection of cellular heterogeneity, with considerable potential for advancing diagnostics and therapies in regenerative medicine, virology and cancer biology.

The relationship between nanoscale genome organization and gene expression in mouse embryonic stem cells during pluripotency transition.

During early development, gene expression is tightly regulated. However, how genome organization controls gene expression during the transition from naïve embryonic stem cells to epiblast stem cells is still poorly understood. Using single-molecule microscopy approaches to reach nanoscale resolution, we show that genome remodeling affects gene transcription during pluripotency transition. Specifically, after exit from the naïve pluripotency state, chromatin becomes less compacted, and the OCT4 transcription factor has lower mobility and is more bound to its cognate sites. In epiblast cells, the active transcription hallmark, H3K9ac, decreases within the Oct4 locus, correlating with reduced accessibility of OCT4 and, in turn, with reduced expression of Oct4 nascent RNAs. Despite the high variability in the distances between active pluripotency genes, distances between Nodal and Oct4 decrease during epiblast specification. In particular, highly expressed Oct4 alleles are closer to nuclear speckles during all stages of the pluripotency transition, while only a distinct group of highly expressed Nodal alleles are in close proximity to Oct4 when associated with a nuclear speckle in epiblast cells. Overall, our results provide new insights into the role of the spatiotemporal genome remodeling during mouse pluripotency transition and its correlation with the expression of key pluripotency genes.

The Wnt-dependent master regulator NKX1-2 controls mouse pre-implantation development.

Embryo size, specification, and homeostasis are regulated by a complex gene regulatory and signaling network. Here we used gene expression signatures of Wnt-activated mouse embryonic stem cell (mESC) clones to reverse engineer an mESC regulatory network. We identify NKX1-2 as a novel master regulator of preimplantation embryo development. We find that Nkx1-2 inhibition reduces nascent RNA synthesis, downregulates genes controlling ribosome biogenesis, RNA translation, and transport, and induces severe alteration of nucleolus structure, resulting in the exclusion of RNA polymerase I from nucleoli. In turn, NKX1-2 loss of function leads to chromosome missegregation in the 2- to 4-cell embryo stages, severe decrease in blastomere numbers, alterations of tight junctions (TJs), and impairment of microlumen coarsening. Overall, these changes impair the blastocoel expansion-collapse cycle and embryo cavitation, leading to altered lineage specification and developmental arrest.

Active transcription and epigenetic reactions synergistically regulate meso-scale genomic organization.

In interphase nuclei, chromatin forms dense domains of characteristic sizes, but the influence of transcription and histone modifications on domain size is not understood. We present a theoretical model exploring this relationship, considering chromatin-chromatin interactions, histone modifications, and chromatin extrusion. We predict that the size of heterochromatic domains is governed by a balance among the diffusive flux of methylated histones sustaining them and the acetylation reactions in the domains and the process of loop extrusion via supercoiling by RNAPII at their periphery, which contributes to size reduction. Super-resolution and nano-imaging of five distinct cell lines confirm the predictions indicating that the absence of transcription leads to larger heterochromatin domains. Furthermore, the model accurately reproduces the findings regarding how transcription-mediated supercoiling loss can mitigate the impacts of excessive cohesin loading. Our findings shed light on the role of transcription in genome organization, offering insights into chromatin dynamics and potential therapeutic targets.

LATS1 controls CTCF chromatin occupancy and hormonal response of 3D-grown breast cancer cells.

The cancer epigenome has been studied in cells cultured in two-dimensional (2D) monolayers, but recent studies highlight the impact of the extracellular matrix and the three-dimensional (3D) environment on multiple cellular functions. Here, we report the physical, biochemical, and genomic differences between T47D breast cancer cells cultured in 2D and as 3D spheroids. Cells within 3D spheroids exhibit a rounder nucleus with less accessible, more compacted chromatin, as well as altered expression of ~2000 genes, the majority of which become repressed. Hi-C analysis reveals that cells in 3D are enriched for regions belonging to the B compartment, have decreased chromatin-bound CTCF and increased fusion of topologically associating domains (TADs). Upregulation of the Hippo pathway in 3D spheroids results in the activation of the LATS1 kinase, which promotes phosphorylation and displacement of CTCF from DNA, thereby likely causing the observed TAD fusions. 3D cells show higher chromatin binding of progesterone receptor (PR), leading to an increase in the number of hormone-regulated genes. This effect is in part mediated by LATS1 activation, which favors cytoplasmic retention of YAP and CTCF removal.

CEBPA phase separation links transcriptional activity and 3D chromatin hubs.

Cell identity is orchestrated through an interplay between transcription factor (TF) action and genome architecture. The mechanisms used by TFs to shape three-dimensional (3D) genome organization remain incompletely understood. Here we present evidence that the lineage-instructive TF CEBPA drives extensive chromatin compartment switching and promotes the formation of long-range chromatin hubs during induced B cell-to-macrophage transdifferentiation. Mechanistically, we find that the intrinsically disordered region (IDR) of CEBPA undergoes in vitro phase separation (PS) dependent on aromatic residues. Both overexpressing B cells and native CEBPA-expressing cell types such as primary granulocyte-macrophage progenitors, liver cells, and trophectoderm cells reveal nuclear CEBPA foci and long-range 3D chromatin hubs at CEBPA-bound regions. In short, we show that CEBPA can undergo PS through its IDR, which may underlie in vivo foci formation and suggest a potential role of PS in regulating CEBPA function.

STORM Microscopy and Cluster Analysis for PcG Studies.

Advanced microscopy techniques (such as STORM, STED, and SIM) have recently allowed the visualization of biological samples beyond the diffraction limit of light. Thanks to this breakthrough, the organization of molecules can be revealed within single cells as never before.Here, we describe the application of STochastic Optical Reconstruction Microscopy (STORM) for the study of polycomb group of proteins (PcG) in the context of chromatin organization. We present a clustering algorithm to quantitatively analyze the spatial distribution of nuclear molecules (e.g., EZH2 or its associated chromatin mark H3K27me3) imaged by 2D STORM. This distance-based analysis uses x-y coordinates of STORM localizations to group them into "clusters." Clusters are classified as singles if isolated or into islands if they form a group of closely associated clusters. For each cluster, the algorithm calculates the number of localizations, the area, and the distance to the closest cluster.This approach can be used for every type of adherent cell line and allows the imaging of every protein for which an antibody is available. It represents a comprehensive strategy to visualize and quantify how PcG proteins and related histone marks organize in the nucleus at nanometric resolution.

The magic of unraveling genome architecture and function.

Over the last decades, technological breakthroughs in super-resolution microscopy have allowed us to reach molecular resolution and design experiments of unprecedented complexity. Investigating how chromatin is folded in 3D, from the nucleosome level up to the entire genome, is becoming possible by "magic" (imaging genomic), i.e., the combination of imaging and genomic approaches. This offers endless opportunities to delve into the relationship between genome structure and function. Here, we review recently achieved objectives and the conceptual and technical challenges the field of genome architecture is currently undertaking. We discuss what we have learned so far and where we are heading. We elucidate how the different super-resolution microscopy approaches and, more specifically, live-cell imaging have contributed to the understanding of genome folding. Moreover, we discuss how future technical developments could address remaining open questions.

Cohesin-independent STAG proteins interact with RNA and R-loops and promote complex loading.

Most studies of cohesin function consider the Stromalin Antigen (STAG/SA) proteins as core complex members given their ubiquitous interaction with the cohesin ring. Here, we provide functional data to support the notion that the SA subunit is not a mere passenger in this structure, but instead plays a key role in the localization of cohesin to diverse biological processes and promotes loading of the complex at these sites. We show that in cells acutely depleted for RAD21, SA proteins remain bound to chromatin, cluster in 3D and interact with CTCF, as well as with a wide range of RNA binding proteins involved in multiple RNA processing mechanisms. Accordingly, SA proteins interact with RNA, and R-loops, even in the absence of cohesin. Our results place SA1 on chromatin upstream of the cohesin ring and reveal a role for SA1 in cohesin loading which is independent of NIPBL, the canonical cohesin loader. We propose that SA1 takes advantage of structural R-loop platforms to link cohesin loading and chromatin structure with diverse functions. Since SA proteins are pan-cancer targets, and R-loops play an increasingly prevalent role in cancer biology, our results have important implications for the mechanistic understanding of SA proteins in cancer and disease.

Women's contribution in understanding how topoisomerases, supercoiling, and transcription control genome organization.

One of the biggest paradoxes in biology is that human genome is roughly 2 m long, while the nucleus containing it is almost one million times smaller. To fit into the nucleus, DNA twists, bends and folds into several hierarchical levels of compaction. Still, DNA has to maintain a high degree of accessibility to be readily replicated and transcribed by proteins. How compaction and accessibility co-exist functionally in human cells is still a matter of debate. Here, we discuss how the torsional stress of the DNA helix acts as a buffer, regulating both chromatin compaction and accessibility. We will focus on chromatin supercoiling and on the emerging role of topoisomerases as pivotal regulators of genome organization. We will mainly highlight the major breakthrough studies led by women, with the intention of celebrating the work of this group that remains a minority within the scientific community.

Cohesin controls X chromosome structure remodeling and X-reactivation during mouse iPSC-reprogramming.

Reactivation of the inactive X chromosome is a hallmark epigenetic event during reprogramming of mouse female somatic cells to induced pluripotent stem cells (iPSCs). This involves global structural remodeling from a condensed, heterochromatic into an open, euchromatic state, thereby changing a transcriptionally inactive into an active chromosome. Despite recent advances, very little is currently known about the molecular players mediating this process and how this relates to iPSC-reprogramming in general. To gain more insight, here we perform a RNAi-based knockdown screen during iPSC-reprogramming of mouse fibroblasts. We discover factors important for X chromosome reactivation (XCR) and iPSC-reprogramming. Among those, we identify the cohesin complex member SMC1a as a key molecule with a specific function in XCR, as its knockdown greatly affects XCR without interfering with iPSC-reprogramming. Using super-resolution microscopy, we find SMC1a to be preferentially enriched on the active compared with the inactive X chromosome and that SMC1a is critical for the decompacted state of the active X. Specifically, depletion of SMC1a leads to contraction of the active X both in differentiated and in pluripotent cells, where it normally is in its most open state. In summary, we reveal cohesin as a key factor for remodeling of the X chromosome from an inactive to an active structure and that this is a critical step for XCR during iPSC-reprogramming.

MiOS, an integrated imaging and computational strategy to model gene folding with nucleosome resolution.

The linear sequence of DNA provides invaluable information about genes and their regulatory elements along chromosomes. However, to fully understand gene function and regulation, we need to dissect how genes physically fold in the three-dimensional nuclear space. Here we describe immuno-OligoSTORM, an imaging strategy that reveals the distribution of nucleosomes within specific genes in super-resolution, through the simultaneous visualization of DNA and histones. We combine immuno-OligoSTORM with restraint-based and coarse-grained modeling approaches to integrate super-resolution imaging data with Hi-C contact frequencies and deconvoluted micrococcal nuclease-sequencing information. The resulting method, called Modeling immuno-OligoSTORM, allows quantitative modeling of genes with nucleosome resolution and provides information about chromatin accessibility for regulatory factors, such as RNA polymerase II. With Modeling immuno-OligoSTORM, we explore intercellular variability, transcriptional-dependent gene conformation, and folding of housekeeping and pluripotency-related genes in human pluripotent and differentiated cells, thereby obtaining the highest degree of data integration achieved so far to our knowledge.

Super resolution microscopy reveals how elongating RNA polymerase II and nascent RNA interact with nucleosome clutches.

Transcription and genome architecture are interdependent, but it is still unclear how nucleosomes in the chromatin fiber interact with nascent RNA, and which is the relative nuclear distribution of these RNAs and elongating RNA polymerase II (RNAP II). Using super-resolution (SR) microscopy, we visualized the nascent transcriptome, in both nucleoplasm and nucleolus, with nanoscale resolution. We found that nascent RNAs organize in structures we termed RNA nanodomains, whose characteristics are independent of the number of transcripts produced over time. Dual-color SR imaging of nascent RNAs, together with elongating RNAP II and H2B, shows the physical relation between nucleosome clutches, RNAP II, and RNA nanodomains. The distance between nucleosome clutches and RNA nanodomains is larger than the distance measured between elongating RNAP II and RNA nanodomains. Elongating RNAP II stands between nascent RNAs and the small, transcriptionally active, nucleosome clutches. Moreover, RNA factories are small and largely formed by few RNAP II. Finally, we describe a novel approach to quantify the transcriptional activity at an individual gene locus. By measuring local nascent RNA accumulation upon transcriptional activation at single alleles, we confirm the measurements made at the global nuclear level.

A protocol to quantify chromatin compaction with confocal and super-resolution microscopy in cultured cells.

Here, we describe three complementary microscopy-based approaches to quantify morphological changes of chromatin organization in cultured adherent cells: the analysis of the coefficient of variation of DNA, the measurement of DNA-free nuclear areas, and the quantification of chromatin-associated proteins at the nuclear edge. These approaches rely on confocal imaging and stochastic optical reconstruction microscopy and allow a fast and robust quantification of chromatin compaction. These approaches circumvent inter-variability between imaging conditions and apply to every type of adherent cells. For complete details on the use and execution of this protocol, please refer to Neguembor et al. (2021).

Transcription-mediated supercoiling regulates genome folding and loop formation.

The chromatin fiber folds into loops, but the mechanisms controlling loop extrusion are still poorly understood. Using super-resolution microscopy, we visualize that loops in intact nuclei are formed by a scaffold of cohesin complexes from which the DNA protrudes. RNA polymerase II decorates the top of the loops and is physically segregated from cohesin. Augmented looping upon increased loading of cohesin on chromosomes causes disruption of Lamin at the nuclear rim and chromatin blending, a homogeneous distribution of chromatin within the nucleus. Altering supercoiling via either transcription or topoisomerase inhibition counteracts chromatin blending, increases chromatin condensation, disrupts loop formation, and leads to altered cohesin distribution and mobility on chromatin. Overall, negative supercoiling generated by transcription is an important regulator of loop formation in vivo.

Mesoscale Modeling and Single-Nucleosome Tracking Reveal Remodeling of Clutch Folding and Dynamics in Stem Cell Differentiation.

Nucleosomes form heterogeneous groups in vivo, named clutches. Clutches are smaller and less dense in mouse embryonic stem cells (ESCs) compared to neural progenitor cells (NPCs). Using coarse-grained modeling of the pluripotency Pou5f1 gene, we show that the genome-wide clutch differences between ESCs and NPCs can be reproduced at a single gene locus. Larger clutch formation in NPCs is associated with changes in the compaction and internucleosome contact probability of the Pou5f1 fiber. Using single-molecule tracking (SMT), we further show that the core histone protein H2B is dynamic, and its local mobility relates to the structural features of the chromatin fiber. H2B is less stable and explores larger areas in ESCs compared to NPCs. The amount of linker histone H1 critically affects local H2B dynamics. Our results have important implications for how nucleosome organization and H2B dynamics contribute to regulate gene activity and cell identity.

From research to rapid response: mass COVID-19 testing by volunteers at the Centre for Genomic Regulation.

The COVID-19 pandemic has posed and is continuously posing enormous societal and health challenges worldwide. The research community has mobilized to develop novel projects to find a cure or a vaccine, as well as to contribute to mass testing, which has been a critical measure to contain the infection in several countries. Through this article, we share our experiences and learnings as a group of volunteers at the Centre for Genomic Regulation (CRG) in Barcelona, Spain. As members of the ORFEU project, an initiative by the Government of Catalonia to achieve mass testing of people at risk and contain the epidemic in Spain, we share our motivations, challenges and the key lessons learnt, which we feel will help better prepare the global society to address similar situations in the future.

The Suv420h histone methyltransferases regulate PPAR-γ and energy expenditure in response to environmental stimuli.

Obesity and its associated metabolic abnormalities have become a global emergency with considerable morbidity and mortality. Epidemiologic and animal model data suggest an epigenetic contribution to obesity. Nevertheless, the cellular and molecular mechanisms through which epigenetics contributes to the development of obesity remain to be elucidated. Suv420h1 and Suv420h2 are histone methyltransferases responsible for chromatin compaction and gene repression. Through in vivo, ex vivo, and in vitro studies, we found that Suv420h1 and Suv420h2 respond to environmental stimuli and regulate metabolism by down-regulating peroxisome proliferator-activated receptor gamma (PPAR-γ), a master transcriptional regulator of lipid storage and glucose metabolism. Accordingly, mice lacking Suv420h proteins activate PPAR-γ target genes in brown adipose tissue to increase mitochondria respiration, improve glucose tolerance, and reduce adipose tissue to fight obesity. We conclude that Suv420h proteins are key epigenetic regulators of PPAR-γ and the pathways controlling metabolism and weight balance in response to environmental stimuli.

Wnt/ß-catenin signaling pathway safeguards epigenetic stability and homeostasis of mouse embryonic stem cells.

Mouse embryonic stem cells (mESCs) are pluripotent and can differentiate into cells belonging to the three germ layers of the embryo. However, mESC pluripotency and genome stability can be compromised in prolonged in vitro culture conditions. Several factors control mESC pluripotency, including Wnt/ß-catenin signaling pathway, which is essential for mESC differentiation and proliferation. Here we show that the activity of the Wnt/ß-catenin signaling pathway safeguards normal DNA methylation of mESCs. The activity of the pathway is progressively silenced during passages in culture and this results into a loss of the DNA methylation at many imprinting control regions (ICRs), loss of recruitment of chromatin repressors, and activation of retrotransposons, resulting into impaired mESC differentiation. Accordingly, sustained Wnt/ß-catenin signaling maintains normal ICR methylation and mESC homeostasis and is a key regulator of genome stability.

A cytosolic Ezh1 isoform modulates a PRC2-Ezh1 epigenetic adaptive response in postmitotic cells.

The evolution of chromatin-based epigenetic cell memory may be driven not only by the necessity for cells to stably maintain transcription programs, but also by the need to recognize signals and allow plastic responses to environmental stimuli. The mechanistic role of the epigenome in adult postmitotic tissues, however, remains largely unknown. In vertebrates, two variants of the Polycomb repressive complex (PRC2-Ezh2 and PRC2-Ezh1) control gene silencing via methylation of histone H3 on Lys27 (H3K27me). Here we describe a reversible mechanism that involves a novel isoform of Ezh1 (Ezh1ß). Ezh1ß lacks the catalytic SET domain and acts in the cytoplasm of skeletal muscle cells to control nuclear PRC2-Ezh1 activity in response to atrophic oxidative stress, by regulating Eed assembly with Suz12 and Ezh1α (the canonical isoform) at their target genes. We report a novel PRC2-Ezh1 function that utilizes Ezh1ß as an adaptive stress sensor in the cytoplasm, thus allowing postmitotic cells to maintain tissue integrity in response to environmental changes.