Analysis of chromatin organization and gene expression in T cells identifies functional genes for rheumatoid arthritis

Cell-Based Therapy and Genome Editing as Emerging Therapeutic Approaches to Treat Rheumatoid Arthritis

Rheumatoid arthritis (RA) is an autoimmune disease characterized by chronic inflammation of the joints. Although much remains unknown about the pathogenesis of RA, there is evidence that impaired immune tolerance and the development of RA are related. And it is precisely the restoration of immune tolerance at the site of the inflammation that is the ultimate goal of the treatment of RA. Over the past few decades, significant progress has been made in the treatment of RA, with higher rates of disease remission and improved long-term outcomes. Unfortunately, despite these successes, the proportion of patients with persistent, difficult-to-treat disease remains high, and the task of improving our understanding of the basic mechanisms of disease development and developing new ways to treat RA remains relevant. This review focuses on describing new treatments for RA, including cell therapies and gene editing technologies that have shown potential in preclinical and early clinical trials. In addition, we discuss the opportunities and limitations associated with the use of these new approaches in the treatment of RA.

Mechanism of EHMT2-mediated genomic imprinting associated with Prader-Willi syndrome

Prader-Willi Syndrome (PWS) is caused by loss of expression of paternally expressed genes in the human 15q11.2-q13 imprinting domain. A set of imprinted genes that are active on the paternal but silenced on the maternal chromosome are intricately regulated by a bipartite imprinting center (PWS-IC) located in the PWS imprinting domain. In past work, we discovered that euchromatic histone lysine N-methyltransferase-2 (EHMT2/G9a) inhibitors were capable of un-silencing PWS-associated genes by restoring their expression from the maternal chromosome. Here, in mice lacking the Ehmt2 gene, we document un-silencing of the imprinted Snrpn/Snhg14 gene on the maternal chromosome in the late embryonic and postnatal brain. Using PWS and Angelman syndrome patient derived cells with either paternal or maternal deletion of 15q11-q13, we have found that chromatin of maternal PWS-IC is closed and has compact 3D folding confirmation. We further show that a new and distinct noncoding RNA preferentially transcribed from upstream of the PWS-IC interacts with EHMT2 and forms a heterochromatin complex to silence gene expression of SNRPN in CIS on maternal chromosome. Taken together, these findings demonstrate that allele-specific recruitment of EHMT2 is required to maintain the maternal imprints. Our findings provide novel mechanistic insights and support a new model for imprinting maintenance of the PWS imprinted domain.

Regulation of 3D genome organization during T cell activation

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

Application and perspective of CRISPR/Cas9 genome editing technology in human diseases modeling and gene therapy

The Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) mediated Cas9 nuclease system has been extensively used for genome editing and gene modification in eukaryotic cells. CRISPR/Cas9 technology holds great potential for various applications, including the correction of genetic defects or mutations within the human genome. The application of CRISPR/Cas9 genome editing system in human disease research is anticipated to solve a multitude of intricate molecular biology challenges encountered in life science research. Here, we review the fundamental principles underlying CRISPR/Cas9 technology and its recent application in neurodegenerative diseases, cardiovascular diseases, autoimmune related diseases, and cancer, focusing on the disease modeling and gene therapy potential of CRISPR/Cas9 in these diseases. Finally, we provide an overview of the limitations and future prospects associated with employing CRISPR/Cas9 technology for diseases study and treatment.

Emerging strategies for nanomedicine in autoimmunity

Autoimmune disorders have risen to be among the most prevalent chronic diseases across the globe, affecting approximately 5-7% of the population. As autoimmune diseases steadily rise in prevalence, so do the number of potential therapeutic strategies to combat them. In recent years, fundamental research investigating autoimmune pathologies has led to the emergence of several cellular targets that provide new therapeutic opportunities. However, key challenges persist in terms of accessing and specifically combating the dysregulated, self-reactive cells while avoiding systemic immune suppression and other off-target effects. Fortunately, the continued advancement of nanomedicines may provide strategies to address these challenges and bring innovative autoimmunity therapies to the clinic. Through precise engineering and rational design, nanomedicines can possess a variety of physicochemical properties, surface modifications, and cargoes, allowing for specific targeting of therapeutics to pathological cell and organ types. These advances in nanomedicine have been demonstrated in cancer therapies and have the broad potential to advance applications in autoimmunity therapies as well. In this review, we focus on leveraging the power of nanomedicine for prevalent autoimmune disorders throughout the body. We expand on three key areas for the development of autoimmunity therapies - avoiding systemic immunosuppression, balancing interactions with the immune system, and elevating current platforms for delivering complex cargoes - and emphasize how nanomedicine-based strategies can overcome these barriers and enable the development of next-generation, clinically relevant autoimmunity therapies.

