Connor Davidson1,2
Matteo Vecellio- 1Wellcome Centre of Human Genetics, University of Oxford, Oxford, United Kingdom
- 2Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, Botnar Research Centre, University of Oxford, Oxford, United Kingdom
- 3MRC WIMM Centre for Computational Biology, MRC Weatherall Institute for Molecular Medicine, University of Oxford, Oxford, United Kingdom
- 4Centro Ricerche Fondazione Italiana Ricerca Sull’Artrite (FIRA), Fondazione Pisana x la Scienza ONLUS, San Giuliano Terme, Italy
Ankylosing Spondylitis (AS) is a chronic inflammatory arthritis of the spine exhibiting a strong genetic background. The mechanistic and functional understanding of the AS-associated genomic loci, identified with Genome Wide Association Studies (GWAS), remains challenging. Chromosome conformation capture (3C) and derivatives are recent techniques which are of great help in elucidating the spatial genome organization and of enormous support in uncover a mechanistic explanation for disease-associated genetic variants. The perturbation of three-dimensional (3D) genome hierarchy may lead to a plethora of human diseases, including rheumatological disorders. Here we illustrate the latest approaches and related findings on the field of genome organization, highlighting how the instability of 3D genome conformation may be among the causes of rheumatological disease phenotypes. We suggest a new perspective on the inclusive potential of a 3C approach to inform GWAS results in rheumatic diseases. 3D genome organization may ultimately lead to a more precise and comprehensive functional interpretation of AS association, which is the starting point for emerging and more specific therapies.
Introduction
The combination of environmental and genetic factors may lead to the development of complex diseases (Hunter, 2005). Ankylosing Spondylitis (AS) is a common form of arthritis primarily affecting the spine, characterised by inflammation at the entheses (Bridgewood et al., 2019) and sacroiliac joints (Brewerton et al., 1973). AS is a highly heritable disease with more than 100 genomic loci found implicated in increasing the risk (International Genetics of Ankylosing Spondylitis Consortium et al., 2013; Ellinghaus et al., 2016). Genome wide association studies (GWAS) have been very successful in polygenic disease as they identified thousands of common genetic variants or single nucleotide polymorphisms (SNPs), which can have a phenotypical individual effect (Huo et al., 2019; Crouch and Bodmer, 2020). The identification of a causal variant from GWAS data may help our understanding of complex traits biology, suggesting new target genes and methods of controlling them. Unfortunately, disease-associated loci often contain multiple genes making the scenario extremely challenging; genetic variants in proximal vicinity tend to be inherited together, in a phenomenon called linkage disequilibrium (LD), making difficult to identify the causal variant underpinning the association (Cano-Gamez and Trynka, 2020). In recent years, large-scale epigenomic projects have mapped hundreds of thousands of potential regulatory sites in the human genome, but only a small proportion of these elements are proximal to transcription start sites (Bagchi and Iyer, 2016). In AS, we and others were able to successfully identify the causal functional SNPs at the Interleukin 23 Receptor (IL23R), endoplasmic reticulum aminopeptidase 1 (ERAP1) and RUNX Family Transcription Factor 3 (RUNX3) genomic loci elucidating their transcriptional regulation (Keidel et al., 2013; Roberts et al., 2016; Vecellio et al., 2016; Vecellio et al., 2018; Vecellio et al., 2021).
The three-dimensional (3D) organization of the genome is essential in facilitating fundamental processes which occur in the cell nucleus including, transcriptional regulation, DNA damage and replication (McCord et al., 2020). Over the last 20 years, chromosome conformation capture (3C) techniques have been widely used to identify and estimate the frequency of interaction of multiple genomic loci in the genome (Dekker et al., 2002; Phillips-Cremins et al., 2013; Hansen et al., 2018; Haws et al., 2022). In 3C methodology, restriction enzyme digestion followed by re-ligation of cross-linked chromatin in the nucleus of a cell, allows to detect the spatial vicinity between DNA sequences (de Wit and de Laat, 2012). 3C experiments have revealed that chromosomes are folded in complex structures emerging at different scales. These structures can be impacted by disease-associated SNPs, as reported extensively (Gorkin et al., 2019; Anania and Lupianez, 2020; Tsuchiya et al., 2021). Recent studies in polygenic disorders (Girdhar et al., 2022) show that using 3D genome architecture investigation has the utility to clarify the role of disease associated SNPs and to link them to specific genes to understand the phenotype and account for biological function (Khunsriraksakul et al., 2022; Zhao et al., 2022).
Original 3C methods are low-throughput and not able to define if multiple regions interact simultaneously or mutually exclusively. For this reason, several technologies deriving from the standard 3C have been developed, including Hi-C (high-throughput chromosome conformation capture) which allows the analysis of spatial genome organization and chromosome folding through sequencing (Dekker et al., 2002; van Steensel and Dekker, 2010; Dekker et al., 2013).
Here, we explore a selection of 3C methods (i.e. Hi-C, Tiling Capture-C) which might facilitate the understanding of 3D organization and chromosomal interactions and their impact in rheumatic diseases pathophysiology. We emphasise the importance of 3C approaches to inform GWAS interpretation, and their possible future application in precision medicine in prioritizing potential drug targets in polygenic rheumatic disorders including AS.
The 3D genome and the nuclear architecture
Chromosomes fold into compartments, often indicated as A and B, which refer to genomic loci with similar transcriptional activity that physically segregate in 3D space. In addition, contact domains, consist in any visible domains corresponding to elevated chromatin interactions regions. Topologically associated domains (TADs) are the hallmarks of genomic organization (see Figure 1) and are defined as local organizational domains taught to be formed primarily by loop extrusion where boundaries are most conserved during cell differentiation (<1 MB scale) (Gorkin et al., 2019; Goel and Hansen, 2021). TAD dysregulation is linked to various diseases, including neurological disorders and tumorigenesis (Medrano-Fernandez and Barco, 2016; Yang et al., 2019). As demonstrated by Lupianez et al. using CRISPR/Cas9 genome editing and 3C methods, the disruption of TADs might lead to a rewiring of long-range regulatory architectures and result in a pathogenic phenotype. Specifically, the authors focused on rare limb malformations and identified several rearrangements in the epha4/pax3 (EPH Receptor A4/Paired Box 3) locus in mice causing disruption of the TAD, chromatin structural changes and the aberrant expression of developmental genes (Lupianez et al., 2015; Lupianez et al., 2016). In addition, Luppino and others have demonstrated the loss of cohesin leads to reduced chromatin mixing thus affecting the topology and transcriptional bursting frequencies of boundary-proximal genes (Luppino et al., 2020).
FIGURE 1. TADs insulation. Simplified model showing how topologically associated domains (TADs) are insulated by borders which are leaky enough to influence and regulate nearby genes expression (Luppino et al., 2020).
Loop boundaries are enriched in CCCTC-binding factor (CTCF) and cohesin, two architectural necessary proteins involved in long-range genome looping. It has been postulated that these two proteins form domains by a loop extrusion process. First proposed in 2015 following mathematical simulations and polymer modelling, the loop extrusion model (see Figure 2) suggests that a chromatin loop of increasing size is extruded by the cohesin complex until it is stalled by a pair of convergently oriented CTCF-bound sites (de Wit et al., 2015; Sanborn et al., 2015). The cohesin complex extrudes a loop uni- or bi-directionally until it faces an occupied CTCF binding site and then the loop is stabilized (Fudenberg et al., 2017). Several studies have explored the interplay between compartmentalization and chromatin looping, showing that the depletion of CTCF has no effect on compartments (Nora et al., 2017), while cohesin (and/or RAD21) removal brings to the loss of domains and TADs making compartmentalization more prominent (Haarhuis et al., 2017; Wutz et al., 2017). In addition to CTCF and cohesin, boundaries frequently colocalize with active transcription start sites along with additional genomic factors such as YY1, RAD21 and ZNF143 which exhibit enrichment at strong boundaries (Hsieh et al., 2015; Bonev et al., 2017). Boundaries which are depleted of CTCF and YY1 are defined as weak boundaries. In addition, chromatin immunoprecipitation (ChIP)-seq experiments have revealed enrichment for ASH2L (ASH2 like, histone lysine methyltransferase complex subunit), H3K4Me3, SP1 (specificity protein 1) among other factors enriched at these boundaries. The role of RNA polymerase II has been investigated by Hsieh and others (Hsieh et al., 2015) demonstrating that active transcription mediated genome folding has a crucial role in the maintenance of the enhancer-promoter and promoter-promoter domains. Following the inhibition of RNA polymerase II, the intensity of those domains is significantly reduced without affecting higher-order chromatin organization (Hsieh et al., 2020). Stable enhancer-promoter interactions have been observed during the formation of Drosophila embryos, suggesting these interactions are important in developmental stages, cell fate decision and limb formation (Ghavi-Helm et al., 2014; Ji et al., 2016). Abnormalities in enhancer-promoter interactions, such as mutation in encoding proteins genes or enhancer-binding proteins lead to disease like Cornelia de Lange syndrome, often referred as enhanceropathy (Olley et al., 2018). The process of transcription termination at the 3’ of a gene requires the recruitment of specific factors, which cross-talk with the initiation and enhancement machinery required for the start of transcription. 3C approach has been crucial in demonstrating the formation of gene loops showing interaction between transcription factors associated with promoter and those linked with transcription termination (Al-Husini et al., 2020).
FIGURE 2. The loop-extrusion model. The model shows the generation of long-range cis-interactions, following extrusion by cohesin, CTCF cognate sites binding and involvement of other accessory proteins including YY1, RAD21 and ZNF143 (Hansen et al., 2018).
The generation of comprehensive and high-resolution 3D genome maps may facilitate the discovery of small alterations (Chakraborty and Ay, 2019) which can lead a cascade of aberrant molecular events and drive disease phenotype.
3C-derived technologies to map the 3D genome
Recent advances in specific methods to investigate the 3D genome architecture include development of high-throughput chromosome conformation capture (Hi-C), Hi-ChIP and Tiled Capture (Mumbach et al., 2016; Galan et al., 2020; Oudelaar, 2022). All three methods capture genomic regions linked via 3D interactions, that are detected by Next-Generation Sequencing (NGS).
Hi-C was introduced in 2009, using biotinylation to enrich for proximity ligated contacts and thus modifying the library amplification process. Specifically, HiC takes advantage of using universal adapters and primers for high throughput sequencing. The unbiased “all versus all” approach clearly has a real advantage in defining all the genomic interactions genome wide (Lieberman-Aiden et al., 2009). It has been demonstrated that during differentiation, transcriptional changes occur when there is alteration in the strength of long-range interactions and the development of cell-type specific enhancer-promoter contacts (Bonev et al., 2017). Further, these interactions occur primarily in the same TAD and are strongly correlated with gene expression demonstrating how TADs constrain enhancer activity (Symmons et al., 2014).
Recently, Micro-C was developed from Hi-C, moving to a capture-fine (∼1 kb) to a nucleosome level of resolution (∼200bp) (Krietenstein et al., 2020). Formaldehyde and disuccinimidyl glutarate are used for fixation and cross-linking steps, while chromatin is digested with micrococcal nuclease (MNase), instead of restriction enzymes, to reach nucleosome-level resolution. Restriction enzymes sites are not equally distributed and not all the DNA is readily accessible as it is affected by nucleosomal accessibility (Szerlong and Hansen, 2011); using micrococcal nuclease digestion takes advantage of local DNA accessibility and facilitates the retention of intact nucleosomes (Voong et al., 2017). Micro-C has revealed two different classes of CTCF loops: those dependent on RNA-binding region (RBR) and those which are not (Hansen et al., 2019).
Hi-ChIP has been developed to delineate promoter-enhancer interactions by leveraging principles of in situ Hi-C (Mumbach et al., 2016), combining long-range contacts investigation with enrichment of specific histone proteins (i.e H3K27Ac) associated with active regions of the genome. Recently, Chandra and others performed Hi-ChIP to provide evidence of non-coding genetic variants having effect on gene expression (i.e. cis-eQTL) and cell-specific gene regulation in five immune cell types (Chandra et al., 2021). Another method coupling ChIP with 3C is ChIA-PET, a chromatin interaction analysis by paired-end tag sequencing (Fullwood et al., 2009). With this approach is possible to detect all chromatin interactions mediated by a specific protein of interest, immunoprecipitated using specific antibodies (Li et al., 2002; Fullwood et al., 2010). In 2013, ChIA-PET was used to described for the first time the genome-wide chromatin interactions of cohesin (DeMare et al., 2013).
Tiling Capture C was designed to identify at high resolution level the interaction between promoter and enhancer within TADs. This method uses a panel of specific capture oligonucleotides tiled across all contiguous restriction fragments within specified genomic regions, typically around 1Mb (Downes et al., 2022), to obtain enrichment for specific interactions and subsequent targeted sequencing. High-resolution maps with low-cells input makes Tiled-C at forefront of the 3C methods. Oudelaar et al. developed the Tiled-C approach to characterize the chromatin architecture of mouse erythroid cells during in vivo differentiation, focusing on the α-globin locus (Oudelaar et al., 2020). Recently, the same group has generated the most detailed genomic local interaction map at base-pair resolution (20bp), using a micrococcal nuclease (MN) based 3C approach. MN is an enzyme digesting the genome largely independent of DNA sequence: the authors benefit of MN to demonstrate the effects of the depletion of two crucial elements, cohesin and CTCF, on chromatin architecture (Hua et al., 2021; Aljahani et al., 2022). A summary of the different chromatin conformation techniques here presented is showed in Table 1.
TABLE 1. A summary of the different techniques routinely used to analyse the conformation of the genome.
The investigation of 3D genome topology has dramatically evolved over the past decade: technological advances have boost the field to an unprecedented level (Bouwman et al., 2022). Ongoing efforts are made in providing more accurate data analysis for a better understanding and interpretation of the functional consequence of these changes.
Investigating the 3D genome in rheumatic diseases and AS
The 3D genome organization has been investigated in autoimmune and rheumatic disease, but more is yet to come. In 2012, 3C was used to identify physical interaction between the chromosome region 16p13, often associated with increased risk for a plethora of autoimmune diseases, including multiple sclerosis (Zuvich et al., 2011), primary biliary cirrhosis (Mells et al., 2011) and systemic lupus erythematosus (SLE) (Gateva et al., 2009), with DEXI (dexamethasone-induced protein), a gene with previous unknown function, revealed as a strong candidate for autoimmune disease (Davison et al., 2012).
Meddens and others performed circular chromosome conformation capture (4C)-seq to analyse chromatin interactions in inflammatory bowel disease (IBD) susceptibility loci and DNA regulatory elements providing novel relevant candidate genes (Meddens et al., 2016). In a recent work by Carini and others, the authors explored the genomic architecture of whole blood obtained from rheumatoid arthritis (RA) patients investigating a possible chromosome conformation signature (including IFNAR1, IL-21R, IL-23, IL-17A and CXCL13), before and after the administration of methotrexate (MTX) treatment. This was important to identify the non-responders to disease modified anti-rheumatic drugs such as MTX and, whether there was an association between the chromatin signature and RA-specific expression quantitative trait loci (eQTL) (Carini et al., 2018). Integrating data from Hi-C with gene expression profiling and disease activity scores has been successful in a recent work on SLE, where the authors established the genomic interaction landscape identifying specific SLE-associated loops (Zhao et al., 2022).
Nearly 90% of disease-associated SNPs are located in non-coding regions (Ricano-Ponce and Wijmenga, 2013). The role of specific GWAS hits can be elucidated via 3D genome analysis, thus clarifying which genes are influenced by which particular SNP through a spatial connection (Li et al., 2020). Our group has recently demonstrated the presence of a chromatin loop between the AS-associated SNP rs4648889 and the distal promoter of the RUNX3 gene, confirming together with other functional experiments previously reported the primacy of this genetic variant in the association with RUNX3 in AS (Cohen et al., 2021). The complexity of the RUNX3 locus is also confirmed by Capture-C experiments showing multiple interactions among different SNPs and the RUNX3 promoter (personal communication).
In 2015 the Orozco group investigated chromatin interactions between disease-associated genetic variants and their functional targets in B and T cells in four autoimmune disease, including RA, type 1 diabetes, psoriatic arthritis, and juvenile idiopathic arthritis. They performed Capture Hi-C and demonstrated that only few looping interactions were common to both cell lines and disease-associated SNPs interact with candidate genes relevant to the disease and located megabases away (Martin et al., 2015).
It is important to dissect the role of chromatin contacts at diverse genomic loci: the usage of 3D genome structure to perform gene prioritization will be very informative to define new drug targets and evaluate if a therapy is working or a more effective therapy is needed.
As previously demonstrated, genetic variations can influence 3D chromatin conformation, together with accessibility and gene expression (Gorkin et al., 2019). Long-distance eQTLs potentially regulate gene expression and spatial gene regulatory interactions are supposed to be the drivers of the heritability of complex traits. Our completed genome-wide study of chromatin interactions and the regulatory effects of AS-associated genetic variants is unprecedented in the field and it will be very informative in identifying genes and cells to prioritize as therapeutic targets (Brown et al., unpublished data). On the same line, promoter capture Hi-C and RNA-sequencing approaches were recently used to link associated variants of systemic sclerosis (a connective tissue immune-mediated disease) with their target genes, especially in CD4+T cells and CD14+ monocytes obtained from 10 patients and five matched healthy controls. The authors identified new potential targets genes and 15 other potential drug targets for repurposing of drugs already in use in other immune-mediated diseases (Gonzalez-Serna et al., 2022).
In RA, a comprehensive genomic map has been recently generated to link risk-associated genetic variants with functional chromatin interactions, active regulatory DNA elements and differential gene expression in fibroblast-like synoviocytes, providing the proof of concept for a causal role of these cells in RA susceptibility (Ge et al., 2021).
The examples provided in this section show once again the importance of a multimodal approach for the identification of cell types and molecular states critically associated with rheumatic diseases. The generation of a comprehensive and high-resolution 3D genome map may yield insights in disease-associated TAD appearance and chromosome loop strengths.
Discussion
Although chromosome conformation capture is a relatively new field of investigation, understanding 3D genome folding and its influence on gene expression has rapidly grown beside the innovation of 3C methods. The relevance of looping formation and genomic organization together with the identification of architectural proteins associated with boundaries might shed light on the functional and mechanistic implications in different diseases, as demonstrated in these seminal works focused on neurological disorders and cancer (Bharadwaj et al., 2014; Weischenfeldt et al., 2017; Zhu et al., 2020). Further, disease-associated genetic variants may disrupt higher-order genomic organization, due to elimination of annotated boundaries (Ibn-Salem et al., 2014). Normal and disease-associated TAD structure data may yield valuable and perhaps diagnostically important information on gene regulation and disease aetiology. Several diseases could be linked to an aberrant chromatin loop dynamic (Mehrjouy et al., 2018). This suggests an increased interest in studying genome looping and how this may affect gene expression and function in diseases (Krumm and Duan, 2019).
3C approaches might be useful in addressing an unmet clinical need of predicting those patients who will not respond to specific treatments and thus facilitating earlier access to more effective therapies and a better quality of life. The integration of 3C methodologies with functional data (i.e. eQTL) may pinpoint individual loci into a gene regulatory network which is critical to our understanding of complex diseases. Single-cell genomic assays might be a promising tool for the quantification of molecular traits (i.e. transcriptomics, chromatin accessibility, transcription factors occupancy and time-course trajectories of cells) at single cell level (Buenrostro et al., 2015; Zheng et al., 2017; Cuomo et al., 2020; Del Priore et al., 2021). Recently it has been elegantly demonstrated that genetic variants associated with differential binding of PU.1, a master transcription factor regulating myeloid development and having a substantial effect on neutrophils function (McKercher et al., 1996; Siwaponanan et al., 2017), are predominantly cell type specific, associated with specific chromatin state, and regulate enhancer-promoter interactions and downstream gene expression, exhibiting association with IBD susceptibility (Watt et al., 2021). Additionally, PU.1 has been found involved in regulating the interaction loop at DDX60L (Probable ATP-dependent RNA helicase DDX60-like) promoter thus inducing overexpression of this gene in CD4+ T Cells from SLE patients (Zhao et al., 2022).
The dynamic interplay between genomic sequence, 3D chromatin structure and a specific pathological process such AS, is crucial in defining a powerful strategy to discover novel genetic regions potential for diagnostic or therapeutic purposes. Nevertheless, it is important to bear in mind that a combination of more genetic variants (i.e. haplotypes) may have a higher risk on the probability of developing a particular disorder, while individually may have a mild influence.
Concluding remarks
The last decade has been a remarkable time for genetic research in AS, from genetic association studies to genetics-driven novel clinical trials (Baeten et al., 2013; van der Heijde et al., 2019; Izana Bioscience, 2021). The identification of few hundreds genomic loci make the overall analysis challenging, considering most of them are enriched not only in immune cell-specific enhancers, but also in osteoclasts or stromal cells, having a crucial involvement in the pathogenesis of AS. Contemplating the results of GWAS suggests that a traditional in vitro approach to uncover the mechanistic contribution of disease-associated SNPs is not sufficient anymore. Times for integration of genetic fine mapping of AS loci with DNA architecture and 3D chromatin interactions, DNA accessibility, single-cell gene expression and gene editing are mature as they occur in related disorders (Al-Mossawi and Coates, 2018; Penkava et al., 2020; Yager et al., 2021), and it appears only a matter of time before they are fully applied in AS.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
Writing – original draft, CD and MV; Writing – Review and Editing, CD, BPW, CJC, JCK and MV, and all authors reviewed and approved the final 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Al-Husini, N., Medler, S., and Ansari, A. (2020). Crosstalk of promoter and terminator during RNA polymerase II transcription cycle. Biochim. Biophys. Acta Gene Regul. Mech. 1863 (12), 194657. doi:10.1016/j.bbagrm.2020.194657
Al-Mossawi, H., and Coates, L. C. (2018). Personalized medicine - a new reality in psoriatic arthritis? Nat. Rev. Rheumatol. 14 (8), 449–451. doi:10.1038/s41584-018-0043-3
Aljahani, A., Hua, P., Karpinska, M. A., Quililan, K., Davies, J. O. J., and Oudelaar, A. M. (2022). Analysis of sub-kilobase chromatin topology reveals nano-scale regulatory interactions with variable dependence on cohesin and CTCF. Nat. Commun. 13 (1), 2139. doi:10.1038/s41467-022-29696-5
Anania, C., and Lupianez, D. G. (2020). Order and disorder: Abnormal 3D chromatin organization in human disease. Brief. Funct. Genomics 19 (2), 128–138. doi:10.1093/bfgp/elz028
Baeten, D., Baraliakos, X., Braun, J., Sieper, J., Emery, P., van der Heijde, D., et al. (2013). Anti-interleukin-17A monoclonal antibody secukinumab in treatment of ankylosing spondylitis: A randomised, double-blind, placebo-controlled trial. Lancet 382 (9906), 1705–1713. doi:10.1016/S0140-6736(13)61134-4
Bagchi, D. N., and Iyer, V. R. (2016). The determinants of directionality in transcriptional initiation. Trends Genet. 32 (6), 322–333. doi:10.1016/j.tig.2016.03.005
Bharadwaj, R., Peter, C. J., Jiang, Y., Roussos, P., Vogel-Ciernia, A., Shen, E. Y., et al. (2014). Conserved higher-order chromatin regulates NMDA receptor gene expression and cognition. Neuron 84 (5), 997–1008. doi:10.1016/j.neuron.2014.10.032
Bonev, B., Mendelson Cohen, N., Szabo, Q., Fritsch, L., Papadopoulos, G. L., Lubling, Y., et al. (2017). Multiscale 3D genome rewiring during mouse neural development. Cell. 171 (3), 557–572. doi:10.1016/j.cell.2017.09.043
Bouwman, B. A. M., Crosetto, N., and Bienko, M. (2022). The era of 3D and spatial genomics. Trends Genet. 38 (10), 1062–1075. doi:10.1016/j.tig.2022.05.010
Brewerton, D. A., Hart, F. D., Nicholls, A., Caffrey, M., James, D. C., and Sturrock, R. D. (1973). Ankylosing spondylitis and HL-A 27. Lancet 1 (7809), 904–907. doi:10.1016/s0140-6736(73)91360-3
Bridgewood, C., Watad, A., Russell, T., Palmer, T. M., Marzo-Ortega, H., Khan, A., et al. (2019). Identification of myeloid cells in the human enthesis as the main source of local IL-23 production. Ann. Rheum. Dis. 78 (7), 929–933. doi:10.1136/annrheumdis-2018-214944
Buenrostro, J. D., Wu, B., Chang, H. Y., and Greenleaf, W. J. (2015). ATAC-seq: A method for assaying chromatin accessibility genome-wide. Curr. Protoc. Mol. Biol. 109 (21), 1–21. doi:10.1002/0471142727.mb2129s109
Cano-Gamez, E., and Trynka, G. (2020). From GWAS to function: Using functional genomics to identify the mechanisms underlying complex diseases. Front. Genet. 11, 424. doi:10.3389/fgene.2020.00424
Carini, C., Hunter, E., Ramadass, A. S., Green, J., Akoulitchev, A., et al. Scottish Early Rheumatoid Arthritis Inception cohort Investigators (2018). Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis. J. Transl. Med. 16 (1), 18. doi:10.1186/s12967-018-1387-9
Chakraborty, A., and Ay, F. (2019). The role of 3D genome organization in disease: From compartments to single nucleotides. Semin. Cell. Dev. Biol. 90, 104–113. doi:10.1016/j.semcdb.2018.07.005
Chandra, V., Bhattacharyya, S., Schmiedel, B. J., Madrigal, A., Gonzalez-Colin, C., Fotsing, S., et al. (2021). Promoter-interacting expression quantitative trait loci are enriched for functional genetic variants. Nat. Genet. 53 (1), 110–119. doi:10.1038/s41588-020-00745-3
Cohen, C. J., Davidson, C., Selmi, C., Bowness, P., and Knight, J. C. (2021). Wordsworth BP and Vecellio disruption of c-MYC binding and chromosomal looping involving genetic variants associated with ankylosing spondylitis upstream of the RUNX3 promoter. Front. Genet. 12, 741867. doi:10.3389/fgene.2021.741867
Crouch, D. J. M., and Bodmer, W. F. (2020). Polygenic inheritance, GWAS, polygenic risk scores, and the search for functional variants. Proc. Natl. Acad. Sci. U. S. A. 117 (32), 18924–18933. doi:10.1073/pnas.2005634117
Cuomo, A. S. E., Seaton, D. D., McCarthy, D. J., Martinez, I., Bonder, M. J., Garcia-Bernardo, J., et al. (2020). Single-cell RNA-sequencing of differentiating iPS cells reveals dynamic genetic effects on gene expression. Nat. Commun. 11 (1), 810. doi:10.1038/s41467-020-14457-z
Davison, L. J., Wallace, C., Cooper, J. D., Cope, N. F., Wilson, N. K., Smyth, D. J., et al. (2012). Long-range DNA looping and gene expression analyses identify DEXI as an autoimmune disease candidate gene. Hum. Mol. Genet. 21 (2), 322–333. doi:10.1093/hmg/ddr468
de Wit, E., and de Laat, W. (2012). A decade of 3C technologies: Insights into nuclear organization. Genes. Dev. 26 (1), 11–24. doi:10.1101/gad.179804.111
de Wit, E., Vos, E. S., Holwerda, S. J., Valdes-Quezada, C., Verstegen, M. J., Teunissen, H., et al. (2015). CTCF binding polarity determines chromatin looping. Mol. Cell. 60 (4), 676–684. doi:10.1016/j.molcel.2015.09.023
Dekker, J., Marti-Renom, M. A., and Mirny, L. A. (2013). Exploring the three-dimensional organization of genomes: Interpreting chromatin interaction data. Nat. Rev. Genet. 14 (6), 390–403. doi:10.1038/nrg3454
Dekker, J., Rippe, K., Dekker, M., and Kleckner, N. (2002). Capturing chromosome conformation. Science 295 (5558), 1306–1311. doi:10.1126/science.1067799
Del Priore, I., Ma, S., Strecker, J., Jacks, T., LaFave, L. M., and Buenrostro, J. D. (2021). Protocol for single-cell ATAC sequencing using combinatorial indexing in mouse lung adenocarcinoma. Star. Protoc. 2 (2), 100583. doi:10.1016/j.xpro.2021.100583
DeMare, L. E., Leng, J., Cotney, J., Reilly, S. K., Yin, J., Sarro, R., et al. (2013). The genomic landscape of cohesin-associated chromatin interactions. Genome Res. 23 (8), 1224–1234. doi:10.1101/gr.156570.113
Downes, D. J., Smith, A. L., Karpinska, M. A., Velychko, T., Rue-Albrecht, K., Sims, D., et al. (2022). Capture-C: A modular and flexible approach for high-resolution chromosome conformation capture. Nat. Protoc. 17 (2), 445–475. doi:10.1038/s41596-021-00651-w
Ellinghaus, D., Jostins, L., Spain, S. L., Cortes, A., Bethune, J., Han, B., et al. (2016). Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat. Genet. 48 (5), 510–518. doi:10.1038/ng.3528
Fudenberg, G., Abdennur, N., Imakaev, M., Goloborodko, A., and Mirny, L. A. (2017). Emerging evidence of chromosome folding by loop extrusion. Cold Spring Harb. Symp. Quant. Biol. 82, 45–55. doi:10.1101/sqb.2017.82.034710
Fullwood, M. J., Han, Y., Wei, C. L., Ruan, X., and Ruan, Y. (2010). Chromatin interaction analysis using paired-end tag sequencing. Curr. Protoc. Mol. Biol. 89 (21), 15.1–21.15. doi:10.1002/0471142727.mb2115s89
Fullwood, M. J., Liu, M. H., Pan, Y. F., Liu, J., Xu, H., Mohamed, Y. B., et al. (2009). An oestrogen-receptor-alpha-bound human chromatin interactome. Nature 462 (7269), 58–64. doi:10.1038/nature08497
Galan, S., Machnik, N., Kruse, K., Diaz, N., Marti-Renom, M. A., and Vaquerizas, J. M. (2020). CHESS enables quantitative comparison of chromatin contact data and automatic feature extraction. Nat. Genet. 52 (11), 1247–1255. doi:10.1038/s41588-020-00712-y
Gateva, V., Sandling, J. K., Hom, G., Taylor, K. E., Chung, S. A., Sun, X., et al. (2009). A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat. Genet. 41 (11), 1228–1233. doi:10.1038/ng.468
Ge, X., Frank-Bertoncelj, M., Klein, K., McGovern, A., Kuret, T., Houtman, M., et al. (2021). Functional genomics atlas of synovial fibroblasts defining rheumatoid arthritis heritability. Genome Biol. 22 (1), 247. doi:10.1186/s13059-021-02460-6
Ghavi-Helm, Y., Klein, F. A., Pakozdi, T., Ciglar, L., Noordermeer, D., Huber, W., et al. (2014). Enhancer loops appear stable during development and are associated with paused polymerase. Nature 512 (7512), 96–100. doi:10.1038/nature13417
Girdhar, K., Hoffman, G. E., Bendl, J., Rahman, S., Dong, P., Liao, W., et al. (2022). Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains. Nat. Neurosci. 25 (4), 474–483. doi:10.1038/s41593-022-01032-6
Goel, V. Y., and Hansen, A. S. (2021). The macro and micro of chromosome conformation capture. Wiley Interdiscip. Rev. Dev. Biol. 10 (6), e395. doi:10.1002/wdev.395
Gonzalez-Serna, D., Shi, C., Kerick, M., Hankinson, J., Ding, J., McGovern, A., et al. (2022). Functional genomics in primary T cells and monocytes identifies mechanisms by which genetic susceptibility loci influence systemic sclerosis risk. Arthritis Rheumatol. doi:10.1002/art.42396
Gorkin, D. U., Qiu, Y., Hu, M., Fletez-Brant, K., Liu, T., Schmitt, A. D., et al. (2019). Common DNA sequence variation influences 3-dimensional conformation of the human genome. Genome Biol. 20 (1), 255. doi:10.1186/s13059-019-1855-4
Haarhuis, J. H. I., van der Weide, R. H., Blomen, V. A., Yanez-Cuna, J. O., Amendola, M., van Ruiten, M. S., et al. (2017). The cohesin release factor WAPL restricts chromatin loop extension. Cell. 169 (4), 693–707. doi:10.1016/j.cell.2017.04.013
Hansen, A. S., Cattoglio, C., Darzacq, X., and Tjian, R. (2018). Recent evidence that TADs and chromatin loops are dynamic structures. Nucleus 9 (1), 20–32. doi:10.1080/19491034.2017.1389365
Hansen, A. S., Hsieh, T. S., Cattoglio, C., Pustova, I., Saldana-Meyer, R., Reinberg, D., et al. (2019). Distinct classes of chromatin loops revealed by deletion of an RNA-binding region in CTCF. Mol. Cell. 76 (3), 395–411. doi:10.1016/j.molcel.2019.07.039
Haws, S. A., Simandi, Z., Barnett, R. J., and Phillips-Cremins, J. E. (2022). 3D genome, on repeat: Higher-order folding principles of the heterochromatinized repetitive genome. Cell. 185 (15), 2690–2707. doi:10.1016/j.cell.2022.06.052
Hsieh, T. H., Weiner, A., Lajoie, B., Dekker, J., Friedman, N., and Rando, O. J. (2015). Mapping nucleosome resolution chromosome folding in yeast by micro-C. Cell. 162 (1), 108–119. doi:10.1016/j.cell.2015.05.048
Hsieh, T. S., Cattoglio, C., Slobodyanyuk, E., Hansen, A. S., Rando, O. J., Tjian, R., et al. (2020). Resolving the 3D landscape of transcription-linked mammalian chromatin folding. Mol. Cell. 78 (3), 539–553. doi:10.1016/j.molcel.2020.03.002
Hua, P., Badat, M., Hanssen, L. L. P., Hentges, L. D., Crump, N., Downes, D. J., et al. (2021). Defining genome architecture at base-pair resolution. Nature 595 (7865), 125–129. doi:10.1038/s41586-021-03639-4
Hunter, D. J. (2005). Gene-environment interactions in human diseases. Nat. Rev. Genet. 6 (4), 287–298. doi:10.1038/nrg1578
Huo, Y., Li, S., Liu, J., Li, X., and Luo, X. J. (2019). Functional genomics reveal gene regulatory mechanisms underlying schizophrenia risk. Nat. Commun. 10 (1), 670. doi:10.1038/s41467-019-08666-4
Ibn-Salem, J., Kohler, S., Love, M. I., Chung, H. R., Huang, N., Hurles, M. E., et al. (2014). Deletions of chromosomal regulatory boundaries are associated with congenital disease. Genome Biol. 15 (9), 423. doi:10.1186/s13059-014-0423-1
International Genetics of Ankylosing Spondylitis Consortium Cortes, A., Hadler, J., Pointon, J. P., Robinson, P. C., Karaderi, T., et al. (2013). Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat. Genet. 45 (7), 730–738. doi:10.1038/ng.2667
Izana Bioscience, L. (2021). Efficacy and safety of namilumab for moderate-to-severe axial spondyloarthritis (NAMASTE). Available at: https://ClinicalTrialsgov/show/NCT03622658.
Ji, X., Dadon, D. B., Powell, B. E., Fan, Z. P., Borges-Rivera, D., Shachar, S., et al. (2016). 3D chromosome regulatory landscape of human pluripotent cells. Cell. Stem Cell. 18 (2), 262–275. doi:10.1016/j.stem.2015.11.007
Keidel, S., Chen, L., Pointon, J., and Wordsworth, P. (2013). ERAP1 and ankylosing spondylitis. Curr. Opin. Immunol. 25 (1), 97–102. doi:10.1016/j.coi.2012.11.002
Khunsriraksakul, C., McGuire, D., Sauteraud, R., Chen, F., Yang, L., Wang, L., et al. (2022). Integrating 3D genomic and epigenomic data to enhance target gene discovery and drug repurposing in transcriptome-wide association studies. Nat. Commun. 13 (1), 3258. doi:10.1038/s41467-022-30956-7
Krietenstein, N., Abraham, S., Venev, S. V., Abdennur, N., Gibcus, J., Hsieh, T. S., et al. (2020). Ultrastructural details of mammalian chromosome architecture. Mol. Cell. 78 (3), 554–565. doi:10.1016/j.molcel.2020.03.003
Krumm, A., and Duan, Z. (2019). Understanding the 3D genome: Emerging impacts on human disease. Semin. Cell. Dev. Biol. 90, 62–77. doi:10.1016/j.semcdb.2018.07.004
Li, Q. L., Ito, K., Sakakura, C., Fukamachi, H., Inoue, K., Chi, X. Z., et al. (2002). Causal relationship between the loss of RUNX3 expression and gastric cancer. Cell. 109 (1), 113–124. doi:10.1016/s0092-8674(02)00690-6
Li, Y., Tao, T., Du, L., and Zhu, X. (2020). Three-dimensional genome: Developmental technologies and applications in precision medicine. J. Hum. Genet. 65 (6), 497–511. doi:10.1038/s10038-020-0737-7
Lieberman-Aiden, E., van Berkum, N. L., Williams, L., Imakaev, M., Ragoczy, T., Telling, A., et al. (2009). Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science 326 (5950), 289–293. doi:10.1126/science.1181369
Lupianez, D. G., Kraft, K., Heinrich, V., Krawitz, P., Brancati, F., Klopocki, E., et al. (2015). Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions. Cell. 161 (5), 1012–1025. doi:10.1016/j.cell.2015.04.004
Lupianez, D. G., Spielmann, M., and Mundlos, S. (2016). Breaking TADs: How alterations of chromatin domains result in disease. Trends Genet. 32 (4), 225–237. doi:10.1016/j.tig.2016.01.003
Luppino, J. M., Park, D. S., Nguyen, S. C., Lan, Y., Xu, Z., Yunker, R., et al. (2020). Cohesin promotes stochastic domain intermingling to ensure proper regulation of boundary-proximal genes. Nat. Genet. 52 (8), 840–848. doi:10.1038/s41588-020-0647-9
Martin, P., McGovern, A., Orozco, G., Duffus, K., Yarwood, A., Schoenfelder, S., et al. (2015). Capture Hi-C reveals novel candidate genes and complex long-range interactions with related autoimmune risk loci. Nat. Commun. 6, 10069. doi:10.1038/ncomms10069
McCord, R. P., Kaplan, N., and Giorgetti, L. (2020). Chromosome conformation capture and beyond: Toward an integrative view of chromosome structure and function. Mol. Cell. 77 (4), 688–708. doi:10.1016/j.molcel.2019.12.021
McKercher, S. R., Torbett, B. E., Anderson, K. L., Henkel, G. W., Vestal, D. J., Baribault, H., et al. (1996). Targeted disruption of the PU.1 gene results in multiple hematopoietic abnormalities. EMBO J. 15 (20), 5647–5658. doi:10.1002/j.1460-2075.1996.tb00949.x
Meddens, C. A., Harakalova, M., van den Dungen, N. A., Foroughi Asl, H., Hijma, H. J., Cuppen, E. P., et al. (2016). Systematic analysis of chromatin interactions at disease associated loci links novel candidate genes to inflammatory bowel disease. Genome Biol. 17 (1), 247. doi:10.1186/s13059-016-1100-3
Medrano-Fernandez, A., and Barco, A. (2016). Nuclear organization and 3D chromatin architecture in cognition and neuropsychiatric disorders. Mol. Brain 9 (1), 83. doi:10.1186/s13041-016-0263-x
Mehrjouy, M. M., Fonseca, A. C. S., Ehmke, N., Paskulin, G., Novelli, A., Benedicenti, F., et al. (2018). Regulatory variants of FOXG1 in the context of its topological domain organisation. Eur. J. Hum. Genet. 26 (2), 186–196. doi:10.1038/s41431-017-0011-4
Mells, G. F., Floyd, J. A., Morley, K. I., Cordell, H. J., Franklin, C. S., Shin, S. Y., et al. (2011). Genome-wide association study identifies 12 new susceptibility loci for primary biliary cirrhosis. Nat. Genet. 43 (4), 329–332. doi:10.1038/ng.789
Mumbach, M. R., Rubin, A. J., Flynn, R. A., Dai, C., Khavari, P. A., Greenleaf, W. J., et al. (2016). HiChIP: Efficient and sensitive analysis of protein-directed genome architecture. Nat. Methods 13 (11), 919–922. doi:10.1038/nmeth.3999
Nora, E. P., Goloborodko, A., Valton, A. L., Gibcus, J. H., Uebersohn, A., Abdennur, N., et al. (2017). Targeted degradation of CTCF decouples local insulation of chromosome domains from genomic compartmentalization. Cell. 169 (5), 930–944. doi:10.1016/j.cell.2017.05.004
Olley, G., Ansari, M., Bengani, H., Grimes, G. R., Rhodes, J., von Kriegsheim, A., et al. (2018). BRD4 interacts with NIPBL and BRD4 is mutated in a Cornelia de Lange-like syndrome. Nat. Genet. 50 (3), 329–332. doi:10.1038/s41588-018-0042-y
Oudelaar, A. M., Beagrie, R. A., Gosden, M., de Ornellas, S., Georgiades, E., Kerry, J., et al. (2020). Dynamics of the 4D genome during in vivo lineage specification and differentiation. Nat. Commun. 11 (1), 2722. doi:10.1038/s41467-020-16598-7
Oudelaar, A. M. (2022). Deciphering the regulatory logic of a chromatin domain boundary. Nat. Genet. 54 (7), 914–915. doi:10.1038/s41588-022-01086-z
Penkava, F., Velasco-Herrera, M. D. C., Young, M. D., Yager, N., Nwosu, L. N., Pratt, A. G., et al. (2020). Single-cell sequencing reveals clonal expansions of pro-inflammatory synovial CD8 T cells expressing tissue-homing receptors in psoriatic arthritis. Nat. Commun. 11 (1), 4767. doi:10.1038/s41467-020-18513-6
Phillips-Cremins, J. E., Sauria, M. E., Sanyal, A., Gerasimova, T. I., Lajoie, B. R., Bell, J. S., et al. (2013). Architectural protein subclasses shape 3D organization of genomes during lineage commitment. Cell. 153 (6), 1281–1295. doi:10.1016/j.cell.2013.04.053
Ricano-Ponce, I., and Wijmenga, C. (2013). Mapping of immune-mediated disease genes. Annu. Rev. Genomics Hum. Genet. 14, 325–353. doi:10.1146/annurev-genom-091212-153450
Roberts, A. R., Vecellio, M., Chen, L., Ridley, A., Cortes, A., Knight, J. C., et al. (2016). An ankylosing spondylitis-associated genetic variant in the IL23R-IL12RB2 intergenic region modulates enhancer activity and is associated with increased Th1-cell differentiation. Ann. Rheum. Dis. 75 (12), 2150–2156. doi:10.1136/annrheumdis-2015-208640
Sanborn, A. L., Rao, S. S., Huang, S. C., Durand, N. C., Huntley, M. H., Jewett, A. I., et al. (2015). Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes. Proc. Natl. Acad. Sci. U. S. A. 112 (47), E6456–E6465. doi:10.1073/pnas.1518552112
Siwaponanan, P., Siegers, J. Y., Ghazali, R., Ng, T., McColl, B., Ng, G. Z., et al. (2017). Reduced PU.1 expression underlies aberrant neutrophil maturation and function in beta-thalassemia mice and patients. Blood 129 (23), 3087–3099. doi:10.1182/blood-2016-07-730135
Symmons, O., Uslu, V. V., Tsujimura, T., Ruf, S., Nassari, S., Schwarzer, W., et al. (2014). Functional and topological characteristics of mammalian regulatory domains. Genome Res. 24 (3), 390–400. doi:10.1101/gr.163519.113
Szerlong, H. J., and Hansen, J. C. (2011). Nucleosome distribution and linker DNA: Connecting nuclear function to dynamic chromatin structure. Biochem. Cell. Biol. 89 (1), 24–34. doi:10.1139/O10-139
Tsuchiya, H., Ota, M., Sumitomo, S., Ishigaki, K., Suzuki, A., Sakata, T., et al. (2021). Parsing multiomics landscape of activated synovial fibroblasts highlights drug targets linked to genetic risk of rheumatoid arthritis. Ann. Rheum. Dis. 80 (4), 440–450. doi:10.1136/annrheumdis-2020-218189
van der Heijde, D., Song, I. H., Pangan, A. L., Deodhar, A., van den Bosch, F., Maksymowych, W. P., et al. (2019). Efficacy and safety of upadacitinib in patients with active ankylosing spondylitis (SELECT-AXIS 1): A multicentre, randomised, double-blind, placebo-controlled, phase 2/3 trial. Lancet 394 (10214), 2108–2117. doi:10.1016/S0140-6736(19)32534-6
van Steensel, B., and Dekker, J. (2010). Genomics tools for unraveling chromosome architecture. Nat. Biotechnol. 28 (10), 1089–1095. doi:10.1038/nbt.1680
Vecellio, M., Chen, L., Cohen, C. J., Cortes, A., Li, Y., Bonham, S., et al. (2021). Functional genomic analysis of a RUNX3 polymorphism associated with ankylosing spondylitis. Arthritis Rheumatol. 73 (6), 980–990. doi:10.1002/art.41628
Vecellio, M., Cortes, A., Roberts, A. R., Ellis, J., Cohen, C. J., Knight, J. C., et al. (2018). Evidence for a second ankylosing spondylitis-associated RUNX3 regulatory polymorphism. RMD Open 4 (1), e000628. doi:10.1136/rmdopen-2017-000628
Vecellio, M., Roberts, A. R., Cohen, C. J., Cortes, A., Knight, J. C., Bowness, P., et al. (2016). The genetic association of RUNX3 with ankylosing spondylitis can be explained by allele-specific effects on IRF4 recruitment that alter gene expression. Ann. Rheum. Dis. 75 (8), 1534–1540. doi:10.1136/annrheumdis-2015-207490
Voong, L. N., Xi, L., Wang, J. P., and Wang, X. (2017). Genome-wide mapping of the nucleosome landscape by micrococcal nuclease and chemical mapping. Trends Genet. 33 (8), 495–507. doi:10.1016/j.tig.2017.05.007
Watt, S., Vasquez, L., Walter, K., Mann, A. L., Kundu, K., Chen, L., et al. (2021). Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease. Nat. Commun. 12 (1), 2298. doi:10.1038/s41467-021-22548-8
Weischenfeldt, J., Dubash, T., Drainas, A. P., Mardin, B. R., Chen, Y., Stutz, A. M., et al. (2017). Pan-cancer analysis of somatic copy-number alterations implicates IRS4 and IGF2 in enhancer hijacking. Nat. Genet. 49 (1), 65–74. doi:10.1038/ng.3722
Wutz, G., Varnai, C., Nagasaka, K., Cisneros, D. A., Stocsits, R. R., Tang, W., et al. (2017). Topologically associating domains and chromatin loops depend on cohesin and are regulated by CTCF, WAPL, and PDS5 proteins. EMBO J. 36 (24), 3573–3599. doi:10.15252/embj.201798004
Yager, N., Cole, S., Lledo Lara, A., Maroof, A., Penkava, F., Knight, J. C., et al. (2021). Ex vivo mass cytometry analysis reveals a profound myeloid proinflammatory signature in psoriatic arthritis synovial fluid. Ann. Rheum. Dis. 80 (12), 1559–1567. doi:10.1136/annrheumdis-2021-220280
Yang, M., Vesterlund, M., Siavelis, I., Moura-Castro, L. H., Castor, A., Fioretos, T., et al. (2019). Proteogenomics and Hi-C reveal transcriptional dysregulation in high hyperdiploid childhood acute lymphoblastic leukemia. Nat. Commun. 10 (1), 1519. doi:10.1038/s41467-019-09469-3
Zhao, M., Feng, D., Hu, L., Liu, L., Wu, J., Hu, Z., et al. (2022). 3D genome alterations in T cells associated with disease activity of systemic lupus erythematosus. Ann. Rheumatic Dis. 2022, 2022-222653. annrheumdis-2022-222653. doi:10.1136/ard-2022-222653
Zheng, G. X., Terry, J. M., Belgrader, P., Ryvkin, P., Bent, Z. W., Wilson, R., et al. (2017). Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 8, 14049. doi:10.1038/ncomms14049
Zhu, H., Uuskula-Reimand, L., Isaev, K., Wadi, L., Alizada, A., Shuai, S., et al. (2020). Candidate cancer driver mutations in distal regulatory elements and long-range chromatin interaction networks. Mol. Cell. 77 (6), 1307–1321. doi:10.1016/j.molcel.2019.12.027
Zuvich, R. L., Bush, W. S., McCauley, J. L., Beecham, A. H., De Jager, P. L., et al. International Multiple Sclerosis Genetics Consortium (2011). Interrogating the complex role of chromosome 16p13.13 in multiple sclerosis susceptibility: Independent genetic signals in the CIITA-CLEC16A-SOCS1 gene complex. Hum. Mol. Genet. 20 (17), 3517–3524. doi:10.1093/hmg/ddr250
Keywords: genomics, three dimensional genome, ankylosing spondylitis, chromosome conformation capture (3C), topologically associated domain (TAD), rheumatic and musculoskeletal disease
Citation: Davidson C, Wordsworth BP, Cohen CJ, Knight JC and Vecellio M (2023) Chromosome conformation capture approaches to investigate 3D genome architecture in Ankylosing Spondylitis. Front. Genet. 14:1129207. doi: 10.3389/fgene.2023.1129207
Received: 21 December 2022; Accepted: 16 January 2023;
Published: 25 January 2023.
Edited by:
Jordi Pérez-Tur, Institute of Biomedicine of Valencia (CSIC), SpainReviewed by:
Athar Ansari, Wayne State University, United StatesCopyright © 2023 Davidson, Wordsworth, Cohen, Knight and Vecellio. 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: Matteo Vecellio, matteo.vecellio@ndorms.ox.ac.uk