The stochastic nature of genome organization and function

The chromosome folding problem and how cells solve it

Every cell must solve the problem of how to fold its genome. We describe how the folded state of chromosomes is the result of the combined activity of multiple conserved mechanisms. Homotypic affinity-driven interactions lead to spatial partitioning of active and inactive loci. Molecular motors fold chromosomes through loop extrusion. Topological features such as supercoiling and entanglements contribute to chromosome folding and its dynamics, and tethering loci to sub-nuclear structures adds additional constraints. Dramatically diverse chromosome conformations observed throughout the cell cycle and across the tree of life can be explained through differential regulation and implementation of these basic mechanisms. We propose that the first functions of chromosome folding are to mediate genome replication, compaction, and segregation and that mechanisms of folding have subsequently been co-opted for other roles, including long-range gene regulation, in different conditions, cell types, and species.

Evolutionary cues of gene fusion and fission in plants

KEY MESSAGE

Here, we systematically analyzed the potential fusion and fission events of neighboring genes in Arabidopsis genome and analyzed the influence on the protein targeting.

Local volume concentration, packing domains, and scaling properties of chromatin

We propose the Self Returning Excluded Volume (SR-EV) model for the structure of chromatin based on stochastic rules and physical interactions. The SR-EV rules of return generate conformationally defined domains observed by single-cell imaging techniques. From nucleosome to chromosome scales, the model captures the overall chromatin organization as a corrugated system, with dense and dilute regions alternating in a manner that resembles the mixing of two disordered bi-continuous phases. This particular organizational topology is a consequence of the multiplicity of interactions and processes occurring in the nuclei, and mimicked by the proposed return rules. Single configuration properties and ensemble averages show a robust agreement between theoretical and experimental results including chromatin volume concentration, contact probability, packing domain identification and size characterization, and packing scaling behavior. Model and experimental results suggest that there is an inherent chromatin organization regardless of the cell character and resistant to an external forcing such as RAD21 degradation.

Hi-C techniques: from genome assemblies to transcription regulation

The invention of chromosome conformation capture (3C) techniques, in particular the key method Hi-C providing genome-wide information about chromatin contacts, revolutionized the way we study the three-dimensional organization of the nuclear genome and how it affects transcription, replication, and DNA repair. Because the frequency of chromatin contacts between pairs of genomic segments predictably relates to the distance in the linear genome, the information obtained by Hi-C has also proved useful for scaffolding genomic sequences. Here, we review recent improvements in experimental procedures of Hi-C and its various derivatives, such as Micro-C, HiChIP, and Capture Hi-C. We assess the advantages and limitations of the techniques, and present examples of their use in recent plant studies. We also report on progress in the development of computational tools used in assembling genome sequences.

Loss of multi-level 3D genome organization during breast cancer progression

Breast cancer entails intricate alterations in genome organization and expression. However, how three-dimensional (3D) chromatin structure changes in the progression from a normal to a breast cancer malignant state remains unknown. To address this, we conducted an analysis combining Hi-C data with lamina-associated domains (LADs), epigenomic marks, and gene expression in an in vitro model of breast cancer progression. Our results reveal that while the fundamental properties of topologically associating domains (TADs) are overall maintained, significant changes occur in the organization of compartments and subcompartments. These changes are closely correlated with alterations in the expression of oncogenic genes. We also observe a restructuring of TAD-TAD interactions, coinciding with a loss of spatial compartmentalization and radial positioning of the 3D genome. Notably, we identify a previously unrecognized interchromosomal insertion event, wherein a locus on chromosome 8 housing the MYC oncogene is inserted into a highly active subcompartment on chromosome 10. This insertion is accompanied by the formation of de novo enhancer contacts and activation of MYC, illustrating how structural genomic variants can alter the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the loss of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

Local Volume Concentration, Packing Domains and Scaling Properties of Chromatin