Genomic Exploration of the Brain in People Infected with HIV-Recent Progress and the Road Ahead

PURPOSE OF REVIEW

The adult human brain harbors billions of microglia and other myeloid and lymphoid cells highly susceptible to HIV infection and retroviral insertion into the nuclear DNA. HIV infection of the brain is important because the brain is a potentially large reservoir site that may be a barrier to HIV cure strategies and because infection can lead to the development of HIV-associated neurocognitive disorder. To better understand both the central nervous system (CNS) reservoir and how it can cause neurologic dysfunction, novel genomic, epigenomic, transcriptomic, and proteomic approaches need to be employed. Several characteristics of the reservoir are important to learn, including where the virus integrates, whether integrated proviruses are intact or defective, whether integrated proviruses can be reactivated from a latent state to seed ongoing infection, and how this all impacts brain function.

RECENT FINDINGS

Here, we discuss similarities and differences of viral integration sites between brain and blood and discuss evidence for and against the hypothesis that in the absence of susceptible T-lymphocytes in the periphery, the virus housing in the infected brain is not able to sustain a systemic infection. Moreover, microglia from HIV + brains across a wide range of disease severity appear to share one type of common alteration, which is defined by downregulated expression, and repressive chromosomal compartmentalization, for microglial genes regulating synaptic connectivity. Therefore, viral infection of the brain, including in immunocompetent cases with near-normal levels of CD4 blood lymphocytes, could be associated with an early disruption in microglia-dependent neuronal support functions, contributing to cognitive and neurological deficits in people living with HIV.

Co-localization of clusters of TCR-regulated genes with TAD rearrangements

BACKGROUND

Gene expression has long been known to be influenced by the relative proximity of DNA regulatory elements. Topologically associating domains (TADs) are self-interacting genomic regions involved in regulating gene expression by controlling the proximity of these elements. Prior studies of TADs and their biological roles have revealed correlations between TAD changes and cellular differentiation. Here, we used Hi-C and RNA-seq data to correlate TCR-induced changes in TAD structure and gene expression in human CD4+ T cells.

RESULTS

We developed a pipeline, Differentially Expressed Gene Enrichment Finder (DEGEF), that identifies regions of differentially expressed gene enrichment. Using DEGEF, we found that TCR-regulated genes cluster non-uniformly across the genome and that these clusters preferentially localized in regions of TAD rearrangement. Interestingly, clusters of upregulated genes preferentially formed new Hi-C contacts compared to downregulated clusters, suggesting that TCR-activated CD4+ T cells may regulate genes by changing stimulatory contacts rather than inhibitory contacts.

CONCLUSIONS

Our observations support a significant relationship between TAD rearrangements and changes in local gene expression. These findings indicate potentially important roles for TAD rearrangements in shaping their local regulatory environments and thus driving differential expression of nearby genes during CD4+ T cell activation. Moreover, they provide new insights into global mechanisms that regulate gene expression.

3D chromatin reprogramming primes human memory TH2 cells for rapid recall and pathogenic dysfunction

Memory T cells provide long-lasting defense responses through their ability to rapidly reactivate, but how they efficiently "recall" an inflammatory transcriptional program remains unclear. Here, we show that human CD4⁺ memory T helper 2 (T_H2) cells carry a chromatin landscape synergistically reprogrammed at both one-dimensional (1D) and 3D levels to accommodate recall responses, which is absent in naive T cells. In memory T_H2 cells, recall genes were epigenetically primed through the maintenance of transcription-permissive chromatin at distal (super)enhancers organized in long-range 3D chromatin hubs. Precise transcriptional control of key recall genes occurred inside dedicated topologically associating domains ("memory TADs"), in which activation-associated promoter-enhancer interactions were preformed and exploited by AP-1 transcription factors to promote rapid transcriptional induction. Resting memory T_H2 cells from patients with asthma showed premature activation of primed recall circuits, linking aberrant transcriptional control of recall responses to chronic inflammation. Together, our results implicate stable multiscale reprogramming of chromatin organization as a key mechanism underlying immunological memory and dysfunction in T cells.

