Immune Regulation in Time and Space: The Role of Local- and Long-Range Genomic Interactions in Regulating Immune Responses
- 1Division of Chemical, Systems, and Synthetic Biology, Department of Integrative Biomedical Sciences, Institute of Infectious Disease & Molecular Medicine, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
- 2Radboud Institute for Molecular Life Sciences (RIMLS), Radboud University Medical Center, Nijmegen, Netherlands
- 3Epigenomics & Single Cell Biophysics Group, Department of Cell Biology, Radboud University, Nijmegen, Netherlands
- 4Department of Human Genetics, Radboud University Medical Center, Nijmegen, Netherlands
Mammals face and overcome an onslaught of endogenous and exogenous challenges in order to survive. Typical immune cells and barrier cells, such as epithelia, must respond rapidly and effectively to encountered pathogens and aberrant cells to prevent invasion and eliminate pathogenic species before they become overgrown and cause harm. On the other hand, inappropriate initiation and failed termination of immune cell effector function in the absence of pathogens or aberrant tissue gives rise to a number of chronic, auto-immune, and neoplastic diseases. Therefore, the fine control of immune effector functions to provide for a rapid, robust response to challenge is essential. Importantly, immune cells are heterogeneous due to various factors relating to cytokine exposure and cell-cell interaction. For instance, tissue-resident macrophages and T cells are phenotypically, transcriptionally, and functionally distinct from their circulating counterparts. Indeed, even the same cell types in the same environment show distinct transcription patterns at the single cell level due to cellular noise, despite being robust in concert. Additionally, immune cells must remain quiescent in a naive state to avoid autoimmunity or chronic inflammatory states but must respond robustly upon activation regardless of their microenvironment or cellular noise. In recent years, accruing evidence from next-generation sequencing, chromatin capture techniques, and high-resolution imaging has shown that local- and long-range genome architecture plays an important role in coordinating rapid and robust transcriptional responses. Here, we discuss the local- and long-range genome architecture of immune cells and the resultant changes upon pathogen or antigen exposure. Furthermore, we argue that genome structures contribute functionally to rapid and robust responses under noisy and distinct cellular environments and propose a model to explain this phenomenon.
Survival requires that organisms withstand the constant assault of external and internal threats. The most fundamental of these is coexistence with potential pathogens with the means to invade and cause harm in myriad ways. Consequently, mammals have evolved a complex network of immune responses. In addition to monitoring and responding to incursions, surveilling the patterns of normal tissue development and function is an essential immune process. Given the unrelenting burden this places on the organism, it follows that the immune cells responsible have developed, what is in many ways, a unique phenotype and regulatory state. Namely, immune cells must recognize and respond to infection immediately or run the risk of being overwhelmed by any number of invasive microbial species. The spectrum of potential pathogens is vast, ever-changing, and their virulence factors and means of invasion and evasion almost boundless. A diverse pool of dangers necessitates an equally diverse responsive potential. Fostering phenotypic and functional diversity is thus an indispensable feature of immune cells. This versatility includes a number of requirements: firstly, discretion - overlooking harmless features of internal and external environments thus preventing (or attenuating) unnecessary immune responses are essential in avoiding immune-mediated conditions. Secondly, immune cells must be poised for a rapid response replete with the functional specificity required for resolution. Here then, we have two salient features of the immune response: heterogeneity, and inductive capacity which is programmable. Without the successful synergy of these factors homeostasis is lost: the organism becomes victim to invasive pathogens or, alternatively, fails to contain internal dysregulation. These features are, therefore, crucial to the functioning of immune cells.
The function of immune cells is regulated through gene expression of immune genes. Gene expression patterns for inducible genes are contingent on developmental stage and environmental change. The process is complex and variegated; comprising multiple individual sequences usually in a predetermined order: Chromatin remodeling, transcription, post-transcriptional modification, translation, post-translational modification, and transportation. Regulation can thus occur at every step. Nucleosome remodeling at promoter regions prior to transcription, for example, is well-described (1). Nucleosomes at transcription start sites are depleted allowing RNA polymerase engagement by ATP-dependent chromatin remodeling agents, such as the SWI/SNF complex (1). The second step is the recruitment of general transcription factors and formation of the pre-initiation complex (PIC) (2). RNA polymerase II (RNAPII) within the PIC initiates transcription by synthesizing short RNAs, approximately 25-60nt in length, before the transcription complex pauses, awaiting further positive signaling. At this stage, the negative elongation factor (NELF) binds RNAPII together with DRB Sensitivity Inducing Factor to temporarily inhibit RNAPII reading (2). This pause is followed by elongation. P-TEFb is recruited to the PIC and, in turn, phosphorylates serine residues at the second and fifth amino acid from N-terminal; phosphorylation dissociates NELF allowing elongation to begin and the full-length transcript is subsequently synthesized (3, 4).
These are the minimum requirements for transcription; without additional potentiating factors baseline expression remains low. Synthesizing the number of transcripts of immune genes required following immune challenge necessitates the accruement of other activating elements to induce efficient transcription. Specific transcription factors (TFs) are one such example. Key to immune cells specifically is NFkB - known to induce most inflammatory genes while suppressing a cohort of metabolic genes. In the nucleus, NFkB regulates gene expression by recruitment of P-TEFb. Interestingly, NFkB also has both potentiating and inhibitory effects on gene expression through the variable recruitment of histone acetyltransferases and histone deacetylases (5–7).
The nucleosome is a central feature of DNA arrangement allowing for both compaction and functional control of genes. A strand of DNA makes 1.68 turns around a histone protein octamer, forming a fiber-like structure 11 nm in width. Post-translational modifications of these histones are well-entrenched regulators of this packing process. A variety of local histone modifications are associated with actively transcribed promoters: H3K4me3 recruits RNAPII to the promoter while H3K27ac recruits BRD4 which in turn recruits elongation factor P-TEFb (8–10). Additionally, p300, a histone acetyltransferase which uses acetyl-CoA as a substrate, performs histone crotonylation in macrophages. Although the mechanism remains unknown, crotonylation is associated with a greater degree of gene expression potentiation compared to acetylation (11).
In addition to local genomic structures, distal elements such as enhancers are indispensable in supporting efficient gene expression. Active enhancers are identified by association of H3K4me1 and H3K27ac and a specialized family of RNAs, enhancer RNAs (eRNAs) which are transcribed from these sites (12). Although the mechanism of enhancer action - for example, the necessity of eRNAs - remains controversial, it is believed that enhancers recruit RNAPII to promoters. Orthogonal evidence that RNAPII is present at high levels at active enhancers and that eRNA expression precedes transcription of protein-coding genes supports this theory (13–16). Dynamic local protein milieus, variable histone modifications, and non-coding elements must all interface in time and space to ensure efficient, inducible transcription. This is enabled by refined chromatin structures.
While it has long been understood that the chromatin architecture is purposeful, the role of this reproducible, predictable structure extends beyond the compaction of DNA. Protein complexes and histone modifications at close-range and various forms of chromatin loops at long-range are all argued to play roles in gene regulation. Having established the core functioning of immune gene regulation, it follows that if chromatin structure is to be shown to influence such regulation it must facilitate robust induction and heterogeneity. Accruing evidence demonstrates that features of chromatin structure at a range of scales are central to effective immune regulation. This includes creating permissive environments for rapid gene induction and modifying chromatin structure to encode epigenetic memory and prompt cell activation. An interface exists between predictable processes and random events which govern chromatin structure creating a potential space for heterogeneity and inducibility. This shows clearly that chromatin structure does parallel the needs of immune regulation, and, therefore, provides another avenue of regulatory leverage which could satisfy the needs of the immune response.
Intrinsic to this tension is the potential for things to go awry. While stochasticity hemmed-in by well-regulated processes creates an environment of diversity and flexibility, hyper-activation or irregular gene expression patterns are a risk. Auto-immunity is a key example of this: misappropriated control mechanisms of the gene regulatory space give rise to seemingly irreversible shifts in expression patterns which are pervasive.
Immune cells are present in virtually all tissues, and both tissue-resident and circulating immune cells play an important role in the immune response. However, these two populations of cells differ in origin, function, and transcriptional program. While circulating macrophages, for example, have a short life-time (up to 1 week) and are continuously replenished from hematopoietic stem cells in bone marrow; microglia - brain-resident macrophages - are derived from the embryonic yolk sac and are replenished locally due to an intrinsic self-renewal capability (17). Due to this, microglia have an extremely long lifespan. It follows that phenomena associated with epigenetic reprogramming, such trained immunity and cell exhaustion, influence the function of these cells for the duration of their lifespan. In general, tissue-resident macrophages display an M2-like phenotype and play an important role in tissue homeostasis (18). For instance; microglia contribute to neuron pruning and neuro-development by expressing the CX3CR1 gene. CX3CR1 knockout mice demonstrate reduced dendritic spine pruning, abnormal synapse maturation, decreased functional connectivity and display behavioral abnormalities (19). As tissue-resident immune cells have unique functions and expression profiles, it is unsurprising that they have unique regulation mechanisms and chromatin structure. PU.1 is a central lineage-determining transcription factor (LDTF) for microglia as well as circulating macrophages. In contrast, IRF8 is a key TF for microglia, but not circulating macrophages (20). In addition to the unique expression of key TFs, nucleosome position and histone modifications in tissue-resident macrophages reflect the unique demands of the surrounding tissue environment. Interestingly, transplantation of macrophages to other tissue sites results in the reprogramming of nucleosome position; suggesting that nucleosome position is determined by local microenvironments likely through cell-cell communication between macrophages and tissue-specific cells (21). Variation in long-range chromatin interactions between tissue-resident and circulating immune cells is largely unknown due to technical limitations of exploring this in vivo. Understanding these differences is essential to developmental epigenetics. For example, both microglia and circulating macrophages are considered macrophages, sharing numerous features, although their origin and developmental pathways are different.
In this review, we will discuss the contribution of chromatin structure to immune regulation by exploring how it enables inductive capacity and heterogeneity. Interrogation at single-cell resolution and improvements in chromatin conformation capture techniques posit an array of mechanisms by which chromatin structures regulate immune genes. Through this, we will address two overarching questions: first, which elements of chromatin architecture regulate immune responses and how this is achieved. And, second, how rapid and versatile immune responses facilitated by chromatin structure underlie immune-mediated conditions, like auto-immunity.
Chromatin Structure as a Property and a Driver of Transcriptional Regulation
Discovery of the staggeringly complex nature of chromatin organization has been fuelled - and limited - by approaches capable of reliably identifying components of 3-dimensional nuclear architecture. Chromatin conformation capture techniques, notably Hi-C and its derivatives, have begun to clarify this area of cellular biology and, in doing so, uncover the role of nuclear organization in controlling gene expression - an area in which there remains active debate. This control takes place at an impressive range of scale and by manifold, seemingly independent, mechanisms. We argue that this regulatory platform created by complex chromatin structures creates a dimension of control which is essential for immune functionality.
The Global Organization of Chromatin
How it is that 2 meters of DNA is folded within a nucleus less than 20 um in diameter is a long-standing question of human biology. Despite being dynamically altered depending on biological events, such as cell-cycle progression and differentiation, genome structure is not random. Stretches of DNA fold in sophisticated arrangements without entanglement. The nucleosome is the most basic scaffolding agent, as described above, with both structural and functional relevance. At a global level, however, the arrangement of nuclear content is becoming less opaque. The nucleus at interphase is divided into compartments comprising gene-rich active regions and gene-scarce inactive regions, termed A and B compartments respectively (22) (Figure 1A). Characterization of A/B compartments has led to the promulgation of liquid-liquid phase separation (LLPS) as the governing force of these large distinct zones which, canonically, represent the largest unit of chromatin association and organization with functional relevance. As technologies have matured, numerous groups have described so-called higher orders of chromatin arrangement. Topologically Associating Domains (TADs) constitute large, looping territories which are stable across cell types (23–25) (Figure 1B). The fundamental difference between chromatin compartments and TADs is size; while compartments are multi-megabase in scale, TADs are generally regarded as being between 0.2 and 0.8 megabases in size. In addition to size, chromatin compartments comprise interacting genomic regions on the same or different chromosomes which are in the same transcriptional state - i.e. either expressed or repressed genes. TADs, on the other hand, interact preferentially with themselves and insulate regions of the genome together with (or apart from) neighboring regulatory elements and functionally associated genes (26). An important further distinction is their conservation across cell-type: TAD boundaries are almost always the same [although debate arises here as well; reviewed by Eres and Gilad (27)], whereas A/B compartments are highly variable.
Figure 1 Diverse chromatin structures across the range of nuclear scale and the main proposed regulatory mechanisms which underlie chromatin structure. (A) CTCF-mediated looping forms TADs. Directional CTCF motifs across the genome are occupied by CTCF in relation to loading/unloading dynamics of nuclear proteins. Similarly, Cohesin binds chromatin (top panel) and extrudes chromatin through its closed ring structure to form a loop (middle panels). Extrusion continues bidirectionally until convergent CTCF occupied sites are encountered. Heterodimerization of CTCFs, Cohesin, and various other proteins form a stable loop anchor complex (bottom panel) and are the basis of most stable TAD structures. (B) RNAPII-mediated looping forms small chromatin loops. TFs bind motifs in the regulatory elements of genes along with RNAPII (top-middle panel). Passive diffusion and random motion which are not ATP-mediated bring these protein structures into proximity with one-another (middle-bottom panel) and, by means of affinity-interactions, chromatin loops form between regulatory elements (bottom panel). These protein-protein interactions may result in liquid-liquid phase separation and the relatively dynamic recruitment of transcriptional machinery, creating a permissive environment for rapid and robust transcription. (C) Nuclear compartments are complex chromatin structures. In interphase, the nucleus is divided into A and B compartments. These domains span multiple megabases and comprise active (gene-dense, euchromatin) and inactive (gene-scarce, heterochromatin) possibly mediated by liquid-liquid phase separation. Green - active or A compartment; red - inactive or B compartment. (D) Multiple different chromatin loop structures are contained within single chromosomes which contribute to A or B compartments. Individual chromosomes can contain contiguous regions of active or inactive chromatin which represent the interactions of various elements of the genome facilitating functional regulation. These chromatin interactions are depicted by a schematic of a Hi-C map demonstrating a small area of a single chromosome in interphase. Triangular regions of varying red intensity represent the likelihood of interactions between chromatin regions in a population of cells. TAD with boundaries occupied by convergent CTCF motifs insulate areas of the genome (left). Similarly, the active compartment contribution of this representative chromosome may house RNAPII-mediate loops mediating interactions between genes, or genes and non-coding regulatory elements (as described in B) (middle). The inactive compartment contribution may be mediated by protein-mediated liquid-liquid phase separation to form repressive compartmental domains in gene-scarce heterochromatin (right). These processes can be either ATP- or non-ATP-mediated.
Small, Variable Domains of Chromatin Structure Facilitate Heterogeneity and Inducibility
The essential component of heterogeneity for immune function suggests is cell-to-cell variation in gene regulation and, similarly, in chromatin structure. Recent single-cell investigations have shown this heterogeneity unequivocally. The knowledge that chromatin moves dynamically and forms heterogeneous structures is not new. Structural changes of the chromosomes throughout the cell cycle, passive diffusion of chromatin at interphase, and the structural differences of the X chromosome have been explored for decades (28, 29). However, this work has been limited - focusing only on specific loci or operating at poor resolution. Recent single-cell analyses have demonstrated genome-wide heterogeneity with chromatin structures that respond to cell need. Single cell sequencing has convincingly shown that nucleosome positioning is not consistent between cells in a population (30–32). Heterochromatin domains comprising phase-separated droplets are now known to be more dynamic than previously assumed (33, 34). Chromatin interactions are similarly variegated: the frequency of DNA loop formation is approximately 30% of alleles at maximum depending on the locus and RNAPII-mediated loops are heterogeneous (35–37). Likewise, there is diversity between TADs regarded to be conserved across cell-type and species (35, 38–40). Evidently, chromatin structure is dynamic and highly heterogeneous across the range of scale. These findings raise important questions: why is chromatin structure so variable and how is this achieved?
Consistent improvement in technologies which reliably interrogate chromatin interaction networks, the recognition that there appear to be smaller sub-TAD regions within largely static TADs, and the immensity and apparent stochasticity of A/B compartments have all strengthened the case for smaller, more variable components of chromatin structure. These components could reasonably contribute both towards gene expression regulation and DNA compaction. The source of heterogenous genome structures may be derived from the dynamic nature of the genome throughout the cell cycle. Unsurprisingly, distinct cell-cycle stages reveal distinct chromatin structures. However, even within the same cell cycle stage, individual cells harbor some degree of distinct chromatin structure variability. Two interacting components govern this variability. Firstly, ATP-dependent directional events which are unaffected by thermodynamic fluctuations, such as nucleosome depletion by chromatin remodeling complexes and partitioning by actin (41, 42). Secondly, ATP-independent non-directional events which include nucleosome deposition by histone chaperone, the binding of proteins to DNA, and passive diffusion of chromatin (28, 43, 44).
Many long-range DNA loops are protein-mediated with a number of important proteins being implicated in their formation. This mechanism represents a confluence of directional and non-directional events. The most extensively studied DNA-loops are CTCF-mediated loops. Generally housed within TAD boundaries, these loops are appreciated by punctate signals at the summits of associating domains in Hi-C maps representing the close association of structural proteins (26, 45). The loop extrusion model is posited as the dominant theory of how these domains are formed describing their active formation by SMC proteins, especially cohesin complexes (Figure 1A) (46–49). In brief, these proteins dock and travel divergently along sequences of DNA, extruding chromatin through their looped structure, until they reach TAD boundary - occupied CTCF sites which are orientated in the correct manner. Thereafter, a stable loop ‘base’ is formed through, amongst other processes, the heterodimerization of CTCF and Cohesin rings which secures the TAD. This is supported by convergent CTCF-binding sites being a prominent feature of TAD boundaries and it is now largely agreed that convergent CTCF-binding motifs are the predominant boundary element required for TAD formation (48, 49). This mechanism explains the non-overlapping property of TADs and the predictable location of aberrant TADs following CTCF deletion (45, 48).
There remains debate as to the intricacies of this mechanism and it is incompletely understood. For example, CTCF sites arranged in tandem likely form loops of a differing structure which could contribute to the topology of a domain. Additionally, the distribution of genes within CTCF-loops is seemingly determined by their being either constitutive or cell-type specific with the implication that genes localized within the loop must have a means of induction when needed (50). There also appears to be a distinct role for Cohesin complexes and CTCF; perturbation of the functional pathway of either results in different, and sometimes contradictory, effects on chromatin structure (26). However, depleting CTCF or attenuating Cohesin loading result in the near-complete loss of TADs and CTCF-loops (51, 52). Surprisingly, this enhances compartment strength without having a major effect on baseline gene expression (52). Nor does cohesin depletion affect histone marks while a small number of active genes are repressed (53). Thus, an important discovery is that the action of Cohesin, and therefore the role of CTCF-loops, may be implicated in responses to stimuli and lineage commitment as opposed to gene expression at steady state. This is plausible as a number of enhancers and the genes they induce are co-localized by CTCF-looping. Indeed, Cohesin was shown to be necessary for the induction of responsive genes and enhancers in myeloid cells and macrophages exposed to inflammatory stimuli (54).
RNAPII is another example of protein-mediated DNA looping, specifically between promoters and enhancers (55). Unlike CTCF-loops extruded by ATP-dependent means, RNAPII facilitates interactions by affinity and, therefore, does not irreversible characteristic of CTCF-loops. Instead, RNAPII-mediated loops are driven by cell activity. RNAPII, by facilitating interactions between genes and enhancers within a CTCF-loop and an active enhancer or promoter at a loop anchor, may drive the expression of tissue-specific inducible genes in a coordinated manner with house-keeping genes (Figure 1B) (50). Proteins within the nucleus, therefore, have an established role in the formation of chromatin structure. The interplay of ATP-dependent mechanisms (CTCF-loops) and affinity-driven mechanisms governed by protein kinetics (RNAPII-loops) creates potential for inducible gene expression and rapid responses. Furthermore, proteins create a distinction in location between active cell-type specific genes, constitutively active non-specific genes, and active enhancers suggesting that chromatin structures at this level may enable versatile responses.
Heterogeneity thus relates to processes of chromatin formation and movement. However, stochasticity is inevitable and influences even relatively stable directional processes. For example; the loop extrusion process by SMC proteins itself is not influenced by thermodynamic instability, however, the timing and start location may be stochastic owing to cohesin docking. The end location is, likewise, stochastic: cohesin continues extruding DNA until blocked by CTCF and released by cohesin-releasing factors. TFs and CTCF - which hamper cohesin’s advancement - repeatedly bind and unbind DNA with CTCF remaining mostly in the unbound state (56, 57). Therefore, cohesin may advance past CTCF-binding motifs resulting in heterogeneous loop lengths and variable boundaries. Indeed, the average size of CTCF-loops increases following CTCF deletion (58–60).
There are further components to the non-directional forces driving heterogeneity within chromatin structure. Just as TADs comprise a combination of small protein-mediated loops which scaffold the nucleus through directional events, so A/B compartments are likewise a composite of smaller domains (Figure 1C). However, the formation of these DNA loops is stochastic. Compartmental domains of self-associating chromatin form by non-directional forces, frequently containing a few (or even a single) gene loci(us) and the associated regulatory sequences (61). The most frequent interactions are those between transcriptional start- and end-sites of transcribed genes (61). This drives an accumulation of transcriptional machinery and TFs and there appears to be a clear link between transcription and compartmental domain formation. In-line with this, inter-compartmental domain associations appear stochastic and their strength and frequency reliant on the concentration and type of associated proteins in the looped region (33). Therefore, the affinity and the interaction of these protein complexes will determine the cooperativity of chromatin regions. Compartmental domains comprise either active or inactive chromatin replete with differing sets of multivalent protein complexes could interact preferentially with one another in a manner representative of A/B chromatin compartments (26). Indeed, phase separation owing to multivalent TF cooperativity could occur (62, 63). Supporting this model, phase separated droplets of chromatin may be formed by protein complexes with intrinsically disordered domains including RNAPII and TFs (64, 65). Notably; hyper-phosphorylated - active RNAPII - has been shown to form droplets (66). It would follow that this process would form domains of active chromatin. On the other hand, heterochromatin is occupied by a starkly different complement of histone modifications. For example, H3K9me3 is recognized by the protein HP1, recruiting the specific histone methyltransferase for H3K9 (SUV39H1) via the DNA methyltransferase DNMT3A and the DNA methylation binding protein MECP2 (67). As SUV39H, in turn, transfers methylation group to its neighboring histones H3K9me3 marks are propagated and form tightly packed, functionally silenced heterochromatin which is likely to associate with other such regions (Figure 1D). While CTCF-loops adhere to TAD boundaries, compartmental domains appear to be less restrained. CTCF-loops appear to insulate compartmental domains of the same or different compartments therefore modulating interactions among compartmental domains at close range by a directional force (26). Long-range contacts between compartmental domains, however, are under the influence of thermodynamic principles.
A balance thus exists between two independent organizational events which suggests a complex and reciprocal relationship between transcription and genome architecture (52, 68). The non-directional formation of compartments adds a degree of complexity when in concert with ATP-dependent, protein-mediated mechanisms. The accumulation of proteins required for transcription in space and time creates a permissive environment for rapid and robust transcription. Additionally, the formation of phase separated structures requires constant fusion and fission of droplets creating a dynamic chromatin structure capable of versatility.
Inducible Immune Responses: Rapid and Refined Gene Expression Encoded by Genome Structure
Rapid Response Genes Are Prepared for Robust Action by Nucleosome Modification in Innate Immune Cells
A functioning repertoire of innate immune cells requires that they are rapidly and efficiently activated to enact effector functions, such as responding to invading pathogens. This is reliant on a change in epigenome to facilitate a transcriptional shift. Nucleosome remodeling functions as a gatekeeper in this process: determining which enhancers and promoters become nucleosome-free and regulate immune gene transcription. Nucleosome position is determined both by cellular signaling and developmental programming (69, 70). Stimulation by external signals induces TF cascades resulting in remodeling; these are termed signal-dependent transcription factors; SDTFs. On the other hand, developmental regulation of nucleosome position is informed by LTDFs (71, 72). PU.1 (SPI1) is a known LDTF in both innate and adaptive immune cells - T cells and macrophages (73). PU.1 is also a pioneer factor able to bind DNA in spite of nucleosome presence at binding sites and recruits chromatin remodeling complexes (74). However, PU.1 is not sufficient to induce transcription, rather, it binds cooperatively with other TFs such as AP-1 and members of the IRF family, which are essential in inducing transcription (75, 76). In this way, PU.1 functions as a safeguard to prevent unwanted transcription in the absence of other essential TFs (77).
To induce signal-dependent immune genes, SDTFs, such as NFkB, play a key role. NFkB is expressed constitutively, forming a complex with the repressor protein IκB and is present in the cytosol in this form (78). Both immune and barrier cells, such as epithelia, recognize pathogens through cell-surface receptors such as Toll-like Receptors (TLRs) and C-type Lectins. This recognition triggers the phosphorylation and subsequent degradation of IκB, prompting NFkB nuclear translocation after which it occupies proximal promoters on immune genes. In turn, transcription is induced by the recruitment of co-activators and P-TEFb (6, 79).
Transcription of rapid response immune genes is detectable within 15 minutes following pathogen recognition (80). Despite the first step in transcription induction for most immune genes being nucleosome depletion at promoters, the promoters of numerous immune genes are depleted in resting, or naïve, states. Indeed, it has been shown that knockdown of the SWI/SNF complex does not negatively influence rapid response gene expression, despite other genes being suppressed (81). Even more striking is that not only is nucleosome depletion a feature of naïve immune gene promoters; pre-loading of RNAPII has also been observed at these sites (82). This is an important observation as nucleosome depletion may take minutes (83). Indeed, it has been shown to take ~30 min to achieve nucleosome depletion at immune genes promoters (84). Some active histone marks such as H3K4me3 are also found at rapid response gene promoters at naïve state (85). Elongation marks such as histone acetylation and P-TEFb are not found at naïve state, but P-TEFb is in turn recruited by the bromodomain-containing protein Brd4, which detects H4K5/8/12Ac inducibly acquired at rapid response gene promoters (86). Given that first response genes become detectable within 15 minutes, the benefits of removing the nucleosome prior to immune challenge is apparent. This is demonstrated in Innate Lymphoid Cell (ILCs) (87). ILCs, including NK cells, which develop from lymphoid progenitors, express cytokine genes similar to innate immune cells. Of the 3 subtypes of ILCs, nucleosomes at regulatory elements of IFNg are only depleted in group 1 cells. This depletion is regulated developmentally and removed in a stepwise manner. In addition to nucleosome depletion, p300 is present at rapid response gene promoters in the naïve state, however, p300 shows stronger binding after stimulation. Although nucleosome depletion at rapid response genes is developmentally regulated, it remains unknown how nucleosomes are depleted from rapid response genes only. Rapid response gene enhancers are identified by H3K4me1 enrichment in the naïve state (88). To prevent transcription in this state, repressive histone marks, such as H3K9me3, H3K27me, and H4K20me3 are deposited at the enhancer together with the co-repressor recruitment complex Bcl6 (89, 90). In another study, transcriptional responses to TLR4 ligation primarily relied on pre-existing H3K4me2 enriched enhancers while de novo enhancers are regulated in a PU.1- and NFkB-dependent manner (91).
Contrary to previous dogma suggesting that innate immune cells lack memory capability, histone modifications, such as the addition of H3K4me3 at gene promoters, have been shown to infer long-term memory to past infectious challenges (92). This phenomenon, referred to as trained immunity, increases the speed and efficiency of rapid response gene transcription. The most prominent example being circulating macrophages. Following detection of BCG or beta-glucan exposure through C-type Lectin Dec1, their transcriptional profiles are durably altered; enabling robust responses to future immune challenge. Simultaneously, their metabolic profiles are altered. Upregulation of mTOR/HIF1a by Dec1 signaling activates glycolysis, glutaminolysis, and cholesterol synthesis pathways while oxygen consumption is reduced (93). This shift towards an aerobic glycolytic predominance is observed even under high-oxygen conditions. The role of this is the synthesis of intermediate metabolites, rather than maximizing energy production. Fumarate, an intermediate metabolite of the TCA cycle, is known to inhibit the histone demethylation activity of KDM5 (94). Due to the higher concentration of fumarate following mTOR/HIF1a upregulation, demethylation of H3K4me3 at some inflammatory gene promoters, such as IL6 and TNFa, is inhibited. This maintains H3K4me3 at these promoters even after stimuli are removed. Additionally, alpha-ketoglutarate is a cofactor of KDM5 and thus the balance between fumarate and alpha-Ketoglutaric acid influences KDM5 activity (94). Likewise, glutamine synthesized by glutaminolysis is converted to citric acid and its influx into the TCA cycle contributes to the accumulation of fumarate. Cholesterol synthesis is another pathway implicated in enhancing trained immunity; mevalonate, a cholesterol intermediate, contributes to the accumulation of H3K27ac; the precise mechanism of how this occurs remains unclear (95).
Strikingly, this epigenetic reprogramming can occur at the hematopoietic stem cell level, giving rise to trained macrophages which persist for at least 3 months (and perhaps even decades) (96, 97). Conversely, the addition of repressive chromatin marks leads to a slower, attenuated immune response. Chronically activated immune cells lose their ability to produce cytokines, proliferate, and lyse pathogens. The inflammatory gene promoters in these cells are not acetylated by the 2nd challenge (98). These phenotypic features revert to the naïve state following exposure to beta-glucan (99). Therefore, the speed and success of the immune response is fine-tuned by the dynamic modification of histone marks.
Enhancer-Promoter Interactions Are Pre-Formed at Rapid Response Genes in The Steady State in Both Innate Cells and T Cells
Successful rapid responses are enabled by long-range genomic interaction as well as local structures, such as nucleosomes. Enhancers-promoter loops at rapid response genes are pre-established in the naïve state (100). These enhancer-promoter interactions are mediated by both sequence-specific binding proteins, such as CTCF, and non-sequence-specific proteins, such as RNAPII. However, it remains unknown how the loops themselves are established. Interestingly, while immune cell activation leads to nucleosome remodeling at non-rapid response genes, most of the enhancer-promoter interactions at these genes are pre-established. Promoter capture Hi-C (PCHi-C) experiments have shown that activated macrophages (both M1 and M2 polarized) display very similar enhancer-promoter interaction patterns to naïve macrophages (M0 macrophages) (101). Despite this apparent stability of long-range interactions, enhancer-promoter loops are known to be changed by activation. In THP-1 cells, for example, gained enhancer-promoter interactions have been identified by PCHi-C (102). The reasons for the discrepancy in findings between studies which suggest long-term stability being a feature of enhancer-promoter loops and those providing evidence of dynamic change associated with the functional state of cells is still a topic of debate. One possible explanation relates to the dynamics of the genomic structure. PCHi-C experiments in macrophages were performed at 24 hours after activation as opposed to at 4 hours in THP-1 cells. If observed structural features of the genome are emergent properties of a self-organizing system, genome structure would be perturbed temporally by stimulation before returning to stable steady state within 24 hours (103). Despite difficulties in ascribing strict rules to enhancer-promoter interactions, it is clear that these patterns are highly cell-type specific and immune activation minimally alters these interactions, at least in the case of rapid response genes. Therefore, the strategy of pre-forming loops during steady or naïve states to prepare for transcriptional responses is not immune cell-specific. Rather, this is the general strategy employed by responsive genes and occurs in a signal-dependent manner. For example, glucocorticoid receptor response genes also form enhancer-promoter interactions at steady state (104).
Correlations between loop formation and cellular activity states does not suggest that long-range interactions regulate the timing of responsive expression. However, forming a loop before activation seems biologically reasonable. Cohesin needs nearly 40 minutes to form 900 Kb loops by extrusion (53). Considering that the average TAD is 1 Mb, reconstruction of the chromatin interaction pattern requires at least 40 minutes. Furthermore, various DNA-binding proteins present in the nucleus form large complexes before being functional. When cohesin encounters these DNA-binding proteins, loop extrusion takes place at a slower rate than naked DNA in vitro (105). If chromatin re-organization facilitates relatively slow cellular responses, such as cell-cycle progression or differentiation, this is acceptable. However, if rapidity is essential to the response in question, as in immune responses, it would be detrimental to dramatically re-wire enhancer-promoter interactions at every challenge to allow for rapid gene responses. Further work is needed to decipher the relationship between enhancer-promoter interactions and the control of expression timing. An approach might be to investigate the influence of depleting chromatin re-organizing factors on the timing of gene expression. For example; Batf is a pioneer factor which recruits CTCF via the Ets protein family (106). If enhancer-promoter interactions do control expression timing, the deletion of Batf would delay the onset of transcription at affected sites.
In summary, rapid response genes are prepared to initiate cellular function in a developmentally-programmed manner by removing nucleosomes, modifying histone tails, and forming enhancer-promoter loops (Figures 2A, B). Additionally, constitutive expression of NFkB in the naïve state appears to be an important factor in accelerating transcriptional initiation as it omits the need to rapidly transcribe the NFkB gene. However, there remain unanswered questions. LDTFs such PU.1 play an important role in depleting nucleosomes during developmental programming and yet not all PU.1 binding sites situate at rapid response genes. It is likely that other TFs which cooperatively bind PU.1, such as AP-1 and the IRF family, are important in selecting rapid response genes and, thus, control the timing of transcription. Likewise, it is unknown how sites for pre-formed loop formation are selected. Finally, not all inflammatory genes are trained by beta-glucan or BCG treatment and it remains unclear how specific genes are selected. Previously our group identified a novel class of lncRNAs called Immune Priming LncRNAs (IPLs). We showed that UMLILO, an IPL located upstream of IL8, is essential for training the IL8 gene (107). Non-coding transcripts, such as eRNAs, likely play an essential role in defining trainable genes.
Figure 2 Immune cells exhibit changes in chromatin structure which facilitate rapid, robust responses and facilitate long-term memory. (A) Rapid-response genes in innate immune cells are primed for rapid and robust responses in the naïve state. Permissive chromatin marks, acetylation and H3K4 trimethylation, at gene promoters, facilitate the pre-loading and stabilization of RNAPII at rapid-response gene promoters (left panel). Following challenge, P-TEFb phosphorylates RNAPII and active transcription takes place within minutes (right panel). (B) Early-response immune genes are primed following immune challenge and maintain epigenetic memory to facilitate rapid and robust responses in future. In the naïve state, early-response immune gene loci are inaccessible chromatin structure (left panel). Following immune challenge, pioneer TFs (PU.1) recruit the SWI/SNF complex (middle-left panel) which, in turn, histone-modifying enzymes and facilitates the recruitment of permissive histone marks and immune-specific TFs (NFkB) initiate the transcription cascade (middle-right panel). This occurs within 30 min and the epigenetic modification of these loci is maintained to facilitate rapid and robust responses in future (right panel). (C) Rapid-response genes in T cells are primed for rapid and robust responses in the naïve state. Similarly to innate immune cells, permissive chromatin marks, acetylation and H3K4 trimethylation, at gene promoters, facilitate the pre-loading and stabilization of RNAPII at rapid-response gene promoters (left panel). This allows for rapid transcription of these genes within minutes following immune challenge even in adaptive cohorts of cells (right panel). (D) T cell effector genes are epigenetically modified over time following immune challenge to facilitate long-term immune memory and more effective responses following future challenge. In the naïve state, effector gene loci are silenced by repressive chromatin marks (left panel). Following immune challenge, the recruitment of immune-specific TFs (NFAT) (middle-left panel) results in histone modification and the recruitment of RNAPII and subsequent expression of effector genes (middle-right panel). In memory states, these specific histone modifications are maintained (right panel).
Genomic Structure and Immune Activation in T Cells
Macrophages and T cells differ markedly in the activation pathways they follow. In the case of macrophages, interaction between pathogen components and pattern recognition receptors is sufficient to initiate effector responses, such as cytokine production. On the other hand, T cells require multiple steps for successful activation. Firstly, pathogen recognition by T cells relies on antigen presentation by antigen presenting cells - such as dendritic cells - followed by secondary stimulation. For example, CD4+ helper T cells use CD28 to recognize CD80 and CD86 on dendritic cells. Once activated, T cells enter a proliferative state undergoing clonal expansion and differentiate into various effector cell repertoires depending on the local cytokine milieu. Th1 and Th2 are the predominant subtypes of activated CD4+ helper T cells and play important roles in host defense and humoral immunity, similar to M1 and M2 polarized macrophages, respectively. In the case of CD8+ cytotoxic T cells, the majority of activated CD8+ cells become short-lived effector cells. A subset of CD4+ and CD8+ T cells become long-lived memory T cells enabling a rapid and specific adaptive response on subsequent challenge. The mechanism of memory T cell development remains controversial (108). Owing to the multiple steps involved in activation, proliferation, and differentiation, a fully-fledged T cell response occurs in the order of days; unlike macrophages which are activated within hours. It follows that T cells have sufficient time and capacity (through their sophisticated differentiation programmes) to restructure their genome, as opposed to innate immune cells. Here, we discuss the alterations in genome structure occurring in T cells following immune challenge. Changes to genome architecture at sites encoding the variable T Cell Receptor (TCR) to allow for rearrangement are well-described and reviewed elsewhere. As such we will focus on non-TCR regions with a specific emphasis on cytokine genes.
T Cell Activation and Changes in Local Genome Structure
Numerous studies have shown chromatin accessibility to be extensively altered following T cell activation. One of the earliest investigations demonstrates this in the promoter region of Il2 - a growth factor influencing T cell proliferation, survival, and differentiation expressed in both CD4+ and CD8+ T cells at various times post-challenge (109). The Il2 promoter and neighboring enhancers are inaccessible in the naïve, pre-challenge state but are remodeled by SWI/SNF chromatin remodeling complex recruited by NFAT and AP1 after activation (110–113). The role of neighboring enhancers is highly cell type-specific with different enhancers being activated to induce il2 expression in CD4+ and CD8+ T cells (114). Despite this, genome-wide investigation showed that most chromatin accessibility gains are shared between activated CD4+ and CD8+ T cells despite their different effector functions (115). Notably, most chromatin remodeling is complete within 24h as differential peaks at 24h and 72h after activation as probed by ATAC-Seq are largely overlapping (115, 116).
However, in a similar manner to early response genes in innate immune cells, some genes are prepared for activation in the naïve state. Nucleosome depletion and RNAP2 loading along with some active histone marks can be observed at the Tnfa promoter in CD8+ T cells. In contrast, this dynamic change occurs only after activation at the Ifng promoter (117, 118). Interestingly, other early response genes such as Ccl3 and Il2ra are enriched with H3K4me2, but not H3K4me3 as might be expected (117). It remains unknown why not all early response gene promoters are enriched with H3K4me3. It is possible that this is essential for the precise regulation of transcriptional timing.
A subset of CD4+ and CD8+ T cells exhibit transcriptional memory to previous pathogenic encounters which engender rapid and robust transcriptional responses on re-challenge. These memory T cells do so not only by rearrangement of the genomic TCR region, but also epigenetic change at cytokine gene loci, similar to immune training in the case of innate immune cells. Additionally, enhancer regions of the memory-associated genes DUSP10, BCL6, and TNFSF10 become accessible after activation in T cells and these structural changes are stably maintained even after the resolution of an immune challenge (119). Likewise, memory T cell enhancers and promoters are in a poised state demonstrated by the accumulation of active histone marks (120). The mechanism of epigenetic change following activation in this cohort of cells remains unclear. It is possible that, similar to that which occurs in trained immunity, intracellular metabolites in activated CD8+ T cells increase flux through key pathways such as glycolysis and glutaminolysis (121). On the hand, chronic infection is known to lead to T cell exhaustion and this epigenetic landscape is distinct from that of memory T cells (122).
There is ample evidence demonstrating that the dynamic remodeling of chromatin structures and histone modifications are an important characteristic of T cell functionality. A wide range of genes necessary for activated T cells undergo such changes, notably exclusions are those implicated in rapid response such as Tnfa. As in the case of trained immunity in the context of the innate immune response, memory T cells prepare for future responses by remodeling chromatin at cytokine genes (Figures 2C, D).
T Cell Activation and Enhancer-Promoter Interactions
Some of the earliest and most seminal work on T cell genome structure details the interaction of Ifng and Il4, cytokine genes characteristic of Th1 and Th2 CD4+ subsets respectively (123–125). Ifng and Il4 loci interact in the nuclei of mouse naïve T cells despite being located on different chromosomes. Following T cell polarization, however, this interaction is lost and Ifng, for one, interacts preferentially with enhancers in Th1 cells (126). The interactions of Ifng and enhancer, at least some of them, are mediated by T-bet (127). Furthermore, genome-wide analysis of the interactome at Ifng and Il4 promoters revealed induced interactions with other Th1 and Th2 specific genes (128). Although the majority of the interactions are stable during activation, de novo interactions mediated by STAT proteins exist between key genes of the CD4+ response (128).
Further examples of interactions taking place in the genome of T cells are those between Il4, Il5 and Il13 gene loci. These cytokines are transcribed in Th2 cells but not naïve or Th1 cells. Under most circumstances and in many cell types these loci are in contact, fibroblasts being a proven example of this. Interactions with the locus control region (LCR) are also observed in T cells (129). Surprisingly, LCR-promoter interactions occur — although to a weaker degree — in naïve T cells and Th1 cells. This suggests that Th2-associated cytokine genes interact with enhancers to facilitate rapid induction in Th2 cells and thus the expression of Th2-associated cytokines are regulated by the subtype-specific activation of enhancers rather than the interaction (129–131).
Recent genome-wide chromatin interaction studies have identified interactions specific to T cell activation. Interestingly, and similar to their chromatin accessibility profiles, activated CD4+ and CD8+ T cells share the same chromatin interactions despite differences in effector function (115).
Cis-regulatory element interaction studies have intriguing results. One such PCHi-C study showed activated and naïve T cells with similar patterns of enhancer-promoter relationships.
On the other hand, H3K27ac HiChIP experiments have demonstrated that naïve T cells show distinct enhancer-promoter interactions compared to Th17 and regulatory T cells (132). One possible explanation for these discrepancies is methodological disparity. PCHi-C captures all interactions with promoters within the nucleus, including the loops which form to facilitate nuclear compaction, while HiChIP captures primarily active regions identified by H3K27ac pull-down. Therefore, interactions mediated by H3K27ac are underrepresented in PCHi-C data. A further explanation is cell type variety. The biology of T cells is complex; it is plausible that while activated Th1 and Th2 populations have similar genome interaction patterns to naïve T cells, Th17 and Treg cells may show entirely different functional patterns. Furthermore, owing to the fact that Th17 and Treg cells are relatively rare, their interaction patterns may have been diluted by a pooled PCHi-C approach as opposed to cell-sorting occurring in HiChIP approaches.
Heterogeneity in the Immune Response: Dynamic and Heterogeneous Genome Structures Influence the Transcription and Phenotype of Immune Cells
Expression patterns in immune response are known to be highly heterogeneous. Even in consistent genetic and environmental conditions - for example LPS stimulations in cultured cell lines - phenotype and transcription patterns are distinct at the single cell level following activation (133). In vivo, epigenetic modifications by past infection and pathogen variability add further layers of complexity to this heterogeneous response (134). The source of heterogeneity in transcription can be classified as extrinsic or intrinsic noise.
The Contributions of Genomic Structure in Robust Gene Expression
Intrinsic noise arises from the random nature of biochemical reactions in the processes of transcription and translation. Transcription is characterized by discontinuous bursts; the interval between bursts produces heterogeneity and cell-cell variability in transcription (135). To limit perturbation by intrinsic noise, cells insulate expression fluctuations by gene regulatory networks (136). For example, the expression of NFkB - a key TF in immune responses - is heterogeneous at steady-state (137). However, absolute expression levels of NFkB at steady state are not a determinant of activation following stimulation; even cells with a dearth of NFkB are activated efficiently. A feedback loop consisting of NFkB and P50 detects fold changes in NFkB concentration in the nucleus, inducing activation rather than absolute NFkB concentrations (137, 138).
There is relatively little cell-cell variation in the expression of genes involved in transcription, such as kinases and TFs, however, cytokine and chemokine expression varies markedly between cells (139). Some cytokine and chemokines are expressed binomially while others are expressed ubiquitously (140). Given the stochastic forces governing chromatin structure and transcription, it is not surprising that the expression of cytokines and chemokines following activation is highly variable. Rather, the intriguing question is how ubiquitously expressed genes, such as IL8, are stably transcribed from an heterogenous chromatin structure. Despite the frequency of enhancer-promoter interaction at maximum approximately 30%, IL8 is expressed in over 80% of innate immune cells following LPS stimulation (35, 133). One explanation is that direct interaction between the enhancer and promoter, which cross-linking based techniques such as ChIA-Drop and Hi-C, detect, is not essential. The majority of TFs have intrinsically disordered regions, enabling the formation of biological condensates (62). Condensates, namely LLPSs, are membrane-less organelles and allow rapid, reversible condensation of specific proteins and nucleic acid into discrete assemblies that dynamically exchange biomolecules with the surrounding environment. Therefore, enhancers and promoters located in the same condensate could interact and induce transcription without direct contact.
Another possible explanation is the presence of multiple enhancers. Frequently, a single gene interacts with multiple enhancers. This is particularly true of robustly expressing genes determining cell identity which interact with multiple enhancers and form so-called super-enhancers (141, 142). In this scenario, only a portion of available enhancers need interact with a promoter and induce transcription. For example, if 10 enhancers can interact with the same target gene at a probability of 30%, approximately 98% of cells will display an interaction with at least one enhancer. In immune responses, NFkB forms such a super-enhancer by exploiting pre-formed chromatin structures; multiple loops are thus important in robustly expressing NFkB (143, 144). Indeed, redundancy of enhancers is a pervasive feature of the mammalian genome and is essential for robustness of expression (145). Indeed, IL8 forms a super-enhancer in macrophages. Importantly, these two features are not mutually exclusive. The LLPS model suggests that condensates are formed at super-enhancers. Therefore, if at least one enhancer is located within a condensate at close range to the gene of interest gene induction can occur. In this situation, stable ATP-dependent loops formed by cohesin and CTCF could be key factors as deletion of CTCF binding sites increases cell-to-cell variation of gene expression (146).
A third possibility for this robustness of immune gene expression can be explained through the lens of dynamics. For example, transcription is known to confine chromatin movement (147, 148). If immune genes are transcribed, movement of chromatin is confined and this facilitates robust transcription. On the other hand, if no transcription takes place, promoters and enhancers diffuse randomly, increasingly their probability of interaction. Therefore, although there is delay, many cells could express immune genes.
Small-scale genome structures such as nucleosome position and histone modification have also found to be associated to heterogeneity in gene expression. Nucleosome position at promoter is varied, especially lowly expressed genes while it is less heterogenous in actively transcribed gene promoter (149). This robustness of nucleosome position is probably due to the stable nucleosome depletion by ATP-dependent chromatin remodeling complex. For histone modification, H3K4me3 signals positively correlate with high gene expression in immune cells and contribute heterogeneity in gene expression (150). It remains unknown how frequently cells (or allele) pose specific histone modification at regulatory elements of rapid response gene precisely and how histone modifications contribute robustness and heterogeneity in gene expression. Ideally, multi-omics data which detects gene expression and histone modification simultaneously should be applied to elucidate it. Although these techniques are recently developed, the inference of data is not straight-forward to fully understand relationship between histone modification and gene expression (151). A lot of histone modification have been identified and studied relationship with gene expression and they are inter-related mutually. For example, H3K9ac and H3K27ac are associated with active gene expression, but it remains unknown whether H3K9ac is always co-localized with H3K27ac, and whether H3K27ac only promoter shows distinct robustness in gene expression from H3K9ac/H3K27ac co-localized promoter. To understand contribution of histone modifications in gene expression robustness, it must consider gene expression probabilities under specific condition represented by histone modification pattern. Therefore, Baysian approach will be useful for inference of data.
The Influence of Genomic Structure on Cell Fate Following Activation
In the steady-state, populations of THP-1 are not transcriptionally distinct. However, following LPS challenge, M1 and M2-like phenotypes develop despite variation in expression levels due to intrinsic noise (133). If phenotypic decision after stimulation is explained solely by cell-cell variation in gene expression arise from intrinsic noise, the ratio of the M1 and M2-like phenotypes are always constant. However, in macrophages, nitric oxide (NO) adjust and limit overall inflammation intensity, influencing decisions on phenotype (152, 153), suggesting that extrinsic noise such as paracrine signaling influence on phenotypic decision making. Analogous to quorum sensing observed in bacteria, cells detect and respond to population density through cell-cell communication. For example, T cells use IL2 as an autoinducer that regulates population density through STAT5 signaling (154–156). In plants, NO affects chromatin structure, suggesting that alteration at steady-state by NO influences phenotypic decisions (157). Likewise, in many biological processes, such as in differentiation, alteration of chromatin structures precedes expression reprogramming (158–160). This implies that even in cells with the same transcriptional pattern, cell fate could be determined by differing chromatin structures. Therefore, it is possible that heterogeneous chromatin structures at steady-state influence phenotype, despite being apparently indistinguishable in gene expression. Of note, chromatin structure is a probabilistic modulator of gene expression as opposed to an absolute determinant. For example, enhancer-promoter interactions regulate the frequency of transcription bursts but do not eliminate stochasticity (161). This principle also applies to the nucleosome and not all genes with nucleosome-free promoter are transcribed (30). Therefore, specific chromatin conformations at steady-state may increase the probability of M1 (or M2) phenotypes following stimulation, but this is not deterministic (Figure 3).
Figure 3 A proposed model of how genome structure influences phenotypic consequence. (A) In the naïve state, all inflammatory cells are almost indistinguishable with regards to their transcriptional profiles and durable epigenetic changes have yet to take place. (B) Due to thermodynamic fluctuations, however, some cells become primed states transiently at various possible loci. Chromatin structural changes by means of external signals, for example NO as a paracrine signal, can favor a primed state with specific characteristics. (C) If cells encounter activation signals in this primed state, there is a high probability they will become a specific phenotype following activation. For example, pro-inflammatory primed cells become pro-inflammatory cells (and vice versa). On the other hand, owing to the intrinsic stochasticity of all transcription, primed cells can become the opposing phenotype (i.e. inflammatory or anti-inflammatory), although this probability is lower. Rectangular insets below cells diagrams is a representation of the chromatin changes which might occur at specific gene loci.
Due to technical limitations in study, the influence of immune cell heterogeneity in phenotype decision making remains unknown. Recently developed methods exploring chromatin structure and transcription simultaneously will hopefully reveal this relationship in immune responses in more in detail. Similarly, live-cell imaging techniques have the potential to demonstrate the link between chromatin structures and cell fate (31, 158, 162–164).
Immune-Mediated Disease: Rapid Response, Heterogeneity, and Non-Coding Elements and Their Interface With Genome Structure
Rapid Response Genes Are a Risk Factor in Potentiating Auto-Immune Disease
The structures enabling inducible gene expression which is rapid and robust, while indispensable, can cause unintended transcription at the steady state. This aberrant transcription of rapid response inflammatory genes can lead to autoimmune diseases. Specifically, cytokine genes form poised chromatin structure by trained immunity in innate immune cells and memory T cells are implicated in autoimmune disease development (165). For example, rheumatoid arthritis (RA) is mediated by cytokines released by activated macrophages and T cells contributing to damage of cartilage and bone. Synovial macrophages are derived from the embryonic yolk sac and are partially replaced by bone marrow-derived macrophages in the perinatal period. However, an increased number of bone marrow-derived synovial macrophages are associated with the development of RA (166). Specific rapid response genes, such as TNFa and IL8, are constitutively expressed by macrophages in RA patients (167, 168). Additionally, activated Th17 cells infiltrate the synovium and express TNF (169). In these cells, genomic structures are relevant to our understanding of diseases such as RA.
Aberrant rapid response gene expression patterns are in part due to alterations in local genome structure. For example, the promoter of TBX5, a TF responsible for IL8 induction, displays H3K4me3 enrichment and histone acetylation in cells derived from RA patients, leading to the overexpression of TBX5 in fibroblasts (170). Overexpressed IL8 recruits neutrophils, a major source of cartilage-degrading enzymes. In addition to local genome structures, long-range contacts are also implicated in inducing aberrant expression in genes responsible for autoimmune conditions (171). For example, risk variants for RA are located between the gene bodies of TNFAIP3 and OLIG3. Interestingly, these variants physically interact with TNFAIP3 and other immune genes - IL20RA and IFNGR1 - in B and T cells derived from RA patients, but not with OLIG3. This is despite OLIG3 being the most proximal and, therefore, most likely to interact (172). The allele containing the risk variant rs6927172 causes increased contact between the enhancer and IL20RA, resulting in greater expression in T cells. It remains unknown how exactly this variant increases interaction frequency. IL20RA, a receptor of the IL20 family proteins including IL19, IL20 and IL24, activates the JAK-STAT pathway by binding with IL20 and IL24 (173). As the JAK-STAT pathway regulates cytokine gene expression, constitutive activation of JAK-STAT induces cytokine expression, leading to the recruitment of other immune cells in this inflammatory milieu. This, in turn, leads to persistent joint inflammation, proliferative synovitis and, ultimately, damage of the underlying cartilage. Further, it is known that cytokines which display super-enhancer architecture are sensitive to JAK inhibitors (174). These results indicate that local and long-range genomic structure influence autoimmune diseases development and attenuating these aberrant responses mitigate the risks inferred.
Alterations in Genome Structure in Tissue-Resident Immune Cells Are a Risk Factor for Autoimmune Disease
Epigenetic reprogramming at regulatory elements of immune genes, especially rapid response genes in tissue-resident immune cells, comes at the risk of precipitating autoimmune diseases. Because tissue-resident immune cells proliferate by self-renewal once their epigenome is reprogrammed, this can have lifelong transcriptional effects for the cell. Multiple sclerosis (MS) is an important example of an autoimmune disease caused by abnormal activation of tissue-resident immune cells. MS, a neuroinflammatory disease, is characterized by chronic demyelination, oligodendrocyte death, and axonal and neuronal loss (175). Importantly, neurons in the central nervous system appear to be affected in regions of activated microglia. There are over 200 identified risk variants for MS over-and-above MHC-II and these are enriched in brain-resident immune cells, especially microglia (176). This suggests that microglia are a major contributor to the pathophysiology of MS. As in RA, MS risk variants may impact chromatin looping events. Variants located in the intron of CLEC16A interact with DEXI and up-regulate the expression of DEXI in monocytes. It remains unknown how DEXI induces MS. However, DEXI is a gene known to be responsive to dexamethasone, a common treatment for inflammatory disorders, and DEXI is likely involved in the regulation of inflammation (177).
Beyond the classic examples of RA and MS, a greater array of diseases are recognized as autoimmune. Diabetes, for example, although regarded as a metabolic disorder, is known to be fueled by the activation of adipose-resident macrophages (178). Again, risk variants associated with the development of Diabetes Mellitus Type 1 have a hand in changing the configuration of chromatin structure (179). Alzheimer’s Disease (AD) is another disease increasingly classified as auto-immune (180). Risk variants for AD are detected in microglia-specific enhancers but not in neuron- or astrocyte-specific enhancers; suggesting that aberrant activation of microglia could be responsible for AD development (181). A conceptual shift is needed to fully appreciate the causative role of inappropriate immune activation - both innate and adaptive - in disease development. Specifically, investigating chromatin structure’s role in immune regulation may reveal a range of previously-overlooked mediators of disease and therapeutic avenues that remain to be exploited.
The Contribution of Non-Coding Elements to Genome Architecture in Autoimmunity
Numerous studies have identified similar risk variants for autoimmune disease, nearly 90% of which are located in non-coding regions with 60% of these located within enhancers. This suggests that autoimmune diseases can be caused by mis-regulation of gene expression in key immune genes, as opposed to alterations in amino acid sequence (101, 182). Indeed, chromatin conformation investigations have shown that enhancers with risk variants interact with responsible immune genes (101, 132). Of note, it is often assumed that enhancers regulate the nearest genes; yet only 14% of enhancers with risk variants are known to do so, demonstrating the importance of considering long-range genome structure in autoimmune pathophysiology (132). The inference of the variants is complex and the mechanisms by which they influence disease development incompletely described. One possibility is that risk variants change the binding affinity of TFs. For example, the Crohn’s Disease risk variant rs17293632 directly disrupts the AP-1 motif, dysregulating the TGF-β–SMAD3 pathway (182). However, only 10-20% of risk variants are known to directly alter recognizable TF binding motifs (182). Although mutations in flanking regions also alter TF binding affinity, it remains unknown the extent to which this may impact TF binding affinity (183). Interestingly, strain-specific TF binding patterns cannot be explained by mutations (184). Likewise, several studies raise important questions about the effect of mutations on the dynamics of TF binding (185). Another mechanism by which risk variants may cause disease is by means of structural change in the genome. As shown in by the example of risk variant rs6927172 in T cells, there is no doubt that long-range genome structures are intimately involved in the causation of autoimmune disease. In addition, local structures also appears to be regulated by risk variants and a single variant is sufficient to alter chromatin accessibility (186). In this scenario, an important question is how these structural changes are controlled. One potential is the influence of TF binding: it is possible that DNA-DNA interactions are not only mediated by architectural protein such as CTCF, but also by an array of TFs through phase separation. A notable example of the influence of risk variants on loop structure is asthma risk variant which affect CTCF mediated loops and regulate ORMDL3 gene, a negative regulator of IL2 in CD4+ T cells (187). However, as above, this may not be sufficient to explain all cases. Another possibility is the role of certain lncRNAs, such as eRNAs. Some lncRNAs are known to mediate enhancer-promoter interactions (188), likely via their secondary structures and protein binding potential (189). Indeed, the strong linkage disequilibrium between certain risk variants and lncRNA loci are implicated in changes to lncRNA structure (190). For example, it is predicted that rs11153299 and rs2038013, which have a strong linkage disequilibrium with inflammatory bowel disease risk variant, disrupt the structure of lncRNA TRAF3IP2-AS1. Mutations, therefore, may either lead to lncRNA expression suppression or deformation of their secondary structures. Further studies are needed to begin interrogating these unknowns.
Epigenetic Modification as a Precipitant of Autoimmunity in the Context of Genome Architecture
As discussed above, in at least some cases, chromatin structure is implicated in autoimmune disease development. However, alteration of chromatin structure by genetic variants alone is insufficient to precipitate autoimmunity. Most autoimmune conditions occur in older individuals with younger individuals never being risk-variant carriers. HSCs may acquire critical causative mutations only late in life (191). However, we assume this to be a limited driver of autoimmunity as these mutations occur randomly and cannot explain familial patterns of autoimmune disease. To explain this, we propose a model in which altered chromatin structure caused by genetic variants creates a scaffold upon which age-induced epigenetic changes induce abnormal inflammation. For example, the histone deacetylation activity of Sirtuin family proteins, such as SIRT1, declines with age (192). As these senescent cells accumulate histone acetylation marks, enhancers become activated. In cells which are risk-variant carriers, chromatin loops mediated by these variants bridge activated enhancers and inflammatory gene promoters resulting in aberrant expression. As rapid response genes form poised structures in the naïve state, these genes are activated easily in older people carrying risk variants. Techniques able to capture enhancer-promoter interaction and histone modification simultaneously - such as HiChIP and ChIA-Drop - will be powerful tools in fully understanding autoimmune disease development (36, 193). Interestingly, histone modification variations between individuals and single cells increase with age, suggesting that heterogeneity of the chromatin structure may play a role in autoimmune disease development (194). Single cell chromatin profiling techniques will be a key to understand autoimmune disease development as well.
We reviewed that immune cells make rapid response possible by preformed “ready-to-respond” and it’s a risk of autoimmune diseases. Also, we argued that genome structure is heterogenous and how it possibly affects to immune gene expression and decision of phenotype after activation. Still, it is largely unknown the role of heterogeneity of the genome structure in immune response. For example, is plasticity of genome structure important for response to diverse pathogens? Also, it is unknown whether or not there is interrelationship between lower scale genome structure such as nucleosome and higher scale such as DNA looping for heterogeneity. As some TFs mediate DNA loops and depletion of nucleosome is essential for TF binding, DNA looping and nucleosome position will show interrelationship and are not completely independent (195). The techniques that detect multiple data simultaneously such as transcription, nucleosome and DNA loops and live-cell imaging of epigenome will bring a novel insight to these questions.
LD, MM, and YN conceived, designed and wrote a paper. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. Tang L, Nogales E, Ciferri C. Structure and Function of SWI/SNF Chromatin Remodeling Complexes and Mechanistic Implications for Transcription. Prog Biophys Mol Biol (2010) 102:122–8. doi: 10.1016/j.pbiomolbio.2010.05.001
3. Yamada T, Yamaguchi Y, Inukai N, Okamoto S, Mura T, Handa H. P-TEFb-Mediated Phosphorylation of hSpt5 C-Terminal Repeats Is Critical for Processive Transcription Elongation. Mol Cell (2006) 21:227–37. doi: 10.1016/j.molcel.2005.11.024
4. Lu X, Zhu X, Li Y, Liu M, Yu B, Wang Y, et al. Multiple P-TEFbs Cooperatively Regulate the Release of Promoter-Proximally Paused RNA Polymerase II. Nucleic Acids Res (2016) 44:6853–67. doi: 10.1093/nar/gkw571
6. Barboric M, Nissen RM, Kanazawa S, Jabrane-Ferrat N, Peterlin BM. NF-κB Binds P-TEFb to Stimulate Transcriptional Elongation by RNA Polymerase II. Mol Cell (2001) 8:327–37. doi: 10.1016/S1097-2765(01)00314-8
7. Ashburner BP, Westerheide SD, Baldwin AS. The p65 (RelA) Subunit of NF-κB Interacts With the Histone Deacetylase (HDAC) Corepressors HDAC1 and HDAC2 to Negatively Regulate Gene Expression. Mol Cell Biol (2001) 21:7065–77. doi: 10.1128/MCB.21.20.7065-7077.2001
8. Jang MK, Mochizuki K, Zhou M, Jeong H-S, Brady JN, Ozato K. The Bromodomain Protein Brd4 is a Positive Regulatory Component of P-TEFb and Stimulates RNA Polymerase II-Dependent Transcription. Mol Cell (2005) 19:523–34. doi: 10.1016/j.molcel.2005.06.027
9. Dey A, Chitsaz F, Abbasi A, Misteli T, Ozato K. The Double Bromodomain Protein Brd4 Binds to Acetylated Chromatin During Interphase and Mitosis. Proc Natl Acad Sci (2003) 100:8758–63. doi: 10.1073/pnas.1433065100
11. Sabari BR, Tang Z, Huang H, Yong-Gonzalez V, Molina H, Kong HE, et al. Intracellular Crotonyl-CoA Stimulates Transcription Through P300-Catalyzed Histone Crotonylation. Mol Cell (2015) 58:203–15. doi: 10.1016/j.molcel.2015.02.029
13. Arner E, Daub CO, Vitting-Seerup K, Andersson R, Lilje B, Drabløs F, et al. Transcribed Enhancers Lead Waves of Coordinated Transcription in Transitioning Mammalian Cells. Science (2015) 347:1010–4. doi: 10.1126/science.1259418
14. Catarino RR, Stark A. Assessing Sufficiency and Necessity of Enhancer Activities for Gene Expression and the Mechanisms of Transcription Activation. Genes Dev (2018) 32:202–23. doi: 10.1101/gad.310367.117
19. Paolicelli RC, Bolasco G, Pagani F, Maggi L, Scianni M, Panzanelli P, et al. Synaptic Pruning by Microglia is Necessary for Normal Brain Development. Science (2011) 333:1456–8. doi: 10.1126/science.1202529
20. Masuda T, Tsuda M, Yoshinaga R, Tozaki-Saitoh H, Ozato K, Tamura T, et al. IRF8 Is a Critical Transcription Factor for Transforming Microglia Into a Reactive Phenotype. Cell Rep (2012) 1:334–40. doi: 10.1016/j.celrep.2012.02.014
21. Lavin Y, Winter D, Blecher-Gonen R, David E, Keren-Shaul H, Merad M, et al. Tissue-Resident Macrophage Enhancer Landscapes are Shaped by the Local Microenvironment. Cell (2014) 159:1312–26. doi: 10.1016/j.cell.2014.11.018
22. Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome. Science (2009) 326:289–93. doi: 10.1126/science.1181369
23. Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological Domains in Mammalian Genomes Identified by Analysis of Chromatin Interactions. Nature (2012) 485:376–80. doi: 10.1038/nature11082
25. Nora EP, Lajoie BR, Schulz EG, Giorgetti L, Okamoto I, Servant N, et al. Spatial Partitioning of the Regulatory Landscape of the X-inactivation Centre. Nature (2012) 485:381–5. doi: 10.1038/nature11049
28. Levi V, Ruan Q, Plutz M, Belmont AS, Gratton E. Chromatin Dynamics in Interphase Cells Revealed by Tracking in a Two-Photon Excitation Microscope. Biophys J (2005) 89:4275–85. doi: 10.1529/biophysj.105.066670
30. Small EC, Xi L, Wang J-P, Widom J, Licht JD. Single-Cell Nucleosome Mapping Reveals the Molecular Basis of Gene Expression Heterogeneity. Proc Natl Acad Sci (2014) 111:E2462–71. doi: 10.1073/pnas.1400517111
31. Cao J, Cusanovich DA, Ramani V, Aghamirzaie D, Pliner HA, Hill AJ, et al. Joint Profiling of Chromatin Accessibility and Gene Expression in Thousands of Single Cells. Science (2018) 361:1380–5. doi: 10.1126/science.aau0730
32. Cusanovich DA, Daza R, Adey A, Pliner HA, Christiansen L, Gunderson KL, et al. Multiplex Single-Cell Profiling of Chromatin Accessibility by Combinatorial Cellular Indexing. Science (2015) 348:910–4. doi: 10.1126/science.aab1601
33. Larson AG, Elnatan D, Keenen MM, Trnka MJ, Johnston JB, Burlingame AL, et al. Liquid Droplet Formation by HP1α Suggests a Role for Phase Separation in Heterochromatin. Nature (2017) 547:236–40. doi: 10.1038/nature22822
35. Finn EH, Pegoraro G, Brandão HB, Valton A-L, Oomen ME, Dekker J, et al. Extensive Heterogeneity and Intrinsic Variation in Spatial Genome Organization. Cell (2019) 176:1502–15.e10. doi: 10.1016/j.cell.2019.01.020
37. Cattoni DI, Cardozo Gizzi AM, Georgieva M, Di Stefano M, Valeri A, Chamousset D, et al. Single-Cell Absolute Contact Probability Detection Reveals Chromosomes are Organized by Multiple Low-Frequency Yet Specific Interactions. Nat Commun (2017) 8:1753. doi: 10.1038/s41467-017-01962-x
38. Nagano T, Lubling Y, Stevens TJ, Schoenfelder S, Yaffe E, Dean W, et al. Single-Cell Hi-C Reveals Cell-to-Cell Variability in Chromosome Structure. Nature (2013) 502:59–64. doi: 10.1038/nature12593
40. Bintu B, Mateo LJ, Su J-H, Sinnott-Armstrong NA, Parker M, Kinrot S, et al. Super-Resolution Chromatin Tracing Reveals Domains and Cooperative Interactions in Single Cells. Science (2018) 362:eaau1783. doi: 10.1126/science.aau1783
44. De Cian A, Praly E, Ding F, Singh V, Lavelle C, Le Cam E, et al. ATP-Independent Cooperative Binding of Yeast Isw1a to Bare and Nucleosomal DNA. PloS One (2012) 7:e31845. doi: 10.1371/journal.pone.0031845
45. Rao SSP, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D Map of the Human Genome At Kilobase Resolution Reveals Principles of Chromatin Looping. Cell (2014) 159:1665–80. doi: 10.1016/j.cell.2014.11.021
48. Sanborn AL, Rao SSP, Huang S-C, Durand NC, Huntley MH, Jewett AI, et al. Chromatin Extrusion Explains Key Features of Loop and Domain Formation in Wild-Type and Engineered Genomes. Proc Natl Acad Sci (2015) 112:E6456–65. doi: 10.1073/pnas.1518552112
50. Tang Z, Luo OJ, Li X, Zheng M, Zhu JJ, Szalaj P, et al. CTCF-Mediated Human 3D Genome Architecture Reveals Chromatin Topology for Transcription. Cell (2015) 163:1611–27. doi: 10.1016/j.cell.2015.11.024
51. Nora EP, Goloborodko A, Valton A-L, Gibcus JH, Uebersohn A, Abdennur N, et al. Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains From Genomic Compartmentalization. Cell (2017) 169:930–44.e22. doi: 10.1016/j.cell.2017.05.004
52. Schwarzer W, Abdennur N, Goloborodko A, Pekowska A, Fudenberg G, Loe-Mie Y, et al. Two Independent Modes of Chromatin Organization Revealed by Cohesin Removal. Nature (2017) 551:51–6. doi: 10.1038/nature24281
54. Cuartero S, Weiss FD, Dharmalingam G, Guo Y, Ing-Simmons E, Masella S, et al. Control of Inducible Gene Expression Links Cohesin to Hematopoietic Progenitor Self-Renewal and Differentiation. Nat Immunol (2018) 19:932–41. doi: 10.1038/s41590-018-0184-1
55. Allepuz-Fuster P, O’Brien MJ, González-Polo N, Pereira B, Dhoondia Z, Ansari A, et al. RNA Polymerase II Plays an Active Role in the Formation of Gene Loops Through the Rpb4 Subunit. Nucleic Acids Res (2019) 47:8975–87. doi: 10.1093/nar/gkz597
58. Wutz G, Várnai C, Nagasaka K, Cisneros DA, Stocsits RR, Tang W, et al. Topologically Associating Domains and Chromatin Loops Depend on Cohesin and are Regulated by CTCF, WAPL, and PDS5 Proteins. EMBO J (2017) 36:3573–99. doi: 10.15252/embj.201798004
59. Haarhuis JHI, van der Weide RH, Blomen VA, Yáñez-Cuna JO, Amendola M, van Ruiten MS, et al. The Cohesin Release Factor Wapl Restricts Chromatin Loop Extension. Cell (2017) 169:693–707.e14. doi: 10.1016/j.cell.2017.04.013
60. Gassler J, Brandão HB, Imakaev M, Flyamer IM, Ladstätter S, Bickmore WA, et al. A Mechanism of Cohesin-Dependent Loop Extrusion Organizes Zygotic Genome Architecture. EMBO J (2017) 36:3600–18. doi: 10.15252/embj.201798083
61. Rowley MJ, Nichols MH, Lyu X, Ando-Kuri M, Rivera ISM, Hermetz K, et al. Evolutionarily Conserved Principles Predict 3D Chromatin Organization. Mol Cell (2017) 67:837–52.e7. doi: 10.1016/j.molcel.2017.07.022
64. Lin Y, Currie SL, Rosen MK. Intrinsically Disordered Sequences Enable Modulation of Protein Phase Separation Through Distributed Tyrosine Motifs. J Biol Chem (2017) 292:19110–20. doi: 10.1074/jbc.M117.800466
65. van der Lee R, Buljan M, Lang B, Weatheritt RJ, Daughdrill GW, Dunker AK, et al. Classification of Intrinsically Disordered Regions and Proteins. Chem Rev (2014) 114:6589–631. doi: 10.1021/cr400525m
66. Lu H, Yu D, Hansen AS, Ganguly S, Liu R, Heckert A, et al. Phase-Separation Mechanism for C-terminal Hyperphosphorylation of RNA Polymerase II. Nature (2018) 558:318–23. doi: 10.1038/s41586-018-0174-3
70. Teif VB, Vainshtein Y, Caudron-Herger M, Mallm J-P, Marth C, Höfer T, et al. Genome-Wide Nucleosome Positioning During Embryonic Stem Cell Development. Nat Struct Mol Biol (2012) 19:1185–92. doi: 10.1038/nsmb.2419
71. Calderon D, Nguyen MLT, Mezger A, Kathiria A, Müller F, Nguyen V, et al. Landscape of Stimulation-Responsive Chromatin Across Diverse Human Immune Cells. Nat Genet (2019) 51:1494–505. doi: 10.1038/s41588-019-0505-9
72. Kong X, Wei G, Chen N, Zhao S, Shen Y, Zhang J, et al. Dynamic Chromatin Accessibility Profiling Reveals Changes in Host Genome Organization in Response to Baculovirus Infection. PloS Pathog (2020) 16:e1008633. doi: 10.1371/journal.ppat.1008633
74. Minderjahn J, Schmidt A, Fuchs A, Schill R, Raithel J, Babina M, et al. Mechanisms Governing the Pioneering and Redistribution Capabilities of the non-Classical Pioneer PU.1. Nat Commun (2020) 11:402. doi: 10.1038/s41467-019-13960-2
75. Tsatsanis C, Androulidaki A, Alissafi T, Charalampopoulos I, Dermitzaki E, Roger T, et al. Corticotropin-Releasing Factor and the Urocortins Induce the Expression of TLR4 in Macrophages Via Activation of the Transcription Factors PU.1 and AP-1. J Immunol (2006) 176:1869–77. doi: 10.4049/jimmunol.176.3.1869
76. Marecki S, Riendeau CJ, Liang MD, Fenton MJ. PU.1 and Multiple IFN Regulatory Factor Proteins Synergize to Mediate Transcriptional Activation of the Human IL-1β Gene. J Immunol (2001) 166:6829–38. doi: 10.4049/jimmunol.166.11.6829
77. Fischer J, Walter C, Tönges A, Aleth H, Jordão MJC, Leddin M, et al. Safeguard Function of PU.1 Shapes the Inflammatory Epigenome of Neutrophils. Nat Immunol (2019) 20:546–58. doi: 10.1038/s41590-019-0343-z
79. Vanden Berghe W, De Bosscher K, Boone E, Plaisance S, Haegeman G. The Nuclear Factor-κB Engages CBP/p300 and Histone Acetyltransferase Activity for Transcriptional Activation of the Interleukin-6 Gene Promoter. J Biol Chem (1999) 274:32091–8. doi: 10.1074/jbc.274.45.32091
80. Zambrano S, Loffreda A, Carelli E, Stefanelli G, Colombo F, Bertrand E, et al. First Responders Shape a Prompt and Sharp NF-κB-Mediated Transcriptional Response to TNF-α. iScience (2020) 23:101529. doi: 10.1016/j.isci.2020.101529
81. Ramirez-Carrozzi VR, Braas D, Bhatt DM, Cheng CS, Hong C, Doty KR, et al. A Unifying Model for the Selective Regulation of Inducible Transcription by CpG Islands and Nucleosome Remodeling. Cell (2009) 138:114–28. doi: 10.1016/j.cell.2009.04.020
84. Diermeier S, Kolovos P, Heizinger L, Schwartz U, Georgomanolis T, Zirkel A, et al. TNFα Signalling Primes Chromatin for NF-κB Binding and Induces Rapid and Widespread Nucleosome Repositioning. Genome Biol (2014) 15:536. doi: 10.1186/s13059-014-0536-6
85. Escoubet-Lozach L, Benner C, Kaikkonen MU, Lozach J, Heinz S, Spann NJ, et al. Mechanisms Establishing TLR4-Responsive Activation States of Inflammatory Response Genes. PloS Genet (2011) 7:e1002401. doi: 10.1371/journal.pgen.1002401
87. Shih H-Y, Sciumè G, Mikami Y, Guo L, Sun H-W, Brooks SR, et al. Developmental Acquisition of Regulomes Underlies Innate Lymphoid Cell Functionality. Cell (2016) 165:1120–33. doi: 10.1016/j.cell.2016.04.029
88. Mercer EM, Lin YC, Benner C, Jhunjhunwala S, Dutkowski J, Flores M, et al. Multilineage Priming of Enhancer Repertoires Precedes Commitment to the B and Myeloid Cell Lineages in Hematopoietic Progenitors. Immunity (2011) 35:413–25. doi: 10.1016/j.immuni.2011.06.013
89. Barish GD, Yu RT, Karunasiri M, Ocampo CB, Dixon J, Benner C, et al. Bcl-6 and NF- B Cistromes Mediate Opposing Regulation of the Innate Immune Response. Genes Dev (2010) 24:2760–5. doi: 10.1101/gad.1998010
91. Kaikkonen MU, Spann NJ, Heinz S, Romanoski CE, Allison KA, Stender JD, et al. Remodeling of the Enhancer Landscape During Macrophage Activation Is Coupled to Enhancer Transcription. Mol Cell (2013) 51:310–25. doi: 10.1016/j.molcel.2013.07.010
92. Netea MG, Domínguez-Andrés J, Barreiro LB, Chavakis T, Divangahi M, Fuchs E, et al. Defining Trained Immunity and its Role in Health and Disease. Nat Rev Immunol (2020) 20:375–88. doi: 10.1038/s41577-020-0285-6
93. Cheng S-C, Quintin J, Cramer RA, Shepardson KM, Saeed S, Kumar V, et al. mTOR- and HIF-1 -Mediated Aerobic Glycolysis as Metabolic Basis for Trained Immunity. Science (2014) 345:1250684–1250684. doi: 10.1126/science.1250684
94. Arts RJW, Novakovic B, ter Horst R, Carvalho A, Bekkering S, Lachmandas E, et al. Glutaminolysis and Fumarate Accumulation Integrate Immunometabolic and Epigenetic Programs in Trained Immunity. Cell Metab (2016) 24:807–19. doi: 10.1016/j.cmet.2016.10.008
95. Bekkering S, Arts RJW, Novakovic B, Kourtzelis I, van der Heijden CDCC, Li Y, et al. Metabolic Induction of Trained Immunity Through the Mevalonate Pathway. Cell (2018) 172:135–46.e9. doi: 10.1016/j.cell.2017.11.025
96. Mitroulis I, Ruppova K, Wang B, Chen L-S, Grzybek M, Grinenko T, et al. Modulation of Myelopoiesis Progenitors Is an Integral Component of Trained Immunity. Cell (2018) 172:147–61.e12. doi: 10.1016/j.cell.2017.11.034
97. Kaufmann E, Sanz J, Dunn JL, Khan N, Mendonça LE, Pacis A, et al. BCG Educates Hematopoietic Stem Cells to Generate Protective Innate Immunity Against Tuberculosis. Cell (2018) 172:176–90.e19. doi: 10.1016/j.cell.2017.12.031
98. Miller BC, Sen DR, Al Abosy R, Bi K, Virkud YV, LaFleur MW, et al. Subsets of Exhausted CD8+ T Cells Differentially Mediate Tumor Control and Respond to Checkpoint Blockade. Nat Immunol (2019) 20:326–36. doi: 10.1038/s41590-019-0312-6
99. Novakovic B, Habibi E, Wang S-Y, Arts RJW, Davar R, Megchelenbrink W, et al. β-Glucan Reverses the Epigenetic State of LPS-Induced Immunological Tolerance. Cell (2016) 167:1354–68.e14. doi: 10.1016/j.cell.2016.09.034
101. Javierre BM, Burren OS, Wilder SP, Kreuzhuber R, Hill SM, Sewitz S, et al. Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters. Cell (2016) 167:1369–84.e19. doi: 10.1016/j.cell.2016.09.037
102. Pradhananga S, Spalinskas R, Poujade F-A, Eriksson P, Sahlén P. Promoter Anchored Interaction Landscape of THP-1 Macrophages Captures Early Immune Response Processes. Cell Immunol (2020) 355:104148. doi: 10.1016/j.cellimm.2020.104148
104. D’Ippolito AM, McDowell IC, Barrera A, Hong LK, Leichter SM, Bartelt LC, et al. Pre-Established Chromatin Interactions Mediate the Genomic Response to Glucocorticoids. Cell Syst (2018) 7:146–60.e7. doi: 10.1016/j.cels.2018.06.007
105. Stigler J, Çamdere GÖ, Koshland DE, Greene EC. Single-Molecule Imaging Reveals a Collapsed Conformational State for DNA-Bound Cohesin. Cell Rep (2016) 15:988–98. doi: 10.1016/j.celrep.2016.04.003
106. Pham D, Moseley CE, Gao M, Savic D, Winstead CJ, Sun M, et al. Batf Pioneers the Reorganization of Chromatin in Developing Effector T Cells Via Ets1-Dependent Recruitment of Ctcf. Cell Rep (2019) 29:1203–20.e7. doi: 10.1016/j.celrep.2019.09.064
107. Fanucchi S, Fok ET, Dalla E, Shibayama Y, Börner K, Chang EY, et al. Immune Genes are Primed for Robust Transcription by Proximal Long Noncoding RNAs Located in Nuclear Compartments. Nat Genet (2019) 51:138–50. doi: 10.1038/s41588-018-0298-2
109. Sojka DK, Bruniquel D, Schwartz RH, Singh NJ. IL-2 Secretion by CD4 + T Cells In Vivo Is Rapid, Transient, and Influenced by TCR-Specific Competition. J Immunol (2004) 172:6136–43. doi: 10.4049/jimmunol.172.10.6136
111. Chen X, Wang J, Woltring D, Gerondakis S, Shannon MF. Histone Dynamics on the Interleukin-2 Gene in Response to T-Cell Activation. Mol Cell Biol (2005) 25:3209–19. doi: 10.1128/MCB.25.8.3209-3219.2005
113. Yukawa M, Jagannathan S, Vallabh S, Kartashov AV, Chen X, Weirauch MT, et al. AP-1 Activity Induced by Co-Stimulation is Required for Chromatin Opening During T Cell Activation. J Exp Med (2020) 217:e20182009. doi: 10.1084/jem.20182009
114. Yui MA, Hernández-Hoyos G, Rothenberg EV. A New Regulatory Region of the IL-2 Locus That Confers Position-Independent Transgene Expression. J Immunol (2001) 166:1730–9. doi: 10.4049/jimmunol.166.3.1730
115. Bediaga NG, Coughlan HD, Johanson TM, Garnham AL, Naselli G, Schröder J, et al. Multi-Level Remodelling of Chromatin Underlying Activation of Human T Cells. Sci Rep (2021) 11:528. doi: 10.1038/s41598-020-80165-9
116. Yang J, McGovern A, Martin P, Duffus K, Ge X, Zarrineh P, et al. Analysis of Chromatin Organization and Gene Expression in T Cells Identifies Functional Genes for Rheumatoid Arthritis. Nat Commun (2020) 11:4402. doi: 10.1038/s41467-020-18180-7
117. Russ BE, Olshanksy M, Smallwood HS, Li J, Denton AE, Prier JE, et al. Distinct Epigenetic Signatures Delineate Transcriptional Programs During Virus-Specific CD8+ T Cell Differentiation. Immunity (2014) 41:853–65. doi: 10.1016/j.immuni.2014.11.001
118. Denton AE, Russ BE, Doherty PC, Rao S, Turner SJ. Differentiation-Dependent Functional and Epigenetic Landscapes for Cytokine Genes in Virus-Specific CD8+ T Cells. Proc Natl Acad Sci (2011) 108:15306–11. doi: 10.1073/pnas.1112520108
119. Tu WJ, Hardy K, Sutton CR, McCuaig R, Li J, Dunn J, et al. Priming of Transcriptional Memory Responses Via the Chromatin Accessibility Landscape in T Cells. Sci Rep (2017) 7:44825. doi: 10.1038/srep44825
123. Schoenborn JR, Dorschner MO, Sekimata M, Santer DM, Shnyreva M, Fitzpatrick DR, et al. Comprehensive Epigenetic Profiling Identifies Multiple Distal Regulatory Elements Directing Transcription of the Gene Encoding Interferon-γ. Nat Immunol (2007) 8:732–42. doi: 10.1038/ni1474
127. Hatton RD, Harrington LE, Luther RJ, Wakefield T, Janowski KM, Oliver JR, et al. A Distal Conserved Sequence Element Controls Ifng Gene Expression by T Cells and NK Cells. Immunity (2006) 25:717–29. doi: 10.1016/j.immuni.2006.09.007
128. Hakim O, Sung M-H, Nakayamada S, Voss TC, Baek S, Hager GL. Spatial Congregation of STAT Binding Directs Selective Nuclear Architecture During T-cell Functional Differentiation. Genome Res (2013) 23:462–72. doi: 10.1101/gr.147652.112
131. Cai S, Lee CC, Kohwi-Shigematsu T. SATB1 Packages Densely Looped, Transcriptionally Active Chromatin for Coordinated Expression of Cytokine Genes. Nat Genet (2006) 38:1278–88. doi: 10.1038/ng1913
132. Mumbach MR, Satpathy AT, Boyle EA, Dai C, Gowen BG, Cho SW, et al. Enhancer Connectome in Primary Human Cells Identifies Target Genes of Disease-Associated DNA Elements. Nat Genet (2017) 49:1602–12. doi: 10.1038/ng.3963
133. Fischer C, Metsger M, Bauch S, Vidal R, Böttcher M, Grote P, et al. Signals Trigger State-Specific Transcriptional Programs to Support Diversity and Homeostasis in Immune Cells. Sci Signal (2019) 12:eaao5820. doi: 10.1126/scisignal.aao5820
134. Avraham R, Haseley N, Brown D, Penaranda C, Jijon HB, Trombetta JJ, et al. Pathogen Cell-to-Cell Variability Drives Heterogeneity in Host Immune Responses. Cell (2015) 162:1309–21. doi: 10.1016/j.cell.2015.08.027
139. Hagai T, Chen X, Miragaia RJ, Rostom R, Gomes T, Kunowska N, et al. Gene Expression Variability Across Cells and Species Shapes Innate Immunity. Nature (2018) 563:197–202. doi: 10.1038/s41586-018-0657-2
140. Shalek AK, Satija R, Adiconis X, Gertner RS, Gaublomme JT, Raychowdhury R, et al. Single-Cell Transcriptomics Reveals Bimodality in Expression and Splicing in Immune Cells. Nature (2013) 498:236–40. doi: 10.1038/nature12172
142. Whyte WA, Orlando DA, Hnisz D, Abraham BJ, Lin CY, Kagey MH, et al. Master Transcription Factors and Mediator Establish Super-Enhancers At Key Cell Identity Genes. Cell (2013) 153:307–19. doi: 10.1016/j.cell.2013.03.035
143. Kolovos P, Georgomanolis T, Koeferle A, Larkin JD, Brant L, Nikolicć M, et al. Binding of Nuclear Factor κB to Noncanonical Consensus Sites Reveals its Multimodal Role During the Early Inflammatory Response. Genome Res (2016) 26:1478–89. doi: 10.1101/gr.210005.116
144. Brown JD, Lin CY, Duan Q, Griffin G, Federation A, Paranal RM, et al. Nf-κB Directs Dynamic Super Enhancer Formation in Inflammation and Atherogenesis. Mol Cell (2014) 56:219–31. doi: 10.1016/j.molcel.2014.08.024
145. Osterwalder M, Barozzi I, Tissières V, Fukuda-Yuzawa Y, Mannion BJ, Afzal SY, et al. Enhancer Redundancy Provides Phenotypic Robustness in Mammalian Development. Nature (2018) 554:239–43. doi: 10.1038/nature25461
146. Ren G, Jin W, Cui K, Rodrigez J, Hu G, Zhang Z, et al. CTCF-Mediated Enhancer-Promoter Interaction Is a Critical Regulator of Cell-to-Cell Variation of Gene Expression. Mol Cell (2017) 67:1049–58.e6. doi: 10.1016/j.molcel.2017.08.026
147. Nagashima R, Hibino K, Ashwin SS, Babokhov M, Fujishiro S, Imai R, et al. Single Nucleosome Imaging Reveals Loose Genome Chromatin Networks Via Active RNA Polymerase II. J Cell Biol (2019) 218:1511–30. doi: 10.1083/jcb.201811090
149. Lai B, Gao W, Cui K, Xie W, Tang Q, Jin W, et al. Principles of Nucleosome Organization Revealed by Single-Cell Micrococcal Nuclease Sequencing. Nature (2018) 562:281–5. doi: 10.1038/s41586-018-0567-3
150. Ku WL, Nakamura K, Gao W, Cui K, Hu G, Tang Q, et al. Single-Cell Chromatin Immunocleavage Sequencing (scChIC-seq) to Profile Histone Modification. Nat Methods (2019) 16:323–5. doi: 10.1038/s41592-019-0361-7
151. Zhu C, Zhang Y, Li YE, Lucero J, Behrens MM, Ren B. Joint Profiling of Histone Modifications and Transcriptome in Single Cells From Mouse Brain. Nat Methods (2021) 18:283–92. doi: 10.1038/s41592-021-01060-3
154. Hart Y, Reich-Zeliger S, Antebi YE, Zaretsky I, Mayo AE, Alon U, et al. Paradoxical Signaling by a Secreted Molecule Leads to Homeostasis of Cell Levels. Cell (2014) 158:1022–32. doi: 10.1016/j.cell.2014.07.033
155. Almeida ARM, Amado IF, Reynolds J, Berges J, Lythe G, Molina-París C, et al. Quorum-Sensing in CD4+ T Cell Homeostasis: A Hypothesis and a Model. Front Immunol (2012) 3:1–15. doi: 10.3389/fimmu.2012.00125
156. Feinerman O, Jentsch G, Tkach KE, Coward JW, Hathorn MM, Sneddon MW, et al. Single-Cell Quantification of IL-2 Response by Effector and Regulatory T Cells Reveals Critical Plasticity in Immune Response. Mol Syst Biol (2010) 6:437. doi: 10.1038/msb.2010.90
157. Ageeva-Kieferle A, Rudolf EE, Lindermayr C. Redox-Dependent Chromatin Remodeling: A New Function of Nitric Oxide as Architect of Chromatin Structure in Plants. Front Plant Sci (2019) 10:625. doi: 10.3389/fpls.2019.00625
158. Ma S, Zhang B, LaFave LM, Earl AS, Chiang Z, Hu Y, et al. Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin. Cell (2020) 183:1103–16.e20. doi: 10.1016/j.cell.2020.09.056
159. Sato Y, Hilbert L, Oda H, Wan Y, Heddleston JM, Chew T-L, et al. Histone H3K27 Acetylation Precedes Active Transcription During Zebrafish Zygotic Genome Activation as Revealed by Live-Cell Analysis. Development (2019) 146:dev179127. doi: 10.1242/dev.179127
166. Tu J, Hong W, Guo Y, Zhang P, Fang Y, Wang X, et al. Ontogeny of Synovial Macrophages and the Roles of Synovial Macrophages From Different Origins in Arthritis. Front Immunol (2019) 10:1146. doi: 10.3389/fimmu.2019.01146
168. Georganas C, Liu H, Perlman H, Hoffmann A, Thimmapaya B, Pope RM. Regulation of IL-6 and IL-8 Expression in Rheumatoid Arthritis Synovial Fibroblasts: The Dominant Role for NF-κB But Not C/EBPβ or C-Jun. J Immunol (2000) 165:7199–206. doi: 10.4049/jimmunol.165.12.7199
169. Pène J, Chevalier S, Preisser L, Vénéreau E, Guilleux M-H, Ghannam S, et al. Chronically Inflamed Human Tissues are Infiltrated by Highly Differentiated Th17 Lymphocytes. J Immunol (2008) 180:7423–30. doi: 10.4049/jimmunol.180.11.7423
170. Karouzakis E, Trenkmann M, Gay RE, Michel BA, Gay S, Neidhart M. Epigenome Analysis Reveals TBX5 as a Novel Transcription Factor Involved in the Activation of Rheumatoid Arthritis Synovial Fibroblasts. J Immunol (2014) 193:4945–51. doi: 10.4049/jimmunol.1400066
172. McGovern A, Schoenfelder S, Martin P, Massey J, Duffus K, Plant D, et al. Capture Hi-C Identifies a Novel Causal Gene, IL20RA, in the Pan-Autoimmune Genetic Susceptibility Region 6q23. Genome Biol (2016) 17:212. doi: 10.1186/s13059-016-1078-x
173. Kragstrup TW, Andersen T, Heftdal LD, Hvid M, Gerwien J, Sivakumar P, et al. The IL-20 Cytokine Family in Rheumatoid Arthritis and Spondyloarthritis. Front Immunol (2018) 9:2226. doi: 10.3389/fimmu.2018.02226
176. International Multiple Sclerosis Genetics Consortium*†. Multiple Sclerosis Genomic Map Implicates Peripheral Immune Cells and Microglia in Susceptibility. Science (2019) 365:eaav7188. doi: 10.1126/science.aav7188
177. Davison LJ, Wallace C, Cooper JD, Cope NF, Wilson NK, Smyth DJ, et al. Long-Range DNA Looping and Gene Expression Analyses Identify DEXI as an Autoimmune Disease Candidate Gene. Hum Mol Genet (2012) 21:322–33. doi: 10.1093/hmg/ddr468
178. Tsalamandris S, Antonopoulos AS, Oikonomou E, Papamikroulis G-A, Vogiatzi G, Papaioannou S, et al. The Role of Inflammation in Diabetes: Current Concepts and Future Perspectives. Eur Cardiol Rev (2019) 14:50–9. doi: 10.15420/ecr.2018.33.1
179. Tomlinson MJ, Pitsillides A, Pickin R, Mika M, Keene KL, Hou X, et al. Fine Mapping and Functional Studies of Risk Variants for Type 1 Diabetes At Chromosome 16p13.13. Diabetes (2014) 63:4360–8. doi: 10.2337/db13-1785
181. Nott A, Holtman IR, Coufal NG, Schlachetzki JCM, Yu M, Hu R, et al. Brain Cell Type–Specific Enhancer–Promoter Interactome Maps and Disease - Risk Association. Science (2019) 366:1134–9. doi: 10.1126/science.aay0793
182. Farh KK-H, Marson A, Zhu J, Kleinewietfeld M, Housley WJ, Beik S, et al. Genetic and Epigenetic Fine Mapping of Causal Autoimmune Disease Variants. Nature (2015) 518:337–43. doi: 10.1038/nature13835
183. Yella VR, Bhimsaria D, Ghoshdastidar D, Rodríguez-Martínez JA, Ansari AZ, Bansal M. Flexibility and Structure of Flanking DNA Impact Transcription Factor Affinity for its Core Motif. Nucleic Acids Res (2018) 46:11883–97. doi: 10.1093/nar/gky1057
184. Link VM, Duttke SH, Chun HB, Holtman IR, Westin E, Hoeksema MA, et al. Analysis of Genetically Diverse Macrophages Reveals Local and Domain-wide Mechanisms That Control Transcription Factor Binding and Function. Cell (2018) 173:1796–809.e17. doi: 10.1016/j.cell.2018.04.018
187. Schmiedel BJ, Seumois G, Samaniego-Castruita D, Cayford J, Schulten V, Chavez L, et al. 17q21 Asthma-Risk Variants Switch CTCF Binding and Regulate IL-2 Production by T Cells. Nat Commun (2016) 7:13426. doi: 10.1038/ncomms13426
190. Ren C, Liu F, Ouyang Z, An G, Zhao C, Shuai J, et al. Functional Annotation of Structural ncRNAs Within Enhancer RNAs in the Human Genome: Implications for Human Disease. Sci Rep (2017) 7:15518. doi: 10.1038/s41598-017-15822-7
193. Mumbach MR, Rubin AJ, Flynn RA, Dai C, Khavari PA, Greenleaf WJ, et al. HiChIP: Efficient and Sensitive Analysis of Protein-Directed Genome Architecture. Nat Methods (2016) 13:919–22. doi: 10.1038/nmeth.3999
194. Cheung P, Vallania F, Warsinske HC, Donato M, Schaffert S, Chang SE, et al. Single-Cell Chromatin Modification Profiling Reveals Increased Epigenetic Variations With Aging. Cell (2018) 173:1385–97.e14. doi: 10.1016/j.cell.2018.03.079
Keywords: genome structure, immune regulation, rapid response, heterogeneity, autoimmune diseases
Citation: Devenish LP, Mhlanga MM and Negishi Y (2021) Immune Regulation in Time and Space: The Role of Local- and Long-Range Genomic Interactions in Regulating Immune Responses. Front. Immunol. 12:662565. doi: 10.3389/fimmu.2021.662565
Received: 01 February 2021; Accepted: 26 April 2021;
Published: 11 May 2021.
Edited by:Amy L. Kenter, University of Illinois at Chicago, United States
Reviewed by:Andrew D. Wells, Children’s Hospital of Philadelphia, United States
Masaki Miyazaki, Kyoto University, Japan
Copyright © 2021 Devenish, Mhlanga and Negishi. 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: Yutaka Negishi, email@example.com