We propose the Self Returning Excluded Volume (SR-EV) model for the structure of chromatin based on stochastic rules and physical interactions. The SR-EV rules of return generate conformationally-defined domains observed by single cell imaging techniques. From nucleosome to chromosome scales, the model captures the overall chromatin organization as a corrugated system, with dense and dilute regions alternating in a manner that resembles the mixing of two disordered bi-continuous phases. This particular organizational topology is a consequence of the multiplicity of interactions and processes occurring in the nuclei, and mimicked by the proposed return rules. Single configuration properties and ensemble averages show a robust agreement between theoretical and experimental results including chromatin volume concentration, contact probability, packing domain identification and size characterization, and packing scaling behavior. Model and experimental results suggest that there is an inherent chromatin organization regardless of the cell character and resistant to an external forcing such as Rad21 degradation.

Three-dimensional chromatin reorganization regulates B cell development during ageing

The contribution of three-dimensional genome organization to physiological ageing is not well known. Here we show that large-scale chromatin reorganization distinguishes young and old bone marrow progenitor (pro-) B cells. These changes result in increased interactions at the compartment level and reduced interactions within topologically associated domains (TADs). The gene encoding Ebf1, a key B cell regulator, switches from compartment A to B with age. Genetically reducing Ebf1 recapitulates some features of old pro-B cells. TADs that are most reduced with age contain genes important for B cell development, including the immunoglobulin heavy chain (Igh) locus. Weaker intra-TAD interactions at Igh correlate with altered variable (V), diversity (D) and joining (J) gene recombination. Our observations implicate three-dimensional chromatin reorganization as a major driver of pro-B cell phenotypes that impair B lymphopoiesis with age.

Stochasticity in genetics and gene regulation

Development from fertilized egg to functioning multi-cellular organism requires precision. There is no precision, and often no survival, without plasticity. Plasticity is conferred partly by stochastic variation, present inherently in all biological systems. Gene expression levels fluctuate ubiquitously through transcription, alternative splicing, translation and turnover. Small differences in gene expression are exploited to trigger early differentiation, conferring distinct function on selected individual cells and setting in motion regulatory interactions. Non-selected cells then acquire new functions along the spatio-temporal developmental trajectory. The differentiation process has many stochastic components. Meiotic segregation, mitochondrial partitioning, X-inactivation and the dynamic DNA binding of transcription factor assemblies-all exhibit randomness. Non-random X-inactivation generally signals deleterious X-linked mutations. Correct neural wiring, such as retina to brain, arises through repeated confirmatory activity of connections made randomly. In immune system development, both B-cell antibody generation and the emergence of balanced T-cell categories begin through stochastic trial and error followed by functional selection. Aberrant selection processes lead to immune dysfunction. DNA sequence variants also arise through stochastic events: some involving environmental fluctuation (radiation or presence of pollutants), or genetic repair system malfunction. The phenotypic outcome of mutations is also fluid. Mutations may be advantageous in some circumstances, deleterious in others. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Oncogenic enhancers prime quiescent metastatic cells to escape NK immune surveillance by eliciting transcriptional memory

Metastasis arises from disseminated tumour cells (DTCs) that are characterized by intrinsic phenotypic plasticity and the capability of seeding to secondary organs. DTCs can remain latent for years before giving rise to symptomatic overt metastasis. In this context, DTCs fluctuate between a quiescent and proliferative state in response to systemic and microenvironmental signals including immune-mediated surveillance. Despite its relevance, how intrinsic mechanisms sustain DTCs plasticity has not been addressed. By interrogating the epigenetic state of metastatic cells, we find that tumour progression is coupled with the activation of oncogenic enhancers that are organized in variable interconnected chromatin domains. This spatial chromatin context leads to the activation of a robust transcriptional response upon repeated exposure to retinoic acid (RA). We show that this adaptive mechanism sustains the quiescence of DTCs through the activation of the master regulator SOX9. Finally, we determine that RA-stimulated transcriptional memory increases the fitness of metastatic cells by supporting the escape of quiescent DTCs from NK-mediated immune surveillance. Overall, these findings highlight the contribution of oncogenic enhancers in establishing transcriptional memories as an adaptive mechanism to reinforce cancer dormancy and immune escape, thus amenable for therapeutic intervention.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication. Specifically, burst frequency and burst size are up-regulated by E-P communication strength, following power laws with positive exponents. Analysis of the scaling exponents further reveals that burst frequency is preferentially regulated. Bursting kinetics are down-regulated by E-P genomic distance with negative power-law exponents, and this negative modulation desensitizes at large distances. The mutual information between burst frequency (or burst size) and E-P spatial distance further reveals essential characteristics of the information transfer from E-P communication to transcriptional bursting kinetics. These findings, which are in agreement with experimental observations, not only reveal fundamental principles of E-P communication in transcriptional bursting but also are essential for understanding cellular decision-making.