3D genome organization and epigenetic regulation in autoimmune diseases

Three-dimensional (3D) genomics is an emerging field of research that investigates the relationship between gene regulatory function and the spatial structure of chromatin. Chromatin folding can be studied using chromosome conformation capture (3C) technology and 3C-based derivative sequencing technologies, including chromosome conformation capture-on-chip (4C), chromosome conformation capture carbon copy (5C), and high-throughput chromosome conformation capture (Hi-C), which allow scientists to capture 3D conformations from a single site to the entire genome. A comprehensive analysis of the relationships between various regulatory components and gene function also requires the integration of multi-omics data such as genomics, transcriptomics, and epigenomics. 3D genome folding is involved in immune cell differentiation, activation, and dysfunction and participates in a wide range of diseases, including autoimmune diseases. We describe hierarchical 3D chromatin organization in this review and conclude with characteristics of C-techniques and multi-omics applications of the 3D genome. In addition, we describe the relationship between 3D genome structure and the differentiation and maturation of immune cells and address how changes in chromosome folding contribute to autoimmune diseases.

Recent advances of osteoimmunology research in rheumatoid arthritis: From single-cell omics approach

ABSTRACT

Cellular immune responses as well as generalized and periarticular bone loss are the key pathogenic features of rheumatoid arthritis (RA). Under the pathological conditions of RA, dysregulated inflammation and immune processes tightly interact with skeletal system, resulting in pathological bone damage via inhibition of bone formation or induction of bone resorption. Single-cell omics technologies are revolutionary tools in the field of modern biological research.They enable the display of the state and function of cells in various environments from a single-cell resolution, thus making it conducive to identify the dysregulated molecular mechanisms of bone destruction in RA as well as the discovery of potential therapeutic targets and biomarkers. Here, we summarize the latest findings of single-cell omics technologies in osteoimmunology research in RA. These results suggest that single-cell omics have made significant contributions to transcriptomics and dynamics of specific cells involved in bone remodeling, providing a new direction for our understanding of cellular heterogeneity in the study of osteoimmunology in RA.

Single-cell detection of primary transcripts, their genomic loci and nuclear factors by 3D immuno-RNA/DNA FISH in T cells

Over the past decades, it has become increasingly clear that higher order chromatin folding and organization within the nucleus is involved in the regulation of genome activity and serves as an additional epigenetic mechanism that modulates cellular functions and gene expression programs in diverse biological processes. In particular, dynamic allelic interactions and nuclear locations can be of functional importance during the process of lymphoid differentiation and the regulation of immune responses. Analyses of the proximity between chromatin and/or nuclear regions can be performed on populations of cells with high-throughput sequencing approaches such as chromatin conformation capture ("3C"-based) or DNA adenine methyltransferase identification (DamID) methods, or, in individual cells, by the simultaneous visualization of genomic loci, their primary transcripts and nuclear compartments within the 3-dimensional nuclear space using Fluorescence In Situ Hybridization (FISH) and immunostaining. Here, we present a detailed protocol to simultaneously detect nascent RNA transcripts (3D RNA FISH), their genomic loci (3D DNA FISH) and/or their chromosome territories (CT paint DNA FISH) combined with the antibody-based detection of various nuclear factors (immunofluorescence). We delineate the application and effectiveness of this robust and reproducible protocol in several murine T lymphocyte subtypes (from differentiating thymic T cells, to activated splenic and peripheral T cells) as well as other murine cells, including embryonic stem cells, B cells, megakaryocytes and macrophages.

HIV integration in the human brain is linked to microglial activation and 3D genome remodeling

