Chromatin Network Analyses: Towards Structure-Function Relationships in Epigenomics
Vera Pancaldi1,2*- 1Centre de Recherches en Cancérologie de Toulouse (CRCT), Institut National de la Santé et de la Recherche Médicale (Inserm) U1037, Centre National de la Recherche Scientifique (CNRS) U5071, Université Paul Sabatier, Toulouse, France
- 2Barcelona Supercomputing Center, Barcelona, Spain
Recent technological advances have allowed us to map chromatin conformation and uncover the genome’s spatial organization of the genome inside the nucleus. These experiments have revealed the complexities of genome folding, characterized by the presence of loops and domains at different scales, which can change across development and in different cell types. There is strong evidence for a relationship between the topological properties of chromatin contacts and cellular phenotype. Chromatin can be represented as a network, in which genomic fragments are the nodes and connections represent experimentally observed spatial proximity of two genomically distant regions in a specific cell type or biological condition. With this approach we can consider a variety of chromatin features in association with the 3D structure, investigating how nuclear chromatin organization can be related to gene regulation, replication, malignancy, phenotypic variability and plasticity. We briefly review the results obtained on genome architecture through network theoretic approaches. As previously observed in protein-protein interaction networks and many types of non-biological networks, external conditions could shape network topology through a yet unidentified structure-function relationship. Similar to scientists studying the brain, we are confronted with a duality between a spatially embedded network of physical contacts, a related network of correlation in the dynamics of network nodes and, finally, an abstract definition of function of this network, related to phenotype. We summarise major developments in the study of networks in other fields, which we think can suggest a path towards better understanding how 3D genome configuration can impact biological function and adaptation to the environment.
Introduction
The 4D Nucleome: Features and Dynamics of Chromatin Contacts
For just over 10 years we have had the tools to explore chromosome conformation inside the nucleus with Hi‐C (Lieberman-Aiden et al., 2009) and related techniques. Understanding the organization of the overall dynamic 3D structure of DNA in the nucleus with associated proteins and nucleic acids [the nucleome (Dekker et al., 2017)] remains an important step towards better understanding the connection between genome and phenotype. Beside the variations of Hi‐C [ChIA-PET (Fullwood et al., 2009), HiChIP (Mumbach et al., 2016), micro-C (Hsieh et al., 2015)], alternative approaches provide an independent picture of nuclear organization. For example, microscopy based techniques allow us to visualize previously tagged chromatin regions for inference of 3D interactions [DNA FISH Oligopaints (Beliveau et al., 2014), Hi-M (Gizzi et al., 2019)]. Moreover, methods that infer 3D proximity by segregating genomic fragments in different locations in the nucleus permit the inference of 3D contact maps independent of proximity ligation [GAM (Beagrie et al., 2017), SPRITE (Quinodoz et al., 2018; Quinodoz et al., 2020)], which is the traditional process to identify chromatin contacts by most Hi‐C and derived methods. Finally, imaging integrated approaches including RNA and chromatin mark quantification in single cells (Takei et al., 2021) are shedding an unprecedented light on single cell nucleome organization.
Given the large regions of the genome that do not contain genes, which are often the focus of our attention, variations of Hi‐C that allow us to study interactions involving a specific subset of chromosomal regions such as Capture Hi‐C have been particularly useful. In particular promoter capture Hi‐C (PCHi‐C) (Schoenfelder et al., 2015a; Mifsud et al., 2015) confirmed the existence and importance of long-range chromatin interactions and identified the role of important regulators of developmental genes, such as the Polycomb complex, in creating a highly interconnected core of genes in mouse embryonic stem cells (Schoenfelder et al., 2015b). Importantly, PCHi‐C allows us to study long-range interactions that can only be captured by Hi‐C using extremely deep sequencing.
The increased availability of chromatin conformation data has enabled detailed analysis of chromatin rearrangements during differentiation (Bonev et al., 2017), progression to malignancy (Malod-Dognin et al., 2020; Vilarrasa-Blasi et al., 2021), senescence (Chandra et al., 2015; Sati et al., 2020) and even during the cell cycle (Nagano et al., 2017; Zhang et al., 2019). Many results have been achieved with extensive use of mathematical models including polymer models, which assume DNA to be a polymer with specific properties that can be related to epigenomic modifications. As summarised in a recent review (Di Stefano et al., 2021), mainly two modelling approaches, a data-driven one and a bottom-up one starting from principles of polymer behaviour, have enabled exploration of dynamics at different scales. Single locus-scale dynamics was modelled in great detail, describing specific processes involved in gene activation and related chromatin changes. On the other hand, coarser models can describe large-scale chromatin rearrangements, compaction and decompaction events and movements to or from the nuclear lamina.
Relating Chromatin Structure to Transcription and Replication
There is some evidence of a tight connection between chromatin contacts and transcription. More specifically, changes in the 3D looping structures of chromatin in neurons following their activation have been shown to affect gene expression on different time scales (Beagan et al., 2020). Also during reprogramming of mouse B cells into pluripotency, specific open chromatin 3D enhancer structures (network hubs) around active genes were shown to form in association to gene expression activation and disappear during gene silencing (Di Stefano et al., 2020). Interestingly, microscopy based techniques have also investigated chromatin dynamics and showed a caging effect that restricts chromatin motion in the proximity of a transcribed gene (Germier et al., 2017), suggesting that chromatin dynamics at various scales is closely associated with transcriptional regulation (Shaban et al., 2020). Whether structure drives the activation of transcription in gene assemblies (transcription factories) (Iborra et al., 1996), or mRNA production reshapes chromatin structure (Hilbert et al., 2021), it is clear that structure and transcription are tightly intertwined. Nevertheless a lack of a clear relation between transcription activation and 3D contacts, at least at single-cell level (Espinola et al., 2021; Ing-Simmons et al., 2021) leaves us with open questions.
On a larger scale, the role of typical structural features, such as Topologically Associating Domains (TADs), in gene expression remains hotly debated (Cavalheiro et al., 2021), with evidence that in some specific loci, perturbing TAD structure can lead to substantial changes in expression patterns, even related to pathology and morphological aberrations (Lupiáñez et al., 2015), while at the genome wide level, expression programs can be resilient to drastic chromosomal rearrangements (Ghavi-Helm et al., 2019).
At a whole nucleome scale, the formation of TAD cliques (highly interconnected groups of TADs) towards the nuclear lamina has been associated with differentiation (Paulsen et al., 2019) and large TAD-level rearrangements have been observed also during reprogramming towards pluripotency (Di Stefano et al., 2020).
Topological analysis of the chromatin network has also led to the identification of drastic rearrangements involved in the progression from healthy to malignant cells. Despite the added complications of considering genomic alterations that are characteristics of cancer cells, clear 3D rearrangements have been observed to accompany oncogenic transformation in B cells, beyond those encountered in normal B cell differentiation (Malod-Dognin et al., 2020; Vilarrasa-Blasi et al., 2021). Exploring nucleome dynamics in differentiation of T cells, it was proposed that some of these rearrangements can ensure the irreversibility and stability of the differentiation status (Hu et al., 2018). Recent results mapping chromatin architecture in neurons suggest the importance of modifications in the TAD structure in neuronal function. More specifically, specific TAD melting events involving particularly long genes important for neuronal functions are up-regulated by local chromatin decompaction (Winick-Ng et al., 2020).
Whereas the concept of a relationship between 3D genome structure and transcription has been controversial and has not been elucidated in the details, it has been known for decades that loops of chromatin are involved in replication and in the formation of replicons, held together by cohesin (Edenberg and Huberman 1975; Guillou et al., 2010). It has also been clear for years that TADs can be seen as units of replication timing (RT), and replication timing domains are clearly connected to phenotypic states and cell types (Ryba et al., 2011; Pope et al., 2014; Boulos et al., 2015), suggesting a strong connection between 3D structure and replication acting on larger scales than those of replicons. Recent work featuring the integration of replication datasets with chromatin contact maps reinforces these findings (Jodkowska et al., 2019; Marchal et al., 2019) and suggests that chromatin might be predominantly organised in structures that ensure specific replication programs, which might become the substrate of transcriptional organization. More specifically, we have found that regions with similar replication timing are close to each other linearly but also in 3D through long range (tens of megabase) interactions (Madrid-Mencía et al., 2020) and we have also observed clustering or chromatin regions containing origins with the same efficiency of activation (effectively the number of cells that use specific origins at any division) across large distances (Jodkowska et al., 2019). Recent reports of the effect of knocking out a main replication regulator, Rif1, have highlighted its possible impact in shaping RT in the entire genome and even an effect on epigenetic marks and organization (Klein et al., 2021). These results have started shedding light on the mechanistic connection between replication, 3D structure and the epigenomic state of the cell.
Variability, Cellular Environment and Epigenomic Rearrangements
Genes that are variable (across single cells) are also more strongly regulated (along time-courses) and more evolvable (across species) (Ciliberti et al., 2007; Lehner 2008; Tirosh et al., 2009; Sigalova et al., 2020). Inter-individual non-genetic variability is of fundamental relevance in medicine and recent studies have highlighted its importance also in healthy individuals (Chen et al., 2016; Ecker et al., 2017). Chromatin is the common substrate that unifies these three different biological processes and we have proposed that the chromatin context, including local and global conformation, could be considered as one of the important factors affecting variability in methylation, gene expression and ultimately phenotype (Ecker et al., 2018). At a local level, it has been shown that the chromatin state can impact gene expression variability. For example, genes repressed by binding of the Polycomb complex, were shown to have more variable expression levels across single cells (Kar et al., 2017). Moreover, promoter features defining specific genomic characteristics of transcriptional start sites (promoter architecture) were identified to be the main determinants of expression variability across single drosophila embryos and across individuals in specific human tissues (Sigalova et al., 2020). Indeed expression variability across individuals was also predictive of differential expression after genetic perturbations or changes in environmental conditions. Promoter features that were found to be important in determining variability include the narrow or broad transcription initiation region, Transcription Factor (TF) occupancy or increased regulatory complexity. Since chromatin properties and TF binding patterns are highly correlated with 3D chromatin interactions, it can be speculated that groups of genes that have similar values for these features would be forming preferential contacts in 3D.
Regions of the genome bound by Polycomb were found to cluster in 3D (Schoenfelder et al., 2015b; Pancaldi et al., 2016), as well as super-enhancer regions (data not shown), suggesting that the folding of chromatin in the nucleus might generate “areas” of higher or lower gene expression variability. This tendency for specific genomic regions to form clusters in 3D can be quantified using Chromatin Assortativity, which measures the significance of correlation between values of a feature on a chromatin region and other regions that interact with it (Pancaldi et al., 2016; Madrid-Mencía et al., 2020) against random expectation. Chromatin assortativity of the promoter features from Sigalova et al. (2020) projected onto a promoter interaction network for human cells derived from chromosome capture experiments indeed shows clustering of promoters of genes with similar sequence and transcriptional features in the 3D genome (Figure 1).
FIGURE 1. A chromatin assortativity analysis performed on the chromatin network shows preferential interactions of genes that have similar promoter characteristics. (A) Schematic describing the chromatin network approach in which nodes are chromatin fragments and edges describe experimentally identified 3D contacts between chromatin fragments, (B) Projection of promoter width from Sigalova et al. (2020) averaged across many tissues onto the combined Promoter Capture Hi‐C network for haematopoietic cells (Javierre et al., 2016). In this network nodes are gene promoters and node colour represents values of “promoter width” with increasing width shown as darker purple. (C) Chromatin assortativity analysis here measures whether promoters with similar values of the promoter features detailed in Sigalova et al. interact preferentially in the network from (A). For each feature a ChAs value is calculated and the values obtained in 100 randomized networks is used to estimate a Z-score, shown on the barplot for each feature. The ChAs calculation is performed using distance preserving randomizations with the ChAseR package (Madrid-Mencía et al., 2020). High values of the Z-score denote preferential 3D distal interactions of promoters with similar values for these features.
So far, we have described how specific regions of the genome can potentially be identified in 3D as having particular characteristic profiles of expression variability. Therefore we can wonder if changes in the external conditions of the cell might affect these nuclear organization patterns to promote adaptation. For example, strong evidence exists for an increase of overall expression and phenotypic variability when cells are grown under stressful environments (López-Maury et al., 2008) or in cancer (Brock et al., 2009) which involves strong regulation of response genes. Looking at a more global level, the mechanisms that control cellular variability in adaptation to external conditions could be reflected as rearrangements in chromatin organization. Indeed, the external environment impacts the organization of chromatin in the nucleus in space and its dynamics in time (the 4D nucleome). This is true across species and contexts, starting from changes of polycomb occupancy during heat-shock in drosophila (Li et al., 2015), to activation of transposable elements associated to chromatin decompaction in heat-stressed plants (Sun et al., 2020), to chromatin architecture rearrangements after osmotic stress in yeast (Amat et al., 2019) and during quiescence (Guidi et al., 2015).
Fibroblasts at late passages also showed telomeres shortening and shifting towards the nuclear center (Bridger et al., 2000), while differentiation, pluripotency and senescence also lead to widespread chromosomal rearrangements and drastic epigenomic changes (Bonev et al., 2017; Sati et al., 2020; Chovanec et al., 2021). These findings reinforce the concept that chromosomal organization is tightly related to cellular phenotype and function, which are strongly connected to the cell’s environment, opening the interesting question of which one between chromatin structure and phenotype is the cause or consequence.
We can therefore hypothesize the existence of a relationship between the global epigenomic characteristics of a cell population, especially 3D chromatin structure, and its plasticity in responding to changing environments or stimulation. To reveal this connection we are in need of a unified framework to represent the epigenome in an integrated fashion and with multiscale resolution.
One of the main challenges in modelling the epigenome is the lack of a complete understanding of the physics of chromatin (intended as DNA and all associated proteins, RNAs and modifications) and the reconciliation of different techniques to generate experimental data which are all affected by specific biases (for example, sequencing biases for Hi‐C related techniques, proximity/contact detection in microscopy based techniques) (Fiorillo et al., 2021). Moreover, we lack tools proposing integrated models covering different time and space scales. Network approaches to chromatin structure can provide a unified multi-scale picture of genome organization, while offering metrics for the overall characteristics of the 3D epigenome, to potentially relate them to phenotypes.
Chromatin Network Approaches
The concept of using network theory to interpret chromatin structure datasets has been around for a decade. As soon as Hi-C contact maps were generated, they were interpreted as distance matrices and easily transformed into adjacency matrices of chromatin networks (Botta et al., 2010; Boulos et al., 2013; Boulos, 2014; Babaei et al., 2015; Morlot et al., 2016; Boulos et al., 2017). With increasing resolution and new chromosome capture techniques more chromatin network approaches were proposed (Pancaldi et al., 2016; Thibodeau et al., 2017; Norton et al., 2018; Huang et al., 2019; Malod-Dognin et al., 2020; Chovanec et al., 2021). The increasing interest in this framework has also prompted the development of 3D genome network visualizers, which allow users to explore the topology (contact patterns) of genomes in a non-linear fashion (Thibodeau et al., 2016; Madrid-Mencía et al., 2020; Ramirez et al., 2020; Chovanec et al., 2021) as well as to visualize other chromatin datasets in a 3D context. This approach can be considered complementary to 3D visualization of constraint based and polymer models of chromatin (Dekker et al., 2013; Lesne et al., 2014; Paulsen et al., 2017; Fiorillo et al., 2021).
In 2012 a pioneering paper by Sandhu et al. (2012) proposed an interesting analysis of the chromatin networks evinced by RNA Polymerase 2 ChIA-PET, a technique which provides 3D contacts between chromatin fragments mediated by polymerase. Despite the incompleteness of the obtained network, they merged different chromosome fragments within specific genomic regions to recover a connected network of promoters from the data set. Detailed analysis of this network’s topology highlighted communities related to specific biological functions and the presence of rich club nodes (highly connected nodes with other high degree nodes), representing important cellular processes, and less connected nodes (spokes), with specific developmental functions and which are enriched for genomic mutations and genetic polymorphisms. These results suggested that evolution has shaped the 3D genome structure.
More recently, similar network topology analyses applied to Hi‐C data targeted at investigating smaller regions within the genome (meso-scale) have identified the presence of core-periphery structures in certain TADs and suggest that the location of single nucleotide variants within these two regions of the TAD can determine the impact the variants have on diseases (Huang et al., 2019).
Network representations are particularly useful when using PCHi‐C datasets. As we have previously shown, promoter capture Hi‐C datasets in mESCs can be easily represented as networks involving promoters and other genomic regions (Other ends or Promoter Interacting Regions) generating networks which feature a main large connected component and various other smaller ones (Pancaldi et al., 2016). Despite not having time-resolved or single-cell chromatin contact networks, we used two statistical network properties (betweenness centrality and bridgeness) to estimate whether nodes with different chromatin marks were highly connected simultaneously with all their partners or with one at a time, akin to the definition of date and party hubs in PPIs (Han et al., 2004). Using these and other network statistics, we established that Polycomb mediated interactions are to some extent more stable and fixed, across cells or time (party hub). On the contrary, RNAPII was predicted to have more variable interactions, changing its partners across cells or time (date hub). The tools that we have produced (the ChAseR package and GARDEN-NET website) allow seamless integration and analysis of any property defined along the genome (gene properties such as transcriptional levels or variability, functional categories or evolutionary characteristics, chromatin features, or any user-defined data) onto the chromatin structures defined by any experimental chromatin contact datasets (Madrid-Mencía et al., 2020).
Importantly, networks allow us to visualize chromatin at different scales. Since TADs are defined similarly to network modules, namely as domains that have more interactions within them than with other domains, we can imagine grouping chromatin fragments within a TAD into a single node and investigating a network of TADs, leading to a hierarchical representation of the genome (Jodkowska et al., 2019). At the same time we can combine different datasets (for example Hi‐C and PCHi‐C) in a single network, to preserve a high number of long-range contacts while also considering the linear proximity of regions along the genome (Malod-Dognin et al., 2020).
Chovanec et al. (2021) presented a network approach to visualize and study the difference in 3D chromatin conformation between naive and primed human pluripotent stem cells (PSCs) using Canvas, a tool to represent chromatin fragment interaction networks in a multi-scale framework. They identified cell type specific clusters displaying coordinated gene expression and corresponding in some cases to functional units. Interestingly, they showed that TAD border insulation (separation between TADs) is stronger in primed compared to naive PSC, suggesting a coupling between phenotype and structure that would see naive cells with a “weaker” TAD structure. They suggested that smaller communities would come together into larger domains during the priming process, which is compatible with the formation of more functionally constrained units relying on long-range (>1 MB) polycomb enriched interactions between developmentally related genes.
Interestingly, the application of network theory analyses to chromatin allows us to peer into the organizing principles that underlie the nucleome in search for a structure-function relationship in chromatin networks.
Chromatin Structure-Function Relationships: Multi-Disciplinary Concepts Towards an Epigenome-Phenotype Connection
Applying network theory approaches to biological networks has provided insight on different aspects of biological processes, from evolution and the structure of protein-protein interaction networks, all the way to mechanistic hypotheses on mode of action of drugs from looking at drug-disease-symptom interaction networks (Sonawane et al., 2019). The common theme through these approaches is the realization that a structure-function relationship exists between the topology and structural characteristics of a given network and the behaviour of the entity represented by that network.
To better examine the concept of structure-function relationships, we can turn to another biological entity that rivals the genome in complexity and mystery: the brain. Two separate approaches can be used to describe and study brain networks (Lynn and Bassett 2019). First, through imaging and tracing of water molecules along white matter, we can define a spatially embedded structural network connecting different regions of the brain. This approach is so far limited to a meso-scale level, involving the grouping of different neurons, but the latest technology can currently capture the dynamics of these meso-scale contacts in humans and primates. Second, we can take into account the functional relationship between different neurons, according to the principle that neural connections are strengthened when the two connected neurons fire together and that the emergence of collective firing behaviour of multiple neurons can give rise to cognitive functions. This approach, which relies on functional Magnetic Resonance Imaging (fMRI) allows us to identify correlation of activities between different nodes of the structural network (normally only possible at the brain region level including tens of thousands of neurons). A certain correlation between structural brain networks and functional brain networks is expected and identified, but the relation is certainly not trivial and the dynamics of the functional correlation might have to be considered (Suárez et al., 2020). These two networks and their relationship are specific to particular developmental stages and altered in mental disease.
A correspondence between genome and brain networks becomes evident, with a parallelism between chromatin and brain structural networks versus co-expression and functional brain networks (Figure 2). Whereas in brain networks the clear connection between the two networks is represented by the signals exchanged between neurons, we are still discovering how the numerous molecular processes taking place inside a cell could provide a connection between 3D genome and phenotype. We could imagine the chromatin structure network and the gene/protein regulatory interaction or co-expression network as two components of a multilayer network (Kivelä et al., 2014), with phenotypes defined as an emergent collective behaviour on the gene co-expression network. A multilayer network has been successfully applied to the C. elegans connectome, showing the potential of such an approach in epigenomics (Bentley et al., 2016). A further important point of convergence between these two network types is their spatial character. More and more we are able to study chromatin structure through microscopy in single cells, providing us with spatial coordinates of genomic fragments. These new approaches allow us to go beyond the topology of the structure and create real maps of the genome inside each nucleus. A spatial perspective on chromatin could suggest exploiting a range of techniques that have been applied to other kinds of spatial networks, which normally focus on spreading phenomena, including mobility, epidemiology and power-grid networks, among many others (Barthélemy 2011). Such an approach in epigenomics would open new perspectives while requiring thorough considerations of the physics of the nucleus. To what extent an equivalence between spatial distance in chromatin and geographical distance is appropriate, together with the importance of transport on the chromatin network [information flow (Juan et al., 2016) or physical flow of molecules (Cortini and Filion 2018)] are likely to become topics of interest. A related issue that could have interesting repercussions in the epigenomics of diseases is the study of failures in these spatial networks and how resilient these networks are to attack (Schneider et al., 2011).
FIGURE 2. Proposing a structure-function relationship linking the epigenome to the phenotype. We propose an analogy between the brain and chromatin represented as multilayer networks. In the brain, a structural network of physical connections between neurons and brain regions can be identified and related to a functional network describing correlation in firing of distant brain regions, which collectively give rise to cognitive functions. Similarly, the physical interactions of genomic regions inside the nucleus (left bottom layer: chromatin network, nodes are gene promoters with contacts retrieved by promoter-capture Hi‐C) can be related to gene-coexpression networks (left top layer: nodes are genes and edges represent correlation of expression between them across different conditions) that ultimately determine cellular phenotypes, which can be interpreted as specific gene expression patterns associated to distinct phenotypes.
To further study this structure-function relationship and relate the cellular network to its environment, we can rely on results connecting network properties to their environment and behaviour. Multiple network metrics and topological features have been employed to describe networks’ response to specific conditions or perturbations, their resilience or plasticity (robustness or sensitivity) (Bar-Yam and Epstein 2004; Schneider et al., 2011; Csermely et al., 2015). Modularity in protein interaction networks has been attributed to the need for optimizing multiple tasks at once (Kashtan and Alon 2005) and was also observed in metabolic networks (Ravasz et al., 2002), again suggesting that structure in biological networks reflects environmental pressure. Particularly when looking at the topology of gene regulatory networks, changes in their architecture can be related to the cell’s environment (Luscombe et al., 2004).
In yeast, stress is associated with a differentiation of the expression programme, which can be reflected into a global rearrangement of the transcriptome associated with a drastic change of the protein interaction network. More specifically, heat stress produces a disaggregation in the interaction network of budding yeast (Mihalik and Csermely 2011) and we were ourselves able to show a similar effect in fission yeast (Lehtinen et al., 2013). We observed a decrease in module interconnectedness after stress treatment and measured a decrease in the overlap between network modules, representing a clear systemic disaggregation of the network topology. Interestingly, other approaches focussed on changes in metabolism before and after stress treatment arrived at similar conclusions using metabolic networks (Gopalacharyulu et al., 2009).
Some of these effects are reminiscent of phenomena that have been observed in social and ecological networks. In complex systems like social networks of traders in the stock market or schools of fish in the ocean, specific dynamics of synchronization have been related to success [prey evasion for fish or large profits for traders (Saavedra et al., 2011)]. Moreover, the structures of messaging networks (in which nodes are traders, edges represent instant messaging interactions across the day) were shown to be drastically affected by the situation of the stock market (Romero et al., 2016). In days of relative calm and predictable stock market fluctuations, traders exchange messages with all other traders, whereas on the days of strongly unpredictable price shocks the trader’s networks display a “turtled up” structure, with predominant connections to the network core. These are clear examples of relationships between network structure and its function in response to the environment.
In line with this view, it has been proposed that specific components in a biological network can have different functions, namely the core would serve to ensure a prompt response to known stimuli, while the periphery would contain elements needed to respond to unexpected situations (Csermely 2018). This interesting perspective implies that, on evolutionary timescales, repeated stimuli would lead to a change in the structure of the network corresponding to learning of the new conditions. Moreover, depending on the level of predictability of external conditions, the network might need to adapt its structure according to external conditions. These principles could govern the chromatin rearrangements that are observed in cells in response to changes in their external environment.
Discussion
In summary, a large body of literature exists regarding system level analyses on protein-protein or gene-coexpression networks as well as networks from other disciplines. The notion that a cell’s phenotype is somehow encoded in its genome suggests that looking at the epigenomic networks under this light could also be relevant. When considering the 3D genome, Csermely’s hypothesis (Csermely 2018) might suggest that, through evolution, a core 3D network of highly conserved genes, involved in basic cellular functions, has been expanded by addition of new genes involved in new functions, specific to cell types and tissues. Importantly, changes in chromatin organization can thus encode changes in the cellular phenotype in response to varying needs of the cell to adapt to its external context. We anticipate that a network science perspective on (epi)genomes will continue to provide insights and new understanding for exploiting the rich datasets that will be produced across cell types, with potential applications in cancer and disease in general.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.embopress.org/action/downloadSupplement?doi=10.15252%2Fmsb.20209539&file=msb209539-sup-0002-Datasets.zip.
Author Contributions
VP conceived and wrote the manuscript.
Funding
VP research is funded by INSERM, Fondation Toulouse Cancer Santé and Pierre Fabre Research Institute as part of the Chair of Bioinformatics in Oncology of the CRCT.
Conflict of Interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We acknowledge Flavien Raynal for help with generating Figure 1 and Alexis Coullomb for critical reading of the manuscript.
References
Amat, R., Böttcher, R., Le Dily, F., Vidal, E., Quilez, J., Cuartero, Y., et al. (2019). Rapid Reversible Changes in Compartments and Local Chromatin Organization Revealed by Hyperosmotic Shock. Genome Res. 29 (1), 18–28. doi:10.1101/gr.238527.118
Babaei, S., Mahfouz, A., Hulsman, M., Lelieveldt, B. P., de Ridder, J., and Reinders, M. (2015). Hi-C Chromatin Interaction Networks Predict Co-expression in the Mouse Cortex. Plos Comput. Biol. 11 (5), e1004221. doi:10.1371/journal.pcbi.1004221
Bar-Yam, Y., and Epstein, I. R. (2004). Response of Complex Networks to Stimuli. Proc. Natl. Acad. Sci. 101 (13), 4341–4345. doi:10.1073/pnas.0400673101
Barthélemy, M. (2011). Spatial Networks. Phys. Rep. 499 (1), 1–101. doi:10.1016/j.physrep.2010.11.002
Beagan, J. A., Pastuzyn, E. D., Fernandez, L. R., Guo, M. H., FengGuo, K., Titus, K. R., et al. (2020). Three-Dimensional Genome Restructuring across Timescales of Activity-Induced Neuronal Gene Expression. Nat. Neurosci. 23 (6), 707–717. doi:10.1038/s41593-020-0634-6
Beagrie, R. A., Scialdone, A., Schueler, M., Kraemer, D. C., Chotalia, M., Xie, S. Q., et al. (2017). Complex Multi-Enhancer Contacts Captured by Genome Architecture Mapping. Nature 543 (7646), 519–524. doi:10.1038/nature21411
Beliveau, B. J., Apostolopoulos, N., and Wu, C. T. (2014). Visualizing Genomes with Oligopaint FISH Probes. Curr. Protoc. Mol. Biol. 105 (January), 14.23. doi:10.1002/0471142727.mb1423s105
Bentley, B., Branicky, R., Barnes, C. L., Chew, Y. L., Yemini, E., Bullmore, E. T., et al. (2016). The Multilayer Connectome of Caenorhabditis Elegans. Plos Comput. Biol. 12 (12), e1005283. doi:10.1371/journal.pcbi.1005283
Bonev, B., Mendelson Cohen, N., Szabo, Q., Fritsch, L., Papadopoulos, G. L., Lubling, Y., et al. (2017). Multiscale 3D Genome Rewiring during Mouse Neural Development. Cell 171 (3), 557–e24. doi:10.1016/j.cell.2017.09.043
Botta, M., Haider, S., Leung, I. X., Lio, P., Mozziconacci, J., and Mozziconacci, J. (2010). Intra- and Inter-chromosomal Interactions Correlate with CTCF Binding Genome Wide. Mol. Syst. Biol. 6, 426. doi:10.1038/msb.2010.79
Boulos, R. E., Arneodo, A., Jensen, P., and Audit, B. (2013). Revealing Long-Range Interconnected Hubs in Human Chromatin Interaction Data Using Graph Theory. Phys. Rev. Lett. 111 (11), 118102. doi:10.1103/PhysRevLett.111.118102
Boulos, R. E., Drillon, G., Argoul, F., Arneodo, A., and Audit, B. (2015). Structural Organization of Human Replication Timing Domains. FEBS Lett. 589 (20PartA), 2944–2957. doi:10.1016/j.febslet.2015.04.015
Boulos, R. E., Tremblay, N., Arneodo, A., Borgnat, P., and Audit, B. (2017). Multi-Scale Structural Community Organisation of the Human Genome. BMC Bioinformatics 18 (1), 209. doi:10.1186/s12859-017-1616-x
Boulos, R. E., Julienne, H., Baker, A., Chen, C.-L., Petryk, N., Kahli, M., et al. (2014). From the Chromatin Interaction Network to the Organization of the Human Genome into Replication N/U-Domains. New J. Phys. 16 (11), 115014. doi:10.1088/1367-2630/16/11/115014
Bridger, J. M., Boyle, S., Kill, I. R., and Bickmore, W. A. (2000). Re-Modelling of Nuclear Architecture in Quiescent and Senescent Human Fibroblasts. Curr. Biol. 10 (3), 149–152. doi:10.1016/S0960-9822(00)00312-2
Brock, A., Chang, H., and Huang, S. (2009). Non-genetic Heterogeneity-Aa Mutation-independent Driving Force for the Somatic Evolution of Tumours. Nat. Rev. Genet. 10 (5), 336–342. doi:10.1038/nrg2556
Cardozo Gizzi, A. M., Cattoni, D. I., Fiche, J. B., Espinola, S. M., Gurgo, J., Messina, O., et al. (2019). Microscopy-Based Chromosome Conformation Capture Enables Simultaneous Visualization of Genome Organization and Transcription in Intact Organisms. Mol. Cel 74 (1), 212–e5. doi:10.1016/j.molcel.2019.01.011
Cavalheiro, G. R., Pollex, T., and Furlong, E. E. (2021). To Loop or Not to Loop: What Is the Role of TADs in Enhancer Function and Gene Regulation? Curr. Opin. Genet. Develop. 67 (April), 119–129. doi:10.1016/j.gde.2020.12.015
Chandra, T., Ewels, P. A., Schoenfelder, S., Furlan-Magaril, M., Wingett, S. W., Kirschner, K., et al. (2015). Global Reorganization of the Nuclear Landscape in Senescent Cells. Cell Rep 10 (4), 471–483. doi:10.1016/j.celrep.2014.12.055
Chen, L., Ge, B., Casale, F. P., Vasquez, L., Kwan, T., Garrido-Martín, D., et al. (2016). Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells. Cell 167 (5), 1398–e24. doi:10.1016/j.cell.2016.10.026
Chovanec, P., Collier, A. J., Krueger, C., Várnai, C., Semprich, C. I., Schoenfelder, S., et al. (2021). Widespread Reorganisation of Pluripotent Factor Binding and Gene Regulatory Interactions between Human Pluripotent States. Nat. Commun. 12 (1), 2098. doi:10.1038/s41467-021-22201-4
Ciliberti, S., Martin, O. C., and Wagner, A. (2007). Robustness Can Evolve Gradually in Complex Regulatory Gene Networks with Varying Topology. Plos Comput. Biol. 3 (2), e15. doi:10.1371/journal.pcbi.0030015
Cortini, R., and Filion, G. J. (2018). Theoretical Principles of Transcription Factor Traffic on Folded Chromatin. Nat. Commun. 9 (1), 1740. doi:10.1038/s41467-018-04130-x
Csermely, P., Hódsági, J., Korcsmáros, T., Módos, D., Perez-LopezPerez-Lopez, Á. R., Szalay, K., et al. (2015). Cancer Stem Cells Display Extremely Large Evolvability: Alternating Plastic and Rigid Networks as a Potential Mechanism: Network Models, Novel Therapeutic Target Strategies, and the Contributions of Hypoxia, Inflammation and Cellular Senescence. Semin. Cancer Biol. 30 (February), 42–51. doi:10.1016/j.semcancer.2013.12.004
Csermely, P. (2018). The Wisdom of Networks: A General Adaptation and Learning Mechanism of Complex Systems: The Network Core Triggers Fast Responses to Known Stimuli; Innovations Require the Slow Network Periphery and Are Encoded by Core-Remodeling. Bioessays 40 (1). doi:10.1002/bies.201700150
Dekker, J., Belmont, A. S., Guttman, M., Leshyk, V. O., Lis, J. T., Lomvardas, S., et al. (2017). The 4D Nucleome Project. Nature 549 (7671), 219–226. doi:10.1038/nature23884
Dekker, J., Marti-Renom, M. A., and Mirny, L. A. (2013). Exploring the Three-Dimensional Organization of Genomes: Interpreting Chromatin Interaction Data. Nat. Rev. Genet. 14 (6), 390–403. doi:10.1038/nrg3454
Di Stefano, M., Stadhouders, R., Farabella, I., Castillo, D., Serra, F., Graf, T., et al. (2020). Transcriptional Activation during Cell Reprogramming Correlates with the Formation of 3D Open Chromatin Hubs. Nat. Commun. 11 (1), 2564. doi:10.1038/s41467-020-16396-1
Di Stefano, M., Paulsen, J., Jost, D., and Marti-Renom, M. A. (2021). 4D Nucleome Modeling. Curr. Opin. Genet. Develop. 67 (April), 25–32. doi:10.1016/j.gde.2020.10.004
Ecker, S., Chen, L., Pancaldi, V., Bagger, F. O., Fernández, J. M., Carrillo de Santa Pau, E., et al. (2017). Genome-Wide Analysis of Differential Transcriptional and Epigenetic Variability across Human Immune Cell Types. Genome Biol. 18 (1), 18. doi:10.1186/s13059-017-1156-8
Ecker, S., Pancaldi, V., Valencia, A., Beck, S., and Paul, D. S. (2018). Epigenetic and Transcriptional Variability Shape Phenotypic Plasticity. BioEssays 40 (2). doi:10.1002/bies.201700148
Edenberg, H. J., and Huberman, J. A. (1975). Eukaryotic Chromosome Replication. Annu. Rev. Genet. 9 (1), 245–284. doi:10.1146/annurev.ge.09.120175.001333
Espinola, S. M., Götz, M., Bellec, M., Messina, O., Fiche, J. B., Houbron, C., et al. (2021). Cis-Regulatory Chromatin Loops Arise before TADs and Gene Activation, and Are Independent of Cell Fate during Early Drosophila Development. Nat. Genet. 53 (4), 477–486. doi:10.1038/s41588-021-00816-z
Fiorillo, L., Musella, F., Conte, M., Kempfer, R., Chiariello, A. M., Bianco, S., et al. (2021). Comparison of the Hi-C, GAM and SPRITE Methods Using Polymer Models of Chromatin. Nat. Methods 18 (5), 482–490. doi:10.1038/s41592-021-01135-1
Fullwood, M. J., Liu, M. H., Liu Pan, Y. F., Pan, J., Xu, H., Mohamed, Y. B., et al. (2009). An Oestrogen-Receptor-Alpha-Bound Human Chromatin Interactome. Nature 462 (7269), 58–64. doi:10.1038/nature08497
Germier, T., Kocanova, S., Walther, N., Bancaud, A., Shaban, H. A., Sellou, H., et al. (2017). Real-Time Imaging of a Single Gene Reveals Transcription-Initiated Local Confinement. Biophys. J. 113 (7), 1383–1394. doi:10.1016/j.bpj.2017.08.014
Ghavi-Helm, Y., Jankowski, A., Meiers, S., Viales, R. R., Korbel, J. O., and Furlong, E. E. M. (2019). Highly Rearranged Chromosomes Reveal Uncoupling between Genome Topology and Gene Expression. Nat. Genet. 51 (8), 1272–1282. doi:10.1038/s41588-019-0462-3
Gopalacharyulu, P. V., Velagapudi, V. R., Lindfors, E., Halperin, E., and Oresic, M. (2009). Dynamic Network Topology Changes in Functional Modules Predict Responses to Oxidative Stress in Yeast. Mol. Biosyst. 5 (3), 276–287. doi:10.1039/b815347g
Guidi, M., Ruault, M., Marbouty, M., Loïodice, I., Cournac, A., Billaudeau, C., et al. (2015). Spatial Reorganization of Telomeres in Long-Lived Quiescent Cells. Genome Biol. 16 (January), 206. doi:10.1186/s13059-015-0766-2
Guillou, E., Ibarra, A., Coulon, V., Casado-Vela, J., Rico, D., Casal, I., et al. (2010). Cohesin Organizes Chromatin Loops at DNA Replication Factories. Genes Dev. 24 (24), 2812–2822. doi:10.1101/gad.608210
Han, J. D., Bertin, N., Hao, T., Goldberg, D. S., Berriz, G. F., Zhang, L. V., et al. (2004). Evidence for Dynamically Organized Modularity in the Yeast Protein-Protein Interaction Network. Nature 430 (6995), 88–93. doi:10.1038/nature02555
Hilbert, L., Sato, Y., Kuznetsova, K., Bianucci, T., Kimura, H., Jülicher, F., et al. (2021). Transcription Organizes Euchromatin via Microphase Separation. Nat. Commun. 12 (1), 1360. doi:10.1038/s41467-021-21589-3
Hsieh, T. H., Weiner, A., Lajoie, B., Dekker, J., Friedman, N., Rando, O. J., et al. (2015). Mapping Nucleosome Resolution Chromosome Folding in Yeast by Micro-C. Cell 162 (1), 108–119. doi:10.1016/j.cell.2015.05.048
Hu, G., Cui, K., Fang, D., Hirose, S., Wang, X., Wangsa, D., et al. (2018). Transformation of Accessible Chromatin and 3D Nucleome Underlies Lineage Commitment of Early T Cells. Immunity 48 (2), 227–e8. doi:10.1016/j.immuni.2018.01.013
Huang, H., Chen, S. T., Titus, K. R., Titus, D. J., Emerson, D. S., and Phillips-Cremins, J. E. (2019). A Subset of Topologically Associating Domains Fold into Mesoscale Core-Periphery Networks. Sci. Rep. 9 (1), 9526. doi:10.1038/s41598-019-45457-9
Iborra, F. J., Pombo, A., Jackson, D. A., and Cook, P. R. (1996). Active RNA Polymerases Are Localized within Discrete Transcription ‘factories' in Human Nuclei. J. Cel Sci 109 (Pt 6), 1427–1436. doi:10.1242/jcs.109.6.1427
Ing-Simmons, E., Vaid, R., Bing, X. Y., Levine, M., Mannervik, M., and Vaquerizas, J. M. (2021). Independence of Chromatin Conformation and Gene Regulation during Drosophila Dorsoventral Patterning. Nat. Genet. 53 (4), 487–499. doi:10.1038/s41588-021-00799-x
Javierre, B. M., Burren, O. S., Wilder, S. P., Kreuzhuber, R., Hill, S. M., Sewitz, S., et al. (2016). Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters. Cell 167 (5), 1369–e19. doi:10.1016/j.cell.2016.09.037
Jodkowska, K., Pancaldi, V., Almeida, R., Rigau, M., Graña-Castro, O., Fernández-JustelFernández-Justel, J. M., et al. (2019). Three-Dimensional Connectivity and Chromatin Environment Mediate the Activation Efficiency of Mammalian DNA Replication Origins. BioRxiv [Preprint]. doi:10.1101/644971
Juan, D., Perner, J., Carrillo de Santa Pau, E., Marsili, S., Ochoa, D., Chung, H. R., et al. (2016). Epigenomic Co-localization and Co-evolution Reveal a Key Role for 5hmC as a Communication Hub in the Chromatin Network of ESCs. Cel Rep 14 (5), 1246–1257. doi:10.1016/j.celrep.2016.01.008
Kar, G., Kim, J. K., Kolodziejczyk, A. A., Natarajan, K. N., Torlai Triglia, E., Mifsud, B., et al. (2017). Flipping between Polycomb Repressed and Active Transcriptional States Introduces Noise in Gene Expression. Nat. Commun. 8 (1), 36. doi:10.1038/s41467-017-00052-2
Kashtan, N., and Alon, U. (2005). Spontaneous Evolution of Modularity and Network Motifs. Proc. Natl. Acad. Sci. U S A. 102 (39), 13773–13778. doi:10.1073/pnas.0503610102
Kivelä, M., Arenas, A., Barthelemy, M., Gleeson, J. P., Moreno, Y., and Porter, M. A. (2014). Multilayer Networks. J. Complex Networks 2 (3), 203–271. doi:10.1093/comnet/cnu016
Klein, K. N., Zhao, P. A., Lyu, X., Sasaki, T., Bartlett, D. A., Singh, A. M., et al. (2021). Replication Timing Maintains the Global Epigenetic State in Human Cells. Science 372 (6540), 371–378. doi:10.1126/science.aba5545
Lehner, B. (2008). Selection to Minimise Noise in Living Systems and its Implications for the Evolution of Gene Expression. Mol. Syst. Biol. 4 (January), 170. doi:10.1038/msb.2008.11
Lehtinen, S., Marsellach, F. X., Codlin, S., Schmidt, A., Clément-Ziza, M., Beyer, A., et al. (2013). Stress Induces Remodelling of Yeast Interaction and Co-expression Networks. Mol. Biosyst. 9 (7), 1697–1707. doi:10.1039/c3mb25548d
Lesne, A., Riposo, J., Roger, P., Cournac, A., and Mozziconacci, J. (2014). 3D Genome Reconstruction from Chromosomal Contacts. Nat. Methods 11 (11), 1141–1143. doi:10.1038/nmeth.3104
Li, L., Lyu, X., Hou, C., Takenaka, N., Huy, Q., Ong, C. T., et al. (2015). Widespread Rearrangement of 3D Chromatin Organization Underlies Polycomb-Mediated Stress-Induced Silencing. Mol. Cel 58 (2), 216–231. doi:10.1016/j.molcel.2015.02.023
Lieberman-Aiden, E., Van Berkum, N. L., Williams, L., Imakaev, M., Ragoczy, T., Telling, A., et al. (2009). Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome. Science 326 (5950), 289–293. doi:10.1126/science.1181369
López-Maury, L., Marguerat, S., and Bähler, J. (2008). Tuning Gene Expression to Changing Environments: From Rapid Responses to Evolutionary Adaptation. Nat. Rev. Genet. 9 (8), 583–593. doi:10.1038/nrg2398
Lupiáñez, D. G., Kraft, K., Heinrich, V., Krawitz, P., Brancati, F., Klopocki, E., et al. (2015). Disruptions of Topological Chromatin Domains Cause Pathogenic Rewiring of Gene-Enhancer Interactions. Cell 161 (5), 1012–1025. doi:10.1016/j.cell.2015.04.004
Luscombe, N. M., Babu, M. M., Yu, H., Snyder, M., Teichmann, S. A., and Gerstein, M. (2004). Genomic Analysis of Regulatory Network Dynamics Reveals Large Topological Changes. Nature 431 (7006), 308–312. doi:10.1038/nature02782
Lynn, C. W., and Bassett, D. S. (2019). The Physics of Brain Network Structure, Function and Control. Nat. Rev. Phys. 1 (5), 318–332. doi:10.1038/s42254-019-0040-8
Madrid-Mencía, M., Raineri, E., Cao, T. B. N., and Pancaldi, V. (2020). Using GARDEN-NET and ChAseR to Explore Human Haematopoietic 3D Chromatin Interaction Networks. Nucleic Acids Res. 48 (8), 4066–4080. doi:10.1093/nar/gkaa159
Malod-Dognin, N., Pancaldi, V., Valencia, A., and Pržulj, N. (2020). Chromatin Network Markers of Leukemia. Bioinformatics 36 (Suppl. l_1), i455–63. doi:10.1093/bioinformatics/btaa445
Marchal, C., Sima, J., and Gilbert, D. M. (2019). Control of DNA Replication Timing in the 3D Genome. Nat. Rev. Mol. Cel Biol 20 (12), 721–737. doi:10.1038/s41580-019-0162-y
Mifsud, B., Tavares-Cadete, F., Young, A. N., Sugar, R., Schoenfelder, S., Ferreira, L., et al. (2015). Mapping Long-Range Promoter Contacts in Human Cells with High-Resolution Capture Hi-C. Nat. Genet. 47 (6), 598–606. doi:10.1038/ng.3286
Mihalik, Á., and Csermely, P., (2011). Heat Shock Partially Dissociates the Overlapping Modules of the Yeast Protein-Protein Interaction Network: A Systems Level Model of Adaptation. Plos Comput. Biol. 7 (10), e1002187. doi:10.1371/journal.pcbi.1002187
Morlot, J.-B., Mozziconacci, J., and Lesne, A. (2016). Network Concepts for Analyzing 3D Genome Structure from Chromosomal Contact Maps. EPJ Nonlinear Biomed. Phys. 4 (1), 1–15. doi:10.1140/epjnbp/s40366-016-0029-5
Mumbach, M. R., Rubin, A. J., Flynn, R. A., Dai, C., Khavari, P. A., Greenleaf, W. J., et al. (2016). HiChIP: Efficient and Sensitive Analysis of Protein-Directed Genome Architecture. Nat. Methods 13 (11), 919–922. doi:10.1038/nmeth.3999
Nagano, T., Lubling, Y., Várnai, C., Dudley, C., Leung, W., Baran, Y., et al. (2017). Cell-cycle Dynamics of Chromosomal Organization at Single-Cell Resolution. Nature 547 (7661), 61–67. doi:10.1038/nature23001
Norton, H. K., Emerson, D. J., Huang, H., Kim, J., Titus, K. R., Gu, S., et al. (2018). Detecting Hierarchical Genome Folding with Network Modularity. Nat. Methods 15 (2), 119–122. doi:10.1038/nmeth.4560
Pancaldi, V., Carrillo-de-Santa-Pau, E., Javierre, B. M., Juan, D., Fraser, P., Spivakov, M., et al. (2016). Integrating Epigenomic Data and 3D Genomic Structure with a New Measure of Chromatin Assortativity. Genome Biol. 17 (1), 152. doi:10.1186/s13059-016-1003-3
Paulsen, J., Liyakat Ali, T. M., Nekrasov, M., Delbarre, E., Baudement, M.-O., Kurscheid, S., et al. (2019). Long-Range Interactions between Topologically Associating Domains Shape the Four-Dimensional Genome during Differentiation. Nat. Genet. 51 (5), 835–843. doi:10.1038/s41588-019-0392-0
Paulsen, J., Sekelja, M., Oldenburg, A. R., Barateau, A., Briand, N., Delbarre, E., et al. (2017). Chrom3D: Three-Dimensional Genome Modeling from Hi-C and Nuclear Lamin-Genome Contacts. Genome Biol. 18 (1), 21. doi:10.1186/s13059-016-1146-2
Pope, B. D., Ryba, T., Dileep, V., Yue, F., Wu, W., Denas, O., et al. (2014). Topologically Associating Domains Are Stable Units of Replication-Timing Regulation. Nature 515 (7527), 402–405. doi:10.1038/nature13986
Quinodoz, S. A., Bhat, P., Ollikainen, N., Jachowicz, J. W., Banerjee, A. K., Chovanec, P., et al. (2020). RNA Promotes the Formation of Spatial Compartments in the Nucleus. BioRxiv [Preprint]. doi:10.1101/2020.08.25.267435
Quinodoz, S. A., Ollikainen, N., Tabak, B., Palla, A., Schmidt, J. M., Detmar, E., et al. (2018). Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus. Cell 174 (3), 744–e24. doi:10.1016/j.cell.2018.05.024
Ramirez, R. N., Bedirian, K., Gray, S. M., and Diallo, A. (2020). DNA Rchitect: an R Based Visualizer for Network Analysis of Chromatin Interaction Data. Bioinformatics 36 (2), 644–646. doi:10.1093/bioinformatics/btz608
Ravasz, E., Somera, A. L., Mongru, D. A., Oltvai, Z. N., and Barabási, A. L. (2002). Hierarchical Organization of Modularity in Metabolic Networks. Science 297 (5586), 1551–1555. doi:10.1126/science.1073374
Romero, D. M., Uzzi, B., and Kleinberg, J. (2016). Social Networks under Stress. ArXiv:1602.00572 [Physics] [Preprint]. doi:10.1145/2872427.2883063l
Ryba, T., Hiratani, I., Sasaki, T., Battaglia, D., Kulik, M., Zhang, J., et al. (2011). Replication Timing: a Fingerprint for Cell Identity and Pluripotency. Plos Comput. Biol. 7 (10), e1002225. doi:10.1371/journal.pcbi.1002225
Saavedra, S., Hagerty, K., and Uzzi, B. (2011). Synchronicity, Instant Messaging, and Performance Among Financial Traders. Proc. Natl. Acad. Sci. U S A. 108 (13), 5296–5301. doi:10.1073/pnas.1018462108
Sandhu, K. S., Li, G., Poh, H. M., Quek, Y. L., Sia, Y. Y., Peh, S. Q., et al. (2012). Large-Scale Functional Organization of Long-Range Chromatin Interaction Networks. Cel Rep 2 (5), 1207–1219. doi:10.1016/j.celrep.2012.09.022
Sati, S., Bonev, B., Szabo, Q., Jost, D., Bensadoun, P., Serra, F., et al. (2020). 4D Genome Rewiring during Oncogene-Induced and Replicative Senescence. Mol. Cel 78 (3), 522–e9. doi:10.1016/j.molcel.2020.03.007
Schneider, C. M., Moreira, A. A., Andrade, J. S., Havlin, S., and Herrmann, H. J. (2011). Mitigation of Malicious Attacks on Networks. Proc. Natl. Acad. Sci. U S A. 108 (10), 3838–3841. doi:10.1073/pnas.1009440108
Schoenfelder, S., Furlan-Magaril, M., Mifsud, B., Tavares-Cadete, F., Sugar, R., Javierre, B. M., et al. (2015a). The Pluripotent Regulatory Circuitry Connecting Promoters to Their Long-Range Interacting Elements. Genome Res. 25 (4), 582–597. doi:10.1101/gr.185272.114
Schoenfelder, S., Sugar, R., Dimond, A., Javierre, B. M., Armstrong, H., Mifsud, B., et al. (2015b). Polycomb Repressive Complex PRC1 Spatially Constrains the Mouse Embryonic Stem Cell Genome. Nat. Genet. 47 (10), 1179–1186. doi:10.1038/ng.3393
Shaban, H. A., Barth, R., and Bystricky., K. (2020). Navigating the Crowd: Visualizing Coordination between Genome Dynamics, Structure, and Transcription. Genome Biol. 21 (1), 278. doi:10.1186/s13059-020-02185-y
Sigalova, O. M., Shaeiri, A., Forneris, M., Furlong, E. E., and Zaugg, J. B. (2020). Predictive Features of Gene Expression Variation Reveal Mechanistic Link with Differential Expression. Mol. Syst. Biol. 16 (8), e9539. doi:10.15252/msb.20209539
Sonawane, A. R., Weiss, S. T., Glass, K., and Sharma, A. (2019). Network Medicine in the Age of Biomedical Big Data. Front. Genet. 10, 294. doi:10.3389/fgene.2019.00294
Suárez, L. E., Markello, R. D., Betzel, R. F., and Misic., B. (2020). Linking Structure and Function in Macroscale Brain Networks. Trends Cogn. Sci. 24 (4), 302–315. doi:10.1016/j.tics.2020.01.008
Sun, L., Jing, Y., Liu, X., Li, Q., Xue, Z., Cheng, Z., et al. (2020). Heat Stress-Induced Transposon Activation Correlates with 3D Chromatin Organization Rearrangement in Arabidopsis. Nat. Commun. 11 (1), 1886. doi:10.1038/s41467-020-15809-5
Takei, Y., Yun, J., Zheng, S., Ollikainen, N., Pierson, N., White, J., et al. (2021). Integrated Spatial Genomics Reveals Global Architecture of Single Nuclei. Nature 590 (7845), 344–350. doi:10.1038/s41586-020-03126-2
Thibodeau, A., Márquez, E. J., Luo, O., Ruan, Y., Menghi, F., Shin, D. G., et al. (2016). QuIN: A Web Server for Querying and Visualizing Chromatin Interaction Networks. Plos Comput. Biol. 12 (6), e1004809. doi:10.1371/journal.pcbi.1004809
Thibodeau, A., Márquez, E. J., Shin, D. G., Vera-Licona, P., and Ucar, D. (2017). Chromatin Interaction Networks Revealed Unique Connectivity Patterns of Broad H3K4me3 Domains and Super Enhancers in 3D Chromatin. Sci. Rep. 7 (1), 14466. doi:10.1038/s41598-017-14389-7
Tirosh, I., Barkai, N., and Verstrepen, K. J. (2009). Promoter Architecture and the Evolvability of Gene Expression. J. Biol. 8 (11), 95. doi:10.1186/jbiol204
Vilarrasa-Blasi, R., Soler-Vila, P., Verdaguer-Dot, N., Russiñol, N., Di Stefano, M., Chapaprieta, V., et al. (2021). Dynamics of Genome Architecture and Chromatin Function during Human B Cell Differentiation and Neoplastic Transformation. Nat. Commun. 12 (1), 651. doi:10.1038/s41467-020-20849-y
Winick-Ng, W., Alexander, K., Harabula, I., Zea Redondo, L., Meijer, M., Serebreni, L., et al. (2020). Cell-Type Specialization in the Brain Is Encoded by Specific Long-Range Chromatin Topologies. BioRxiv [Preprint]. doi:10.1101/2020.04.02.020990
Keywords: chromatin networks, nucleome, epigenomics, variability, structure-function, complex networks, Hi-C
Citation: Pancaldi V (2021) Chromatin Network Analyses: Towards Structure-Function Relationships in Epigenomics. Front. Bioinform. 1:742216. doi: 10.3389/fbinf.2021.742216
Received: 15 July 2021; Accepted: 04 October 2021;
Published: 27 October 2021.
Edited by:
Baldomero Oliva, Pompeu Fabra University, SpainReviewed by:
Arda Halu, Brigham and Women’s Hospital and Harvard Medical School, United StatesYu Xia, McGill University, Canada
Copyright © 2021 Pancaldi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vera Pancaldi, vera.pancaldi@inserm.fr