Fire Burn and Cauldron Bubble: What Is in Your Genome Editing Brew?

Genome editing is a rapidly evolving biotechnology with the potential to transform many sectors of industry such as agriculture, biomanufacturing, and medicine. This technology is enabled by an ever-growing portfolio of biomolecular reagents that span the central dogma, from DNA to RNA to protein. In this paper, we draw from our unique perspective as the National Metrology Institute of the United States to bring attention to the importance of understanding and reporting genome editing formulations accurately and promoting concepts to verify successful delivery into cells. Achieving the correct understanding may be hindered by the way units, quantities, and stoichiometries are reported in the field. We highlight the variability in how editing formulations are reported in the literature and examine how a reference molecule could be used to verify the delivery of a reagent into cells. We provide recommendations on how more accurate reporting of editing formulations and more careful verification of the steps in an editing experiment can help set baseline expectations of reagent performance, toward the aim of enabling genome editing studies to be more reproducible. We conclude with a future outlook on technologies that can further our control and enable our understanding of genome editing outcomes at the single-cell level.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance. Interestingly, we find that long-range E-P interactions can induce bimodal and trimodal mRNA distributions. The 4D nucleome equation also allows for model selection and parameter inference. When this equation is applied to the mouse embryonic stem cell smRNA-FISH data and the E-P genomic-distance data, the predicted E-P contact probability and mRNA distribution are in good agreement with experimental results. Further statistical inference indicates that the E-P interactions prefer to modulate the mRNA level by controlling promoter activation and transcription initiation rates. Our model and results provide quantitative insights into both spatiotemporal gene-expression determinants (i.e., long-range E-P interactions) and cellular fates during development.

Local Volume Concentration, Packing Domains and Scaling Properties of Chromatin

We propose the Self Returning Excluded Volume (SR-EV) model for the structure of chromatin based on stochastic rules and physical interactions that is able to capture the observed behavior across imaging and sequencing based measures of chromatin organization. The SR-EV model takes the return rules of the Self Returning Random Walk, incorporates excluded volume interactions, chain connectivity and expands the length scales range from 10 nm to over 1 micron. The model is computationally fast and we created thousands of configurations that we grouped in twelve different ensembles according to the two main parameters of the model. The analysis of the configurations was done in a way completely analogous to the experimental treatments used to determine chromatin volume concentration, contact probability, packing domain identification and size characterization, and packing scaling behavior. We find a robust agreement between the theoretical and experimental results. The overall organization of the model chromatin is corrugated, with dense packing domains alternating with a very dilute regions in a manner that resembles the mixing of two disordered bi-continuous phases. The return rules combined with excluded volume interactions lead to the formation of packing domains. We observed a transition from a short scale regime to a long scale regime occurring at genomic separations of ~ 4 × 104 base pairs or ~ 100 nm in distance. The contact probability reflects this transition with a change in the scaling exponent from larger than -1 to approximately -1. The analysis of the pair correlation function reveals that chromatin organizes following a power law scaling with exponent D ∈ { 2 , 3 } in the transition region between the short and long distance regimes.

Chromosome and protein folding: In search for unified principles

Structural biology has traditionally focused on the structures of proteins, short nucleic acids, small molecules, and their complexes. However, it is now widely recognized that the 3D organization of chromosomes should also be included in this list, despite significant differences in scale and complexity of organization. Here we highlight some notable similarities between the folding processes that shape proteins and chromosomes. Both biomolecules are folded by two types of processes: the affinity-mediated interactions, and by active (ATP-dependent) processes. Both chromosome and proteins in vivo can have partially unstructured and non-equilibrium ensembles with yet to be understood functional roles. By analyzing these biological systems in parallel, we can uncover universal principles of biomolecular organization that transcend specific biopolymers.