To explore genome organization and function in the HIV-infected brain, we applied single-nuclei transcriptomics, cell-type-specific chromosomal conformation mapping, and viral integration site sequencing (IS-seq) to frontal cortex from individuals with encephalitis (HIVE) and without (HIV+). Derepressive changes in 3D genomic compartment structures in HIVE microglia were linked to the transcriptional activation of interferon (IFN) signaling and cell migratory pathways, while transcriptional downregulation and repressive compartmentalization of neuronal health and signaling genes occurred in both HIVE and HIV+ microglia. IS-seq recovered 1,221 brain integration sites showing distinct genomic patterns compared with peripheral lymphocytes, with enrichment for sequences newly mobilized into a permissive chromatin environment after infection. Viral transcription occurred in a subset of highly activated microglia comprising 0.33% of all nuclei in HIVE brain. Our findings point to disrupted microglia-neuronal interactions in HIV and link retroviral integration to remodeling of the microglial 3D genome during infection.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Phenotypic selection during laboratory evolution of yeast populations leads to a genome-wide sustainable chromatin compaction shift

In a previous study, we have shown how microbial evolution has resulted in a persistent reduction in expression after repeatedly selecting for the lowest PGAL1-YFP-expressing cells. Applying the ATAC-seq assay on samples collected from this 28-day evolution experiment, here we show how genome-wide chromatin compaction changes during evolution under selection pressure. We found that the chromatin compaction was altered not only on GAL network genes directly impacted by the selection pressure, showing an example of selection-induced non-genetic memory, but also at the whole-genome level. The GAL network genes experienced chromatin compaction accompanying the reduction in PGAL1-YFP reporter expression. Strikingly, the fraction of global genes with differentially compacted chromatin states accounted for about a quarter of the total genome. Moreover, some of the ATAC-seq peaks followed well-defined temporal dynamics. Comparing peak intensity changes on consecutive days, we found most of the differential compaction to occur between days 0 and 3 when the selection pressure was first applied, and between days 7 and 10 when the pressure was lifted. Among the gene sets enriched for the differential compaction events, some had increased chromatin availability once selection pressure was applied and decreased availability after the pressure was lifted (or vice versa). These results intriguingly show that, despite the lack of targeted selection, transcriptional availability of a large fraction of the genome changes in a very diverse manner during evolution, and these changes can occur in a relatively short number of generations.

Rearrangement of T Cell genome architecture regulates GVHD

WAPL, cohesin's DNA release factor, regulates three-dimensional (3D) chromatin architecture. The 3D chromatin structure and its relevance to mature T cell functions is not well understood. We show that in vivo lymphopenic expansion, and alloantigen-driven proliferation, alters the 3D structure and function of the genome in mature T cells. Conditional deletion of WAPL, cohesin's DNA release factor, in T cells reduced long-range genomic interactions and altered chromatin A/B compartments and interactions within topologically associating domains (TADs) of the chromatin in T cells at baseline. WAPL deficiency in T cells reduced loop extensions, changed expression of cell cycling genes and reduced proliferation following in vitro and in vivo stimulation, and reduced severity of graft-versus-host disease (GVHD) following experimental allogeneic hematopoietic stem cell transplantation. These data collectively characterize 3D genomic architecture of T cells in vivo and demonstrate biological and clinical implications for its disruption by cohesin release factor WAPL.

Multi-omics analysis of naïve B cells of patients harboring the C104R mutation in TACI

Common variable immunodeficiency (CVID) is the most prevalent form of symptomatic primary immunodeficiency in humans. The genetic cause of CVID is still unknown in about 70% of cases. Ten percent of CVID patients carry heterozygous mutations in the tumor necrosis factor receptor superfamily member 13B gene (TNFRSF13B), encoding TACI. Mutations in TNFRSF13B alone may not be sufficient for the development of CVID, as 1% of the healthy population carry these mutations. The common hypothesis is that TACI mutations are not fully penetrant and additional factors contribute to the development of CVID. To determine these additional factors, we investigated the perturbations of transcription factor (TF) binding and the transcriptome profiles in unstimulated and CD40L/IL21-stimulated naïve B cells from CVID patients harboring the C104R mutation in TNFRSF13B and compared them to their healthy relatives with the same mutation. In addition, the proteome of stimulated naïve B cells was investigated. For functional validation, intracellular protein concentrations were measured by flow cytometry. Our analysis revealed 8% less accessible chromatin in unstimulated naïve B cells and 25% less accessible chromatin in class-switched memory B cells from affected and unaffected TACI mutation carriers compared to healthy donors. The most enriched TF binding motifs in TACI mutation carriers involved members from the ETS, IRF, and NF-κB TF families. Validation experiments supported dysregulation of the NF-κB and MAPK pathways. In steady state, naïve B cells had increased cell death pathways and reduced cell metabolism pathways, while after stimulation, enhanced immune responses and decreased cell survival were detected. Using a multi-omics approach, our findings provide valuable insights into the impaired biology of naïve B cells from TACI mutation carriers.

3D Genome Organization as an Epigenetic Determinant of Transcription Regulation in T Cells

In the heart of innate and adaptive immunity lies the proper spatiotemporal development of several immune cell lineages. Multiple studies have highlighted the necessity of epigenetic and transcriptional regulation in cell lineage specification. This mode of regulation is mediated by transcription factors and chromatin remodelers, controlling developmentally essential gene sets. The core of transcription and epigenetic regulation is formulated by different epigenetic modifications determining gene expression. Apart from “classic” epigenetic modifications, 3D chromatin architecture is also purported to exert fundamental roles in gene regulation. Chromatin conformation both facilitates cell-specific factor binding at specified regions and is in turn modified as such, acting synergistically. The interplay between global and tissue-specific protein factors dictates the epigenetic landscape of T and innate lymphoid cell (ILC) lineages. The expression of global genome organizers such as CTCF, YY1, and the cohesin complexes, closely cooperate with tissue-specific factors to exert cell type-specific gene regulation. Special AT-rich binding protein 1 (SATB1) is an important tissue-specific genome organizer and regulator controlling both long- and short-range chromatin interactions. Recent indications point to SATB1’s cooperation with the aforementioned factors, linking global to tissue-specific gene regulation. Changes in 3D genome organization are of vital importance for proper cell development and function, while disruption of this mechanism can lead to severe immuno-developmental defects. Newly emerging data have inextricably linked chromatin architecture deregulation to tissue-specific pathophysiological phenotypes. The combination of these findings may shed light on the mechanisms behind pathological conditions.

Cohesin-Mediated Chromatin Interactions and Autoimmunity

Proper physiological functioning of any cell type requires ordered chromatin organization. In this context, cohesin complex performs important functions preventing premature separation of sister chromatids after DNA replication. In partnership with CCCTC-binding factor, it ensures insulator activity to organize enhancers and promoters within regulatory chromatin. Homozygous mutations and dysfunction of individual cohesin proteins are embryonically lethal in humans and mice, which limits in vivo research work to embryonic stem cells and progenitors. Conditional alleles of cohesin complex proteins have been generated to investigate their functional roles in greater detail at later developmental stages. Thus, genome regulation enabled by action of cohesin proteins is potentially crucial in lineage cell development, including immune homeostasis. In this review, we provide current knowledge on the role of cohesin complex in leukocyte maturation and adaptive immunity. Conditional knockout and shRNA-mediated inhibition of individual cohesin proteins in mice demonstrated their importance in haematopoiesis, adipogenesis and inflammation. Notably, these effects occur rather through changes in transcriptional gene regulation than through expected cell cycle defects. This positions cohesin at the crossroad of immune pathways including NF-kB, IL-6, and IFNγ signaling. Cohesin proteins emerged as vital regulators at early developmental stages of thymocytes and B cells and after antigen challenge. Human genome-wide association studies are remarkably concordant with these findings and present associations between cohesin and rheumatoid arthritis, multiple sclerosis and HLA-B27 related chronic inflammatory conditions. Furthermore, bioinformatic prediction based on protein-protein interactions reveal a tight connection between the cohesin complex and immune relevant processes supporting the notion that cohesin will unearth new clues in regulation of autoimmunity.

Enhancer-silencer transitions in the human genome