Beyond SAHF: An integrative view of chromatin compartmentalization during senescence

Cellular senescence, a persistent form of cell cycle arrest, has been linked to the formation of heterochromatic foci, accompanied by additional concentric epigenetic layers. However, senescence is a highly heterogeneous phenotype, and the formation of these structures is context dependent. Recent developments in the understanding of the high-order chromatin organization have opened new avenues for contextualizing the nuclear and chromatin phenotypes of senescence. Oncogene-induced senescence displays prominent foci and typically exhibits increased chromatin compartmentalization, based on the chromosome conformation assays, as marked by increased transcompaction and segregation of the heterochromatin and euchromatin. However, other types of senescence (e.g., replicative senescence) exhibit comparatively lower levels of compartmentalization. Thus, a more integrative view of the global rearrangement of the chromatin architecture that occurs during senescence is emerging, with potential functional implications for the heterogeneity of the senescence phenotype.

Long interspersed nuclear element 1 and B1/Alu repeats blueprint genome compartmentalization

The organization of the genome into euchromatin and heterochromatin has been known for almost 100 years [1]. More than 50% of mammalian genomes contain repetitive sequences [2,3]. Recently, a functional link between the genome and its folding has been identified [4,5]. Homotypic clustering of long interspersed nuclear element 1 (LINE1 or L1) and B1/Alu retrotransposons forms grossly exclusive nuclear domains that characterize and predict heterochromatin and euchromatin, respectively. The spatial segregation of L1 and B1/Alu-rich compartments is conserved in mammalian cells and can be rebuilt during the cell cycle and established de novo in early embryogenesis. Inhibition of L1 RNA drastically weakened homotypic repeat contacts and compartmental segregation, indicating that L1 plays a more significant role than just being a compartmental marker. This simple and inclusive genetic coding model of L1 and B1/Alu in shaping the macroscopic structure of the genome provides a plausible explanation for the remarkable conservation and robustness of its folding in mammalian cells. It also proposes a conserved core structure on which subsequent dynamic regulation takes place.

Methods for mapping 3D-chromosome architecture around nucleoli

The nucleolus is the largest subcompartment of the nucleus, known to be the place of ribosome biogenesis. Emerging evidence has started to implicate the nucleolus in the organization of chromosomes in the nucleus. Genomic domains contacting the nucleolus are defined as nucleolar associated domains (NADs) and are generally characterized by repressive chromatin states. However, the role of the nucleolus in genome architecture remains still not fully understood mainly because the lack of a membrane has challenged the establishment of methods for accurate identification of NADs. Here, we will discuss recent advances on methods to identify and characterize NADs, discuss their improvements relative to old methods, and highlight future perspectives.

Homologous chromosome associations in domains before meiosis could facilitate chromosome recognition and pairing in wheat

The increasing human population demands an increase in crop yields that must be implemented through breeding programmes to ensure a more efficient and sustainable production of agro-food products. In the framework of breeding, genetic crosses are developed between cultivated species such as wheat and their relative species that are used as genetic donors to transfer desirable agronomic traits into the crop. Unfortunately, interspecific associations between chromosomes from the donor species and the cultivar are rare during meiosis, the process to produce gametes in organisms with sexual reproduction, hampering the transfer of genetic variability into wheat. In addition, little is known about how homologous (equivalent) chromosomes initiate interaction and recognition within the cell nucleus to enter meiosis. In this context, we aim to get insight into wheat chromatin structure, particularly the distribution of homologous chromosomes within the cell nucleus and their putative interactions in premeiotic stages to facilitate chromosome associations and recombination at the beginning of meiosis. Cytogenetics allows the study of both the structure and the behaviour of chromosomes during meiosis and is key in plant breeding. In this study we visualized an extra pair of barley homologous chromosomes in a wheat genetic background to study the spatial distribution, arrangements and interactions occurring exclusively between this pair of homologous chromosomes during premeiosis using fluorescence in situ hybridization (FISH). Our results suggest that homologous chromosomes can initiate interactions in premeiotic stages that could facilitate the processes of specific chromosome recognition and association occurring at the onset of meiosis.