Dual-function regulatory elements (REs), acting as enhancers in some cellular contexts and as silencers in others, have been reported to facilitate the precise gene regulatory response to developmental signals in Drosophila melanogaster. However, with few isolated examples detected, dual-function REs in mammals have yet to be systematically studied. We herein investigated this class of REs in the human genome and profiled their activity across multiple cell types. Focusing on enhancer–silencer transitions specific to the development of T cells, we built an accurate deep learning classifier of REs and identified about 12,000 silencers active in primary peripheral blood T cells that act as enhancers in embryonic stem cells. Compared with regular silencers, these dual-function REs are evolving under stronger purifying selection and are enriched for mutations associated with disease phenotypes and altered gene expression. In addition, they are enriched in the loci of transcriptional regulators, such as transcription factors (TFs) and chromatin remodeling genes. Dual-function REs consist of two intertwined but largely distinct sets of binding sites bound by either activating or repressing TFs, depending on the type of RE function in a given cell line. This indicates the recruitment of different TFs for different regulatory modes and a complex DNA sequence composition of these REs with dual activating and repressive encoding. With an estimated >6% of cell type–specific human silencers acting as dual-function REs, this overlooked class of REs requires a specific investigation on how their inherent functional plasticity might be a contributing factor to human diseases.

Reduced chromatin accessibility to CD4 T cell super-enhancers encompassing susceptibility loci of rheumatoid arthritis

Background

Rheumatoid arthritis (RA) is an inflammatory disease that manifests as a preclinical stage of systemic autoimmunity followed by chronic progressive synovitis. Disease-associated genetic SNP variants predominantly map to non-coding, regulatory regions of functional importance in CD4 T cells, implicating these cells as key regulators. A better understanding of the epigenome of CD4 T cells holds the promise of providing information on the interaction between genetic susceptibility and exogenous factors.

Methods

We mapped regions of chromatin accessibility using ATAC-seq in peripheral CD4 T cell subsets of patients with RA (n=18) and compared them to T cells from patients with psoriatic arthritis (n=11) and age-matched healthy controls (n=10). Transcripts of selected genes were quantified using qPCR.

Findings

RA-associated epigenetic signatures were identified that in part overlapped between central and effector memory CD4 T cells and that were to a lesser extent already present in naïve cells. Sites more accessible in RA were highly enriched for the motif of the transcription factor (TF) CTCF suggesting differences in the three-dimensional chromatin structure. Unexpectedly, sites with reduced chromatin accessibility were enriched for motifs of TFs pertinent for T cell function. Most strikingly, super-enhancers encompassing RA-associated SNPs were less accessible. Analysis of selected transcripts and published DNA methylation patterns were consistent with this finding. The preferential loss in accessibility at these super-enhancers was seen in patients with high and low disease activity and on a variety of immunosuppressive treatment modalities.

Interpretation

Disease-associated genes are epigenetically less poised to respond in CD4 T cells from patients with established RA.

Funding

This work was supported by I01 BX001669 from the Veterans Administration.

Genome Editing Using CRISPR-Cas9 and Autoimmune Diseases: A Comprehensive Review

Autoimmune diseases are disorders that destruct or disrupt the body's own tissues by its own immune system. Several studies have revealed that polymorphisms of multiple genes are involved in autoimmune diseases. Meanwhile, gene therapy has become a promising approach in autoimmune diseases, and clustered regularly interspaced palindromic repeats and CRISPR-associated protein 9 (CRISPR-Cas9) has become one of the most prominent methods. It has been shown that CRISPR-Cas9 can be applied to knock out proprotein convertase subtilisin/kexin type 9 (PCSK9) or block PCSK9, resulting in lowering low-density lipoprotein cholesterol. In other studies, it can be used to treat rare diseases such as ornithine transcarbamylase (OTC) deficiency and hereditary tyrosinemia. However, few studies on the treatment of autoimmune disease using CRISPR-Cas9 have been reported so far. In this review, we highlight the current and potential use of CRISPR-Cas9 in the management of autoimmune diseases. We summarize the potential target genes for immunomodulation using CRISPR-Cas9 in autoimmune diseases including rheumatoid arthritis (RA), inflammatory bowel diseases (IBD), systemic lupus erythematosus (SLE), multiple sclerosis (MS), type 1 diabetes mellitus (DM), psoriasis, and type 1 coeliac disease. This article will give a new perspective on understanding the use of CRISPR-Cas9 in autoimmune diseases not only through animal models but also in human models. Emerging approaches to investigate the potential target genes for CRISPR-Cas9 treatment may be promising for the tailored immunomodulation of some autoimmune diseases in the near future.

Fine mapping with epigenetic information and 3D structure

Since 2005, thousands of genome-wide association studies (GWAS) have been published, identifying hundreds of thousands of genetic variants that increase risk of complex traits such as autoimmune diseases. This wealth of data has the potential to improve patient care, through personalized medicine and the identification of novel drug targets. However, the potential of GWAS for clinical translation has not been fully achieved yet, due to the fact that the functional interpretation of risk variants and the identification of causal variants and genes are challenging. The past decade has seen the development of great advances that are facilitating the overcoming of these limitations, by utilizing a plethora of genomics and epigenomics tools to map and characterize regulatory elements and chromatin interactions, which can be used to fine map GWAS loci, and advance our understanding of the biological mechanisms that cause disease.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

Association between CTLA-4 gene polymorphism and risk of rheumatoid arthritis: a meta-analysis

Cytotoxic T lymphocyte-associated protein 4 (CTLA-4) gene polymorphisms may be involved in the risk of Rheumatoid arthritis (RA). However, evidence for the association remains controversial. Therefore, we performed a meta-analysis to confirm the relationship between CTLA-4 gene polymorphisms and RA. The pooled odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to assess the strength of association. Stratified analysis was conducted by ethnicity. In total, 66 case-control studies including 21681 cases and 23457 controls were obtained. For rs3087243 polymorphism, significant association was detected in Asians (A vs. G: OR=0.77, 95%CI=0.65-0.90, P=0.001; AA vs. GG: OR=0.67, 95%CI=0.48-0.94, P=0.02) and Caucasians (A vs. G: OR=0.89, 95%CI=0.86-0.93, P<0.00001; AA vs. GG: OR=0.81, 95%CI=0.75-0.88, P<0.00001). For rs231775 polymorphism, significant association was observed in the overall (G vs. A: OR =1.16, 95%CI=1.08-1.25, P<0.0001; GG vs. AA: OR=1.29, 95%CI=1.12-1.50, P=0.0006), and in Asians (G vs. A: OR=1.27, 95%CI=1.10-1.47, P=0.001; GG vs. AA: OR=1.58, 95%CI=1.24-2.01, P=0.0002), but not in Caucasians. However, there was no association between rs5742909 polymorphism and RA. This meta-analysis confirmed that rs3087243 and rs231775 polymorphism were associated with the risk of RA in both overall population and ethnic-specific analysis, but there was no association between rs5742909 polymorphism and RA risk.

Genome editing to define the function of risk loci and variants in rheumatic disease

Discoveries in human genetic studies have revolutionized our understanding of complex rheumatic and autoimmune diseases, including the identification of hundreds of genetic loci and single nucleotide polymorphisms that potentially predispose individuals to disease. However, in most cases, the exact disease-causing variants and their mechanisms of action remain unresolved. Functional follow-up of these findings is most challenging for genomic variants that are in non-coding genomic regions, where the large majority of common disease-associated variants are located, and/or that probably affect disease progression via cell type-specific gene regulation. To deliver on the therapeutic promise of human genetic studies, defining the mechanisms of action of these alleles is essential. Genome editing technology, such as CRISPR-Cas, has created a vast toolbox for targeted genetic and epigenetic modifications that presents unprecedented opportunities to decipher disease-causing loci, genes and variants in autoimmunity. In this Review, we discuss the past 5-10 years of progress in resolving the mechanisms underlying rheumatic disease-associated alleles, with an emphasis on how genomic editing techniques can enable targeted dissection and mechanistic studies of causal autoimmune risk variants.

Mapping the landscape of chromatin dynamics during naïve CD4+ T-cell activation

T-cell activation induces context-specific gene expression programs that promote energy generation and biosynthesis, progression through the cell cycle and ultimately cell differentiation. The aim of this study was to apply the omni ATAC-seq method to characterize the landscape of chromatin changes induced by T-cell activation in mature naïve CD4+ T-cells. Using a well-established ex vivo protocol of canonical T-cell receptor signaling, we generated genome-wide chromatin maps of naïve T-cells from pediatric donors in quiescent or recently activated states. We identified thousands of individual chromatin accessibility peaks that are associated with T-cell activation, the majority of which were annotated intronic and intergenic enhancer regions. A core set of 3268 gene promoters underwent chromatin remodeling and concomitant changes in gene expression in response to activation, and were enriched in multiple pathways controlling cell cycle regulation, metabolism, inflammatory response genes and cell survival. Leukemia inhibitory factor (LIF) was among those factors that gained the highest accessibility and expression, in addition to IL2-STAT5 dependent chromatin remodeling in the T-cell activation response. Using publicly available data we found the chromatin response was far more dynamic at 24-h compared with 72-h post-activation. In total 546 associations were reproduced at both time-points with similar strength of evidence and directionality of effect. At the pathways level, the IL2-STAT5, KRAS signalling and UV response pathways were replicable at both time-points, although differentially modulated from 24 to 72 h post-activation.

Deciphering the Complexity of 3D Chromatin Organization Driving Lymphopoiesis and Lymphoid Malignancies

Proper lymphopoiesis and immune responses depend on the spatiotemporal control of multiple processes, including gene expression, DNA recombination and cell fate decisions. High-order 3D chromatin organization is increasingly appreciated as an important regulator of these processes and dysregulation of genomic architecture has been linked to various immune disorders, including lymphoid malignancies. In this review, we present the general principles of the 3D chromatin topology and its dynamic reorganization during various steps of B and T lymphocyte development and activation. We also discuss functional interconnections between architectural, epigenetic and transcriptional changes and introduce major key players of genomic organization in B/T lymphocytes. Finally, we present how alterations in architectural factors and/or 3D genome organization are linked to dysregulation of the lymphopoietic transcriptional program and ultimately to hematological malignancies.

Immune Regulation in Time and Space: The Role of Local- and Long-Range Genomic Interactions in Regulating Immune Responses

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.

Perspectives on the Genetic Associations of Ankylosing Spondylitis

Ankylosing spondylitis (AS) is a common form of inflammatory spinal arthritis with a complex polygenic aetiology. Genome-wide association studies have identified more than 100 loci, including some involved in antigen presentation (HLA-B27, ERAP1, and ERAP2), some in Th17 responses (IL6R, IL23R, TYK2, and STAT3), and others in macrophages and T-cells (IL7R, CSF2, RUNX3, and GPR65). Such observations have already helped identify potential new therapies targeting IL-17 and GM-CSF. Most AS genetic associations are not in protein-coding sequences but lie in intergenic regions where their direct relationship to particular genes is difficult to assess. They most likely reflect functional polymorphisms concerned with cell type-specific regulation of gene expression. Clarifying the nature of these associations should help to understand the pathogenic pathways involved in AS better and suggest potential cellular and molecular targets for drug therapy. However, even identifying the precise mechanisms behind the extremely strong HLA-B27 association with AS has so far proved elusive. Polygenic risk scores (using all the known genetic associations with AS) can be effective for the diagnosis of AS, particularly where there is a relatively high pre-test probability of AS. Genetic prediction of disease outcomes and response to biologics is not currently practicable.

Characterisation of CD4+ T-cell subtypes using single cell RNA sequencing and the impact of cell number and sequencing depth

CD4+ T-cells represent a heterogeneous collection of specialised sub-types and are a key cell type in the pathogenesis of many diseases due to their role in the adaptive immune system. By investigating CD4+ T-cells at the single cell level, using RNA sequencing (scRNA-seq), there is the potential to identify specific cell states driving disease or treatment response. However, the impact of sequencing depth and cell numbers, two important factors in scRNA-seq, has not been determined for a complex cell population such as CD4+ T-cells. We therefore generated a high depth, high cell number dataset to determine the effect of reduced sequencing depth and cell number on the ability to accurately identify CD4+ T-cell subtypes. Furthermore, we investigated T-cell signatures under resting and stimulated conditions to assess cluster specific effects of stimulation. We found that firstly, cell number has a much more profound effect than sequencing depth on the ability to classify cells; secondly, this effect is greater when cells are unstimulated and finally, resting and stimulated samples can be combined to leverage additional power whilst still allowing differences between samples to be observed. While based on one individual, these results could inform future scRNA-seq studies to ensure the most efficient experimental design.